友達からの問題。
を[tex:0
という問題。
一回の試行の期待値がだから2くらいじゃね?と思ってシミュレーションしてみる。
# 和が初めて1を超えるときを返す関数 f<- function(vec){ return(head(which(vec>1),1)) } pp<- seq(1,10000,by=100) x<- rep(0,length(pp)) for(p in pp){ a<- matrix(runif(12*p),nc=12) asum<- t(apply(a,1,cumsum)) trial<- apply(asum,1,f) x[which(pp==p)]<- sum(trial)/p } plot(pp,x,type="l")