Bu veri kümesinde bazı belirli zıtlıklar ardından R tekrarlanan önlemler Anova çalıştırmak için çalışıyorum. Bence doğru yaklaşım Anova()
araç paketinden kullanmak olacaktır
.
Verilerimi verileri ?Anova
kullanarak
alınan örnekle OBrienKaiser
açıklayalım (Not: Örnekte cinsiyet faktörünü göz ardı ettim):
Denekler faktörü, tedavi (3 seviye: kontrol, A, B) ve tekrarlanan 2 arasında bir tasarıma sahibiz -Ölçümler (denekler içinde) faktörleri, faz (3 seviye: ön test, son test, takip) ve saat (5 seviye: 1 ila 5).
Standart ANOVA tablosu (örnek (Anova) 'dan farklı olarak Tip 3 Kareler Toplamı'na geçtim, alanımın istediği bu):
require(car)
phase <- factor(rep(c("pretest", "posttest", "followup"), c(5, 5, 5)),
levels=c("pretest", "posttest", "followup"))
hour <- ordered(rep(1:5, 3))
idata <- data.frame(phase, hour)
mod.ok <- lm(cbind(pre.1, pre.2, pre.3, pre.4, pre.5, post.1, post.2, post.3, post.4, post.5, fup.1, fup.2, fup.3, fup.4, fup.5) ~ treatment, data=OBrienKaiser)
av.ok <- Anova(mod.ok, idata=idata, idesign=~phase*hour, type = 3)
summary(av.ok, multivariate=FALSE)
Şimdi, en yüksek dereceli etkileşimin önemli olacağını düşünün (durum böyle değildir) ve aşağıdaki kontrastlarla daha fazla araştırmak istiyoruz:
1 ve 2 saatleri ile saat 3 (kontrast 1) ile saat 1 ve 2 saatleri arasında bir fark var mı? tedavi koşullarında saat 4 ve 5'e (kontrast 2) karşı mı (birlikte A&B)?
Başka bir deyişle, bu karşıtlıkları nasıl belirleyebilirim:
((treatment %in% c("A", "B")) & (hour %in% 1:2))
karşı((treatment %in% c("A", "B")) & (hour %in% 3))
((treatment %in% c("A", "B")) & (hour %in% 1:2))
karşı((treatment %in% c("A", "B")) & (hour %in% 4:5))
Benim fikrim, gerekli olmayan tedavi koşulu (kontrol) ile uyumlu başka bir ANOVA çalıştırmak olacaktır:
mod2 <- lm(cbind(pre.1, pre.2, pre.3, pre.4, pre.5, post.1, post.2, post.3, post.4, post.5, fup.1, fup.2, fup.3, fup.4, fup.5) ~ treatment, data=OBrienKaiser, subset = treatment != "control")
av2 <- Anova(mod2, idata=idata, idesign=~phase*hour, type = 3)
summary(av2, multivariate=FALSE)
Bununla birlikte, hala 1 ve 2 saatleri 3 ve 1 ve 2 ve 4 ve 5 ile karşılaştırarak uygun konu içi kontrast matrisinin nasıl kurulacağına dair hiçbir fikrim yok. Gereksiz tedavi grubunun ihmal edilmesinin, genel hata terimini değiştirdiği için gerçekten iyi bir fikir olup olmadığından emin değilim.
Gitmeden önce de gitmeyi Anova()
düşünüyordum lme
. Ancak, ANOVA ders kitabı anove(lme)
ile standart ANOVA'daki olası negatif varyanslardan (bunlara izin verilmezlme
) kaynaklanan geri dönüşler arasında F ve p değerleri arasında küçük farklılıklar vardır . Buna bağlı olarak, biri bana gls
tekrarlanan ölçümler ANOVA'nın takılmasına izin veren bir işaret etti , ancak kontrast argümanı yok.
Açıklığa kavuşturmak için: İstenen kontrastların anlamlı olup olmadığını yanıtlayan bir F veya t testi (tip III toplam kareler kullanarak) istiyorum.
Güncelleme:
Zaten cevap gelmedi, R-yardımına çok benzer bir soru soruldu .
Benzer sorular R-help'de de bir süre önce sorulmuştu. Ancak, cevaplar da sorunu çözmedi.
Güncelleme (2015):
Bu soru hala bir aktivite ürettiğinden, tezlerin ve temelde diğer tüm kontrastların belirlenmesi , afeks vinyetinde tarif edildiği gibi paketle afex
birlikte paketle nispeten kolay bir şekilde yapılabilir .lsmeans
treatment
, 3) her kişi ortalamasıprePostFup
için, 4) seviyeleri üzerinden her kişi için ortalama saat 1, 2 (= veri grubu 1) saatin yanı sıra 3,4 (= veri grubu 2), 5) 2 bağımlı grup için t-testi uygular. Maxwell ve Delaney (2004) ve Kirk (1995) iç tasarımlarda toplanmış bir hata terimiyle tezat oluşturmayı caydırdığından, bu basit bir alternatif olabilir.