Gelmaneando
Hoy, gelmaneo así:
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 |
bar <- function(n, reps = 1e4){ foo <- function(n){ x <- rnorm(n) tmp <- t.test(x) c(tmp$p.value, abs(mean(x))) } res <- replicate(reps, foo(n)) tmp <- t(res) tmp <- tmp[tmp[,1] < 0.05,] tmp[,2] } res <- lapply(c(3, 10, 20, 50, 100), bar) sapply(res, mean) #[1] 0.8662636 0.6583157 0.4934551 0.3240322 0.2337086 |
Resumo: Fabrico un montón de errores de tipo I. Recuérdese: error de tipo I implica artículo… más