Ограниченный выход кубического сплайна в Rms-пакете после cph

Question

Я разрабатываю модель регрессии ЦОГ в R.
Модель, которую я сейчас использую, выглядит следующим образом:

fh <- cph(S ~ rcs(MPV,4) + rcs(age,3) + BMI + smoking + hyperten + gender +
    rcs(FVCPP,3) + TLcoPP, x=TRUE, y=TRUE, surv=TRUE, time.inc=2*52) 

Если я тогда хочу посмотреть на это с

print(fh, latex = TRUE) 

Я получаю 3 coefs/SE/Wald и т. д. для MPV (MVP, MVP' and MVP'') и 2 для age (age, age'),
Может кто-нибудь объяснить мне, что это за результаты? то есть я полагаю, что они имеют отношение к ограниченным кубическим сплайнам, которые я добавил.

Answer 1

Когда ты пишешь rcs(MPV,4)вы определяете количество узлов, используемых в сплайне; в этом случае 4. Аналогично rcs(age,3) определяет сплайн с 3 узлами. Из-за ограничений по идентификации вычитается 1 узел из каждого сплайна. Вы можете думать об этом как об определении перехвата для каждого сплайна. Так rcs(Age,3) является линейной комбинацией 2 нелинейных базисных функций и перехвата, в то время как rcs(MPV,4) является линейной комбинацией 3 нелинейных базисных функций и точки пересечения, т. е.

а также

В приведенных выше обозначениях то, что вы получаете из оператора печати, это коэффициенты регрессии а также с соответствующими стандартными ошибками, p-значениями и т. д. Перехват а также обычно устанавливаются в ноль, но они важны, потому что без них подпрограмма подбора модели не имеет представления о том, где на оси у ограничивать сплайны.

И последнее замечание: на самом деле вы можете быть более заинтересованы в summary(fh),