สารฆ่าแมลงกลุ่ม neonicotinoid ส่งผลที่ไม่ได้ตั้งใจหรือไม่?
ในแบบฝึกหัดสุดท้ายด้านการทดสอบสมมติฐานก่อนจะนำทุกอย่างมารวมกันในกรณีศึกษาของบทถัดไป เราจะศึกษาผลกระทบของสารฆ่าแมลงกลุ่ม neonicotinoid ที่มีต่อการสืบพันธุ์ของผึ้ง สารฆ่าแมลงเหล่านี้ถูกใช้อย่างแพร่หลายในสหรัฐอเมริกาเพื่อกำจัดเพลี้ยและแมลงศัตรูพืชชนิดอื่นๆ
ในการศึกษาล่าสุด Straub และคณะ (Proc. Roy. Soc. B, 2016) ได้ศึกษาผลกระทบของ neonicotinoids ต่ออสุจิของผึ้งที่ช่วยผสมเกสร ในแบบฝึกหัดนี้และแบบฝึกหัดถัดไป เราจะวิเคราะห์ว่าการรับสารฆ่าแมลงส่งผลต่อจำนวนอสุจิที่มีชีวิตต่อสารละลายครึ่งมิลลิลิตรอย่างไร
เริ่มต้นด้วยการทำ EDA ตามปกติ ให้พล็อต ECDF ของจำนวนอสุจิที่มีชีวิตสำหรับผึ้งที่ไม่ได้รับสาร (เก็บไว้ในอาร์เรย์ NumPy ชื่อ control) และผึ้งที่ได้รับสารฆ่าแมลง (เก็บไว้ในอาร์เรย์ NumPy ชื่อ treated)
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Statistical Thinking in 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()