Lanzamiento de una moneda (distribución geométrica)

Se lanza una moneda equilibrada hasta que se obtiene cara (éxito). Se trata de simular un gran número de veces esta experiencia y calcular, por simulación y teóricamente (mediante la distribución geométrica), la probabilidad de obtener el primer éxito en una determinada repetición. El script pide en primer lugar en qué repetición se desea el primer éxito y, después, cuántas simulaciones se van a realizar.

NOTA: El siguiente script se debe ejecutar paso a paso.

#INICIO -------------------------
rm(list=ls(all=TRUE))

n<-as.numeric(readline("Número de lanzamiento en el que aparece por primera vez éxito: "))
n2<-as.numeric(readline("Introduce el número total de repeticiones del experimento: "))
moneda<-0:1 #éxito = 1
cont<-0 #contador para ver las veces que se tiene éxito en la tirada n

for (i in 1:n2)
{
        X<-0 #v.a "número de lanzamiento en el que aparece por primera vez éxito"
        while(sample(moneda,1)!=1){
                X<-X+1
        }
        if(X+1==n){
                cont<-cont+1
        }
}

cat("Resultado de la simulación: ",cont/n2,"\n")
cat("Resultado teórico:",dgeom(n-1,1/2),"\n")
#FIN -------------------------

No hay comentarios:

Publicar un comentario