### 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.0116  -1.0649  -0.5783   1.0511   1.7122  

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 0.00000503 ***
---
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 = 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.81480964346477 -0.8001456 0.62931735 0.9181320  0.4492636   1.876329
pnde -0.14191048 0.20585321 -0.6893770 0.49058602767169 -0.5453753 0.26155440 0.8676989  0.5796242   1.298948
tnie -0.03461486 0.05570722 -0.6213712 0.53435543702666 -0.1437990 0.07456928 0.9659774  0.8660618   1.077420
tnde -0.14540466 0.20754994 -0.7005767 0.48356720628084 -0.5521951 0.26138575 0.8646723  0.5756848   1.298729
pnie -0.03112068 0.04924827 -0.6319141 0.52744303630556 -0.1276455 0.06540416 0.9693586  0.8801653   1.067590
te   -0.17652533 0.21379178 -0.8256881 0.40898103742048 -0.5955495 0.24249886 0.8381775  0.5512596   1.274430
pm    0.18243073 0.03357056  5.4342479 0.00000005502803  0.1166336 0.24822781        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.81480964346477 -0.8001456 0.62931735 0.9181320  0.4492636   1.876329
pnde -0.14877687 0.21060693 -0.7064196 0.47992722511692 -0.5615589 0.26400514 0.8617614  0.5703193   1.302135
tnie -0.03790634 0.06047761 -0.6267830 0.53080146847813 -0.1564403 0.08062760 0.9628031  0.8551826   1.083967
tnde -0.15284616 0.21656735 -0.7057673 0.48033284669676 -0.5773104 0.27161805 0.8582617  0.5614063   1.312086
pnie -0.03383705 0.05295419 -0.6389873 0.52283122376196 -0.1376254 0.06995126 0.9667290  0.8714251   1.072456
te   -0.18668321 0.22193166 -0.8411743 0.40025028833116 -0.6216613 0.24829485 0.8297065  0.5370515   1.281838
pm    0.18823296 0.03457526  5.4441517 0.00000005205285  0.1204667 0.25599922        NA         NA         NA
