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

Я разрабатываю модель регрессии ЦОГ в 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').
Может кто-нибудь объяснить мне, что это за выходные данные?то есть я полагаю, что они имеют отношение к ограниченным кубическим сплайнам, которые я добавил.

1 Ответ

0 голосов
/ 02 марта 2019

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

и

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

В качестве заключительного замечания вы можете быть на самом делебольше интересует вывод summary(fh).

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...