วิธีการทดสอบค่ามัธยฐานของอารมณ์ในอาร์
การทดสอบค่ามัธยฐานของอารมณ์ ใช้เพื่อเปรียบเทียบค่ามัธยฐานของกลุ่มอิสระตั้งแต่ 2 กลุ่มขึ้นไป
ฟังก์ชัน median_test ของคลัง เหรียญ สามารถใช้เพื่อทำการทดสอบนี้ใน R ซึ่งใช้ไวยากรณ์ต่อไปนี้:
median_test (ตอบกลับ ~ กลุ่มข้อมูล)
ทอง:
- การตอบสนอง: เวกเตอร์ของค่าการตอบสนอง
- group: เวกเตอร์ของค่าการจัดกลุ่ม
- ข้อมูล: กรอบข้อมูลที่มีการตอบสนองและเวกเตอร์กลุ่ม
ตัวอย่างต่อไปนี้แสดงให้เห็นถึงวิธีการใช้ฟังก์ชันนี้เพื่อทำการทดสอบอารมณ์ค่ามัธยฐานใน R
ตัวอย่าง: การทดสอบค่ามัธยฐานของอารมณ์ใน R
สมมติว่าครูต้องการทราบว่าวิธีการศึกษาสองวิธีที่แตกต่างกันทำให้คะแนนสอบของนักเรียนในชั้นเรียนแตกต่างกันหรือไม่ เพื่อทดสอบสิ่งนี้ เธอสุ่มให้นักเรียน 10 คนใช้วิธีการเรียนแบบหนึ่ง และนักเรียนอีก 10 คนใช้วิธีการเรียนอีกวิธีหนึ่ง หลังจากผ่านไปสองสัปดาห์ นักเรียนแต่ละคนจะสอบแบบเดียวกัน
เธอตัดสินใจใช้การทดสอบค่ามัธยฐานของ Mood’s เพื่อพิจารณาว่าคะแนนสอบมัธยฐานของทั้งสองกลุ่มแตกต่างกันหรือไม่
ขั้นตอนที่ 1: สร้างกรอบข้อมูล
#createdata method = rep(c('method1', 'method2'), each=10) score = c(75, 77, 78, 83, 83, 85, 89, 90, 91, 97, 77, 80, 84, 84, 85, 90, 92, 92, 94, 95) examData = data.frame(method, score) #viewdata examData method score 1 method1 75 2 method1 77 3 method1 78 4 method1 83 5 method1 83 6 method1 85 7 method1 89 8 method1 90 9 method1 91 10 method1 97 11 method2 77 12 method2 80 13 method2 84 14 method2 84 15 method2 85 16 method2 90 17 method2 92 18 method2 92 19 method2 94 20 method2 95
ขั้นตอนที่ 2: ทำการทดสอบอารมณ์แบบมัธยฐาน
#load the coin library library(corner) #perform Mood's Median Test median_test(score~method, data = examData) #output Asymptotic Two-Sample Brown-Mood Median Test data: score by method (method1, method2) Z = -0.43809, p-value = 0.6613 alternative hypothesis: true mu is not equal to 0
ค่า p ของการทดสอบคือ 0.6613 เนื่องจากค่านี้ไม่น้อยกว่า 0.05 เราจึงไม่สามารถปฏิเสธสมมติฐานว่างได้ เราไม่มีหลักฐานเพียงพอที่จะกล่าวว่ามีความแตกต่างที่มีนัยสำคัญทางสถิติในคะแนนสอบมัธยฐานระหว่างทั้งสองกลุ่ม
ตามค่าเริ่มต้น ฟังก์ชันนี้จะกำหนดคะแนนเป็น 0 ให้กับการสังเกตที่เท่ากับค่ามัธยฐานทุกประการ อย่างไรก็ตาม คุณสามารถระบุค่านี้เป็น 0.5 หรือ 1 ได้โดยใช้อาร์กิวเมนต์ mid.score
ตัวอย่างเช่น โค้ดต่อไปนี้ดำเนินการทดสอบอารมณ์ค่ามัธยฐานเดียวกันทุกประการ แต่จะกำหนดค่า 0.5 ให้กับการสังเกตที่เท่ากับค่ามัธยฐาน:
#perform Mood's Median Test median_test(score~method, mid.score="0.5" , data = examData) #output Asymptotic Two-Sample Brown-Mood Median Test data: score by method (method1, method2) Z = -0.45947, p-value = 0.6459 alternative hypothesis: true mu is not equal to 00
ค่า p ทดสอบกลายเป็น 0.6459 ซึ่งต่ำกว่าค่า p ก่อนหน้าเล็กน้อยที่ 0.6613 อย่างไรก็ตาม บทสรุปของการทดสอบยังคงเหมือนเดิม: เราไม่มีหลักฐานเพียงพอที่จะบอกว่าค่ามัธยฐานของการสอบระหว่างทั้งสองกลุ่มมีความแตกต่างกันอย่างมีนัยสำคัญ