1สูตรและวิธีทำ
สถิติสำหรับการปฏิบัติที่อิงหลักฐาน
- ความไวและความจำเพาะ
- Se = a/(a + c), Sp = d/(b + d)
- ค่าพยากรณ์
- PPV = a/(a + b), NPV = d/(c + d)
- ความแม่นยำ
- (a + d)/n
- อัตราส่วนความควรจะเป็น (likelihood ratio)
- LR+ = Se/(1 − Sp), LR− = (1 − Se)/Sp
- PPV ที่ความชุก p
- Se·p / (Se·p + (1 − Sp)(1 − p))
- ออดส์
- ออดส์ = p/(1 − p); p = ออดส์/(1 + ออดส์); ออดส์หลังตรวจ = ออดส์ก่อนตรวจ × LR
- ตัววัดความเสี่ยง
- RR = EER/CER, ARR = CER − EER, RRR = ARR/CER
- จำนวนผู้ป่วยที่ต้องรักษาเพื่อให้ได้ประโยชน์ (NNT)
- NNT = 1/ARR ปัดขึ้น; NNH = 1/(EER − CER)
- อัตราส่วนออดส์
- OR = ad/(bc)
- อัตราอุบัติการณ์
- ผู้ป่วยรายใหม่ / person-time
วิธีการ
- วาดตาราง 2 × 2 และระบุ a, b, c และ d ก่อนคำนวณ
- อ่านความไวและความจำเพาะตามคอลัมน์ อ่านค่าพยากรณ์ตามแถว
- สำหรับความชุกอื่น ให้นึกภาพ 10 000 คนและเติมตารางด้วย Se, Sp และความชุก
- สำหรับอัตราส่วนความควรจะเป็น ให้แปลงความน่าจะเป็นเป็นออดส์ คูณ แล้วแปลงกลับ
- สำหรับการทดลอง ให้หา EER และ CER ก่อน แล้ว RR = EER/CER, ARR = CER − EER และ NNT = 1/ARR ปัดขึ้น
สถิติเชิงพรรณนา
- ค่าเฉลี่ย
- x̄ = Σx/n
- ส่วนเบี่ยงเบนมาตรฐานตัวอย่าง
- s = √(Σ(x − x̄)²/(n − 1))
- มัธยฐาน
- ค่ากลางของข้อมูลที่เรียงลำดับแล้ว
- คะแนนมาตรฐาน z
- z = (x − μ)/σ
- การแจกแจงปกติ
- 68% อยู่ภายใน ±1σ, 95% อยู่ภายใน ±2σ
ช่วงความเชื่อมั่น
- ทราบ σ
- x̄ ± z·σ/√n
- ไม่ทราบ σ
- x̄ ± t·s/√n, df = n − 1
- ค่าวิกฤต z
- 90%: 1.645, 95%: 1.96, 99%: 2.576
- ขนาดตัวอย่าง
- n = (z·σ/E)², rounded up
วิธีการ
- ตัดสินใจระหว่าง z (ทราบ σ ของประชากร) กับ t (มีเพียง s ของตัวอย่าง โดยองศาเสรี n − 1)
- หาค่าวิกฤตสำหรับระดับความเชื่อมั่น
- ค่าความคลาดเคลื่อน: E = ค่าวิกฤต × ส่วนเบี่ยงเบนมาตรฐาน/√n
- ช่วง: x̄ ± E
การทดสอบสมมติฐานและช่วงความเชื่อมั่น
- การทดสอบ z สำหรับค่าเฉลี่ย
- z = (x̄ − μ0)/(σ/√n)
- การทดสอบ t สำหรับค่าเฉลี่ย
- t = (x̄ − μ0)/(s/√n), df = n − 1
- สองทางที่ระดับ 5%
- ปฏิเสธ H0 ถ้า |z| > 1.96
- ช่วงความเชื่อมั่น
- x̄ ± z·σ/√n
วิธีการ
- ระบุ H₀ และ H₁ และตัดสินว่าการทดสอบเป็นทางเดียวหรือสองทาง
- คำนวณตัวสถิติทดสอบ: (x̄ − μ₀)/(σ/√n) หรือใช้ s และ t
- เปรียบเทียบกับค่าวิกฤต หรือหาค่า p
- ปฏิเสธ H₀ ถ้าตัวสถิติเลยค่าวิกฤต (p ต่ำกว่า α) และระบุข้อสรุปเป็นคำพูด
2โจทย์ฝึก
1.สถิติสำหรับการปฏิบัติที่อิงหลักฐานปานกลาง
แบบฝึกหัด (ค่าสมมติ ไม่ใช่คำแนะนำทางคลินิก) ในการศึกษาแบบ cohort มีการติดตามกลุ่มคนเป็นระยะเวลาต่างกัน (ดูตาราง: จำนวนคนและจำนวนปีของการติดตามแต่ละกลุ่ม) ระหว่างการติดตามมีผู้ป่วยรายใหม่ 11 ราย จงหาอัตราอุบัติการณ์ IR ต่อ 1000 person-years ปัดเป็นทศนิยมหนึ่งตำแหน่ง
50 × 2.5 ปี60 × 5 ปี11 รายใหม่
IR = ต่อ 1000 person-years
2.สถิติสำหรับการปฏิบัติที่อิงหลักฐานยาก
แบบฝึกหัด (ค่าสมมติ ไม่ใช่คำแนะนำทางคลินิก) ผู้ป่วยมีความน่าจะเป็นก่อนตรวจที่จะเป็นโรค 10% การตรวจที่มีความไว 84% และความจำเพาะ 85% ให้ผลเป็น บวก จงหาความน่าจะเป็นหลังตรวจ p′ ที่จะเป็นโรคเป็นเปอร์เซ็นต์ ปัดเป็นทศนิยมหนึ่งตำแหน่ง
ก่อนตรวจ 10%Se = 84%Sp = 85%บวก
p′ = %
3.สถิติเชิงพรรณนาปานกลาง
หาส่วนเบี่ยงเบนมาตรฐานตัวอย่าง s ของข้อมูล: 54, 57, 33, 48, 54
n = 5
s =
4.สถิติเชิงพรรณนายาก
ปริมาณหนึ่งมีการแจกแจงปกติ ค่าเฉลี่ย 60 และส่วนเบี่ยงเบนมาตรฐาน 8 สัดส่วน P ของค่าที่ต่ำกว่า 76 คือเท่าใด
μ = 60σ = 8x = 76
P =
5.ช่วงความเชื่อมั่นปานกลาง
ตัวอย่าง 8 ชิ้นมีค่าเฉลี่ย 11 และส่วนเบี่ยงเบนมาตรฐาน s = 1 ประชากรมีการแจกแจงใกล้เคียงปกติ หาค่าความคลาดเคลื่อน E ของช่วงความเชื่อมั่น 90% สำหรับค่าเฉลี่ย
n = 8x̄ = 11s = 190%t(7) = 1.895
E =
6.ช่วงความเชื่อมั่นยาก
เพลาที่กลึง 11 ชิ้นมีส่วนเบี่ยงเบนมาตรฐานของตัวอย่าง 10.1 µm สมมติว่าเส้นผ่านศูนย์กลางมีการแจกแจงปกติ จงหา ขีดจำกัดล่าง σ ของช่วงความเชื่อมั่น 95% สำหรับส่วนเบี่ยงเบนมาตรฐานของประชากร เมื่อ df = 10 การแจกแจง χ² 2.5% อยู่ต่ำกว่า 3.247 และ 2.5% อยู่สูงกว่า 20.483
n = 11s = 10.1 µmχ² = 3.247, 20.483
σ = µm
7.การทดสอบสมมติฐานและช่วงความเชื่อมั่นปานกลาง
พนักงาน 16 คนทำแบบทดสอบการพิมพ์ก่อนและหลังหลักสูตรอบรม บันทึกผลต่าง d = หลัง − ก่อน (คำต่อนาที) ของแต่ละคน ค่าเฉลี่ยของผลต่างคือ 2.1 และส่วนเบี่ยงเบนมาตรฐานของผลต่างคือ 8.4 จงหาค่าสถิติทดสอบแบบจับคู่ t สำหรับ H₀: μd = 0
n = 16d̄ = 2.1sd = 8.4
t =
8.การทดสอบสมมติฐานและช่วงความเชื่อมั่นยาก
ตัวอย่างการวัด 16 ค่ามีค่าเฉลี่ย 102.3 ส่วนเบี่ยงเบนมาตรฐานของประชากรคือ 5 หาขีดจำกัดล่าง L ของช่วงความเชื่อมั่น 90% สำหรับค่าเฉลี่ย
x̄ = 102.3σ = 5n = 1690%
L =
3สารละลาย
1.สถิติสำหรับการปฏิบัติที่อิงหลักฐาน
คำตอบ: IR = 25.9 ต่อ 1000 person-years
อัตราส่วนความควรจะเป็นและตัววัดความเสี่ยง: LR+ = ความไว/(1 − ความจำเพาะ), LR− = (1 − ความไว)/ความจำเพาะ ในการทดลอง EER และ CER คืออัตราการเกิดเหตุการณ์เมื่อมีและไม่มีการรักษา: RR = EER/CER, ARR = CER − EER, RRR = ARR/CER และ NNT = 1/ARR ปัดขึ้น ในการศึกษาแบบ case-control ใช้ได้เฉพาะอัตราส่วนออดส์ OR = ad/(bc) อัตราอุบัติการณ์ = ผู้ป่วยรายใหม่/person-time
อัตราอุบัติการณ์
IR = 11/425·1000 = 25.9 ต่อ 1000 person-years
IR = 25.9 ต่อ 1000 person-years
2.สถิติสำหรับการปฏิบัติที่อิงหลักฐาน
คำตอบ: p′ = 38.4 %
ค่าพยากรณ์ขึ้นกับความชุก: PPV = Se·p/(Se·p + (1 − Sp)(1 − p)) และ NPV = Sp(1 − p)/(Sp(1 − p) + (1 − Se)p) เมื่อใช้อัตราส่วนความควรจะเป็น: ออดส์ก่อนตรวจ = p/(1 − p), ออดส์หลังตรวจ = ออดส์ก่อนตรวจ × LR, ความน่าจะเป็น = ออดส์/(1 + ออดส์) สำหรับประชากร N บวกจริง = N·p·Se และบวกเท็จ = N·(1 − p)·(1 − Sp) NNT = 1/(CER × RRR) ปัดขึ้น
อัตราส่วนความควรจะเป็น
LR+ = 0.84/(1 − 0.85) = 5.6
ความน่าจะเป็นหลังตรวจ
p′ = 0.6222/(1 + 0.6222) = 0.3836 ≈ 38.4%
3.สถิติเชิงพรรณนา
คำตอบ: s = 9.628
สถิติเชิงพรรณนา: ค่าเฉลี่ยคือผลรวมหารด้วย n มัธยฐานคือค่ากลางของข้อมูลที่เรียงลำดับแล้ว (ค่าเฉลี่ยของสองค่ากลางเมื่อ n เป็นจำนวนคู่) ส่วนเบี่ยงเบนมาตรฐานตัวอย่างคือ s = √(Σ(x − x̄)²/(n − 1)) คะแนนมาตรฐาน z = (x − μ)/σ บอกว่าค่าหนึ่งอยู่ห่างจากค่าเฉลี่ยกี่เท่าของส่วนเบี่ยงเบนมาตรฐาน
กำลังสองของความเบี่ยงเบน
Σ(x − x̄)² = 4.8² + 7.8² + (−16.2)² + (−1.2)² + 4.8² = 370.8
4.สถิติเชิงพรรณนา
คำตอบ: P = 0.9772
สถิติเชิงพรรณนา: ค่าเฉลี่ยคือผลรวมหารด้วย n มัธยฐานคือค่ากลางของข้อมูลที่เรียงลำดับแล้ว (ค่าเฉลี่ยของสองค่ากลางเมื่อ n เป็นจำนวนคู่) ส่วนเบี่ยงเบนมาตรฐานตัวอย่างคือ s = √(Σ(x − x̄)²/(n − 1)) คะแนนมาตรฐาน z = (x − μ)/σ บอกว่าค่าหนึ่งอยู่ห่างจากค่าเฉลี่ยกี่เท่าของส่วนเบี่ยงเบนมาตรฐาน
ตารางการแจกแจงปกติมาตรฐาน
5.ช่วงความเชื่อมั่น
คำตอบ: E = 0.67
ช่วงความเชื่อมั่นของค่าเฉลี่ย: x̄ ± ค่าความคลาดเคลื่อน เมื่อทราบ σ ของประชากร ค่าความคลาดเคลื่อนคือ z·σ/√n เมื่อมีเพียงส่วนเบี่ยงเบนมาตรฐานตัวอย่าง s ให้ใช้ t ที่องศาเสรี n − 1: t·s/√n ค่าวิกฤต: z = 1.645 (90%), 1.96 (95%), 2.576 (99%) เพื่อให้ได้ค่าความคลาดเคลื่อน E หรือน้อยกว่า ต้องมี n ≥ (z·σ/E)² ปัดขึ้น
ค่าวิกฤต
df = n − 1 = 7, t = 1.895
6.ช่วงความเชื่อมั่น
คำตอบ: σ = 7.057 µm
ช่วงสำหรับความแปรปรวน: (n − 1)s²/σ² มีการแจกแจง χ² ที่มีองศาเสรี n − 1 ช่วง 95% สำหรับ σ² ไปจาก (n − 1)s² หารด้วยค่า χ² ด้านบน (2.5% อยู่เหนือค่านั้น) ถึง (n − 1)s² หารด้วยค่า χ² ด้านล่าง (2.5% อยู่ต่ำกว่าค่านั้น) ช่วงสำหรับ σ คือรากที่สองของทั้งสองขีดจำกัด และไม่สมมาตรรอบ s
ผลรวมกำลังสอง
(n − 1)s² = 10·10.1² = 1020
ขีดจำกัดของ σ²
1020/20.483 = 49.8
1020/3.247 = 314.2
ช่วง 95% สำหรับ σ: [7.057, 17.72] µm
7.การทดสอบสมมติฐานและช่วงความเชื่อมั่น
คำตอบ: t = 1
การทดสอบค่าเฉลี่ย: ตัวสถิติทดสอบวัดว่าค่าเฉลี่ยตัวอย่างอยู่ห่างจากค่าที่อ้างใน H₀ กี่เท่าของความคลาดเคลื่อนมาตรฐาน: z = (x̄ − μ₀)/(σ/√n) เมื่อทราบ σ และ t = (x̄ − μ₀)/(s/√n) เมื่อประมาณจากตัวอย่าง ช่วงความเชื่อมั่นคือ x̄ ± (ค่าวิกฤต)·(ความคลาดเคลื่อนมาตรฐาน)
ความคลาดเคลื่อนมาตรฐานของผลต่างเฉลี่ย
เมื่อ df = 15 ค่าวิกฤตสองด้านที่ 5% คือ 2.131 |t| น้อยกว่า จึงไม่ปฏิเสธ H₀
8.การทดสอบสมมติฐานและช่วงความเชื่อมั่น
คำตอบ: L = 100.2
การทดสอบค่าเฉลี่ย: ตัวสถิติทดสอบวัดว่าค่าเฉลี่ยตัวอย่างอยู่ห่างจากค่าที่อ้างใน H₀ กี่เท่าของความคลาดเคลื่อนมาตรฐาน: z = (x̄ − μ₀)/(σ/√n) เมื่อทราบ σ และ t = (x̄ − μ₀)/(s/√n) เมื่อประมาณจากตัวอย่าง ช่วงความเชื่อมั่นคือ x̄ ± (ค่าวิกฤต)·(ความคลาดเคลื่อนมาตรฐาน)
ขีดจำกัดของช่วง
L = 102.3 − 2.056 = 100.2