Перевірка гіпотези: чи стали дзьоби глибшими у 2012 році?
Ваш графік ECDF і визначений довірчий інтервал досить чітко показують, що дзьоби G. scandens на острові Дафна-Мейджор стали глибшими. Але чи можливо, що цей ефект виник випадково? Іншими словами, яка ймовірність отримати спостережену різницю середніх глибин дзьоба, якщо справжні середні однакові?
Будьте уважні! Ми перевіряємо гіпотезу не про те, що глибини дзьоба походять з одного й того самого розподілу. Для цього доречний перестановочний тест. Гіпотеза полягає в тому, що середні рівні. Щоб виконати цю перевірку, потрібно зсунути обидва набори даних так, щоб вони мали однакове середнє, а потім застосувати бутстреп-вибірку, щоб обчислити різницю середніх.
Ця вправа є частиною курсу
Статистичне мислення в Python (Частина 2)
Інструкції до вправи
- Створіть об'єднаний масив глибин дзьоба за 1975 і 2012 роки та обчисліть і збережіть його середнє значення.
- Зсунути
bd_1975іbd_2012так, щоб їхні середні дорівнювали середньому, щойно обчисленому для об'єднаного набору даних. - Візьміть по 10 000 бутстреп-реплікатів середнього для глибин дзьоба 1975 і 2012 років.
- Відніміть реплікати за 1975 рік від реплікатів за 2012 рік, щоб отримати бутстреп-реплікати різниці.
- Обчисліть і виведіть p-значення. Спостережена різниця середніх, яку ви обчислили в попередній вправі, усе ще є у вашому просторі імен як
mean_diff.
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
# Compute mean of combined data set: combined_mean
combined_mean = ____(____((bd_1975, bd_2012)))
# Shift the samples
bd_1975_shifted = ____
bd_2012_shifted = ____
# Get bootstrap replicates of shifted data sets
bs_replicates_1975 = ____
bs_replicates_2012 = ____
# Compute replicates of difference of means: bs_diff_replicates
bs_diff_replicates = ____
# Compute the p-value
p = np.sum(____ >= ____) / len(____)
# Print p-value
print('p =', p)