Ben binomial cevap değişkeni ve kategorik bir yordayıcı ile bir glmm çalıştırıyorum. Rasgele efekt, veri toplama için kullanılan iç içe tasarım tarafından verilir. Veriler şöyle görünür:
m.gen1$treatment
[1] sucrose control protein control no_injection .....
Levels: no_injection control sucrose protein
m.gen1$emergence
[1] 1 0 0 1 0 1 1 1 1 1 1 0 0....
> m.gen1$nest
[1] 1 1 1 2 2 3 3 3 3 4 4 4 .....
Levels: 1 2 3 4 5 6 8 10 11 13 15 16 17 18 20 22 24
Çalıştığım ilk model şöyle görünüyor
m.glmm.em.<-glmer(emergence~treatment + (1|nest),family=binomial,data=m.gen1)
Bunun gibi görünen iki uyarı alıyorum:
Warning messages:
1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 0.0240654 (tol = 0.001, component 4)
2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model is nearly unidentifiable: large eigenvalue ratio
- Rescale variables?
Model özeti, tedavilerden birinin burada görebileceğiniz alışılmadık derecede büyük bir standart hataya sahip olduğunu gösterir:
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 2.565 1.038 2.472 0.0134 *
treatmentcontrol -1.718 1.246 -1.378 0.1681
treatmentsucrose 16.863 2048.000 0.008 0.9934
treatmentprotein -1.718 1.246 -1.378 0.1681
Glmer kontrolünden farklı optimize edicileri ve diğer paketlerden işlevleri denedim ve benzer bir çıktı alıyorum. Rastgele efekti görmezden gelerek modeli çalıştırdım ve sorun devam ediyor. Verileri araştırırken yüksek bir Std ile tedavi olduğunu fark ettim. hata yalnızca yanıt değişkeninde başarılı olur. Sadece soruna neden olup olmadığını kontrol etmek için o tedavi için bir "başarısızlık" ile sahte bir veri noktası ekledim ve model sorunsuz çalışır ve makul standart hata verir. Bunu burada görebilirsiniz:
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.4090 1.6712 2.040 0.0414 *
treatmentcontrol -1.8405 1.4290 -1.288 0.1978
treatmentsucrose -0.2582 1.6263 -0.159 0.8738
treatmentprotein -2.6530 1.5904 -1.668 0.0953 .
Sezgimin iyi bir tahminin önlenmesi için bu tedavinin başarısızlığı konusunda doğru olup olmadığını ve bu soruna nasıl çözüm bulabileceğimi merak ediyordum.
Şimdiden teşekkürler!