CommencerCommencer gratuitement

Générer des réplicats par permutation

Comme expliqué dans la vidéo, un réplicat par permutation est une valeur unique d’une statistique calculée à partir d’un échantillon par permutation. Tout comme la fonction draw_bs_reps() que vous avez écrite au chapitre 2 vous est utile pour générer des réplicats bootstrap, il est utile d’avoir une fonction similaire, draw_perm_reps(), pour générer des réplicats par permutation. Vous allez écrire cette fonction utile 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.

Cet exercice fait partie du cours

Réflexion statistique en Python (Partie 2)

Afficher le cours

Instructions

  • Définissez une fonction avec la signature suivante : draw_perm_reps(data_1, data_2, func, size=1).
    • Initialisez un tableau destiné à contenir les réplicats par permutation avec np.empty().
    • Écrivez une boucle for pour :
      • Calculer un échantillon par permutation en utilisant votre fonction permutation_sample()
      • Passer ces échantillons à func() pour calculer le réplicat et stocker le résultat dans votre tableau de réplicats.
    • Retournez le tableau de réplicats.

Exercice interactif pratique

Essayez cet exercice en complétant cet exemple de code.

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