Flexsurvreg-应用协变量

我正在使用

flexsurv

软件包中的

flexsurvreg

，以使Gompertz模型适合生存数据。每个参数上有一个分类和一个连续协变量的例子如下：

library(survival)
library(flexsurv)
Call:
flexsurvreg(formula = Surv(time, status) ~ sex + shape(ph.ecog), 
data = lung, dist = "gompertz")

Estimates: 
                data mean  est        L95%       U95%       se         exp(est)     L95%       U95%     
shape                  NA   0.000659  -0.000210   0.001527   0.000443         NA         NA         NA
rate                   NA   0.003355   0.002042   0.005512   0.000850         NA         NA         NA
sex              1.396476  -0.525441  -0.852719  -0.198163   0.166982   0.591294   0.426254   0.820236
shape(ph.ecog)   0.951542   0.000937   0.000375   0.001498   0.000287   1.000937   1.000375   1.001500

N = 227,  Events: 164,  Censored: 63
Total time at risk: 69522
Log-likelihood = -1139.932, df = 4
AIC = 2287.864


par$coefficients
shape          rate          sex            shape(ph.ecog) 
0.0006588492  -5.6973226550  -0.5254411364   0.0009368314 

请注意，摘要中的利率术语与

par$coefficients

调用之间的差异在于后者显示为前者的自然对数

根据渐晕图，Gompertz危险函数定义为——

Gompertz distribution with shape parameter a and rate parameter b has hazard function H(x: a, b) = b.e^{ax}

我的问题是，如何将协变量应用于形状和速率参数，以导出新的危害函数

新的形状参数是

shape\u new=shape+ph.ecog\u value（介于1-5之间）*shape（ph.ecog）=0.0006588492+（一些值介于1和5之间）*0.0009368314

，新的比率是

rate\u new=rate+sex=-5.6973226550-0.5254411364

（根据我前面示例中的变量）

注-之前的问题是相关的，但是，它不涉及形状/速率参数的协变量-