CommencerCommencer gratuitement

Un test d’hypothèse bootstrap à un échantillon

Une autre jeune grenouille a été étudiée, la grenouille C, et vous souhaitez vérifier si les grenouilles B et C ont des forces d’impact similaires. Malheureusement, vous ne disposez pas des forces d’impact de la grenouille C, mais vous savez que leur moyenne est de 0,55 N. Comme vous n’avez pas les données brutes, vous ne pouvez pas faire de test de permutation et vous ne pouvez pas évaluer l’hypothèse selon laquelle les forces des grenouilles B et C proviennent de la même distribution. Vous allez donc tester une autre hypothèse, moins restrictive : la force moyenne de frappe de la grenouille B est égale à celle de la grenouille C.

Pour mettre en place le test d’hypothèse par bootstrap, vous prendrez la moyenne comme statistique de test. Rappelez-vous : votre objectif est de calculer la probabilité d’obtenir une force d’impact moyenne inférieure ou égale à celle observée pour la grenouille B si l’hypothèse que la vraie moyenne des forces d’impact de la grenouille B est égale à celle de la grenouille C est vraie. Vous commencez par translater toutes les données de la grenouille B de sorte que la moyenne soit 0,55 N. Cela revient à ajouter la moyenne de la force de la grenouille C et à soustraire la moyenne de la force de la grenouille B à chaque mesure de la grenouille B. Les autres propriétés de la distribution de la grenouille B, comme la variance, restent inchangées.

Cet exercice fait partie du cours

Réflexion statistique en Python (Partie 2)

Afficher le cours

Instructions

  • Translatez les forces d’impact de la grenouille B de sorte que sa moyenne soit 0,55 N.
  • Utilisez votre fonction draw_bs_reps() pour effectuer 10 000 réplicats bootstrap de la moyenne de vos forces translatées.
  • Calculez la p-valeur en trouvant la fraction de vos réplicats bootstrap qui sont inférieurs à la force d’impact moyenne observée de la grenouille B. Notez que la variable d’intérêt ici est force_b.
  • Affichez votre p-valeur.

Exercice interactif pratique

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

# Make an array of translated impact forces: translated_force_b
translated_force_b = ____

# Take bootstrap replicates of Frog B's translated impact forces: bs_replicates
bs_replicates = draw_bs_reps(____, ____, 10000)

# Compute fraction of replicates that are less than the observed Frog B force: p
p = np.sum(____ <= np.mean(____)) / 10000

# Print the p-value
print('p = ', ____)
Modifier et exécuter le code