Kom igångKom igång gratis

Permutationstest: vildtyp mot heterozygot

Testa hypotesen att heterozygotens och vildtypens boutlängder är identiskt fördelade med hjälp av ett permutationstest.

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

Fallstudier i statistiskt tänkande

Visa kurs

Övningsinstruktioner

  • Beräkna skillnaden i medelvärde (heterozygot minus vildtyp för boutlängder) för de faktiska datamängderna och lagra resultatet i variabeln diff_means_exp. numpy-arrayerna bout_lengths_wt och bout_lengths_het finns redan i ditt namnrum.
  • Dra 10 000 permutationsreplikat av skillnaden i medelvärde med hjälp av dcst.draw_perm_reps(). Du kan även använda funktionen dcst.diff_of_means() och lagra resultatet i perm_reps.
  • Beräkna p-värdet, där "minst lika extremt som" definieras som att skillnaden i medelvärde under nollhypotesen är större än eller lika med det experimentellt observerade värdet.
  • Skriv ut p-värdet på skärmen.

Interaktiv övning med praktiskt arbete

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

# Compute the difference of means: diff_means_exp
diff_means_exp = ____ - ____

# Draw permutation replicates: perm_reps
perm_reps = ____(____, ____, 
                               ____, size=____)

# Compute the p-value: p-val
p_val = ____(____ >= ____) / len(____)

# Print the result
print('p =', p_val)
Redigera och kör kod