### 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:
glm(formula = spiders ~ trt + bili_bin + age + male + stage, 
    family = binomial(link = "logit"), data = data)

Coefficients:
             Estimate Std. Error z value Pr(>|z|)    
(Intercept) -3.637453   1.040043  -3.497 0.000470 ***
trt         -0.213652   0.299018  -0.715 0.474910    
bili_bin     1.099999   0.306126   3.593 0.000327 ***
age         -0.006266   0.014950  -0.419 0.675139    
male        -1.629057   0.588212  -2.770 0.005614 ** 
stage        0.917829   0.215159   4.266 1.99e-05 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 332.32  on 275  degrees of freedom
Residual deviance: 278.77  on 270  degrees of freedom
AIC: 290.77

Number of Fisher Scoring iterations: 5

### Mediation analysis 
             est         se          Z         p      lower      upper
cde  -0.25638256 0.35882155 -0.7145127 0.4749102 -0.9596599 0.44689474
pnde -0.25638256 0.35882155 -0.7145127 0.4749102 -0.9596599 0.44689474
tnie -0.04659889 0.07309229 -0.6375350 0.5237764 -0.1898571 0.09665937
tnde -0.25638256 0.35882155 -0.7145127 0.4749102 -0.9596599 0.44689474
pnie -0.04659889 0.07309229 -0.6375350 0.5237764 -0.1898571 0.09665937
te   -0.30298146 0.36608377 -0.8276287 0.4078808 -1.0204925 0.41452956
pm    0.13479265 0.25743161  0.5236056 0.6005529 -0.3697640 0.63934933
      exp(est) exp(lower) exp(upper)
cde  0.7738459  0.3830231   1.563450
pnde 0.7738459  0.3830231   1.563450
tnie 0.9544702  0.8270773   1.101485
tnde 0.7738459  0.3830231   1.563450
pnie 0.9544702  0.8270773   1.101485
te   0.7386128  0.3604174   1.513658
pm          NA         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.2563826 0.35882155 -0.7145127 0.4749102 -0.9596599 0.4468947 0.7738459
pnde -0.2563826 0.35882155 -0.7145127 0.4749102 -0.9596599 0.4468947 0.7738459
tnie -0.0523180 0.08080341 -0.6474726 0.5173261 -0.2106898 0.1060538 0.9490270
tnde -0.2563826 0.35882155 -0.7145127 0.4749102 -0.9596599 0.4468947 0.7738459
pnie -0.0523180 0.08080341 -0.6474726 0.5173261 -0.2106898 0.1060538 0.9490270
te   -0.3087006 0.36767629 -0.8395988 0.4011334 -1.0293329 0.4119317 0.7344006
pm    0.1485140 0.27672463  0.5366851 0.5914851 -0.3938563 0.6908843        NA
     exp(lower) exp(upper)
cde   0.3830231   1.563450
pnde  0.3830231   1.563450
tnie  0.8100253   1.111882
tnde  0.3830231   1.563450
pnie  0.8100253   1.111882
te    0.3572452   1.509731
pm           NA         NA