1

lrmパッケージからを使用して、最初にロジスティック回帰を実行しようとしていますRMS。私のモデルは では正常に動作しますが、 では動作しglmませんlrm。

model1 <- lrm( Outcome30Days ~ ISS1 + ISS2 + as.factor(GCSgr)+ 
as.factor(Gender)*as.factor(agegr),data=sub2, x=T, y=T, se.fit=T)

ISS1とISS2が削除された場合 、モデルは実行されますが、これら 2 つの変数では実行されません。エラーメッセージ：

Unable to fit model using “lrm.fit”

lrmブートストラップを使用するパッケージvalidateは (明らかに) でのみ動作するため、で実行する必要がありますlrm。どんな助けでも大歓迎です。

4

1 に答える 1

1

lrm は、glm よりも独立変数間の相関に対する許容度が低くなります。モデルが glm で実行され、いくつかの変数を削除したときに lrm で実行される場合、これがおそらく問題です。幸いなことに、 tol 引数を使用して許容範囲を調整できます。デフォルトでは tol=1e-7 です。tol=1e-9 に変更してみてください。コードは次のようになります。

model1 <- lrm( Outcome30Days ~ ISS1 + ISS2 + as.factor(GCSgr)+ 
as.factor(Gender)*as.factor(agegr),data=sub2, x=T, y=T, se.fit=T, tol=1e-9)

ペナルティを変更すると対数尤度が変化し、結果に影響を与える可能性があるため、これはペナルティをいじるよりも優れています。

于 2015-11-20T01:44:38.260 に答える