Showing the Law of Large Numbers - Book Example Continued - Jan 22, 2021

N=1000
prob=0
for (n in 1:N){
  freq=0
  for (i in 1:n){
    trial =sample(1:6, size = 6, replace = TRUE)
    x=trial==1:6
    if (any(x)) freq=freq+1
  }
  prob=append(prob,freq/n)
  
}
prob=prob[-1]  

plot(prob,xlab="n",ylab="freq/n", type="l",col="blue", main="Fraction of Experiments having at least one match")
abline(h=.665)