『どっちを選ぶ』のためのメモ
- セミナーのための道具
A <- function(){ r <- runif(1) if(r<pa){ plot(c(1),cex=40) }else{ plot(c(1),cex=60,pch=20) } } B <- function(){ r <- runif(1) if(r<pb){ plot(c(1),cex=40) }else{ plot(c(1),cex=60,pch=20) } } pas <- runif(10,0.2,0.8) pbs <- runif(10,0.2,0.8) s <- sample(1:10,1) pa <- pas[s] pb <- pbs[s] A() B()