Модель Кокса и проверка тренда

Question

Если у меня есть модель пропорциональных рисков Кокса с предиктором события категориальной переменной (и другими ковариатами, такими как возраст и т. Д.), Скажем, например, интересующая меня категориальная переменная - это размер опухоли, который может быть 0-10, 10-20, 20-30, например, и я вижу тенденцию к увеличению ЧСС смерти с увеличением размера опухоли, как я могу вычислить в r и получить a p?

Answer 1

Существует некоторая опасность в принятии вашей неполной спецификации, так как у вас может быть классификация размеров, которая испортит это, например: 0-10, 10-20, 20-30, 30-80, 80-120, 120-240. Если вы не тщательно продумали свой фактор, чтобы правильно упорядочить возрастающие уровни, то, что я собираюсь предложить как ярлык, потерпит неудачу.

 survmdl <- coxph(Surv(tie,event) ~ as.numeric(fact), data=dat)

Вы получите "тест тренда", который будет интерпретироваться как увеличение в журнале для каждой категории (Опасность) для увеличения размера, и это будет один коэффициент. Так что опубликуйте свои фактические уровни факторов, если вы хотите получить более продуманный ответ.