Bu daha kavramsal bir soru, ama kullandığım R
zaman içindeki paketlerden bahsedeceğim R
. Amaç, öngörme amacıyla doğrusal bir modele uymak ve daha sonra rasgele etkilerin bulunamayacağı yerlerde tahminler yapmak ise, karışık efektler modeli kullanmanın herhangi bir faydası var mı yoksa bunun yerine sabit efekt modeli mi kullanılmalı?
Örneğin, ağırlık ve boy ile ilgili başka bilgilerle ilgili verilerim varsa ve aşağıdaki modeli kullanarak lme4
, öznenin bir faktör olduğu durumlarda seviyeleri ():
mod1 <- lmer(weight ~ height + age + (1|subject), data=df, REML=F)
Daha sonra yeni boy ve yaş verilerini kullanarak modelden ağırlığı tahmin etmek istiyorum. Açıkçası, orijinal verilerdeki denek sapması modelde yakalanmaktadır, ancak bu bilgiyi tahminde kullanmak mümkün müdür? Diyelim ki bazı yeni boy ve yaş verilerim var ve ağırlığı tahmin etmek istiyorum, bunu şu şekilde yapabilirim:
predict(mod1,newdata=newdf) # newdf columns for height, age, subject
Bu kullanılır predict.merMod
ve ben (yeni) konular için bir sütun dahil edebilir newdf
veya set re.form =~0
. İlk örnekte, modelin 'yeni' konu faktörleri ile ne yaptığı açık değildir ve ikinci örnekte modelde yakalanan özne sapması, tahmin için dikkate alınmayacak mı (ortalama üzerinden)?
Her iki durumda da, sabit etkili doğrusal modelin daha uygun olabileceği bana benziyor. Gerçekten, eğer anlayışım doğruysa, o zaman sabit bir etki modeli, tahminte rastgele etki kullanılmıyorsa, karma modelle aynı değerleri tahmin etmelidir. Durum böyle mi olmalı? İçinde R
değil, örneğin:
mod1 <- lmer(weight ~ height + age + (1|subject), data=df, REML=F)
predict(mod1,newdata=newdf, re.form=~0) # newdf columns for height, age, subject
aşağıdakilere farklı sonuçlar verir:
mod2 <- lm(weight ~ height + age, data=df)
predict(mod2,newdata=newdf) # newdf columns for height, age