### 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.09095    1.02743  -1.062 0.288315    
trt         -0.22689    0.32051  -0.708 0.479009    
age         -0.02366    0.01586  -1.492 0.135744    
male         1.77461    0.47253   3.756 0.000173 ***
stage        0.66646    0.19644   3.393 0.000692 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 264.30  on 195  degrees of freedom
Residual deviance: 234.79  on 191  degrees of freedom
AIC: 244.79

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.13924532 0.20578471 -0.6766553 0.4986247 -0.5425759 0.26408530
tnie -0.04368292 0.06632251 -0.6586441 0.5101244 -0.1736727 0.08630680
tnde -0.14356613 0.20649001 -0.6952691 0.4868867 -0.5482791 0.26114685
pnie -0.03936212 0.05868296 -0.6707589 0.5023741 -0.1543786 0.07565437
te   -0.18292824 0.21580295 -0.8476633 0.3966255 -0.6058943 0.24003777
pm    0.22244537 0.37667763  0.5905457 0.5548248 -0.5158292 0.96071997
      exp(est) exp(lower) exp(upper)
cde  0.9181320  0.4492636   1.876329
pnde 0.8700146  0.5812491   1.302239
tnie 0.9572574  0.8405720   1.090141
tnde 0.8662635  0.5779435   1.298418
pnie 0.9614025  0.8569475   1.078590
te   0.8328279  0.5455863   1.271297
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.15476207 0.22007602 -0.7032210 0.4819181 -0.5861032 0.27657901
tnie -0.05132540 0.07678422 -0.6684368 0.5038548 -0.2018197 0.09916891
tnde -0.16058470 0.23343498 -0.6879205 0.4915029 -0.6181089 0.29693946
pnie -0.04550277 0.06636576 -0.6856362 0.4929425 -0.1755773 0.08457173
te   -0.20608747 0.23971029 -0.8597356 0.3899348 -0.6759110 0.26373607
pm    0.23011979 0.35727806  0.6440916 0.5195160 -0.4701323 0.93037193
      exp(est) exp(lower) exp(upper)
cde  0.9181320  0.4492636   1.876329
pnde 0.8566190  0.5564916   1.318611
tnie 0.9499695  0.8172423   1.104253
tnde 0.8516457  0.5389627   1.345734
pnie 0.9555170  0.8389726   1.088251
te   0.8137619  0.5086928   1.301785
pm          NA         NA         NA