Spatial Survival 区域数据

问：

如何在 R 中模拟空间区域生存数据？请帮我举个例子。

以下代码模拟来自此模型的数据

num_studies <- 50
 num_patients <- 200
 tot_patients <- num_studies * num_patients
 cov <- data.frame(id = 1:tot_patients,
 study = rep(1:num_studies, each = num_patients),
 treat = rbinom(tot_patients, 1, 0.5))

 pop_treat_effect <- -0.5
 study_treat_effect <- pop_treat_effect + rnorm(num_studies, 0, 0.5)
pars <- data.frame(treat = rep(study_treat_effect, each = num_patients))
 set.seed(908070)
dat <- simsurv(dist = "gompertz", lambdas = 0.1, gammas = 0.05,
 x = cov, betas = pars)
 dat <- merge(cov, dat)

但我想从以下模型进行模拟