我目前正在使用R apckagesurvreg中的函数估计具有对数逻辑分布的生存模型(特别是加速故障时间模型)。survival我想模拟符合 King 等人的预期生存时间。(2001),但我不确定从survreg回归输出计算对数逻辑分布的预期生存时间所需的链接函数。我在下面包含了一个最小的工作示例:
library(survival)
data(kidney)
survreg(formula = Surv(time, status) ~
age +
cluster(id),
data = kidney,
dist = "loglogistic",
robust = TRUE)
## Value Std. Err (Naive SE) z p
## (Intercept) 4.38127 0.6783 0.5338 6.46 1.05e-10
## age -0.00298 0.0135 0.0114 -0.22 8.26e-01
## Log(scale) -0.23009 0.0732 0.1038 -3.14 1.67e-03
## Scale= 0.794
## Log logistic distribution
## Loglik(model)= -342 Loglik(intercept only)= -342
## Chisq= 0.07 on 1 degrees of freedom, p= 0.79
## (Loglikelihood assumes independent observations)
## Number of Newton-Raphson Iterations: 3
## n= 76
我只是想知道如何根据survreg输出的估计参数计算预期生存时间。
参考文献: King, G.、Tomz, M. 和 Wittenberg, J. (2000)。充分利用统计分析:改进解释和演示。美国政治学杂志,347-361。