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

Есть ли у неоникотиноидных инсектицидов непредвиденные последствия?

В качестве финального упражнения по проверке гипотез — прежде чем мы объединим всё вместе в итоговом разборе в следующей главе — вы изучите влияние неоникотиноидных инсектицидов на размножение пчёл. Эти инсектициды широко применяются в Соединённых Штатах для борьбы с тлёй и другими вредителями, наносящими ущерб растениям.

В недавнем исследовании Straub и соавт. (Proc. Roy. Soc. B, 2016) изучили влияние неоникотиноидов на сперму пчёл-опылителей. В этом и следующем упражнении вы проанализируете, как обработка пестицидом повлияла на количество живых сперматозоидов в половине миллилитра семенной жидкости.

Для начала, как обычно, проведём разведочный анализ данных. Постройте ECDF количества живых сперматозоидов для необработанных пчёл (хранятся в массиве NumPy control) и пчёл, обработанных пестицидом (хранятся в массиве NumPy treated).

Это упражнение является частью курса

Статистическое мышление на Python (часть 2)

Посмотреть курс

Инструкции к упражнению

  • С помощью функции ecdf() получите значения x,y из массивов control и treated для построения ECDF.
  • Постройте оба ECDF на одном графике.
  • Отступы, легенда и подписи осей уже заданы за вас. Нажмите «Отправить ответ», чтобы увидеть результат!

Интерактивное практическое упражнение

Попробуйте выполнить это упражнение, дополнив этот пример кода.

# 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()
Редактировать и запускать код