anova.rms問題RCS（）項與有效值包方差分析功能

Question

I'm有問題：

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) 

-everything工作好了，但是當我嘗試將模型用於比較的影響年齡非線性

anova(mod1,mod2) 

我得到

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

這是什麼意思？我能做些什麼來規避它？

//中號

Answer 1

您應該能夠使用單因素方差分析（模2）的輸出作爲一個方法來評估的意義，但最好的答案會是比較-2 *日誌（可能性）的統計數據。 anova.rms函數的設計不適合採用兩種模型。假定第二個和後來的未命名參數是模型中的術語名稱，而不是適合的對象。 （請注意，與rcs條款，你不會看到個別條款的總和等於完整的模型卡方值。我已經問了Harrell關於這一點，他說，「手工」進行交叉模型比較）。

Answer 2

這種比較是使用lrtest（每米莎的評論）來完成。