Чи мають неонікотиноїдні інсектициди небажані наслідки?
На завершення розділу з перевірки гіпотез, перш ніж об'єднати все в кейс-стаді у наступному розділі, ви дослідите вплив неонікотиноїдних інсектицидів на розмноження бджіл. Ці інсектициди дуже широко застосовують у США для боротьби з попелицями та іншими шкідниками, що пошкоджують рослини.
У нещодавньому дослідженні Straub та ін. (Proc. Roy. Soc. B, 2016) вивчали вплив неонікотиноїдів на сперму запилювальних бджіл. У цій та наступній вправах ви дослідите, як обробка пестицидом вплинула на кількість живих сперматозоїдів на півмілілітра еякуляту.
Спочатку, як завжди, виконаємо EDA. Побудуйте ECDF для кількості живих сперматозоїдів у необроблених бджіл (збережено в масиві NumPy control) і бджіл, оброблених пестицидом (збережено в масиві NumPy treated).
Ця вправа є частиною курсу
Статистичне мислення в Python (Частина 2)
Інструкції до вправи
- Використайте вашу функцію
ecdf()для генерування пар значеньx,yз масивівcontrolіtreatedдля побудови ECDF. - Побудуйте обидві ECDF на одному графіку.
- Поля, легенду та підписи осей уже налаштовано. Натисніть Submit Answer (Надіслати відповідь), щоб побачити результат!
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
# Compute x,y values for ECDFs
x_control, y_control = ____
x_treated, y_treated = ____
# Plot the ECDFs
plt.plot(____, ____, marker='.', linestyle='none')
plt.plot(____, ____, marker='.', linestyle='none')
# Set the margins
plt.margins(0.02)
# Add a legend
plt.legend(('control', 'treated'), loc='lower right')
# Label axes and show plot
plt.xlabel('millions of alive sperm per mL')
plt.ylabel('ECDF')
plt.show()