Estimation jackknife de base — moyenne
Le rééchantillonnage jackknife est une méthode plus ancienne, moins utilisée que le bootstrap. Elle reste toutefois utile à connaître pour réaliser une estimation jackknife de base. Dans cet exercice, nous allons calculer l'estimation jackknife de la moyenne. Retournons à l'usine de clés anglaises.
Vous êtes propriétaire d'une usine de clés anglaises et vous souhaitez mesurer la longueur moyenne des clés pour vous assurer qu'elles respectent certaines spécifications. Votre usine produit des milliers de clés chaque jour, mais il est irréaliste de mesurer la longueur de chacune. Vous avez cependant accès à un échantillon représentatif de 100 clés. Utilisons l'estimation jackknife pour obtenir la longueur moyenne.
Examinez la variable wrench_lengths dans le Shell.
Cette activité fait partie du cours
Simulation statistique en Python
Instructions de l’exercice
- Obtenez un échantillon jackknife en retirant itérativement une observation de
wrench_lengthset en l'assignant àjk_sample. - Calculez la moyenne de
jk_sampleet ajoutez-la àmean_lengths. - Enfin, calculez l'estimation jackknife
mean_lengths_jkcomme la moyenne du tableaumean_lengths.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# Leave one observation out from wrench_lengths to get the jackknife sample and store the mean length
mean_lengths, n = [], len(wrench_lengths)
index = np.arange(n)
for i in range(n):
jk_sample = ____[index != i]
mean_lengths.append(____)
# The jackknife estimate is the mean of the mean lengths from each sample
mean_lengths_jk = ____(np.array(mean_lengths))
print("Jackknife estimate of the mean = {}".format(mean_lengths_jk))