Adâncimea ciocului este heritabilă la G. scandens?
Heritabilitatea adâncimii ciocului la G. scandens pare scăzută. Este posibil ca această heritabilitate observată să fi apărut doar întâmplător și că adâncimea ciocului să nu fie de fapt heritabilă la această specie. Vei testa această ipoteză în exercițiul de față, utilizând un test de permutare pe perechi.
Acest exercițiu face parte din cursul
Gândire statistică în Python (Partea 2)
Instrucțiuni pentru exercițiu
- Inițializează array-ul de replici ale heritabilității. Vom genera 10.000 de perechi de replici de permutare.
- Scrie o buclă
forpentru a genera replicile.- Permută array-ul
bd_parent_scandensfolosindnp.random.permutation(). - Calculează heritabilitatea dintre array-ul permutat și array-ul
bd_offspring_scandensfolosind funcțiaheritability()pe care ai scris-o în exercițiul anterior. Stochează rezultatul în array-ul de replici.
- Permută array-ul
- Calculează p-valoarea ca numărul de replici mai mari decât valoarea observată
heritability_scandenscalculată în exercițiul anterior.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
# 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)