### Mediator model

Call:
lm(formula = bili ~ trt + age + male + stage, data = data)

Residuals:
    Min      1Q  Median      3Q     Max 
-4.1854 -2.5979 -1.2638  0.5727 24.2813 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) -2.54577    1.76354  -1.444 0.150018    
trt          0.60723    0.54646   1.111 0.267464    
age          0.02419    0.02706   0.894 0.372025    
male        -0.61583    0.84467  -0.729 0.466584    
stage        1.26164    0.32279   3.909 0.000117 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 4.472 on 271 degrees of freedom
Multiple R-squared:  0.06893,	Adjusted R-squared:  0.05518 
F-statistic: 5.016 on 4 and 271 DF,  p-value: 0.0006483

### Outcome model

Call:
survival::survreg(formula = Surv(time, status) ~ trt + bili + 
    trt:bili + age + male + stage, data = data, dist = "weibull")
               Value Std. Error     z       p
(Intercept)  9.97519    0.45912 21.73 < 2e-16
trt          0.15705    0.15386  1.02   0.307
bili        -0.07941    0.03454 -2.30   0.021
age         -0.00564    0.00634 -0.89   0.374
male        -0.34136    0.16281 -2.10   0.036
stage       -0.45524    0.08900 -5.12 3.1e-07
trt:bili    -0.00866    0.01931 -0.45   0.654
Log(scale)  -0.40002    0.07290 -5.49 4.1e-08

Scale= 0.67 

Weibull distribution
Loglik(model)= -1144.2   Loglik(intercept only)= -1203.8
	Chisq= 119.08 on 6 degrees of freedom, p= 2.5e-23 
Number of Newton-Raphson Iterations: 6 
n= 276 

### Mediation analysis 
             est         se          Z         p      lower      upper
cde   0.17391646 0.16646529  1.0447612 0.2961334 -0.1523495 0.50018243
pnde  0.19515050 0.19458469  1.0029078 0.3159054 -0.1862285 0.57652949
tnie -0.07237770 0.06598648 -1.0968566 0.2727041 -0.2017088 0.05695342
tnde  0.18757866 0.18386189  1.0202150 0.3076265 -0.1727840 0.54794134
pnie -0.06480587 0.05936083 -1.0917277 0.2749528 -0.1811510 0.05153923
te    0.12277279 0.19904265  0.6168165 0.5373558 -0.2673436 0.51288922
pm   -0.64968243 1.28051585 -0.5073599 0.6119023 -3.1594474 1.86008252
      exp(est) exp(lower) exp(upper)
cde  1.1899562  0.8586881   1.649022
pnde 1.2154939  0.8300839   1.779851
tnie 0.9301795  0.8173329   1.058607
tnde 1.2063251  0.8413193   1.729689
pnie 0.9372494  0.8343094   1.052890
te   1.1306275  0.7654100   1.670110
pm          NA         NA         NA

Evaluated at:
avar: trt
 a1 (intervened value of avar) = 2.3
 a0 (reference value of avar)  = 1.1
mvar: bili
 m_cde (intervend value of mvar for cde) = 1.4
cvar: age male stage
 c_cond (covariate vector value) = 50 1 2

Note that effect estimates can vary over m_cde and c_cond values when interaction = TRUE.

### Re-evaluation at c_cond = cmean
             est         se          Z         p      lower      upper
cde   0.17391646 0.16646529  1.0447612 0.2961334 -0.1523495 0.50018243
pnde  0.17595734 0.16898757  1.0412443 0.2977622 -0.1552522 0.50716689
tnie -0.07237770 0.06598648 -1.0968566 0.2727041 -0.2017088 0.05695342
tnde  0.16838550 0.16104372  1.0455887 0.2957510 -0.1472544 0.48402541
pnie -0.06480587 0.05936083 -1.0917277 0.2749528 -0.1811510 0.05153923
te    0.10357964 0.17583167  0.5890841 0.5558049 -0.2410441 0.44820337
pm   -0.76285103 1.59473315 -0.4783565 0.6323965 -3.8884706 2.36276851
      exp(est) exp(lower) exp(upper)
cde  1.1899562  0.8586881   1.649022
pnde 1.1923872  0.8561992   1.660580
tnie 0.9301795  0.8173329   1.058607
tnde 1.1833927  0.8630744   1.622593
pnie 0.9372494  0.8343094   1.052890
te   1.1091341  0.7858070   1.565497
pm          NA         NA         NA