EDA: พล็อต ECDF ของความยาว active bout
Active bout คือช่วงเวลาที่ปลาเคลื่อนไหวต่อเนื่องโดยไม่หยุด ให้พล็อต ECDF ของความยาว active bout สำหรับปลาสายพันธุ์กลายพันธุ์และสายพันธุ์ปกติในคืนที่ 7 ของชีวิต ชุดข้อมูลอยู่ใน numpy array ชื่อ bout_lengths_wt และ bout_lengths_mut โดยความยาว bout มีหน่วยเป็นนาที
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
กรณีศึกษาด้านการคิดเชิงสถิติ
คำแนะนำการฝึกหัด
- นำเข้าโมดูล
dc_stat_thinkโดยใช้ชื่อแทนว่าdcstเพื่อให้ใช้ฟังก์ชันต่าง ๆ ได้ - สร้างค่า x และ y สำหรับพล็อต ECDF ของปลาสายพันธุ์ปกติ (
bout_lengths_wt) โดยใช้dcst.ecdf()แล้วเก็บผลลัพธ์ไว้ในnumpyarray ชื่อx_wtและy_wt - ทำเช่นเดียวกันกับปลาสายพันธุ์กลายพันธุ์ (
bout_lengths_mut) โดยเก็บผลลัพธ์ไว้ในnumpyarray ชื่อx_mutและy_mut - ใช้
plt.plot()พล็อต ECDF ทั้งสองเส้นในรูปแบบจุดบนกราฟเดียวกัน โดยต้องระบุ keyword argumentsmarker='.'และlinestyle='none' - แสดงกราฟโดยใช้
plt.show()
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Import the dc_stat_think module as dcst
____
# Generate x and y values for plotting ECDFs
x_wt, y_wt = ____
____, ____ = ____
# Plot the ECDFs
_ = ____(____, ____, ____='.', linestyle='____')
_ = ____(____, ____, ____='.', linestyle='____')
# Make a legend, label axes, and show plot
_ = plt.legend(('wt', 'mut'))
_ = plt.xlabel('active bout length (min)')
_ = plt.ylabel('ECDF')
____