Czy neonikotynoidy mają niezamierzone skutki uboczne?
To ostatnie ćwiczenie z testowania hipotez przed przejściem do studium przypadku w następnym rozdziale. Zbadasz wpływ insektycydów neonikotynoidowych na rozmnażanie pszczół. Środki te są szeroko stosowane w Stanach Zjednoczonych do zwalczania mszyc i innych szkodników niszczących rośliny.
W niedawnym badaniu Straub i in. (Proc. Roy. Soc. B, 2016) zbadali wpływ neonikotynoidów na plemniki pszczół zapylających. W tym i kolejnym ćwiczeniu przeanalizujesz, jak zastosowanie pestycydu wpłynęło na liczbę żywych plemników w połowie mililitra nasienia.
Rozpocznij od EDA, jak zwykle. Narysuj wykresy ECDF dla liczby żywych plemników u pszczół nieleczonych (zapisanych w tablicy NumPy control) oraz u pszczół poddanych działaniu pestycydu (zapisanych w tablicy NumPy treated).
To ćwiczenie jest częścią kursu
Myślenie statystyczne w Pythonie (część 2)
Instrukcje do ćwiczenia
- Użyj funkcji
ecdf(), aby wygenerować wartościx,yz tabliccontrolitreatedpotrzebne do narysowania wykresów ECDF. - Narysuj oba wykresy ECDF na tym samym wykresie.
- Marginesy, legenda i etykiety osi zostały już za ciebie ustawione. Kliknij Prześlij odpowiedź, aby zobaczyć wynik!
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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()