Permutationstest: Wildtyp versus Heterozygot
Teste mit einem Permutationstest die Hypothese, dass die Bout-Längen von Heterozygoten und Wildtyp identisch verteilt sind.
Diese Übung ist Teil des Kurses
Fallstudien zum statistischen Denken
Anleitung zur Übung
- Berechne die Differenz der Mittelwerte (Heterozygot minus Wildtyp-Bout-Längen) der tatsächlichen Datensätze und speichere das Ergebnis in der Variable
diff_means_exp. Dienumpy-Arraysbout_lengths_wtundbout_lengths_hetsind bereits in deinem Namespace vorhanden. - Ziehe 10.000 Permutationsreplikate der Differenz der Mittelwerte mit
dcst.draw_perm_reps(). Du kannst auch die Funktiondcst.diff_of_means()verwenden und dein Ergebnis inperm_repsspeichern. - Berechne den p-Wert, wobei "mindestens so extrem wie" so definiert ist, dass die Differenz der Mittelwerte unter der Nullhypothese größer ist als oder gleich der experimentell beobachteten.
- Gib den p-Wert auf dem Bildschirm aus.
Interaktive Übung
Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.
# 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)