กลับไปหน้าบทความ

อ่าน 5 นาที

การทดสอบอัตราส่วนความน่าจะเป็นแบบลำดับ

การทดสอบอัตราส่วนความน่าจะเป็นแบบลำดับ (SPRT) เป็นการทดสอบสมมติฐานแบบลำดับ เฉพาะ ที่พัฒนาโดยAbraham Wald และต่อมาได้รับการพิสูจน์ว่าเหมาะสมที่สุดโดย Wald และJacob Wolfowitz...

การทดสอบอัตราส่วนความน่าจะเป็นแบบลำดับ

การทดสอบอัตราส่วนความน่าจะเป็นแบบลำดับ (SPRT) เป็นการทดสอบสมมติฐานแบบลำดับ เฉพาะ ที่พัฒนาโดยAbraham Wald [ 1 ]และต่อมาได้รับการพิสูจน์ว่าเหมาะสมที่สุดโดย Wald และJacob Wolfowitz [ 2 ] ผลลัพธ์ของ Neyman และ Pearson ในปี 1933เป็นแรงบันดาลใจให้ Wald ปรับปรุงใหม่เป็นปัญหาการวิเคราะห์แบบลำดับ ในทางตรงกันข้าม บทพิสูจน์ของ Neyman-Pearson เสนอหลักเกณฑ์คร่าวๆสำหรับกรณีที่รวบรวมข้อมูลทั้งหมดแล้ว (และทราบอัตราส่วนความน่าจะเป็น)

แม้ว่าเดิมทีจะพัฒนาขึ้นเพื่อใช้ใน การศึกษา การควบคุมคุณภาพในด้านการผลิต แต่ SPRT ได้รับการกำหนดสูตรเพื่อใช้ในการทดสอบด้วยคอมพิวเตอร์ของผู้เข้ารับการทดสอบที่เป็นมนุษย์ในฐานะเกณฑ์การยุติ[ 3 ] [ 4 ] [ 5 ]

ทฤษฎี

เช่นเดียวกับการทดสอบสมมติฐาน แบบคลาสสิก SPRT เริ่มต้นด้วยสมมติฐานสองข้อ เช่นชม0{\displaystyle H_{0}}และชม1{\displaystyle H_{1}}สำหรับสมมติฐานว่างและสมมติฐานทางเลือกตามลำดับ จะต้องระบุให้ชัดเจนดังนี้:

ชม0:พี=พี0{\displaystyle H_{0}:p=p_{0}}
ชม1:พี=พี1{\displaystyle H_{1}:p=p_{1}}

ขั้นตอนต่อไปคือการคำนวณผลรวมสะสมของอัตราส่วนลอการิทึมความน่าจะเป็นบันทึกΛฉัน{\displaystyle \log \แลมบ์ดา _{i}}เมื่อมีข้อมูลใหม่เข้ามา: ด้วยเอส0=0{\displaystyle S_{0}=0}ดังนั้นสำหรับฉัน{\displaystyle i}=1,2,...,

เอสฉัน=เอสฉัน1+บันทึกΛฉัน{\displaystyle S_{i}=S_{i-1}+\log \Lambda _{i}}

กฎการหยุดคือระบบการกำหนดเกณฑ์อย่างง่าย:

  • เอ<เอสฉัน<{\displaystyle a<S_{i}<b}: ติดตามสถานการณ์อย่างต่อเนื่อง ( ความไม่เท่าเทียมที่สำคัญ )
  • เอสฉัน{\displaystyle S_{i}\geq b}: ยอมรับชม1{\displaystyle H_{1}}
  • เอสฉันเอ{\displaystyle S_{i}\leq a}: ยอมรับชม0{\displaystyle H_{0}}

