4

rms パッケージからのロジスティック回帰を confint() で動作させることができません。以下に例を示します。

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() モデルには confint が実装されていないと思います。

私の質問

way に代わる便利な方法はありますか? rms パッケージ用に作成された他の代替手段はありますか?

4

2 に答える 2

4

confint.default(fit)動作するようです。...を生成するより正確なプロファイル信頼区間ではなく、Wald信頼区間を構築していることに注意してくださいconfint.glm()...

class(fit); methods(class="lrm"); methods(class="rms")明らかな代替品を提案しないでください...

あなたはbootcov()ブートストラップの信頼区間を調べたり、友達を探したりするかもしれません(しかし、私はまだこれらを機能させていません...)

于 2012-12-25T22:15:51.743 に答える