bazı veriler Ar GLM ile% 95 güven aralığı (..) Get İşte
dat = data.frame(y = c(9,7,7,7,5,6,4,6,3,5,1,5), x = c(1,1,2,2,3,3,4,4,5,5,6,6), color = rep(c('a','b'),6))
ve bu verilerin arsa fonksiyonu ile bir model çalıştırırken
require(ggplot)
ggplot(dat, aes(x=x,y=y, color=color)) + geom_point() + geom_smooth(method='lm')
isterlerse vardır MCMCglmm()
...
require(MCMCglmm)
summary(MCMCglmm(fixed = y~x/color, data=dat))
ben bana izin tahmin için alt ve üst% 95 aralığını alırsanız iki yansımalar es (renk = a ve renk = b) önemli ölçüde farklıdır.
bu çıkışa bakarak ...
summary(glm(y~x/color, data=dat))
... Ben güven aralığını göremiyorum!
Sorum şu:
fonksiyonunu glm()
kullanırken nasıl tahminler için bu alt ve üst% 95 aralığında bir güven var?
** MASS ** paketinde 'confint' işlevi var, inanıyorum. – joran
... Ayrıca, güven aralıkları karşılaştırarak istatistiksel önemi değerlendirerek dikkatli olun. Bazı önemli sonuçları kaçırırsınız. – joran
Bu durumda, sadece "x: colorb" katsayısı için güven sınırlarının 0'a eşit olup olmadığına bakın. –