### 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 + trt:bili_bin + age + 
    male + stage, family = poisson(link = "log"), data = data)

Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)  -3.063320   0.781858  -3.918 8.93e-05 ***
trt          -0.173563   0.356264  -0.487   0.6261    
bili_bin      0.596285   0.632082   0.943   0.3455    
age          -0.004541   0.009065  -0.501   0.6165    
male         -1.119448   0.457165  -2.449   0.0143 *  
stage         0.607366   0.138421   4.388 1.14e-05 ***
trt:bili_bin  0.073132   0.402330   0.182   0.8558    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for poisson family taken to be 1)

    Null deviance: 198.14  on 275  degrees of freedom
Residual deviance: 161.32  on 269  degrees of freedom
AIC: 335.32

Number of Fisher Scoring iterations: 6

### Mediation analysis 
             est         se          Z         p      lower      upper
cde  -0.08541412 0.36466561 -0.2342259 0.8148096 -0.8001456 0.62931735
pnde -0.14191048 0.20585321 -0.6893770 0.4905860 -0.5453753 0.26155440
tnie -0.03461486 0.05570722 -0.6213712 0.5343554 -0.1437990 0.07456928
tnde -0.14540466 0.20754994 -0.7005767 0.4835672 -0.5521951 0.26138575
pnie -0.03112068 0.04924827 -0.6319141 0.5274430 -0.1276455 0.06540416
te   -0.17652533 0.21379178 -0.8256881 0.4089810 -0.5955495 0.24249886
pm    0.18243073 0.32926771  0.5540499 0.5795447 -0.4629221 0.82778358
      exp(est) exp(lower) exp(upper)
cde  0.9181320  0.4492636   1.876329
pnde 0.8676989  0.5796242   1.298948
tnie 0.9659774  0.8660618   1.077420
tnde 0.8646723  0.5756848   1.298729
pnie 0.9693586  0.8801653   1.067590
te   0.8381775  0.5512596   1.274430
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 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.08541412 0.36466561 -0.2342259 0.8148096 -0.8001456 0.62931735
pnde -0.14877687 0.21060693 -0.7064196 0.4799272 -0.5615589 0.26400514
tnie -0.03790634 0.06047761 -0.6267830 0.5308015 -0.1564403 0.08062760
tnde -0.15284616 0.21656735 -0.7057673 0.4803328 -0.5773104 0.27161805
pnie -0.03383705 0.05295419 -0.6389873 0.5228312 -0.1376254 0.06995126
te   -0.18668321 0.22193166 -0.8411743 0.4002503 -0.6216613 0.24829485
pm    0.18823296 0.32554942  0.5782009 0.5631285 -0.4498322 0.82629810
      exp(est) exp(lower) exp(upper)
cde  0.9181320  0.4492636   1.876329
pnde 0.8617614  0.5703193   1.302135
tnie 0.9628031  0.8551826   1.083967
tnde 0.8582617  0.5614063   1.312086
pnie 0.9667290  0.8714251   1.072456
te   0.8297065  0.5370515   1.281838
pm          NA         NA         NA