График Взвешенная частота Матрицы

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

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)

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

Спасибо.