### 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: MASS::glm.nb(formula = platelet ~ trt + bili + trt:bili + age + male + stage, data = data, init.theta = 8.026508119, link = log) Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 5.951502 0.146535 40.615 < 2e-16 *** trt 0.029041 0.054966 0.528 0.59726 bili -0.004408 0.018344 -0.240 0.81009 age -0.003110 0.002174 -1.431 0.15254 male -0.071854 0.067831 -1.059 0.28946 stage -0.086204 0.026925 -3.202 0.00137 ** trt:bili 0.001582 0.010347 0.153 0.87846 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 (Dispersion parameter for Negative Binomial(8.0265) family taken to be 1) Null deviance: 301.23 on 275 degrees of freedom Residual deviance: 282.26 on 269 degrees of freedom AIC: 3280.4 Number of Fisher Scoring iterations: 1 Theta: 8.027 Std. Err.: 0.694 2 x log-likelihood: -3264.431 ### Mediation analysis est se Z p lower upper cde 0.0375082129 0.057204512 0.65568626 0.5120260 -0.07461057 0.149626997 pnde 0.0371379246 0.058205894 0.63804405 0.5234450 -0.07694353 0.151219381 tnie -0.0005602178 0.005717730 -0.09797906 0.9219489 -0.01176676 0.010646326 tnde 0.0385215759 0.054979681 0.70065114 0.4835208 -0.06923662 0.146279770 pnie -0.0019438691 0.006062439 -0.32064145 0.7484821 -0.01382603 0.009938292 te 0.0365777068 0.056360413 0.64899643 0.5163407 -0.07388667 0.147042086 pm -0.0156020123 0.154411406 -0.10104184 0.9195172 -0.31824281 0.287038783 exp(est) exp(lower) exp(upper) cde 1.0382205 0.9281048 1.161401 pnde 1.0378362 0.9259421 1.163252 tnie 0.9994399 0.9883022 1.010703 tnde 1.0392732 0.9331059 1.157520 pnie 0.9980580 0.9862691 1.009988 te 1.0372549 0.9287769 1.158403 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.0375082129 0.057204512 0.65568626 0.5120260 -0.07461057 0.149626997 pnde 0.0406452170 0.052692261 0.77136977 0.4404878 -0.06262972 0.143920150 tnie -0.0005602178 0.005717730 -0.09797906 0.9219489 -0.01176676 0.010646326 tnde 0.0420288683 0.053160311 0.79060615 0.4291739 -0.06216343 0.146221162 pnie -0.0019438691 0.006062439 -0.32064145 0.7484821 -0.01382603 0.009938292 te 0.0400849992 0.052626063 0.76169482 0.4462422 -0.06306019 0.143230187 pm -0.0142617213 0.147527623 -0.09667153 0.9229872 -0.30341055 0.274887106 exp(est) exp(lower) exp(upper) cde 1.0382205 0.9281048 1.161401 pnde 1.0414825 0.9392912 1.154792 tnie 0.9994399 0.9883022 1.010703 tnde 1.0429246 0.9397293 1.157452 pnie 0.9980580 0.9862691 1.009988 te 1.0408992 0.9388870 1.153995 pm NA NA NA