### 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 = "exponential")
               Value Std. Error     z       p
(Intercept) 10.85310    0.65341 16.61 < 2e-16
trt          0.16870    0.23098  0.73   0.465
bili        -0.11558    0.05129 -2.25   0.024
age         -0.00885    0.00929 -0.95   0.340
male        -0.48473    0.24212 -2.00   0.045
stage       -0.59061    0.12635 -4.67 2.9e-06
trt:bili    -0.00046    0.02904 -0.02   0.987

Scale fixed at 1 

Exponential distribution
Loglik(model)= -1156.8   Loglik(intercept only)= -1206.3
	Chisq= 98.9 on 6 degrees of freedom, p= 4.3e-19 
Number of Newton-Raphson Iterations: 5 
n= 276 

### Mediation analysis 
             est         se          Z         p      lower      upper
cde   0.20166160 0.24935165  0.8087438 0.4186625 -0.2870587 0.69038186
pnde  0.20303585 0.30250303  0.6711862 0.5021019 -0.3898592 0.79593089
tnie -0.08499397 0.07810734 -1.0881688 0.2765206 -0.2380815 0.06809361
tnde  0.20263339 0.28524234  0.7103903 0.4774621 -0.3564313 0.76169811
pnie -0.08459151 0.07779868 -1.0873130 0.2768985 -0.2370741 0.06789110
te    0.11804188 0.30295462  0.3896355 0.6968061 -0.4757383 0.71182203
pm   -0.79674486 2.16591355 -0.3678563 0.7129804 -5.0418574 3.44836769
      exp(est) exp(lower) exp(upper)
cde  1.2234339  0.7504677   1.994477
pnde 1.2251164  0.6771522   2.216503
tnie 0.9185178  0.7881384   1.070466
tnde 1.2246234  0.7001706   2.141910
pnie 0.9188876  0.7889328   1.070249
te   1.1252912  0.6214261   2.037701
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.20166160 0.24935165  0.8087438 0.4186625 -0.2870587 0.69038186
pnde  0.20201570 0.26134450  0.7729862 0.4395305 -0.3102101 0.71424150
tnie -0.08499397 0.07810734 -1.0881688 0.2765206 -0.2380815 0.06809361
tnde  0.20161324 0.24784059  0.8134795 0.4159432 -0.2841454 0.68737186
pnie -0.08459151 0.07779868 -1.0873130 0.2768985 -0.2370741 0.06789110
te    0.11702173 0.26492462  0.4417171 0.6586940 -0.4022210 0.63626444
pm   -0.80328878 2.01019061 -0.3996083 0.6894451 -4.7431900 3.13661242
      exp(est) exp(lower) exp(upper)
cde  1.2234339  0.7504677   1.994477
pnde 1.2238672  0.7332929   2.042637
tnie 0.9185178  0.7881384   1.070466
tnde 1.2233748  0.7526572   1.988483
pnie 0.9188876  0.7889328   1.070249
te   1.1241439  0.6688329   1.889410
pm          NA         NA         NA