EDR: Compute the empirical detection rate for Type I errors and...

Description Usage Arguments References See Also Examples

View source: R/summary_functions.R

Description

Computes the detection rate for determining empirical Type I error and power rates using information from p-values.

Usage

1
EDR(p, alpha = 0.05)

Arguments

p

a numeric vector or matrix/data.frame of p-values from the desired statistical estimator. If a matrix, each statistic must be organized by column, where the number of rows is equal to the number of replications

alpha

the nominal detection rate to be studied (typical values are .10, .05, and .01). Default is .05

References

Sigal, M. J., & Chalmers, R. P. (2016). Play it again: Teaching statistics with Monte Carlo simulation. Journal of Statistics Education, 24(3), 136-156. http://www.tandfonline.com/doi/full/10.1080/10691898.2016.1246953

See Also

ECR

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
rates <- numeric(100)
for(i in 1:100){
   dat <- rnorm(100)
   rates[i] <- t.test(dat)$p.value
}

EDR(rates)
EDR(rates, alpha = .01)

# multiple rates at once
rates <- cbind(runif(1000), runif(1000))
EDR(rates)

SimDesign documentation built on Aug. 15, 2017, 9:02 a.m.

Related to EDR in SimDesign...