范例程式码

data(anorexia, package = "MASS")
anorex.1 <- glm(Postwt ~ Prewt + Treat + offset(Prewt),
            family = gaussian, data = anorexia)
confint(anorex.1)


这产生

                 2.5 %     97.5 %
(Intercept) 23.5253133 76.0169047
Prewt       -0.8814505 -0.2496272
TreatCont   -7.8082428 -0.3858882
TreatFT      0.3818011  8.7443242


我现在想访问这些置信区间,例如,创建4个变量。首选的方法是什么?我知道我可以用

Intercept <- c(confint(anorex.1)[1],confint(anorex.1)[5])
Prewt <- c(confint(anorex.1)[2],confint(anorex.1)[6])


等等。。。但是这似乎很笨拙。

最佳答案

将结果矩阵分配给一个名称,然后按名称提取行:

resCI <- confint(anorex.1)
IntCI <- resCI[ "(Intercept)", ]
PrewtCI <- resCI[ "Prewt", ]

关于r - 访问由confint()返回的对象的元素,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/29904524/

10-09 08:55