### Mediator model

Call:
glm(formula = bili_bin ~ trt + age + male + stage, family = binomial(link = "logit"), 
    data = data)

Coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept) -1.53024    0.85116  -1.798  0.07220 .  
trt         -0.17117    0.25982  -0.659  0.51003    
age         -0.01386    0.01299  -1.067  0.28610    
male         1.33046    0.43911   3.030  0.00245 ** 
stage        0.74640    0.16356   4.563 5.03e-06 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 382.49  on 275  degrees of freedom
Residual deviance: 349.60  on 271  degrees of freedom
AIC: 359.6

Number of Fisher Scoring iterations: 4

### Outcome model

Call:
survival::survreg(formula = Surv(time, status) ~ trt + bili_bin + 
    age + male + stage, data = data, dist = "weibull")
                Value Std. Error     z       p
(Intercept) 10.577076   0.476246 22.21 < 2e-16
trt          0.108623   0.128565  0.84   0.398
bili_bin    -1.097700   0.159144 -6.90 5.3e-12
age         -0.013558   0.006467 -2.10   0.036
male         0.000915   0.172044  0.01   0.996
stage       -0.408418   0.093846 -4.35 1.3e-05
Log(scale)  -0.340606   0.071754 -4.75 2.1e-06

Scale= 0.711 

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

### Mediation analysis 
            est         se         Z         p      lower     upper exp(est)
cde  0.13034747 0.15427801 0.8448869 0.3981740 -0.1720319 0.4327268 1.139224
pnde 0.13034747 0.15427801 0.8448869 0.3981740 -0.1720319 0.4327268 1.139224
tnie 0.05422319 0.08206429 0.6607404 0.5087788 -0.1066199 0.2150662 1.055720
tnde 0.13034747 0.15427801 0.8448869 0.3981740 -0.1720319 0.4327268 1.139224
pnie 0.05422319 0.08206429 0.6607404 0.5087788 -0.1066199 0.2150662 1.055720
te   0.18457066 0.17449117 1.0577651 0.2901626 -0.1574257 0.5265671 1.202702
pm   0.31315831 0.41135947 0.7612765 0.4464919 -0.4930914 1.1194080       NA
     exp(lower) exp(upper)
cde   0.8419523   1.541455
pnde  0.8419523   1.541455
tnie  0.8988673   1.239944
tnde  0.8419523   1.541455
pnie  0.8988673   1.239944
te    0.8543403   1.693110
pm           NA         NA

Evaluated at:
avar: trt
 a1 (intervened value of avar) = 2.3
 a0 (reference value of avar)  = 1.1
mvar: bili_bin
 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 do not vary over m_cde and c_cond values when interaction = FALSE.

### Re-evaluation at c_cond = cmean
            est         se         Z         p       lower     upper exp(est)
cde  0.13034747 0.15427801 0.8448869 0.3981740 -0.17203187 0.4327268 1.139224
pnde 0.13034747 0.15427801 0.8448869 0.3981740 -0.17203187 0.4327268 1.139224
tnie 0.05021844 0.07568389 0.6635288 0.5069920 -0.09811926 0.1985561 1.051501
tnde 0.13034747 0.15427801 0.8448869 0.3981740 -0.17203187 0.4327268 1.139224
pnie 0.05021844 0.07568389 0.6635288 0.5069920 -0.09811926 0.1985561 1.051501
te   0.18056591 0.17162465 1.0520977 0.2927547 -0.15581223 0.5169440 1.197895
pm   0.29647484 0.39721335 0.7463869 0.4554337 -0.48204903 1.0749987       NA
     exp(lower) exp(upper)
cde   0.8419523   1.541455
pnde  0.8419523   1.541455
tnie  0.9065408   1.219640
tnde  0.8419523   1.541455
pnie  0.9065408   1.219640
te    0.8557199   1.676895
pm           NA         NA