Bir karma efekt modeli kullanarak bazı verileri analiz etmeye çalışıyorum. Topladığım veriler zaman içinde farklı genotipteki bazı genç hayvanların ağırlığını temsil ediyor.
Burada önerilen yaklaşımı kullanıyorum: https://gribblelab.wordpress.com/2009/03/09/repeated-measures-anova-using-r/
Özellikle 2 numaralı çözümü kullanıyorum
Bu yüzden benim gibi bir şey var
require(nlme)
model <- lme(weight ~ time * Genotype, random = ~1|Animal/time,
data=weights)
av <- anova(model)
Şimdi, bazı çoklu karşılaştırmalar yapmak istiyorum. Kullanarak multcomp
yapabilirim:
require(multcomp)
comp.geno <- glht(model, linfct=mcp(Genotype="Tukey"))
print(summary(comp.geno))
Ve tabii ki, zamanla aynı şeyi yapabilirim.
İki sorum var:
mcp
Zaman ve Genotip arasındaki etkileşimi görmek için nasıl kullanırım ?Koşarken
glht
bu uyarıyı alıyorum:covariate interactions found -- default contrast might be inappropriate
Bunun anlamı ne? Güvenle görmezden gelebilir miyim? Ya da önlemek için ne yapmalıyım?
EDIT: Bu PDF'i şöyle derdim :
Bu durumda ilgi parametrelerini otomatik olarak belirlemek mümkün olmadığından, çoklu hesaptaki mcp () varsayılan olarak yalnızca temel etkiler için karşılaştırmalar oluşturacak, değişkenleri ve etkileşimleri dikkate almayacaktır . 1.1-2 sürümünden bu yana, sırasıyla argümanları interaction_average = TRUE ve covariate_average = TRUE argümanlarını kullanarak etkileşimli terimler ve değişkenler üzerinde ortalama belirleyebilir, oysaki 1.0-0'den eski sürümler etkileşimli terimler üzerinden otomatik olarak ortalaması alınır. Ancak, kullanıcılara istedikleri kontrast kümesini manuel olarak yazmalarını öneririz.Bunu, genellikle yüksek dereceli etkileşim terimleri olan modellerde meydana gelen, varsayılan kontrastın ne olduğuna dair bir şüphe olduğunda bunu yapmanız gerekir. Bu konuda daha fazla tartışma ve örnekler için Hsu (1996), Bölüm ~ 7 ve Searle (1971), Bölüm ~ 7.3'e bakınız.
Bu kitaplara erişimim yok, ama belki burada birileri var?