Hebben neonicotinoïde insecticiden onbedoelde gevolgen?
Als laatste oefening in hypothesetoetsing voordat we alles samenbrengen in onze casestudy in het volgende hoofdstuk, onderzoek je de effecten van neonicotinoïde insecticiden op de voortplanting van bijen. Deze insecticiden worden in de Verenigde Staten zeer veel gebruikt om bladluizen en andere plagen die planten beschadigen te bestrijden.
In een recente studie onderzochten Straub et al. (Proc. Roy. Soc. B, 2016) de effecten van neonicotinoïden op het sperma van bestuivende bijen. In deze en de volgende oefening bekijk je hoe de pesticidebehandeling invloed had op het aantal levende zaadcellen per halve milliliter sperma.
We beginnen, zoals gebruikelijk, met EDA. Plot ECDF’s van het aantal levende zaadcellen voor onbehandelde bijen (opgeslagen in de NumPy-array control) en bijen behandeld met pesticide (opgeslagen in de NumPy-array treated).
Deze oefening maakt deel uit van de cursus
Statistical Thinking in Python (deel 2)
Oefeninstructies
- Gebruik je
ecdf()-functie omx,y-waarden te genereren uit de arrayscontrolentreatedom de ECDF’s te plotten. - Plot de ECDF’s in dezelfde figuur.
- De marges zijn al voor je ingesteld, net als de legenda en aslabels. Klik op Antwoord verzenden om het resultaat te zien!
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# 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()