cosa (version 1.2.0)

bcrd3r2: Blocked Cluster-level Regression Discontinuity (Three-level Design, Discontinuity at Level 2)

Description

Use mdes.bcrd3r2() to calculate minimum detectable effect size, power.bcrd3r2() to calculate statistical power, and cosa.bcrd3r2() for constrained optimal sample allocation.

Usage

cosa.bcrd3r2(cn1 = 0, cn2 = 0, cn3 = 0, cost = NULL,
             n1 = NULL, n2 = NULL, n3 = NULL, p = NULL,
             n0 = c(10, 3, 100 + g3), p0 = .499, constrain = "power",
             round = TRUE, max.power = FALSE,
             local.solver = c("LBFGS", "SLSQP", "MMA", "COBYLA"),
             rhots = NULL, k1 = -6, k2 = 6, dists = "normal",
             power = .80, es = .25, alpha = .05, two.tailed = TRUE,
             rho2, rho3, omega3, g3 = 0, r21 = 0, r22 = 0, r2t3 = 0)
mdes.bcrd3r2(power = .80, alpha = .05, two.tailed = TRUE,
             rhots = NULL, k1 = -6, k2 = 6, dists = "normal",
             rho2, rho3, omega3, r21 = 0, r22 = 0, r2t3 = 0,
             g3 = 0, p = .50, n1, n2, n3)
power.bcrd3r2(es = .25, alpha = .05, two.tailed = TRUE,
              rhots = NULL, k1 = -6, k2 = 6, dists = "normal",
              rho2, rho3, omega3, r21 = 0, r22 = 0, r2t3 = 0,
              g3 = 0, p = .50, n1, n2, n3)

Arguments

cn1

marginal cost per level 1 unit in treatment and control conditions.

cn2

marginal cost per level 2 unit in treatment and control conditions.

cn3

marginal cost per level 3 unit.

cost

total cost or budget.

n1

average number of level 1 units per level 2 unit.

n2

average number of level 2 units per level 3 unit.

n3

number of level 3 units.

p

proportion of level 2 units in treatment condition.

n0

vector of starting values for n1, n2, n3 (positional). Starting values are replaced with averages when sample sizes are constrained by bounds.

p0

starting value for p when rhots = 0 and p = NULL. Starting value is replaced with average when p is constrained by bounds.

power

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

es

effect size (Cohen's d).

alpha

probability of type I error (\(\alpha\)).

two.tailed

logical; TRUE for two-tailed hypothesis testing.

constrain

character; "cost", "power", or "mdes".

round

logical; TRUE for rounded COSA solution.

max.power

logical; TRUE for maximizing power instead of minimizing variance.

local.solver

subset of c("LBFGS", "SLSQP", "MMA", "COBYLA")

rhots

correlation between the treatment and the scoring variable. Specify rhots = 0 to obtain results equivalent to random assignment designs.

k1

left truncation point (in standard deviation units from full normal distribution mean), ignored when rhots is not NULL or dists = "uniform".

k2

right truncation point (in standard deviation units from full normal distribution mean), ignored when rhots is not NULL or dists = "uniform".

dists

distribution of the scoring variable; "normal" or "uniform". By default, dists = "normal" specification implies a truncated normal distribution with k1 = -6 and k2 = 6.

rho2

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

rho3

proportion of variance in the outcome between level 3 units (unconditional ICC3).

omega3

ratio of the treatment effect variance between level 3 units to the variance in the outcome between level 3 units.

g3

number of covariates at level 3.

r21

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

r22

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

r2t3

proportion of treatment effect variance between level 3 units explained by level 3 covariates.

Value

parms

list of parameters used in the function.

df

degrees of freedom.

sse

standardized standard error.

cosa

constrained optimal sample allocation.

mdes

minimum detectable effect size and (1 - \(\alpha\))% confidence limits.

power

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

Examples

# NOT RUN {
# cost constrained - optimize n2 and n3
cosa.bcrd3r2(constrain = "cost", cost = 20000,
             cn1 = c(5, 3), cn2 = 10, cn3 = 30,
             es = .20, rho2 = .20, rho3 = .10, omega3 = .10,
             r21 = .20, r22 = .20, r2t3 = .20,
             g3 = 1, p = .50, n1 = 25)

# minimum detectable effect size
mdes.bcrd3r2(power = .80, rho2 = .20, rho3 = .10, omega3 = .10,
             r21 = .20, r22 = .30, r2t3 = .40,
             g3 = 1, p = .50, n1 = 13.39, n2 = 3, n3 = 277)

# statistical power
power.bcrd3r2(es = .20, rho2 = .20, rho3 = .10, omega3 = .10,
              r21 = .20, r22 = .30, r2t3 = .40,
              g3 = 1, p = .50, n1 = 13.39, n2 = 3, n3 = 277)
# }