glm yaptım ve sadece her bir katsayı standart hatalarını ayıklamak istiyorum. İnternette se.coef() işlevini gördüm, ancak çalışmıyor, "Error: could not find function "se.coef"".Standart hataları ayıklama glm
Bulunduğunuz bilgiler, summary() tarafından döndürülen coefficients nesnesinde depolanır. Sen thusly ayıklayabileceğiniz:
#Example from ?glm
counts <- c(18,17,15,20,10,20,25,13,12)
outcome <- gl(3,1,9)
treatment <- gl(3,3)
print(d.AD <- data.frame(treatment, outcome, counts))
glm.D93 <- glm(counts ~ outcome + treatment, family=poisson())
#coefficients has the data of interest
> summary(glm.D93)$coefficients
Estimate Std. Error z value Pr(>|z|)
(Intercept) 3.044522e+00 0.1708987 1.781478e+01 5.426767e-71
outcome2 -4.542553e-01 0.2021708 -2.246889e+00 2.464711e-02
outcome3 -2.929871e-01 0.1927423 -1.520097e+00 1.284865e-01
treatment2 1.337909e-15 0.2000000 6.689547e-15 1.000000e+00
treatment3 1.421085e-15 0.2000000 7.105427e-15 1.000000e+00
#So extract the second column
> summary(glm.D93)$coefficients[, 2]
(Intercept) outcome2 outcome3 treatment2 treatment3
0.1708987 0.2021708 0.1927423 0.2000000 0.2000000
summary(glm.D93)$coefficients[, 2] döndürülen her şeyin hızlı bir inceleme için names(summary(glm.D93)) bir göz atın. Eğer summary.glm kontrol ederek daha fazla detay bulabilirsiniz, bu özel bir hesaplama görmek istiyorsanız, bu ayrıntı seviyesi muhtemelen her zaman gerekmiyor olsa da, < 3 istatistikleri sürece.
'glm.D93' nesnesinde saklanan standart hatalar var mı? Str() 'yi kullanarak göz atamadım. Veya 'summary()' hatalarını açıkça hesaplıyor mu? –
@ mindless.panda - AFAIK doğrudan "summary.glm" ile hesaplanır. Fonksiyonu sans '()' konsoluna yazarsanız ve daha sonra yaklaşık 25 satır aşağı kaydırırsanız, nerede hesaplandığını görürsünüz. – Chase
se.coef() aslında çalışır. Ancak bu temel pakette değil: {arm} paketinde: http://www.inside-r.org/packages/cran/arm/docs/se.ranef
Bazı verileri ve örnek kodu koymak için yardım edin. – screechOwl