Generar réplicas por permutación
Como se comentó en el vídeo, una réplica por permutación es un único valor de una estadística calculada a partir de una muestra por permutación. Igual que la función draw_bs_reps() que escribiste en el capítulo 2 te sirve para generar réplicas bootstrap, también es útil tener una función similar, draw_perm_reps(), para generar réplicas por permutación. En este ejercicio vas a escribir esta función tan útil.
La firma de la función es draw_perm_reps(data_1, data_2, func, size=1). Es importante que func sea una función que reciba dos arrays como argumentos. En la mayoría de los casos, func será una función que tú mismo escribas.
Este ejercicio forma parte del curso
Pensamiento estadístico en Python (Parte 2)
Instrucciones del ejercicio
- Define una función con esta firma:
draw_perm_reps(data_1, data_2, func, size=1).- Inicializa un array para almacenar las réplicas por permutación usando
np.empty(). - Escribe un bucle
forpara:- Calcular una muestra por permutación usando tu función
permutation_sample() - Pasar las muestras a
func()para calcular la réplica y guardar el resultado en tu array de réplicas.
- Calcular una muestra por permutación usando tu función
- Devuelve el array de réplicas.
- Inicializa un array para almacenar las réplicas por permutación usando
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
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