¿La profundidad del pico es heredable en G. scandens?
La heredabilidad de la profundidad del pico en G. scandens parece baja. Podría ser que esta heredabilidad observada se deba al azar y que, en realidad, la profundidad del pico no sea heredable en la especie. Aquí pondrás a prueba esa hipótesis. Para hacerlo, realizarás una prueba de permutación por parejas.
Este ejercicio forma parte del curso
Pensamiento estadístico en Python (Parte 2)
Instrucciones del ejercicio
- Inicializa tu array de réplicas de heredabilidad. Tomaremos 10.000 réplicas de permutación por parejas.
- Escribe un bucle
forpara generar tus réplicas.- Permuta el array
bd_parent_scandensusandonp.random.permutation(). - Calcula la heredabilidad entre el array permutado y el array
bd_offspring_scandensusando la funciónheritability()que escribiste en el último ejercicio. Guarda el resultado en el array de réplicas.
- Permuta el array
- Calcula el valor p como el número de réplicas que son mayores que la
heritability_scandensobservada que calculaste en el último ejercicio.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Initialize array of replicates: perm_replicates
perm_replicates = ____
# Draw replicates
for i in range(10000):
# Permute parent beak depths
bd_parent_permuted = ____
perm_replicates[i] = ____
# Compute p-value: p
p = np.sum(____ >= ____) / len(____)
# Print the p-value
print('p-val =', p)