Generare repliche di permutazione
Come discusso nel video, una replica di permutazione è un singolo valore di una statistica calcolata a partire da un campione di permutazione. Così come la funzione draw_bs_reps() che hai scritto nel capitolo 2 è utile per generare repliche bootstrap, è utile avere una funzione simile, draw_perm_reps(), per generare repliche di permutazione. In questo esercizio scriverai questa funzione utile.
La funzione ha la seguente firma di chiamata: draw_perm_reps(data_1, data_2, func, size=1). È importante che func sia una funzione che accetta due array come argomenti. Nella maggior parte dei casi, func sarà una funzione che scrivi tu.
Questo esercizio fa parte del corso
Pensiero statistico in Python (Parte 2)
Istruzioni dell'esercizio
- Definisci una funzione con questa firma:
draw_perm_reps(data_1, data_2, func, size=1).- Inizializza un array per contenere le repliche di permutazione usando
np.empty(). - Scrivi un ciclo
forper:- Calcolare un campione di permutazione usando la tua funzione
permutation_sample() - Passare i campioni a
func()per calcolare la replica e salvare il risultato nel tuo array di repliche.
- Calcolare un campione di permutazione usando la tua funzione
- Restituisci l'array di repliche.
- Inizializza un array per contenere le repliche di permutazione usando
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
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