ที่ไหนเอ{\displaystyle a}และ{\displaystyle b}(เอ<0<<{\displaystyle a<0<b<\infty }ขึ้นอยู่กับข้อผิดพลาดประเภทที่ 1 และประเภทที่ 2ที่ ต้องการα{\displaystyle \alpha }และเบต้า{\displaystyle \beta }สามารถเลือกได้ดังต่อไปนี้:

เอบันทึกเบต้า1α{\displaystyle a\approx \log {\frac {\beta }{1-\alpha }}}และบันทึก1เบต้าα{\displaystyle b\approx \log {\frac {1-\beta }{\alpha }}}

กล่าวอีกนัยหนึ่งคือα{\displaystyle \alpha }และเบต้า{\displaystyle \beta }ต้องตัดสินใจล่วงหน้าเพื่อตั้งค่าเกณฑ์ให้เหมาะสม ค่าตัวเลขจะขึ้นอยู่กับการใช้งาน เหตุผลที่เป็นเพียงค่าประมาณก็คือ ในกรณีแบบไม่ต่อเนื่อง สัญญาณอาจข้ามเกณฑ์ระหว่างตัวอย่าง ดังนั้น ขึ้นอยู่กับค่าปรับของการเกิดข้อผิดพลาดและความถี่ในการสุ่มตัวอย่างอาจตั้งค่าเกณฑ์ให้เข้มงวดมากขึ้นได้ ขอบเขตที่แน่นอนนั้นถูกต้องในกรณีแบบต่อเนื่อง

ตัวอย่าง

ตัวอย่างในตำราเรียนคือการประมาณค่าพารามิเตอร์ของฟังก์ชันการแจกแจงความน่าจะเป็น พิจารณาการแจกแจงแบบเอกซ์โปเนนเชียล :

เอฟθ(x)=θ1อีxθ,x,θ>0{\displaystyle f_{\theta }(x)=\theta ^{-1}e^{-{\frac {x}{\theta }}},\qquad x,\theta >0}

สมมติฐานคือ

{ชม0:θ=θ0ชม1:θ=θ1θ1>θ0.{\displaystyle {\begin{cases}H_{0}:\theta =\theta _{0}\\H_{1}:\theta =\theta _{1}\end{cases}}\qquad \theta _{1}>\theta _{0}.}

ดังนั้น ฟังก์ชันความน่าจะเป็นล็อก (LLF) สำหรับตัวอย่างหนึ่งคือ

บันทึกΛ(x)=บันทึก(θ11อีxθ1θ01อีxθ0)=บันทึก(θ0θ1อีxθ0xθ1)=บันทึก(θ0θ1)+บันทึก(อีxθ0xθ1)=บันทึก(θ1θ0)+(xθ0xθ1)=บันทึก(θ1θ0)+(θ1θ0θ0θ1)x{\displaystyle {\begin{aligned}\log \Lambda (x)&=\log \left({\frac {\theta _{1}^{-1}e^{-{\frac {x}{\theta _{1}}}}}{\theta _{0}^{-1}e^{-{\frac {x}{\theta _{0}}}}}}\right)\\&=\log \left({\frac {\theta _{0}}{\theta _{1}}}e^{{\frac {x}{\theta _{0}}}-{\frac {x}{\theta _{1}}}}\right)\\&=\log \left({\frac {\theta _{0}}{\theta _{1}}}\right)+\log \left(e^{{\frac {x}{\theta _{0}}}-{\frac {x}{\theta _{1}}}}\right)\\&=-\log \left({\frac {\theta _{1}}{\theta _{0}}}\right)+\left({\frac {x}{\theta _{0}}}-{\frac {x}{\theta _{1}}}\right)\\&=-\log \left({\frac {\theta _{1}}{\theta _{0}}}\right)+\left({\frac {x}{\theta _{0}}}-{\frac {x}{\theta _{1}}}\right)\\&=-\log \left({\frac {\theta _{1}}{\theta _{0}}}\right)+\left({\frac {\theta _{1}-\theta _{0}}{\theta _{0}\theta _{1}}}\right)x\end{aligned}}}

ผลรวมสะสมของ LLF สำหรับทุกxคือ

เอสn=ฉัน=1nบันทึกΛ(xฉัน)=nบันทึก(θ1θ0)+(θ1θ0θ0θ1)ฉัน=1nxฉัน{\displaystyle S_{n}=\sum _{i=1}^{n}\log \Lambda (x_{i})=-n\log \left({\frac {\theta _{1}}{\theta _{0}}}\right)+\left({\frac {\theta _{1}-\theta _{0}}{\theta _{0}\theta _{1}}}\right)\sum _{i=1}^{n}x_{i}}

ดังนั้น กฎการหยุดจึงเป็นดังนี้:

เอ<nบันทึก(θ1θ0)+(θ1θ0θ0θ1)ฉัน=1nxฉัน<{\displaystyle a<-n\log \left({\frac {\theta _{1}}{\theta _{0}}}\right)+\left({\frac {\theta _{1}-\theta _{0}}{\theta _{0}\theta _{1}}}\right)\sum _{i=1}^{n}x_{i}<b}

หลังจากจัดเรียงใหม่ ในที่สุดเราก็พบ

เอ+nบันทึก(θ1θ0)<(θ1θ0θ0θ1)ฉัน=1nxฉัน<+nบันทึก(θ1θ0){\displaystyle a+n\log \left({\frac {\theta _{1}}{\theta _{0}}}\right)<\left({\frac {\theta _{1}-\theta _{0}}{\theta _{0}\theta _{1}}}\right)\sum _{i=1}^{n}x_{i}<b+n\log \left({\frac {\theta _{1}}{\theta _{0}}}\right)}

เส้นขอบเขตนั้นเป็นเพียงเส้นขนานสองเส้นที่มีความชันบันทึก(θ1/θ0){\displaystyle \log(\theta _{1}/\theta _{0})}ควรหยุดการสุ่มตัวอย่างเมื่อผลรวมของตัวอย่างมีค่าเกินขอบเขตที่กำหนดให้ทำการสุ่มตัวอย่างต่อไป

แอปพลิเคชัน

การผลิต

การทดสอบนี้ทำบนตัวชี้วัดสัดส่วน และทดสอบว่าตัวแปรpเท่ากับค่าใดค่าหนึ่งจากสองค่าที่ต้องการ คือ p₂บริเวณระหว่างสองค่านี้เรียกว่าบริเวณที่ไม่แตกต่าง (Indifference Regionหรือ IR) ตัวอย่างเช่น สมมติว่าคุณกำลังทำการศึกษาการควบคุมคุณภาพในล็อตสินค้าของโรงงาน ฝ่ายบริหารต้องการให้ล็อตสินค้ามีสินค้าชำรุดไม่เกิน 3% แต่ 1% หรือน้อยกว่านั้นคือล็อตในอุดมคติที่จะผ่านการทดสอบได้อย่างราบรื่น ในตัวอย่างนี้p₁ = 0.01 p₂ = 0.03 และบริเวณระหว่าง ทั้งสองนี้คือ IR เพราะฝ่ายบริหารพิจารณาว่าล็อตเหล่านี้อยู่ในระดับก้ำกึ่งและยอมรับได้ไม่ว่าจะเป็นอย่างไรก็ตาม จะมีการสุ่มตัวอย่างสินค้าทีละชิ้นจากล็อต (การวิเคราะห์ตามลำดับ) จนกว่าการทดสอบจะระบุได้ภายในระดับความคลาดเคลื่อนที่ยอมรับได้ว่าล็อตนั้นเหมาะสมหรือควรถูกปฏิเสธ

การทดสอบผู้เข้ารับการทดสอบที่เป็นมนุษย์

ปัจจุบัน SPRT เป็นวิธีการหลักในการจำแนกผู้เข้าสอบในการทดสอบการจำแนกประเภทด้วยคอมพิวเตอร์ (CCT) ที่มีความยาวแปรผันได้สองตัว p1 และ p2 จะถูกกำหนดโดยการกำหนดคะแนนตัด (เกณฑ์) สำหรับผู้เข้าสอบในตัวชี้วัดสัดส่วนความถูกต้อง และเลือกจุดที่สูงกว่าและต่ำกว่าคะแนนตัดนั้น ตัวอย่างเช่น สมมติว่าคะแนนตัดถูกกำหนดไว้ที่ 70% สำหรับการทดสอบ เราสามารถเลือกp1 = 0.65 p2 = 0.75 จากนั้นการทดสอบจะประเมินความเป็นไปได้ที่คะแนนที่แท้จริงของผู้เข้าสอบในตัวชี้วัดนั้นจะเท่ากับหนึ่งในสองจุดนั้น หากผู้เข้าสอบมีคะแนน 75% พวกเขาจะผ่าน และพวกเขาจะสอบไม่ผ่านหากมีคะแนน 65%

จุดเหล่านี้ไม่ได้ถูกกำหนดขึ้นโดยพลการอย่างสิ้นเชิง ควรมีการกำหนดเกณฑ์คะแนนตัดด้วยวิธีการที่ถูกต้องตามกฎหมายเสมอ เช่นวิธีการ Angoff ที่ปรับปรุงแล้วอีกครั้ง พื้นที่ความไม่แยแสแสดงถึงช่วงคะแนนที่ผู้ออกแบบข้อสอบยอมรับได้ไม่ว่าจะเป็นผ่านหรือตก พารามิเตอร์บนp คือระดับสูงสุดที่ผู้ออกแบบข้อสอบยินดีรับสำหรับการสอบตก (เพราะทุกคนที่ต่ำกว่าระดับนี้มีโอกาสสอบตกสูง) และพารามิเตอร์ล่างp คือระดับต่ำสุดที่ผู้ออกแบบข้อสอบยินดีรับสำหรับการสอบผ่าน (เพราะทุกคนที่อยู่เหนือกว่าระดับนี้มีโอกาสสอบผ่านพอสมควร) แม้ว่าคำจำกัดความนี้อาจดูเหมือนเป็นภาระเล็กน้อย แต่ลองพิจารณากรณีที่มีความสำคัญสูงของการสอบใบอนุญาตประกอบวิชาชีพแพทย์: เราควรพิจารณาว่าบุคคลนั้นอยู่ในระดับใดระดับหนึ่งในสองระดับนี้เมื่อใด?

แม้ว่า SPRT จะถูกนำมาใช้ในการทดสอบครั้งแรกในยุคทฤษฎีการทดสอบแบบคลาสสิก ดังที่ ได้กล่าวไว้ในย่อหน้าก่อนหน้า แต่ Reckase (1983) เสนอแนะให้ใช้ทฤษฎีการตอบสนองต่อข้อสอบ เพื่อกำหนดพารามิเตอร์ p1และโดยกำหนดเกณฑ์คะแนนตัดและขอบเขตความไม่แตกต่างบนเมตริกความสามารถแฝง (theta) และแปลงเป็นเมตริกสัดส่วนเพื่อการคำนวณ การวิจัยเกี่ยวกับ CCT ตั้งแต่นั้นมาได้นำวิธีการนี้มาใช้ด้วยเหตุผลหลายประการ

  1. โดยทั่วไปแล้วคลังข้อสอบขนาดใหญ่จะถูกปรับเทียบด้วย IRT
  2. วิธีนี้ช่วยให้สามารถกำหนดพารามิเตอร์ได้อย่างแม่นยำยิ่งขึ้น
  3. โดยการใช้ฟังก์ชันการตอบสนองต่อแต่ละรายการ พารามิเตอร์ต่างๆ สามารถเปลี่ยนแปลงได้ง่ายระหว่างรายการต่างๆ

การตรวจจับผลลัพธ์ทางการแพทย์ที่ผิดปกติ

Spiegelhalterและคณะ[ 6 ]ได้แสดงให้เห็นว่า SPRT สามารถใช้ในการตรวจสอบประสิทธิภาพของแพทย์ ศัลยแพทย์ และผู้ประกอบวิชาชีพทางการแพทย์อื่นๆ ในลักษณะที่ให้คำเตือนล่วงหน้าเกี่ยวกับผลลัพธ์ที่อาจผิดปกติ พวกเขาแสดงให้เห็นว่ามันสามารถช่วยระบุHarold Shipmanว่าเป็นฆาตกรได้ก่อนที่จะมีการระบุตัวตนจริง[ 6 ]

ส่วนขยาย

แม็กซ์สเปรต

เมื่อไม่นานมานี้ ในปี 2011 ได้มีการนำเสนอส่วนขยายของวิธีการ SPRT ที่เรียกว่า Maximized Sequential Probability Ratio Test (MaxSPRT) [ 7 ]คุณลักษณะเด่นของ MaxSPRT คือการอนุญาตให้มีสมมติฐานทางเลือกแบบผสมด้านเดียว และการแนะนำขอบเขตการหยุดด้านบน วิธีการนี้ถูกนำไปใช้ในการศึกษาวิจัยทางการแพทย์หลายครั้ง[ 8 ]

ดูเพิ่มเติม

อ่านเพิ่มเติม

  • Ghosh, Bhaskar Kumar (1970). การทดสอบสมมติฐานทางสถิติแบบลำดับ . สำนักพิมพ์: Addison-Wesley .
  • Holger Wilker: Sequential-Statistik ใน der Praxis , BoD, Norderstedt 2012, ISBN 978-3848232529.
  • การทดสอบอัตราส่วนความน่าจะเป็นแบบลำดับของ WaldสำหรับRโดย Stéphane Bottine
  • การทดสอบอัตราส่วนความน่าจะเป็นแบบลำดับของ WaldสำหรับPythonโดย Zhenning Yu
ดึงข้อมูลมาจาก " https://en.wikipedia.org/w/index.php?title=Sequential_probability_ratio_test&oldid=1321902819 "

สรุปเนื้อหา

ข้อมูลสำคัญจากบทความ

ข้อมูลสำคัญเกี่ยวกับ การทดสอบอัตราส่วนความน่าจะเป็นแบบลำดับ

การทดสอบอัตราส่วนความน่าจะเป็นแบบลำดับ (SPRT) เป็นการทดสอบสมมติฐานแบบลำดับ เฉพาะ ที่พัฒนาโดยAbraham Wald และต่อมาได้รับการพิสูจน์ว่าเหมาะสมที่สุดโดย Wald และJacob Wolfowitz...

ทฤษฎี

เช่นเดียวกับ การทดสอบสมมติฐาน แบบคลาสสิก SPRT เริ่มต้นด้วยสมมติฐานสองข้อ เช่น ชม 0 {\displaystyle H_{0}} และ ชม 1 {\displaystyle H_{1}} สำหรับ สมมติฐานว่าง และ สมมติฐานทางเลือก ตามลำดับ จะต้องระบุให้ชัดเจนดังนี้:

ตัวอย่าง

ตัวอย่างในตำราเรียนคือ การประมาณค่าพารามิเตอร์ ของ ฟังก์ชันการแจกแจงความน่าจะ เป็น พิจารณา การแจกแจงแบบเอกซ์โปเนนเชียล :

การผลิต

การทดสอบนี้ทำบนตัวชี้วัดสัดส่วน และทดสอบว่าตัวแปร p เท่ากับค่าใดค่าหนึ่งจากสองค่าที่ต้องการ คือ p₂ บริเวณ ระหว่าง สองค่านี้เรียกว่า บริเวณที่ไม่แตกต่าง (Indifference Region หรือ IR) ตัวอย่างเช่น สมมติว่าคุณกำลังทำการศึกษาการควบคุมคุณภาพในล็อตสินค้าของโรงงาน...