La profondità del becco è ereditaria in G. scandens?
L'ereditarietà della profondità del becco in G. scandens sembra bassa. Potrebbe essere che questa ereditarietà osservata sia dovuta al caso e che, in realtà, la profondità del becco non sia davvero ereditaria in questa specie. Qui metterai alla prova questa ipotesi. Per farlo, eseguirai un test di permutazione a coppie.
Questo esercizio fa parte del corso
Pensiero statistico in Python (Parte 2)
Istruzioni dell'esercizio
- Inizializza il tuo array di replicati dell'ereditarietà. Prenderemo 10.000 replicati di permutazione a coppie.
- Scrivi un ciclo
forper generare i replicati.- Permuta l'array
bd_parent_scandensusandonp.random.permutation(). - Calcola l'ereditarietà tra l'array permutato e l'array
bd_offspring_scandensusando la funzioneheritability()che hai scritto nell'ultimo esercizio. Salva il risultato nell'array dei replicati.
- Permuta l'array
- Calcola il p-value come il numero di replicati che sono maggiori dell'
heritability_scandensosservato che hai calcolato nell'ultimo esercizio.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# 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)