Если у меня есть модель пропорциональных рисков Кокса с предиктором события категориальной переменной (и другими ковариантами, такими как возраст и т. д.), скажем, например, категориальная переменная, которая меня интересует, представляет собой размер опухоли, который может быть 0-10, 10- 20, 20-30, например, и я вижу тенденцию к более высокому ЧСС смерти с увеличением размера опухоли, как я могу вычислить r и получить p?
Модель опасностей Кокса и тест на тренд
Ответы (1)
Принятие неполной спецификации сопряжено с некоторой опасностью, поскольку у вас может быть классификация размеров, которая испортит ее, например: 0-10, 10-20, 20-30, 30-80, 80-120, 120- 240. Если вы тщательно не построили свой фактор так, чтобы правильно расположить восходящие уровни, то, что я собираюсь предложить в качестве упрощения, не удастся.
survmdl <- coxph(Surv(tie,event) ~ as.numeric(fact), data=dat)
Вы получите «проверку тенденции», которая будет интерпретироваться как увеличение журнала (опасности) для каждой категории для увеличения размера, и это будет один коэффициент. Так что опубликуйте свои фактические уровни фактора, если вам нужен более продуманный ответ.
person
IRTFM
schedule
16.04.2013
Большое спасибо за оперативный ответ! Я разделил переменную (уровень сывороточного гормона) на квартили, а затем также попытался разделить ее на категории по 10 нг/мл (и в этом последнем случае я получил 7 уровней)...
- person Gloria Dalla Costa; 16.04.2013
Из вашего комментария мне неясно, было ли достаточно информации, и я отмечаю, что вы не приняли ответ, поэтому нужны ли дополнительные комментарии? Если это так, вам нужно отредактировать свой вопрос, чтобы включить гораздо больше деталей (и код R).
- person IRTFM; 16.04.2013
