вопрос связан с двумя разными вопросами, которые я задавал ранее:

1)Воспроизвести частотный матричный график

2)Добавить 95% доверительные пределы к совокупному графику

Я хочу воспроизвести этот сюжет в R:

Я получил это далеко, используя код под графикой:

#Set the number of bets and number of trials and % lines
numbet <- 36 
numtri <- 1000 
#Fill a matrix where the rows are the cumulative bets and the columns are the trials
xcum <- matrix(NA, nrow=numbet, ncol=numtri)
for (i in 1:numtri) {
x <- sample(c(0,1), numbet, prob=c(5/6,1/6), replace = TRUE)
xcum[,i] <- cumsum(x)/(1:numbet)
}
#Plot the trials as transparent lines so you can see the build up
matplot(xcum, type="l", xlab="Number of Trials", ylab="Relative Frequency", main="", col=rgb(0.01, 0.01, 0.01, 0.02), las=1)

Мой вопрос:Как я могу воспроизвести верхний график за один проход, не нанося несколько образцов?

Благодарю.