A profundidade do bico é herdável em G. scandens?
A herdabilidade da profundidade do bico em G. scandens parece baixa. Pode ser que essa herdabilidade observada tenha ocorrido por acaso e que a profundidade do bico não seja realmente herdável na espécie. Você vai testar essa hipótese aqui. Para isso, você fará um teste de permutação por pares.
Este exercício faz parte do curso
Pensamento Estatístico em Python (Parte 2)
Instruções do exercício
- Inicialize seu array de réplicas de herdabilidade. Vamos gerar 10.000 réplicas de permutação por pares.
- Escreva um loop
forpara gerar suas réplicas.- Permute o array
bd_parent_scandensusandonp.random.permutation(). - Calcule a herdabilidade entre o array permutado e o array
bd_offspring_scandensusando a funçãoheritability()que você escreveu no último exercício. Armazene o resultado no array de réplicas.
- Permute o array
- Calcule o p-valor como o número de réplicas que são maiores do que a
heritability_scandensobservada que você calculou no último exercício.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# 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)