### Mediator model

Call:
lm(formula = bili ~ trt + age + male + stage, data = data)

Residuals:
    Min      1Q  Median      3Q     Max 
-4.1854 -2.5979 -1.2638  0.5727 24.2813 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)    
(Intercept) -2.54577    1.76354  -1.444 0.150018    
trt          0.60723    0.54646   1.111 0.267464    
age          0.02419    0.02706   0.894 0.372025    
male        -0.61583    0.84467  -0.729 0.466584    
stage        1.26164    0.32279   3.909 0.000117 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 4.472 on 271 degrees of freedom
Multiple R-squared:  0.06893,	Adjusted R-squared:  0.05518 
F-statistic: 5.016 on 4 and 271 DF,  p-value: 0.0006483

### Outcome model

Call:
glm(formula = spiders ~ trt + bili + trt:bili + age + male + 
    stage, family = poisson(link = "log"), data = data)

Coefficients:
             Estimate Std. Error z value Pr(>|z|)    
(Intercept) -2.872069   0.714776  -4.018 5.87e-05 ***
trt         -0.140866   0.236246  -0.596   0.5510    
bili         0.053733   0.042179   1.274   0.2027    
age         -0.008401   0.009348  -0.899   0.3688    
male        -0.930092   0.455616  -2.041   0.0412 *  
stage        0.661619   0.133304   4.963 6.93e-07 ***
trt:bili    -0.007025   0.023813  -0.295   0.7680    
---
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: 165.31  on 269  degrees of freedom
AIC: 339.31

Number of Fisher Scoring iterations: 5

### Mediation analysis 
             est         se          Z         p       lower      upper
cde  -0.18084137 0.25681081 -0.7041813 0.4813199 -0.68418132 0.32249857
pnde -0.18653375 0.24560138 -0.7594980 0.4475547 -0.66790360 0.29483610
tnie  0.02738003 0.02796929  0.9789320 0.3276136 -0.02743876 0.08219883
tnde -0.19267657 0.23431396 -0.8223009 0.4109057 -0.65192349 0.26657034
pnie  0.03352285 0.03300454  1.0157043 0.3097702 -0.03116486 0.09821056
te   -0.15915372 0.24016903 -0.6626738 0.5075395 -0.62987637 0.31156892
pm   -0.15655508 0.28115067 -0.5568369 0.5776388 -0.70760026 0.39449011
      exp(est) exp(lower) exp(upper)
cde  0.8345677  0.5045031   1.380573
pnde 0.8298306  0.5127824   1.342906
tnie 1.0277583  0.9729343   1.085672
tnde 0.8247487  0.5210426   1.305479
pnie 1.0340911  0.9693158   1.103195
te   0.8528652  0.5326577   1.365566
pm          NA         NA         NA

Evaluated at:
avar: trt
 a1 (intervened value of avar) = 2.3
 a0 (reference value of avar)  = 1.1
mvar: bili
 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.18084137 0.25681081 -0.7041813 0.4813199 -0.68418132 0.32249857
pnde -0.20210463 0.21952306 -0.9206533 0.3572315 -0.63236192 0.22815266
tnie  0.02738003 0.02796929  0.9789320 0.3276136 -0.02743876 0.08219883
tnde -0.20824745 0.21204972 -0.9820690 0.3260659 -0.62385727 0.20736237
pnie  0.03352285 0.03300454  1.0157043 0.3097702 -0.03116486 0.09821056
te   -0.17472460 0.21602781 -0.8088060 0.4186267 -0.59813133 0.24868213
pm   -0.14146685 0.23093505 -0.6125829 0.5401522 -0.59409123 0.31115752
      exp(est) exp(lower) exp(upper)
cde  0.8345677  0.5045031   1.380573
pnde 0.8170094  0.5313353   1.256277
tnie 1.0277583  0.9729343   1.085672
tnde 0.8120061  0.5358734   1.230428
pnie 1.0340911  0.9693158   1.103195
te   0.8396882  0.5498381   1.282334
pm          NA         NA         NA