เปลี่ยนค่าส่วนเบี่ยงเบนมาตรฐานของการแจกแจงแบบปกติ
เราจะสำรวจข้อมูลส่วนสูงของผู้ชายอเมริกันวัยผู้ใหญ่ต่อไป ซึ่งทราบแล้วว่าเป็นการแจกแจงแบบปกติที่มีค่าเฉลี่ย 177 เซนติเมตร และส่วนเบี่ยงเบนมาตรฐาน 8 เซนติเมตร
ในแบบฝึกหัดนี้ จะทำการสุ่มตัวอย่างจากการแจกแจงแบบปกติและคำนวณช่วงความเชื่อมั่น 95% ของค่าเฉลี่ยส่วนสูง แต่คราวนี้จะเปลี่ยนค่าส่วนเบี่ยงเบนมาตรฐานเป็น 15 โดยไม่เปลี่ยนค่าเฉลี่ย แล้วสังเกตว่าค่าเฉลี่ยและช่วงความเชื่อมั่นของส่วนสูงจะเปลี่ยนแปลงไปอย่างไร
มีการนำเข้าไลบรารีต่อไปนี้ให้แล้ว: random, NumPy เป็น np และโมดูล stats ของ SciPy เป็น st
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Monte Carlo Simulations ใน Python
คำแนะนำการฝึกหัด
- สุ่มตัวอย่าง 1,000 ครั้งจากการแจกแจงแบบปกติที่มีค่าเฉลี่ย 177 และส่วนเบี่ยงเบนมาตรฐาน 15 แล้วเก็บผลลัพธ์ไว้ใน
heights_177_15
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
random.seed(1231)
heights_177_8 = st.norm.rvs(loc=177, scale=8, size=1000)
print(np.mean(heights_177_8))
upper = np.quantile(heights_177_8, 0.975)
lower = np.quantile(heights_177_8, 0.025)
print([lower, upper])
# Sample 1,000 times from the normal distribution where the standard deviation is 15
heights_177_15 = ____
print(np.mean(heights_177_15))
upper = np.quantile(heights_177_15, 0.975)
lower = np.quantile(heights_177_15, 0.025)
print([lower, upper])