恣意的なカテゴリー分けで有意差を捏造する（黒い統計家シリーズ） - r-statistics-fanの日記恣意的なカテゴリー分け 連続変数で有意にならなかったからカテゴリ変数にして有意にしてる論文を読んだという話を聞いた これは面白い。つか、よくアクセプトされたな。 黒い統計家になった前提で、どの程度、この技が使えるのかやってみよう。 ### 黒い統計家 nsim <- 1000 res <- matrix(0, nrow = nsim, ncol = 3) for (j in 1:nsim) { n <- 100 y <- rnorm(n, 10, 5) age <- rnorm(n, 50, 15) p.glm <- summary(glm(y ~ 1 + age, family = gaussian))$coefficients[2, 4] y2 <- y[order(age)] p.kuro <- numeric(n - 3) for (i in 2:(n - 2)) { p.kuro
