Pairs bootstrap ของข้อมูลการรู้หนังสือ/อัตราการเจริญพันธุ์
นำฟังก์ชันที่เพิ่งเขียนไปใช้ทำ pairs bootstrap เพื่อพล็อตฮิสโตแกรมแสดงค่าประมาณของความชันจากข้อมูลการไม่รู้หนังสือ/อัตราการเจริญพันธุ์ พร้อมทั้งรายงานช่วงความเชื่อมั่น 95% ของค่าความชัน โดยข้อมูลอยู่ใน NumPy arrays ชื่อ illiteracy และ fertility
เพื่อเป็นการทบทวน draw_bs_pairs_linreg() มี function signature เป็น draw_bs_pairs_linreg(x, y, size=1) และคืนค่า 2 ค่า ได้แก่ bs_slope_reps และ bs_intercept_reps
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Statistical Thinking in Python (ตอนที่ 2)
คำแนะนำการฝึกหัด
- ใช้ฟังก์ชัน
draw_bs_pairs_linreg()เพื่อสร้าง bootstrap replicate ของค่าความชันและจุดตัดแกนจำนวน1000ครั้ง โดยใช้illiteracyเป็นข้อมูลแกน x และfertilityเป็นข้อมูลแกน y - คำนวณและแสดงช่วงความเชื่อมั่น bootstrap 95% ของค่าความชัน
- พล็อตและแสดงฮิสโตแกรมของ replicate ค่าความชัน อย่าลืมใส่ชื่อแกนด้วย ส่วนนี้ได้เตรียมไว้ให้แล้ว กด ส่งคำตอบ เพื่อดูฮิสโตแกรมของคุณได้เลย
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Generate replicates of slope and intercept using pairs bootstrap
bs_slope_reps, bs_intercept_reps = ____
# Compute and print 95% CI for slope
print(np.percentile(____, ____))
# Plot the histogram
_ = plt.hist(bs_slope_reps, bins=50, normed=True)
_ = plt.xlabel('slope')
_ = plt.ylabel('PDF')
plt.show()