データの層化 交絡因子と曝露率に関係があったほうが層化の威力が分かりやすいと思ったので、e=c+0.25として計算してみた。 N <- 1000 e <- 0.5 b <- 0.05 c.max <- 0.5 alpha <- 0.1 n <- 5 myfunc <- function(){ deta <- matrix(nrow=N,ncol=3) deta[,1] …