ira: Individual-level Random Assignment Designs

In PowerUpR: Power Analysis Tools for Multilevel Randomized Experiments

Description

For an individual-level random assignment design (simple randomized controlled trial) use mdes.ira() to calculate the minimum detectable effect size, power.ira() to calculate the statistical power, and mrss.ira() to calculate the minimum required sample size (number of subjects).

If fixed block effects exist (randomized block designs), that is, treatment effect varies from one block to another but it is not random, use mdes.bira2f1() to calculate the minimum detectable effect size, power.bira2f1() to calculate the statistical power, and mrss.bira2f1() to calculate the minimum required sample size (number of subjects per block).

If block effect is constant, that is, treatment effect does not change from one block to another, use mdes.bira2c1() to calculate the minimum detectable effect size, power.bira2c1() to calculate the statistical power, and mrss.bira2c1() to calculate the minimum required sample size (number of subjects per block).

For partially nested individual-level random assignment designs (simple randomized controlled trial with intervention clusters) use mdes.ira_pn() to calculate the minimum detectable effect size, power.ira_pn() to calculate the statistical power, and mrss.ira_pn() to calculate the minimum required sample size (number of subjects).

Usage

mdes.ira(power=.80, alpha=.05, two.tailed=TRUE,
         p=.50, g1=0, r21=0, n)

power.ira(es=.25, alpha=.05, two.tailed=TRUE,
          p=.50, g1=0, r21=0, n)

mrss.ira(es=.25, power=.80, alpha=.05, two.tailed=TRUE,
         n0=10, tol=.10,
         p=.50, g1=0, r21=0)

mdes.bira2f1(power=.80, alpha=.05, two.tailed=TRUE,
             p=.50, g1=0, r21=0, n, J)

power.bira2f1(es=.25, alpha=.05, two.tailed=TRUE,
              p=.50, g1=0, r21=0, n, J)

mrss.bira2f1(es=.25, power=.80, alpha=.05, two.tailed=TRUE,
             J, n0=10, tol=.10,
             p=.50, g1=0, r21=0)

mdes.bira2c1(power=.80, alpha=.05, two.tailed=TRUE,
             p=.50, g1=0, r21=0,
             n, J)

power.bira2c1(es=.25, alpha=.05, two.tailed=TRUE,
              p=.50, g1=0, r21=0,
              n, J)

mrss.bira2c1(es=.25, power=.80, alpha=.05, two.tailed=TRUE,
             J, n0=10, tol=.10,
             p=.50, g1=0, r21=0)

mdes.ira_pn(power=.80, alpha=.05, two.tailed=TRUE,
            df=NULL, ratio_tc_var=1,
            rho_ic=.20, p=.50, r21=0, n, ic_size=1)

power.ira_pn(es=.25,alpha=.05, two.tailed=TRUE,
             df=NULL, ratio_tc_var=1,
             rho_ic=.20, p=.50, r21=0, n, ic_size=1)

mrss.ira_pn(es=.25, power=.80, alpha=.05, two.tailed=TRUE,
            ratio_tc_var=1, z.test=FALSE,
            rho_ic=.20, p=.50, r21=0, ic_size=1, n0=500, tol=.10)

Arguments

power	statistical power (1-β).
es	effect size.
alpha	probability of type I error.
two.tailed	logical; TRUE for two-tailed hypothesis testing, FALSE for one-tailed hypothesis testing.
df	degrees of freedom.
ratio_tc_var	ratio of the treatment group variance to the control group variance (needed for Satterthwaite approximation to df). Ignored when df argument is specified.
rho_ic	proportion of variance in the outcome (for treatment group) that is between intervention clusters.
p	proportion of units randomly assigned to the treatment condition.
g1	number of covariates. g also works.
r21	proportion of variance in the outcome explained by covariates (and fixed blocks if exists). r2 also works.
n	sample size (number of subjects per block, if exist, calculated using simple average or harmonic mean). For partially nested design n is total number of subjects (in all intervention clusters and control group).
J	level 2 sample size (number of blocks).
ic_size	sample size in each intervention cluster on average.
n0	starting value for n.
tol	tolerance to end iterative process for finding n or J.
z.test	logical; TRUE for z-test.

Value

fun	function name.
parms	list of parameters used in power calculation.
df	degrees of freedom.
ncp	noncentrality parameter.
power	statistical power (1-β).
mdes	minimum detectable effect size.
n	sample size (number of subjects).

References

Dong, N., & Maynard, R. (2013). PowerUp!: A tool for calculating minimum detectable effect sizes and minimum required sample sizes for experimental and quasi-experimental design studies. Journal of Research on Educational Effectiveness, 6(1), 24-67. doi: 10.1080/19345747.2012.673143

Lohr, S., Schochet, P. Z., & Sanders, E. (2014). Partially Nested Randomized Controlled Trials in Education Research: A Guide to Design and Analysis. NCER 2014-2000. National Center for Education Research. https://ies.ed.gov/ncer/pubs/20142000/pdf/20142000.pdf

Examples

# cross-checks
mdes.ira(n=250)
power.ira(es=.356, n=250)
mrss.ira(es=.356)

# cross-checks
mdes.bira2c1(n=15, J=20)
power.bira2c1(es=.325, n=15, J=20)
mrss.bira2c1(es=.325, J=20)

# cross-checks
mdes.bira2f1(n=15, J=20)
power.bira2f1(es=.325, n=15, J=20)
mrss.bira2f1(es=.325, J=20)

# cross-checks
mdes.ira_pn(n=250)
power.ira_pn(es=.377, n=250)
mrss.ira_pn(es=.377, n=250)

PowerUpR documentation built on Oct. 25, 2021, 5:06 p.m.