Я реализовал рекуррентную регрессионную модель Кокса с условным интервалом (aka PWP) в R. Однако вместо базовой функции риска (delta_0 (t)), действующей как константа, я хотел бы включить распределение вероятностей на основепродолжительность времени события (это то, что я написал сам), чтобы гарантировать, что базовая функция опасности взвешена по продолжительности, с которой человек подвергался определенному количеству событий.
Например.Когда я рассматриваю 5 повторяющихся событий, базовая функция опасности сбрасывается, и у меня будет некоторый риск, связанный с определенными ковариатами.например, в коде R:
coxph(Surv(tstart,tstop,status) ~ drug + cluster(id) + strata(enum), data=medicalDF)
В моей работе это коллекция медицинских событий.Чем дольше человек имеет повторяющиеся события, тем меньше риск, что они будут иметь новое событие (потому что они проходят лечение).Например, если у вас 10-я головная боль, вероятность того, что ваша головная боль будет продолжительной или очень сильной, намного меньше, чем у вас 5-ая головная боль.
Я потратил много времени на изучение документации R.Есть ли способ изменить базовую функцию базовой опасности?