### Mediator model

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

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.0072  -0.9198  -0.6511   1.1623   1.8965  

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)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.2343  -0.6478  -0.4723   0.3789   1.7449  

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