PowerUpR (version 1.0.4)

bira2r1: Two-Level Blocked Individual-level Random Assignment Design, Treatment at Level 1

Description

Use mdes.bira2r1() to calculate the minimum detectable effect size, power.bira2r1() to calculate the statistical power, and mrss.bira2r1() to calculate the minimum required sample size.

Usage

mdes.bira2r1(power=.80, alpha=.05, two.tailed=TRUE,
             rho2,  omega2, p=.50, g2=0, r21=0, r2t2=0,
             n, J)
power.bira2r1(es=.25, alpha=.05, two.tailed=TRUE,
              rho2,  omega2, g2=0, p=.50, r21=0, r2t2=0,
              n, J)
mrss.bira2r1(es=.25, power=.80, alpha=.05, two.tailed=TRUE,
             n, J0=10, tol=.10,
             rho2, omega2, g2=0, p=.50, r21=0, r2t2=0)

Arguments

power

statistical power \((1-\beta)\).

es

effect size.

alpha

probability of type I error.

two.tailed

logical; TRUE for two-tailed hypothesis testing, FALSE for one-tailed hypothesis testing.

rho2

proportion of variance in the outcome between level 2 units (unconditional ICC2).

omega2

treatment effect heterogeneity as ratio of treatment effect variance among level 2 units to the residual variance at level 2.

p

average proportion of level 1 units randomly assigned to treatment within level 2 units.

g2

number of covariates at level 2.

r21

proportion of level 1 variance in the outcome explained by level 1 covariates.

r2t2

proportion of treatment effect variance among level 2 units explained by level 2 covariates.

n

harmonic mean of level 1 units across level 2 units (or simple average).

J

level 2 sample size.

J0

starting value for J.

tol

tolerance to end iterative process for finding J.

Value

fun

function name.

parms

list of parameters used in power calculation.

df

degrees of freedom.

ncp

noncentrality parameter.

power

statistical power \((1-\beta)\).

mdes

minimum detectable effect size.

J

number of level 2 units.

See Also

cosa.bird2r1

Examples

# NOT RUN {
# cross-checks
mdes.bira2r1(rho2=.17, omega2=.50, n=15, J=20)
power.bira2r1(es=.366, rho2=.17, omega2=.50, n=15, J=20)
mrss.bira2r1(es=.366, rho2=.17, omega2=.50, n=15)
# }