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

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.9999  -0.7220  -0.1242   0.5472   2.6973  

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 exp(est) exp(lower) exp(upper)
cde  0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566 1.031938  0.9314723   1.143239
pnde 0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566 1.031938  0.9314723   1.143239
tnie 0.004488714 0.007229627 0.6208777 0.5346801 -0.009681094 0.01865852 1.004499  0.9903656   1.018834
tnde 0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566 1.031938  0.9314723   1.143239
pnie 0.004488714 0.007229627 0.6208777 0.5346801 -0.009681094 0.01865852 1.004499  0.9903656   1.018834
te   0.035927117 0.052663602 0.6822001 0.4951124 -0.067291646 0.13914588 1.036580  0.9349225   1.149292
pm   0.126912192 0.256351942 0.4950701 0.6205506 -0.375528382 0.62935277       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 exp(est) exp(lower) exp(upper)
cde  0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566 1.031938  0.9314723   1.143239
pnde 0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566 1.031938  0.9314723   1.143239
tnie 0.004550181 0.007261454 0.6266212 0.5309075 -0.009682007 0.01878237 1.004561  0.9903647   1.018960
tnde 0.031438403 0.052259767 0.6015795 0.5474541 -0.070988858 0.13386566 1.031938  0.9314723   1.143239
pnie 0.004550181 0.007261454 0.6266212 0.5309075 -0.009682007 0.01878237 1.004561  0.9903647   1.018960
te   0.035988584 0.052667506 0.6833166 0.4944068 -0.067237832 0.13921500 1.036644  0.9349728   1.149371
pm   0.128430337 0.257804051 0.4981704 0.6183640 -0.376856318 0.63371699       NA         NA         NA
