Aan de slagGa gratis aan de slag

Permutatietoets: wildtype versus heterozygoot

Toets met een permutatietoets de hypothese dat de boutlengtes van heterozygoten en wildtype identiek zijn verdeeld.

Deze oefening maakt deel uit van de cursus

Casestudies in statistisch denken

Cursus bekijken

Oefeninstructies

  • Bereken het verschil van de gemiddelden (heterozygoot min wildtype boutlengtes) van de echte datasets en sla het resultaat op in de variabele diff_means_exp. De numpy-arrays bout_lengths_wt en bout_lengths_het staan al in je namespace.
  • Trek 10.000 permutatiereplicaten van het verschil van de gemiddelden met dcst.draw_perm_reps(). Je kunt ook de functie dcst.diff_of_means() gebruiken; sla je resultaat op in perm_reps.
  • Bereken de p-waarde, waarbij "minstens zo extreem als" betekent dat het verschil van de gemiddelden onder de nulhypothese groter dan of gelijk aan is aan wat experimenteel werd waargenomen.
  • Print de p-waarde op het scherm.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

# 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)
Code bewerken en uitvoeren