Buròcrates, Administradors (Semantic MediaWiki), Curadors (Semantic MediaWiki), Editors (Semantic MediaWiki), Administradors
10.883
modificacions
→El ejemplo de la función ''montehall''
| Línia 756: | Línia 756: | ||
Esta función reproduce la paradoja de [http://math.ucsd.edu/~anistat/chi-an/MonteHallParadox.html Monte Hall]. Se trata de un concursante que está en un concurso donde hay tres puertas. Detrás de una de ellas hay un coche y detrás de las otras dos una cabra (una en cada una de las otras dos puertas). Se supone que al concursante le hace más ilusión ganar el coche que la cabra. El presentador le hace escoger una puerta, pero antes de abrirla lo que hace es abrir una de las otras dos puertas (teniendo en cuenta que el presentador sabe dónde está el coche y abre siempre una puerta con una cabra). Luego le pregunta al concursante qué quiere hacer, si quedarse en la puerta que ha escogido o cambiar a la otra puerta que queda cerrada. | Esta función reproduce la paradoja de [http://math.ucsd.edu/~anistat/chi-an/MonteHallParadox.html Monte Hall]. Se trata de un concursante que está en un concurso donde hay tres puertas. Detrás de una de ellas hay un coche y detrás de las otras dos una cabra (una en cada una de las otras dos puertas). Se supone que al concursante le hace más ilusión ganar el coche que la cabra. El presentador le hace escoger una puerta, pero antes de abrirla lo que hace es abrir una de las otras dos puertas (teniendo en cuenta que el presentador sabe dónde está el coche y abre siempre una puerta con una cabra). Luego le pregunta al concursante qué quiere hacer, si quedarse en la puerta que ha escogido o cambiar a la otra puerta que queda cerrada. | ||
Intuitivamente la gente suele decir que es indiferente pues hay un 50% de probabilidad de acertar pues quedan dos puertas. Pero se demuestra que si se cambia de puerta, se gana en el 66,6% de los casos. El siguiente código lo demuestra, y además ilustra el uso de loops (función '''for''') en R: | Intuitivamente la gente suele decir que es indiferente pues hay un 50% de probabilidad de acertar pues quedan dos puertas. Pero se demuestra que si se cambia de puerta, se gana en el 66,6% de los casos. El siguiente código lo demuestra, y además ilustra el uso de loops (función '''for''') y funciones de generación de número aleatorios en R: | ||
montehall<-function(repeticions,canvi) | |||
{ | |||
cotxe<-trunc(runif(repeticions,1,3.99999999999999999999999999999999)) | |||
eleccio<-trunc(runif(repeticions,1,3.99999999999999999999999999999999)) | |||
res<-0 | |||
for(i in 1:repeticions) | |||
{ | |||
if(cotxe[i]==eleccio[i]) | |||
{ | |||
if(cotxe[i]==1) | |||
{ | |||
obrir<-trunc(runif(1,2,3.99999999999999999999999999999999)) | |||
} | |||
if(cotxe[i]==2) | |||
{ | |||
obrir<-trunc(runif(1,1,2.99999999999999999999999999999999)) | |||
if(obrir==2) | |||
{ | |||
obrir=3 | |||
} | |||
} | |||
if(cotxe[i]==3) | |||
{ | |||
obrir<-trunc(runif(1,1,2.99999999999999999999999999999999)) | |||
} | |||
} | |||
else | |||
{ | |||
trampa<-seq(1:3) | |||
obrir<-subset(trampa,trampa!=cotxe[i] & trampa!=eleccio[i]) | |||
} | |||
if(canvi==T) | |||
{ | |||
trampa2<-seq(1:3) | |||
eleccio2<-subset(trampa2,trampa2!=obrir & trampa2!=eleccio[i]) | |||
} | |||
else | |||
{ | |||
eleccio2<-eleccio[i] | |||
} | |||
if(cotxe[i]==eleccio2) | |||
{ | |||
res<-res+1 | |||
} | |||
} | |||
return(round(100*res/repeticions,digits=2)) | |||
} | |||
Para realizar el experimento 1000 veces, cambiando siempre de puerta: | |||
montehall (1000,T) | |||
El resultado debería aproximarse a 66,6. | |||
=== Documentación de funciones === | === Documentación de funciones === | ||