symulacje

Monty Hall Paradox – symulacja gry

wynik_bz=0
wynik_zm=0
for (i in 1:1000)
{
nagroda=ceiling(runif(1)*3)
wybor1=ceiling(runif(1)*3)
# wariant bez zmiany
if (nagroda==wybor1) wynik_bz=wynik_bz+1
# wariant zmiany
if (wybor1!=nagroda) wynik_zm=wynik_zm+1
else
 {
 if (wybor1==1) {if (runif(1)<0.5) wybor2=2    else wybor2=3}
 if (wybor1==2) {if (runif(1)<0.5) wybor2=1    else wybor2=3}
 if (wybor1==3) {if (runif(1)<0.5) wybor2=1    else wybor2=2}
if (wybor2==nagroda) wynik_zm=wynik_zm+1
 }
}
# Oceny prawdopodobieństw wygranych
wynik_zm/1000
wynik_bz/1000