Generera permutationsreplikat
Som beskrevs i videon är ett permutationsreplikat ett enskilt värde av en statistika beräknad från ett permutationsurval. Precis som funktionen draw_bs_reps() som du skrev i kapitel 2 är användbar för att generera bootstrap-replikat, är det praktiskt att ha en liknande funktion, draw_perm_reps(), för att generera permutationsreplikat. Den funktionen skriver du i den här övningen.
Funktionen har signaturen draw_perm_reps(data_1, data_2, func, size=1). Viktigt att notera: func måste vara en funktion som tar två arrayer som argument. I de flesta fall är func en funktion du skriver själv.
Den här övningen är en del av kursen
Statistiskt tänkande i Python (del 2)
Övningsinstruktioner
- Definiera en funktion med signaturen:
draw_perm_reps(data_1, data_2, func, size=1).- Initialisera en array för att lagra permutationsreplikaten med hjälp av
np.empty(). - Skriv en
for-slinga som:- Beräknar ett permutationsurval med hjälp av din
permutation_sample()-funktion - Skickar urvalen till
func()för att beräkna replikatet och lagrar resultatet i din array med replikat.
- Beräknar ett permutationsurval med hjälp av din
- Returnera arrayen med replikat.
- Initialisera en array för att lagra permutationsreplikaten med hjälp av
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
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