r - 计算生存回归模型的风险比的置信区间

Question

我试图获得survreg()功能风险比的置信区间，但它不像 coxph 那样简单。

Call:
survreg(formula = Surv(survival, DIED) ~ AGE + GENDER + PLATE + 
NEUTRO + NIH, data = IMRAWandIST, dist = "exponential")
           Value Std. Error     z         p
PLATE        0.00236   0.000367  6.42  1.39e-10
EUTRO      -0.02726   0.016695 -1.63  1.03e-01

比例固定为 1

指数分布 Loglik(model)= -4628.6 Loglik(intercept only)= -4736.1 Chisq= 215 on 5 个自由度，p= 0 Newton-Raphson 迭代次数：5 n=917（由于缺失而删除了 28 个观测值）

#estimate of beta
a <- c(coefficients(summary(fitexp)))

print(coef <- (a * -1 * 1 / 1)

#estimate of HR

print(HR <- exp(coef))

结果没有CI，只有标准差，所以我的问题是如何将AFT系数的SE转换为PH系数，然后计算HR的CI。

我有点卡在这里。有人可以帮忙吗？