r - 在 R 中拟合具有时变协变量的全参数比例风险模型

Question

我需要用时变协变量拟合参数 PH 模型（因此，不是 Cox 模型）。我们可以在 R 中做到这一点吗？我听说 survreg 函数无法处理随时间变化的协变量。我徒劳地寻找可以解决这个问题的软件包。

Answer 1

正如@adibender 所写，您可以轻松地估计一个具有恒定基线的模型，该poisson系列的对数时间offset为glm. 这是一个例子

> # Input parameters
> n <- 100         # Number of individuals
> t_max <- 5       # max number of period per individual
> beta <- c(-1, 1) # true coefficient
> 
> # Simulate data
> set.seed(47261114)
> sim_dat <- replicate(
+   n, 
+   {
+     out <- data.frame(
+       tstart = rep(NA_integer_, t_max), 
+       tstop  = rep(NA_integer_, t_max),
+       event  = rep(NA, t_max),
+       x      = rnorm(t_max))
+     
+     for(i in 1:t_max){
+       rate     <- exp(beta %*% c(1, out$x[i]))
+       tstop <- min(rexp(1, rate), 1)
+       out[i, ] <- list(i - 1, i - (1 - tstop),  tstop < 1, out$x[i])
+       if(out$event[i])
+         break
+     }
+     out[!is.na(out$tstart), ]
+   }, simplify = FALSE)
> 
> sim_dat <- do.call(rbind, sim_dat)
> head(sim_dat) # show final data
  tstart     tstop event          x
1      0 0.3018182  TRUE  0.7095841
2      0 0.6724803  TRUE  1.5152877
3      0 1.0000000 FALSE  0.1036868
4      1 2.0000000 FALSE -0.5214508
5      2 2.4831577  TRUE  1.0101403
6      0 1.0000000 FALSE  0.1437594
> 
> # Fit with glm
> glm(event ~ x + offset(log(tstop - tstart)), sim_dat, family = poisson())

Call:  glm(formula = event ~ x + offset(log(tstop - tstart)), family = poisson(), 
    data = sim_dat)

Coefficients:
(Intercept)            x  
    -0.9053       0.9714  

Degrees of Freedom: 248 Total (i.e. Null);  247 Residual
Null Deviance:      382.5 
Residual Deviance: 306.4    AIC: 498.4

flexsurv对于其他发行版，该软件包似乎可以。看到这个帖子。