r - 使用 rms lrm 对象获得限制

Question

我无法从 rms 包中获得逻辑回归以使用 confint()，这是一个示例:

library(rms)
data(mtcars)
dd <- datadist(mtcars)
options(datadist = "dd")
fit <- lrm(am ~ gear + mpg, data=mtcars)
confint(fit)

这给出了错误:

Error: $ operator is invalid for atomic vectors

而 traceback() 给出:

4: profile.glm(object, which = parm, alpha = (1 - level)/4, trace = trace)
3: profile(object, which = parm, alpha = (1 - level)/4, trace = trace)
2: confint.glm(fit)
1: confint(fit)

我想 lrm() 模型没有实现 confint。

我的问题

有没有方便的替代方式？是否为 rms 包创建了其他替代方案？

Answer 1

confint.default(fit)似乎工作。请注意，它正在构建 Wald 置信区间，而不是 confint.glm() 更准确的配置文件置信区间。产生...
class(fit) ; methods(class="lrm") ; methods(class="rms")不要建议任何明显的替代品......

您可以查看 bootcov()和 friend 们的引导置信区间(但我还没有让这些工作......)