Kom igångKom igång gratis

Är näbbdjupet ärftligt hos G. scandens?

Ärftligheten för näbbdjup hos G. scandens verkar låg. Det är möjligt att den observerade ärftligheten uppstod av en slump och att näbbdjupet egentligen inte är ärftligt hos arten. Här ska du testa den hypotesen med ett parpermutationstest.

Den här övningen är en del av kursen

Statistiskt tänkande i Python (del 2)

Visa kurs

Övningsinstruktioner

  • Initialisera din array med replikat av ärftlighet. Vi ska ta 10 000 parpermutationsreplikat.
  • Skriv en for-loop för att generera replikaten.
    • Permutera arrayen bd_parent_scandens med np.random.permutation().
    • Beräkna ärftligheten mellan den permuterade arrayen och arrayen bd_offspring_scandens med hjälp av funktionen heritability() som du skrev i förra övningen. Lagra resultatet i replikatarrayen.
  • Beräkna p-värdet som antalet replikat som är större än det observerade heritability_scandens som du beräknade i förra övningen.

Interaktiv övning med praktiskt arbete

Testa den här övningen genom att slutföra den här exempelkoden.

# 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)
Redigera och kör kod