范例程式码
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/