http://d.hatena.ne.jp/foo22222/20110710/1309865390
トリッキーなプログラミングで短くしてみた
x <- 1:10
s <- x*3
n <- x*10
ans <- array(unlist(mapply(binom.confint, s, n, prior.shape1=1, prior.shape2=1)), dim=c(11, 6, 10))[c(2, 3, 5, 9), 4:6,]
matplot(t(apply(ans, 3, "+")), type="l", col=c("black","red","blue","green"), lty=1)
※コメント投稿者のブログIDはブログ作成者のみに通知されます