Una prueba de hipótesis bootstrap de una muestra
Se estudió otra rana juvenil, la Rana C, y quieres comprobar si la Rana B y la Rana C tienen fuerzas de impacto similares. Por desgracia, no dispones de las fuerzas de impacto de la Rana C, pero sabes que su media es 0.55 N. Como no tienes los datos originales, no puedes hacer una prueba por permutación ni evaluar la hipótesis de que las fuerzas de las Ranas B y C provienen de la misma distribución. Por tanto, probarás otra hipótesis, menos restrictiva: la fuerza media del golpe de la Rana B es igual a la de la Rana C.
Para configurar la prueba de hipótesis con bootstrap, tomarás la media como estadístico de prueba. Recuerda: tu objetivo es calcular la probabilidad de obtener una fuerza de impacto media menor o igual que la observada para la Rana B si fuera cierta la hipótesis de que la media real de las fuerzas de impacto de la Rana B es igual a la de la Rana C. Primero trasladas todos los datos de la Rana B de forma que su media sea 0.55 N. Esto implica sumar la media de la Rana C y restar la media de la Rana B a cada medida de la Rana B. Así, otras propiedades de la distribución de la Rana B, como la varianza, permanecen sin cambios.
Este ejercicio forma parte del curso
Pensamiento estadístico en Python (Parte 2)
Instrucciones del ejercicio
- Traslada las fuerzas de impacto de la Rana B de modo que su media sea 0.55 N.
- Usa tu función
draw_bs_reps()para tomar 10,000 réplicas bootstrap de la media de tus fuerzas trasladadas. - Calcula el valor p encontrando la fracción de tus réplicas bootstrap que son menores que la media de impacto observada de la Rana B. Ten en cuenta que la variable de interés aquí es
force_b. - Imprime tu valor p.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Make an array of translated impact forces: translated_force_b
translated_force_b = ____
# Take bootstrap replicates of Frog B's translated impact forces: bs_replicates
bs_replicates = draw_bs_reps(____, ____, 10000)
# Compute fraction of replicates that are less than the observed Frog B force: p
p = np.sum(____ <= np.mean(____)) / 10000
# Print the p-value
print('p = ', ____)