У меня есть следующая модель:Как построить модель Кокса опасности со шлицами
coxph(Surv(fulength, mortality == 1) ~ pspline(predictor))
где fulength является продолжительность наблюдения (включая смертность), предсказатель является предиктором смертности.
Выход выше команды заключается в следующем:
coef se(coef) se2 Chisq DF p
pspline(predictor), line 0.174 0.0563 0.0562 9.52 1.00 0.002
pspline(predictor), nonl 4.74 3.09 0.200
Как я могу построить эту модель, так что я получаю хорошую извилистую линию с 95% доверительными полосами и отношением рисков по оси у? Что я стремлюсь что-то похожее на это:
HTTP : //stat.ethz.ch/R-manual/R-devel/library/graphics/html/curve.html http://www.r-bloggers.com/plotting-95-confidence-bands-in-r- 2/ – efrem
Я использую пакеты RMS/Hmisc Фрэнка Харрелла, которые, вероятно, способны доставлять что-то очень похожее на выход, хотя я не знаю о правильном сюжете. Я статистически оскорблен составлением нормального распределения под результатами мужчин и женщин. Я не знаю, поддерживает ли rms psplines, потому что Франк предпочитает ограниченные кубические сплайны, но если вы разместите некоторые данные, я с удовольствием попробую. –
Спасибо BondedDust. как вы установили эти пакеты? Когда я пытаюсь установить, я получаю сообщения об ошибках, такие как эта установка пакета «TH.data» имела ненулевой статус выхода – Oposum