Test permutacyjny: typ dziki a heterozygota
Sprawdź hipotezę, że długości epizodów u heterozygoty i typu dzikiego mają identyczny rozkład – użyj do tego testu permutacyjnego.
To ćwiczenie jest częścią kursu
Studia przypadków w myśleniu statystycznym
Instrukcje do ćwiczenia
- Oblicz różnicę średnich (długości epizodów heterozygoty minus typ dziki) na rzeczywistych danych i zapisz wynik w zmiennej
diff_means_exp. Tablicenumpyo nazwachbout_lengths_wtibout_lengths_hetsą już dostępne w przestrzeni nazw. - Wygeneruj 10 000 replik permutacyjnych różnicy średnich, korzystając z funkcji
dcst.draw_perm_reps(). Możesz użyć funkcjidcst.diff_of_means()– wynik zapisz wperm_reps. - Oblicz wartość p, przyjmując za „co najmniej tak ekstremalną" taką różnicę średnich pod hipotezą zerową, która jest większa od zaobserwowanej eksperymentalnie lub jej równa.
- Wyświetl wartość p na ekranie.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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)