Я хочу сгенерировать время выживания из модели пропорциональных рисков Кокса, которая содержит зависящий от времени ковариат. Модель
где генерируется из бинома (1,0.5) и m i ( t ) = β 0 + β 1 X i + β 2 X i t .
Истинные значения параметров используются как
Для не зависящего от времени ковариата (то есть я сгенерирован следующим образом
#For time independent case
# h_0(t) = 1
gamma <- -1
u <- runif(n=100,min=0,max=1)
Xi <- rbinom(n=100,size=1,prob=0.5)
T <- -log(u)/exp(gamma*Xi)
Может ли кто-нибудь помочь мне сгенерировать данные о выживании с изменяющейся во времени ковариатой.
Ответы:
ОК из вашего кода R вы предполагаете экспоненциальное распределение (постоянная опасность) для вашей базовой опасности. Следовательно, ваши опасные функции:
Затем они дают нам функции выживания:
источник