ssMediation.VSMc.logistic: Sample size for testing mediation effect in logistic regression based on Vittinghoff, Sen and McCulloch's (2009) method

Description

Calculate sample size for testing mediation effect in logistic regression based on Vittinghoff, Sen and McCulloch's (2009) method.

Usage

ssMediation.VSMc.logistic(power, 
                          b2, 
                          sigma.m, 
                          p, 
                          corr.xm, 
                          n.lower = 1, 
                          n.upper = 1e+30, 
                          alpha = 0.05, 
                          verbose = TRUE)

Arguments

power

power for testing $b_{2} = 0$ for the logistic regression $\log (p_{i} / (1 - p_{i})) = b 0 + b 1 x_{i} + b 2 m_{i}$ .

regression coefficient for the mediator $m$ in the logistic regression $\log (p_{i} / (1 - p_{i})) = b 0 + b 1 x_{i} + b 2 m_{i}$ .

sigma.m

standard deviation of the mediator.

the marginal prevalence of the outcome.

corr.xm

correlation between the predictor $x$ and the mediator $m$ .

n.lower

lower bound for the sample size.

n.upper

upper bound for the sample size.

alpha

type I error rate.

verbose

logical. TRUE means printing sample size; FALSE means not printing sample size.

Value

sample size.

res.uniroot

results of optimization to find the optimal sample size.

Details

The test is for testing the null hypothesis $b_{2} = 0$ versus the alternative hypothesis $b_{2} \neq 0$ for the logistic regressions: $\log (p_{i} / (1 - p_{i})) = b_{0} + b_{1} x_{i} + b_{2} m_{i}$

Vittinghoff et al. (2009) showed that for the above logistic regression, testing the mediation effect is equivalent to testing the null hypothesis $H_{0} : b_{2} = 0$ versus the alternative hypothesis $H_{a} : b_{2} \neq 0$ .

The full model is $\log (p_{i} / (1 - p_{i})) = b_{0} + b_{1} x_{i} + b_{2} m_{i}$

The reduced model is $\log (p_{i} / (1 - p_{i})) = b_{0} + b_{1} x_{i}$

Vittinghoff et al. (2009) mentioned that if confounders need to be included in both the full and reduced models, the sample size/power calculation formula could be accommodated by redefining corr.xm as the multiple correlation of the mediator with the confounders as well as the predictor.

References

Vittinghoff, E. and Sen, S. and McCulloch, C.E.. Sample size calculations for evaluating mediation. Statistics In Medicine. 2009;28:541-557.

Examples

Run this code

# NOT RUN {
  # example in section 4 (page 545) of Vittinghoff et al. (2009).
  # n=255

  ssMediation.VSMc.logistic(power = 0.80, b2 = log(1.5), sigma.m = 1, p = 0.5, 
    corr.xm = 0.5, alpha = 0.05, verbose = TRUE)

# }

Run the code above in your browser using DataLab

Get 50% off unlimited learning