Esta pregunta está relacionada con dos preguntas diferentes que hice anteriormente:

1)Reproducción de diagrama de matriz de frecuencia

2)Agregue límites de confianza del 95% a la trama acumulativa

Deseo reproducir esta trama en R:

He llegado hasta aquí, usando el código debajo del gráfico:

#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)

Mi pregunta es: ¿Cómo puedo reproducir la trama superior en una pasada, sin trazar varias muestras?

Gracias