おとなりアンテナ | おすすめページ
n <- 1300 p <- 300 # データの生成 set.seed(314) x <- rnorm(n * p, mean = 0, sd = sqrt(1/n)) X <- matrix(x, nrow = n, ncol = p) beta <- matrix(rep(c(10, -10, 0), each = p/3)) prob <- plogis(X %*% beta) y <-
