rcs（）条款的anova.rms问题

Question

我在rms软件包中的anova函数遇到问题：

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

这是什么意思？ 我该怎么做才能规避它？

//中号

Answer 1

您应该能够使用anova（mod2）的输出作为评估显着性的一种方法，但是最好的答案是比较-2 * log（likelihood）统计信息。 anova.rms函数并非旨在采用两种模型拟合。 假定第二个及后续未命名的参数是模型中术语的名称，而不是适合的对象。

（请注意，使用rcs项时，您将看不到各个项的总和等于完整模型卡方值。我已经询问过Harrell，他说要“手动”进行跨模型比较。）

Answer 2

使用lrtest进行此比较（根据Misha的评论）。