とりあえずこちらの3つめのやり方でかいてみた。 library(MCMCpack) batter <- rbind(c(0.71,0.20,0.06,0.01,0.02), # これをいろいろ振ってみる c(0.71,0.20,0.06,0.01,0.02), c(0.71,0.20,0.06,0.01,0.02), c(0.71,0.20,0.06,0.01,0.02), c(0.71,0.20,0.06…