วิธีทำการทดสอบความพอดีของไคสแควร์ใน python
การ ทดสอบความดีพอดีของไคสแควร์ ใช้เพื่อพิจารณาว่าตัวแปรเชิงหมวดหมู่เป็นไปตามการแจกแจงเชิงสมมุติหรือไม่
บทช่วยสอนนี้จะอธิบายวิธีดำเนินการทดสอบไคสแควร์ Goodness-of-Fit ใน Python
ตัวอย่าง: การทดสอบความพอดีของไคสแควร์ใน Python
เจ้าของร้านบอกว่ามีลูกค้ามาที่ร้านในจำนวนเท่ากันทุกวันในสัปดาห์ เพื่อทดสอบสมมติฐานนี้ นักวิจัยจะบันทึกจำนวนลูกค้าที่เข้าร้านในสัปดาห์ที่กำหนดและพบสิ่งต่อไปนี้:
- วันจันทร์: ลูกค้า 50 คน
- วันอังคาร: ลูกค้า 60 คน
- วันพุธ: ลูกค้า 40 คน
- พฤหัสบดี: ลูกค้า 47 คน
- วันศุกร์: ลูกค้า 53 คน
ใช้ขั้นตอนต่อไปนี้เพื่อทำการทดสอบความพอดีของไคสแควร์ใน Python เพื่อตรวจสอบว่าข้อมูลสอดคล้องกับคำกล่าวอ้างของเจ้าของร้านค้าหรือไม่
ขั้นตอนที่ 1: สร้างข้อมูล
ขั้นแรก เราจะสร้างตารางสองตารางเพื่อระบุจำนวนลูกค้าที่สังเกตและคาดหวังในแต่ละวัน:
expected = [50, 50, 50, 50, 50] observed = [50, 60, 40, 47, 53]
ขั้นตอนที่ 2: ทำการทดสอบความพอดีของไคสแควร์
ต่อไป เราสามารถทำการทดสอบความพอดีของไคสแควร์ได้โดยใช้ ฟังก์ชันไคสแควร์ จากไลบรารี SciPy ซึ่งใช้ไวยากรณ์ต่อไปนี้:
ไคสแควร์ (f_obs, f_exp)
ทอง:
- f_obs: อาร์เรย์ของการนับที่สังเกตได้
- f_exp: อาร์เรย์ของการนับที่คาดหวัง โดยค่าเริ่มต้น แต่ละหมวดหมู่จะถือว่ามีแนวโน้มเท่ากัน
รหัสต่อไปนี้แสดงวิธีใช้ฟังก์ชันนี้ในตัวอย่างเฉพาะของเรา:
import scipy.stats as stats #perform Chi-Square Goodness of Fit Test stats.chisquare(f_obs=observed, f_exp=expected) (statistic=4.36, pvalue=0.35947)
สถิติการทดสอบไคสแควร์คือ 4.36 และค่า p ที่สอดคล้องกันคือ 0.35947
โปรดทราบว่าค่า p สอดคล้องกับค่าไคสแควร์ที่มีดีกรีอิสระ n-1 (dof) โดยที่ n คือจำนวนหมวดหมู่ต่างๆ ในกรณีนี้ dof = 5-1 = 4 คุณสามารถใช้ เครื่องคำนวณไคสแควร์เป็นค่า P เพื่อยืนยันว่าค่า p ที่สอดคล้องกับ X 2 = 4.36 โดยมี dof = 4 เท่ากับ 0.35947
โปรดจำไว้ว่าการทดสอบความดีพอดีของไคสแควร์ใช้สมมติฐานว่างและทางเลือกต่อไปนี้:
- H 0 : (สมมติฐานว่าง) ตัวแปรตามหลังการแจกแจงสมมุติ
- H 1 : (สมมติฐานทางเลือก) ตัวแปรไม่เป็นไปตามการแจกแจงสมมุติฐาน
เนื่องจากค่า p (0.35947) ไม่น้อยกว่า 0.05 เราจึงไม่สามารถปฏิเสธสมมติฐานว่างได้ ซึ่งหมายความว่าเราไม่มีหลักฐานเพียงพอที่จะบอกว่าการกระจายตัวของลูกค้าที่แท้จริงนั้นแตกต่างจากที่เจ้าของร้านค้ารายงาน