Frank Harrellの を使用して、を使用rms packageして予測モデルを構築しましたlrm function。
このモデルが別の (lrm-) モデルと比較して、二項イベントで有意に優れた予測値を持っているかどうかを比較したいと思います。
anova(model1, model2)または のpR2 functionようなさまざまな関数を使用pscl libraryして疑似 R^2 を比較しましたが、それらはすべて lrm ベースのモデルでは機能しません。
新しいモデルが以前のモデルより有意に優れているかどうかを確認するにはどうすればよいですか?
更新：サイズまたはステージ（他の変数に加えて）が最良のモデルを提供するかどうかを確認するための例（骨転移の可能性を予測したい場合）を次に示します。
library(rms)
getHdata(prostate)
ddd <- datadist(prostate)
options( datadist = "ddd" )
mod1 = lrm(as.factor(bm) ~ age + sz + rx, data=prostate, x=TRUE, y=TRUE)
mod2 = lrm(as.factor(bm) ~ age + stage + rx, data=prostate, x=TRUE, y=TRUE)