У меня есть следующая модель Кокса PH:
(Время, Событие) ~ X + Y + Z
Я хотел бы получить оценки опасности скорости (я говорю о скоростях опасности НЕ опасности коэффициентов) даны конкретные значения X
, Y
, Z
. Я знаю, что пакет muhaz R может рассчитывать наблюдаемые уровни опасности, но мне интересна предсказанная модель.
Есть ли способ сделать это в R?
Ответы:
Функция в пакете выживания R, чтобы получить базовый уровень опасности
basehaz
.Затем вам нужно умножить его на различные значения чтобы получить конкретную степень опасности с учетом найденных вами коэффициентов.еβ
Простой пример может помочь:
копировании и вставке из итоговых результатов регрессии Кокса, , подразумевая , что женщины в данных имеют показатели опасности ниже , чем у мужчин.ехр( - 0,1929 ) = 0,8245 18 %
summary(reg_fit)
источник
Функция
basehaz
(из предыдущего ответа) обеспечивает совокупную опасность , а не функцию опасности (коэффициент). Я считаю, что этот вопрос был о функции опасности. Оценка функции опасности потребует указания типа сглаживания (например, при оценке плотности). ПакетMuhaz
R может сделать это для одного образца данных. Я не знаю функцию, которая может сделать это для базовой опасности в модели Кокса. Мне тоже это нужно. Я думаю, что мне нужно будет самостоятельно сгладить накопленную опасность.источник