### 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 + age + male + stage, data = data, init.theta = 8.025871972, link = log) Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 5.945393 0.141829 41.919 < 2e-16 *** trt 0.034108 0.043892 0.777 0.43711 bili -0.001686 0.004871 -0.346 0.72929 age -0.003119 0.002172 -1.436 0.15098 male -0.072023 0.067830 -1.062 0.28832 stage -0.086824 0.026574 -3.267 0.00109 ** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 (Dispersion parameter for Negative Binomial(8.0259) family taken to be 1) Null deviance: 301.20 on 275 degrees of freedom Residual deviance: 282.26 on 270 degrees of freedom AIC: 3278.5 Number of Fisher Scoring iterations: 1 Theta: 8.026 Std. Err.: 0.694 2 x log-likelihood: -3264.454 ### Mediation analysis est se Z p lower upper cde 0.040929835 0.05267065 0.7770900 0.4371057 -0.062302736 0.144162406 pnde 0.040929835 0.05267065 0.7770900 0.4371057 -0.062302736 0.144162406 tnie -0.001228225 0.00371725 -0.3304123 0.7410885 -0.008513901 0.006057451 tnde 0.040929835 0.05267065 0.7770900 0.4371057 -0.062302736 0.144162406 pnie -0.001228225 0.00371725 -0.3304123 0.7410885 -0.008513901 0.006057451 te 0.039701610 0.05256235 0.7553242 0.4500545 -0.063318695 0.142721915 pm -0.031573964 0.10426175 -0.3028336 0.7620166 -0.235923233 0.172775306 exp(est) exp(lower) exp(upper) cde 1.0417790 0.9395984 1.155072 pnde 1.0417790 0.9395984 1.155072 tnie 0.9987725 0.9915222 1.006076 tnde 1.0417790 0.9395984 1.155072 pnie 0.9987725 0.9915222 1.006076 te 1.0405003 0.9386443 1.153409 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 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.040929835 0.05267065 0.7770900 0.4371057 -0.062302736 0.144162406 pnde 0.040929835 0.05267065 0.7770900 0.4371057 -0.062302736 0.144162406 tnie -0.001228225 0.00371725 -0.3304123 0.7410885 -0.008513901 0.006057451 tnde 0.040929835 0.05267065 0.7770900 0.4371057 -0.062302736 0.144162406 pnie -0.001228225 0.00371725 -0.3304123 0.7410885 -0.008513901 0.006057451 te 0.039701610 0.05256235 0.7553242 0.4500545 -0.063318695 0.142721915 pm -0.031573964 0.10426175 -0.3028336 0.7620166 -0.235923233 0.172775306 exp(est) exp(lower) exp(upper) cde 1.0417790 0.9395984 1.155072 pnde 1.0417790 0.9395984 1.155072 tnie 0.9987725 0.9915222 1.006076 tnde 1.0417790 0.9395984 1.155072 pnie 0.9987725 0.9915222 1.006076 te 1.0405003 0.9386443 1.153409 pm NA NA NA