ComenzarEmpieza gratis

Prueba de permutación con datos de ranas

La fuerza media del golpe de la Rana A fue de 0,71 newtons (N), y la de la Rana B fue de 0,42 N, para una diferencia de 0,29 N. Es posible que ambas ranas golpeen con la misma fuerza y que esta diferencia observada sea fruto del azar. Vas a calcular la probabilidad de obtener al menos 0,29 N de diferencia en la fuerza media del golpe bajo la hipótesis de que las distribuciones de fuerzas de ambas ranas son idénticas. Usamos una prueba de permutación con una estadística de prueba basada en la diferencia de medias para contrastar esta hipótesis.

Para tu comodidad, los datos se han guardado en los arrays force_a y force_b.

Este ejercicio forma parte del curso

Pensamiento estadístico en Python (Parte 2)

Ver curso

Instrucciones del ejercicio

  • Define una función con firma diff_of_means(data_1, data_2) que devuelva la diferencia de medias entre dos conjuntos de datos: media de data_1 menos media de data_2.
  • Usa esta función para calcular la diferencia empírica de medias observada en las ranas.
  • Genera 10.000 réplicas por permutación de la diferencia de medias.
  • Calcula el valor p.
  • Imprime el valor p.

Ejercicio interactivo práctico

Prueba este ejercicio y completa el código de muestra.

def diff_of_means(data_1, data_2):
    """Difference in means of two arrays."""

    # The difference of means of data_1, data_2: diff
    diff = ____

    return diff

# Compute difference of mean impact force from experiment: empirical_diff_means
empirical_diff_means = ____

# Draw 10,000 permutation replicates: perm_replicates
perm_replicates = draw_perm_reps(____, ____,
                                 ____, size=10000)

# Compute p-value: p
p = np.sum(____ >= ____) / len(____)

# Print the result
print('p-value =', p)
Editar y ejecutar código