▽ほくそ笑む
●11/09 21:24
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 <- rbinom(n, 1, prob) # ﾛｼﾞｽﾃｨｯｸ回帰ﾓﾃﾞﾙの適用 fit <- glm(y ~ X, family = binomial, control = list(maxit = 50))