ssMediation.Sobel {powerMediation} | R Documentation |
Sample size for testing mediation effectd (Sobel's test)
Description
Calculate sample size for testing mediation effect based on Sobel's test.
Usage
ssMediation.Sobel(power,
theta.1a,
lambda.a,
sigma.x,
sigma.m,
sigma.epsilon,
n.lower = 1,
n.upper = 1e+30,
alpha = 0.05,
verbose = TRUE)
Arguments
power |
power of the test. |
theta.1a |
regression coefficient for the predictor in the linear regression linking
the predictor |
lambda.a |
regression coefficient for the mediator in the linear regression linking
the predictor |
sigma.x |
standard deviation of the predictor. |
sigma.m |
standard deviation of the mediator. |
sigma.epsilon |
standard deviation of the random error term
in the linear regression linking
the predictor |
n.lower |
lower bound of the sample size. |
n.upper |
upper bound of the sample size. |
alpha |
type I error rate. |
verbose |
logical. |
Details
The sample size is for testing the null hypothesis
versus the alternative hypothesis
for the linear regressions:
Test statistic is based on Sobel's (1982) test:
where is the estimated standard deviation
of the estimate
using multivariate
delta method:
and
is the
variance
of the estimate
, and
is the variance
of the estimate
,
is the variance of the
mediator
.
From the linear regression , we have the
relationship
. Hence, we can simply
the variance
to
Value
n |
sample size. |
res.uniroot |
results of optimization to find the optimal sample size. |
Note
The test is a two-sided test. For one-sided tests, please double the
significance level. For example, you can set alpha=0.10
to obtain one-sided test at 5% significance level.
Author(s)
Weiliang Qiu stwxq@channing.harvard.edu
References
Sobel, M. E. Asymptotic confidence intervals for indirect effects in structural equation models. Sociological Methodology. 1982;13:290-312.
See Also
powerMediation.Sobel
,
testMediation.Sobel
Examples
ssMediation.Sobel(power=0.8, theta.1a=0.1701, lambda.a=0.1998,
sigma.x=0.57, sigma.m=0.61, sigma.epsilon=0.2,
alpha = 0.05, verbose = TRUE)