проблема anova.rms с терминами rcs ()

Question

У меня проблема с функцией anova в пакете 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

Answer 1

Вы должны иметь возможность использовать вывод anova (mod2) в качестве одного из способов оценки значимости, но лучшим ответом будет сравнение статистики -2 * log (правдоподобия).Функция anova.rms не предназначена для двух подгонок модели.Предполагается, что второй и последующие безымянные аргументы являются именами терминов в модели, а не соответствуют объектам.

(Обратите внимание, что при rcs терминах вы не увидите сумму отдельных терминов, равную полной моделиквадратные значения. Я спросил об этом Харрелла, и он говорит, что нужно проводить кросс-модельные сравнения «вручную».)

Answer 2

Это сравнение выполняется с использованием lrtest (согласно комментарию Миши).