2

rmsパッケージの分散分析機能に問題があります。

require(rms)
getHdata(prostate)

mod1<-cph(Surv(dtime,status!="Alive")~stage+rx+age+wt,data=prostate,x=T,y=T)
mod2<-cph(Surv(dtime,status!="Alive")~stage+rx+rcs(age,4)+wt,data=prostate,x=T,y=T)
anova(mod1)
anova(mod2)

-すべてが正常に機能しますが、非線形性が年齢に与える影響についてモデルを比較しようとすると、

anova(mod1,mod2)

私は得る

Error in anova.rms(mod1, mod2) : factor names not in design:  mod2

これは何を意味するのでしょうか？それを回避するために何ができますか？

// M

4

2 に答える 2

1

有意性を評価する 1 つの方法として anova(mod2) の出力を使用できるはずですが、最善の答えは -2*log(尤度) 統計を比較することです。関数 anova.rms は、2 つのモデル近似を行うようには設計されていません。2 番目以降の名前のない引数は、近似オブジェクトではなく、モデル内の項の名前であると見なされます。

(項を使用すると、個々の項の合計が完全なモデルのカイ 2 乗値と等しくならないことに注意してrcsください。これについて Harrell に尋ねたところ、モデル間の比較は「手作業で」行うようにとのことでした。)

于 2011-09-17T15:44:50.100 に答える