Les insecticides néonicotinoïdes ont-ils des effets indésirables ?
Pour un dernier exercice de test d’hypothèse avant de tout rassembler dans notre étude de cas au chapitre suivant, vous allez examiner les effets des insecticides néonicotinoïdes sur la reproduction des abeilles. Ces insecticides sont très largement utilisés aux États-Unis pour lutter contre les pucerons et autres ravageurs qui endommagent les plantes.
Dans une étude récente, Straub et al. (Proc. Roy. Soc. B, 2016) ont étudié les effets des néonicotinoïdes sur les spermatozoïdes des abeilles pollinisatrices. Dans cet exercice et le suivant, vous allez analyser l’effet du traitement pesticide sur le nombre de spermatozoïdes vivants par demi-millilitre de sperme.
Commençons par une analyse exploratoire, comme d’habitude. Tracez les ECDF du nombre de spermatozoïdes vivants pour les abeilles non traitées (stockées dans le tableau NumPy control) et pour les abeilles traitées par pesticide (stockées dans le tableau NumPy treated).
Cet exercice fait partie du cours
Réflexion statistique en Python (Partie 2)
Instructions
- Utilisez votre fonction
ecdf()pour générer les valeursx,yà partir des tableauxcontrolettreatedafin de tracer les ECDF. - Tracez les ECDF sur le même graphique.
- Les marges ont été définies pour vous, ainsi que la légende et les étiquettes des axes. Cliquez sur Soumettre la réponse pour voir le résultat !
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# 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()