Me estoy volviendo un poco tonto con las probabilidades. - página 7

Nuevo comentario
 
Mathemat:
No, no son equivalentes. Todavía tienes que agruparlos en grupos de cuatro para ver cuántos seises hay.


En centuriones, o legiones, lo entendería. ¿Por qué cuatro de cada uno?
 

Lo entiendo, soy un poco torpe. Hora de dormir :)

 
Porque el evento elemental (el resultado de la prueba) es "al menos un seis en cuatro tiradas". Sólo voy a lanzar un programa.
 
Mathemat:
Porque el evento elemental (el resultado de la prueba) es "caer al menos un seis en cuatro tiradas".

Sorpresa. No lo es.
 
tara:

Svetlana, discúlpame - hemos tenido una pequeña charla con el tocayo. ¿Qué haces esta mañana?

¿Quieres llevarme al cine? :)
 
Swetten:
¿Quieres invitarme al cine? :)

Sí, por una taza de té...
 

Jackass. :)

 

Gracias :)

 
tara:

Te voy a sorprender. No lo es.
¿Qué pasa con eso?
 

Mira, tocayo, esto es una simulación de un juego de matemáticas (4 dados), cien millones de partidas:

#property show_inputs

extern int MAX = 32768;
extern int SERIES = 100000000;

int start( )
{
   int st = GetTickCount( );
   MathSrand( GetTickCount( ) );   
   
   int success = 0;
   for( int i = 0; i < SERIES; i ++ )
   {
      /// В этом маленьком цикле моделируется одна игра (бросок 4 костей). Как только получаем шестерку, игру прекращаем и записываем ее результат как "успех".
      for( int j = 0 ; j < 4; j ++ )
         if( genUniform( ) == 6 )          { success ++ ;   break; }

      ///if( i % 1000000 == 0 )        Comment( i / 1000000 + " mln." );
   }
   
   Print( "success rate = " + ( success + 0.0 ) / SERIES );
   int gone = ( GetTickCount( ) - st ) / 1000.0 ;
   Print( "Total time = " + gone + " sec." );
   return( 0 );
}//+------------------------------------------------------------------+


      int genUniform( )
      {
         int rand = MathRand( );
         return( 1 + 6 * rand / MAX );
      }//+------------------------------------------------------------------+

Resultado:


La simulación de una distribución uniforme de 1 a 6 no es muy precisa, pero el error es pequeño, no más de 0,001.

La C.S. de la desviación de la frecuencia respecto a la probabilidad es MathSqrt( npq ) / n ~ 1/20000, por lo que tampoco en este caso tienes ninguna posibilidad de acercarte a p=2/3.

El valor exacto de la probabilidad (o... er... frecuencia m.o.) es 1 - (5/6)^4 ~ 0,517747.

1234567891011
Nuevo comentario