Как сравнить коэффициенты риска различных уровней в регрессии Кокса с выходом R?
Я запустил модель регрессии Кокса с категориальным предиктором, и я не уверен, могу ли я сравнить разные уровни друг с другом или просто с эталонной категорией?
Например, в наборе данных легких:
lung$ph.ecog=factor(lung$ph.ecog)
fit1=coxph(Surv(lung$time, lung$status) ~ ph.ecog, data=lung)
Я знаю, что вы никогда не установите ph.ecog как фактор, это просто демонстративно. Выход этой модели (ссылка 0):
coef exp(coef) se(coef) z Pr(>|z|)
ph.ecog1 0.3688 1.4461 0.1987 1.857 0.0634 .
ph.ecog2 0.9164 2.5002 0.2245 4.081 4.48e-05 ***
ph.ecog3 2.2080 9.0973 1.0258 2.152 0.0314 *
Итак, я знаю, что ph.ecog=1
имеет коэффициент опасности 1.44
это означает, что это в 1,4 раза чаще, чем ph.ecog=0
и уровень 2 в 2,5 раза чаще, чем 0
, но могу ли я сказать, что ph.ecog
2 является 1.728926 (= 2.5002 / 1.4461) более вероятным, чем ph.ecog=1?
Верно ли это и для других регрессионных моделей, таких как отношения шансов?
Большое спасибо