CommencezCommencez gratuitement

Générer des répliques par permutation

Comme on l'a vu dans la vidéo, une réplique par permutation est une seule valeur d'une statistique calculée à partir d'un échantillon permuté. Tout comme la fonction draw_bs_reps() écrite au chapitre 2 vous est utile pour générer des répliques bootstrap, il est utile d'avoir une fonction similaire, draw_perm_reps(), pour générer des répliques par permutation. Vous allez écrire cette fonction dans cet exercice.

La fonction a la signature d'appel draw_perm_reps(data_1, data_2, func, size=1). Point important : func doit être une fonction qui prend deux tableaux en arguments. Dans la plupart des cas, func sera une fonction que vous aurez écrite vous-même.

Cette activité fait partie du cours

Réflexion statistique en Python (partie 2)

Voir le cours

Instructions de l’exercice

  • Définissez une fonction avec cette signature : draw_perm_reps(data_1, data_2, func, size=1).
    • Initialisez un tableau pour stocker les répliques par permutation à l'aide de np.empty().
    • Écrivez une boucle for pour :
      • Calculer un échantillon par permutation à l'aide de votre fonction permutation_sample()
      • Transmettre les échantillons à func() pour calculer la réplique et stocker le résultat dans votre tableau de répliques.
    • Retournez le tableau de répliques.

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

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
Modifier et exécuter le code