### 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:
lm(formula = alk_phos ~ trt + bili + age + male + stage, data = data)

Residuals:
    Min      1Q  Median      3Q     Max 
-2147.6  -994.2  -638.2   -20.9 11930.7 

Coefficients:
            Estimate Std. Error t value Pr(>|t|)   
(Intercept)  2381.21     835.79   2.849  0.00472 **
trt           -98.06     258.58  -0.379  0.70483   
bili           67.86      28.68   2.366  0.01868 * 
age            -8.01      12.79  -0.626  0.53176   
male          204.49     399.18   0.512  0.60887   
stage         -29.38     156.63  -0.188  0.85136   
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Residual standard error: 2111 on 270 degrees of freedom
Multiple R-squared:  0.02191,	Adjusted R-squared:  0.003798 
F-statistic:  1.21 on 5 and 270 DF,  p-value: 0.3048

### Mediation analysis 
              est         se          Z         p       lower      upper
cde  -117.6684743 310.293779 -0.3792163 0.7045272 -725.833107 490.496158
pnde -117.6684743 310.293779 -0.3792163 0.7045272 -725.833107 490.496158
tnie   49.4484500  49.162356  1.0058194 0.3145025  -46.907997 145.804897
tnde -117.6684743 310.293779 -0.3792163 0.7045272 -725.833107 490.496158
pnie   49.4484500  49.162356  1.0058194 0.3145025  -46.907997 145.804897
te    -68.2200243 312.771078 -0.2181149 0.8273396 -681.240073 544.800025
pm     -0.7248378   3.489952 -0.2076928 0.8354689   -7.565018   6.115342

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 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  -117.6684743 310.293779 -0.3792163 0.7045272 -725.833107 490.496158
pnde -117.6684743 310.293779 -0.3792163 0.7045272 -725.833107 490.496158
tnie   49.4484500  49.162356  1.0058194 0.3145025  -46.907997 145.804897
tnde -117.6684743 310.293779 -0.3792163 0.7045272 -725.833107 490.496158
pnie   49.4484500  49.162356  1.0058194 0.3145025  -46.907997 145.804897
te    -68.2200243 312.771078 -0.2181149 0.8273396 -681.240073 544.800025
pm     -0.7248378   3.489952 -0.2076928 0.8354689   -7.565018   6.115342