LoslegenKostenlos loslegen

Is beak depth heritable at all in G. scandens?

The heritability of beak depth in G. scandens seems low. It could be that this observed heritability was just achieved by chance and beak depth is actually not really heritable in the species. You will test that hypothesis here. To do this, you will do a pairs permutation test.

Diese Übung ist Teil des Kurses

Statistical Thinking in Python (Part 2)

Kurs anzeigen

Anleitung zur Übung

  • Initialize your array of replicates of heritability. We will take 10,000 pairs permutation replicates.
  • Write a for loop to generate your replicates.
    • Permute the bd_parent_scandens array using np.random.permutation().
    • Compute the heritability between the permuted array and the bd_offspring_scandens array using the heritability() function you wrote in the last exercise. Store the result in the replicates array.
  • Compute the p-value as the number of replicates that are greater than the observed heritability_scandens you computed in the last exercise.

Interaktive Übung

Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.

# 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)
Code bearbeiten und ausführen