Prueba de permutación: tipo salvaje frente a heterocigoto
Contrasta la hipótesis de que las duraciones de los intervalos de tipo salvaje y de heterocigoto tienen la misma distribución usando una prueba de permutación.
Este ejercicio forma parte del curso
Casos prácticos de pensamiento estadístico
Instrucciones del ejercicio
- Calcula la diferencia de medias (heterocigoto menos tipo salvaje) de los conjuntos de datos reales y guarda el resultado en la variable
diff_means_exp. Los arrays denumpybout_lengths_wtybout_lengths_hetya están en tu espacio de nombres. - Genera 10.000 réplicas de permutación de la diferencia de medias usando
dcst.draw_perm_reps(). Puedes usar también la funcióndcst.diff_of_means(), y guarda el resultado enperm_reps. - Calcula el valor p, definiendo «al menos tan extremo como» como que la diferencia de medias bajo la hipótesis nula sea mayor o igual que la observada experimentalmente.
- Imprime el valor p en pantalla.
ejercicio interactivo práctico
Prueba este ejercicio completando este código de ejemplo.
# Compute the difference of means: diff_means_exp
diff_means_exp = ____ - ____
# Draw permutation replicates: perm_reps
perm_reps = ____(____, ____,
____, size=____)
# Compute the p-value: p-val
p_val = ____(____ >= ____) / len(____)
# Print the result
print('p =', p_val)