限制使用rms lrm对象

Question

我无法从rms包中获得逻辑回归来处理logistic ()，下面是一个示例：

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

这就产生了错误：

错误：$操作符对原子向量无效

而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()模型实现的。

我的问题

有没有一种方便的替代方法？有没有为rms软件包创建的其他替代方案？

Answer 1

试试这个：

> summary(fit)[ , c("Lower 0.95", "Upper 0.95")]
               Lower 0.95   Upper 0.95
gear        -6.148918e+01 9.756505e+01
 Odds Ratio  1.975094e-27 2.354854e+42
mpg         -1.063706e+00 6.028564e+00
 Odds Ratio  3.451743e-01 4.151185e+02

(这些结果表明完全分离或其他模拟病理。)

Answer 2

confint.default(fit)似乎起作用了。注意，它正在构造Wald置信区间，而不是confint.glm()生成的更精确的剖面置信区间。

class(fit)；methods(class="lrm")；methods(class="rms")不建议任何明显的替换.

您可能会查看bootcov()和朋友的引导置信区间(但我还没有开始工作……)