### 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:
MASS::glm.nb(formula = platelet ~ trt + bili_bin + age + male + 
    stage, data = data, init.theta = 8.132711228, link = log)

Coefficients:
             Estimate Std. Error z value Pr(>|z|)    
(Intercept)  5.971661   0.140575  42.480   <2e-16 ***
trt          0.026199   0.043550   0.602   0.5475    
bili_bin    -0.089038   0.045520  -1.956   0.0505 .  
age         -0.003495   0.002159  -1.618   0.1056    
male        -0.042146   0.068578  -0.615   0.5388    
stage       -0.074391   0.026797  -2.776   0.0055 ** 
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for Negative Binomial(8.1327) family taken to be 1)

    Null deviance: 305.07  on 275  degrees of freedom
Residual deviance: 282.21  on 270  degrees of freedom
AIC: 3274.8

Number of Fisher Scoring iterations: 1


              Theta:  8.133 
          Std. Err.:  0.703 

 2 x log-likelihood:  -3260.806 
### Mediation analysis 
             est          se         Z         p        lower      upper
cde  0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566
pnde 0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566
tnie 0.004488714 0.007229627 0.6208777 0.5346801 -0.009681094 0.01865852
tnde 0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566
pnie 0.004488714 0.007229627 0.6208777 0.5346801 -0.009681094 0.01865852
te   0.035927117 0.052663602 0.6822001 0.4951124 -0.067291646 0.13914588
pm   0.126912192 0.256351942 0.4950701 0.6205506 -0.375528382 0.62935277
     exp(est) exp(lower) exp(upper)
cde  1.031938  0.9314723   1.143239
pnde 1.031938  0.9314723   1.143239
tnie 1.004499  0.9903656   1.018834
tnde 1.031938  0.9314723   1.143239
pnie 1.004499  0.9903656   1.018834
te   1.036580  0.9349225   1.149292
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 do not vary over m_cde and c_cond values when interaction = FALSE.

### Re-evaluation at c_cond = cmean
             est          se         Z         p        lower      upper
cde  0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566
pnde 0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566
tnie 0.004550181 0.007261454 0.6266212 0.5309075 -0.009682007 0.01878237
tnde 0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566
pnie 0.004550181 0.007261454 0.6266212 0.5309075 -0.009682007 0.01878237
te   0.035988584 0.052667506 0.6833166 0.4944068 -0.067237832 0.13921500
pm   0.128430337 0.257804051 0.4981704 0.6183640 -0.376856318 0.63371699
     exp(est) exp(lower) exp(upper)
cde  1.031938  0.9314723   1.143239
pnde 1.031938  0.9314723   1.143239
tnie 1.004561  0.9903647   1.018960
tnde 1.031938  0.9314723   1.143239
pnie 1.004561  0.9903647   1.018960
te   1.036644  0.9349728   1.149371
pm         NA         NA         NA