НачатьНачать бесплатно

Перестановочный тест на данных о лягушках

Средняя сила удара лягушки A составила 0,71 ньютона (Н), а лягушки B — 0,42 Н, то есть разница равна 0,29 Н. Возможно, обе лягушки бьют с одинаковой силой, а наблюдаемое различие объясняется случайностью. Вы вычислите вероятность получить разность средних значений силы удара не менее 0,29 Н при условии, что распределения силы удара у обеих лягушек идентичны. Для проверки этой гипотезы используется перестановочный тест, в котором тестовой статистикой служит разность средних.

Для удобства данные сохранены в массивах 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)
Редактировать и запускать код