Перестановочний тест на даних про жаб
Середня сила удару Жаби A становила 0,71 ньютона (N), а Жаби B — 0,42 N; різниця дорівнює 0,29 N. Можливо, що жаби б'ють із однаковою силою, а спостережувану різницю спричинив випадок. Ви обчислите ймовірність отримати принаймні 0,29 N різниці в середній силі удару за гіпотези, що розподіли сил удару для обох жаб ідентичні. Ми використаємо перестановочний тест із статистикою тесту «різниця середніх», щоб перевірити цю гіпотезу.
Для зручності дані збережено в масивах force_a і force_b.
Ця вправа є частиною курсу
Статистичне мислення в Python (Частина 2)
Інструкції до вправи
- Визначте функцію з сигнатурою виклику
diff_of_means(data_1, data_2), яка повертає різницю середніх між двома наборами даних: середнєdata_1мінус середнєdata_2. - Використайте цю функцію, щоб обчислити емпіричну різницю середніх, спостережену для жаб.
- Згенеруйте 10 000 перестановочних реплік різниці середніх.
- Обчисліть p-значення.
- Виведіть p-значення.
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
def diff_of_means(data_1, data_2):
"""Difference in means of two arrays."""
# The difference of means of data_1, data_2: diff
diff = ____
return diff
# Compute difference of mean impact force from experiment: empirical_diff_means
empirical_diff_means = ____
# Draw 10,000 permutation replicates: perm_replicates
perm_replicates = draw_perm_reps(____, ____,
____, size=10000)
# Compute p-value: p
p = np.sum(____ >= ____) / len(____)
# Print the result
print('p-value =', p)