20/12/15 21:51:03.17 QEfzM0Mb.net
４人に拡張
sim4 <- function(init=1){
i=init # ジャンケン回数
# j4 : 4人でジャンケンしたときに勝者の人数
j4=sample(0:3,1,prob=c(39,12,18,12))
if(j4==1) return(i) # 勝者１人ならi(=1)を返す
while(j4==0){ # あいこが続けば繰り返す
i=i+1
j4=sample(0:3,1,prob=c(39,12,18,12))
}
if(j4==1) return(i) # 勝者1人ならジャンケン回数iを返す
if(j4==2) return(sim2(i+1)) # 勝者2人ならsim2
if(j4==3) return(sim3(i+1)) # 勝者3人ならsim3
}
i=replicate(1e7,sim4())
mean(i)
BEST::plotPost(i)