Generowanie replikacji permutacyjnych
Jak omówiono w filmie, replikacja permutacyjna to pojedyncza wartość statystyki obliczona na podstawie próby permutacyjnej. Podobnie jak funkcja draw_bs_reps(), którą napisałeś w rozdziale 2, jest przydatna do generowania replikacji bootstrapowych, warto mieć analogiczną funkcję – draw_perm_reps() – do generowania replikacji permutacyjnych. W tym ćwiczeniu napiszesz właśnie tę funkcję.
Funkcja ma sygnaturę draw_perm_reps(data_1, data_2, func, size=1). Istotne jest, że func musi być funkcją przyjmującą dwie tablice jako argumenty. W większości przypadków func będzie funkcją napisaną przez ciebie.
To ćwiczenie jest częścią kursu
Myślenie statystyczne w Pythonie (część 2)
Instrukcje do ćwiczenia
- Zdefiniuj funkcję o następującej sygnaturze:
draw_perm_reps(data_1, data_2, func, size=1).- Zainicjuj tablicę do przechowywania replikacji permutacyjnych, używając
np.empty(). - Napisz pętlę
for, która:- Oblicza próbę permutacyjną przy użyciu funkcji
permutation_sample() - Przekazuje próby do
func(), aby obliczyć replikację, i zapisuje wynik w tablicy replikacji.
- Oblicza próbę permutacyjną przy użyciu funkcji
- Zwróć tablicę replikacji.
- Zainicjuj tablicę do przechowywania replikacji permutacyjnych, używając
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
def draw_perm_reps(data_1, data_2, func, size=1):
"""Generate multiple permutation replicates."""
# Initialize array of replicates: perm_replicates
perm_replicates = ____
for i in ____:
# Generate permutation sample
perm_sample_1, perm_sample_2 = ____
# Compute the test statistic
perm_replicates[i] = ____
return perm_replicates