บทสรุปของ Neyman-Pearson สามารถนำไปใช้กับกรณีที่โมฆะอย่างง่ายและทางเลือกไม่ได้อยู่ในตระกูลเดียวกันของการแจกแจงหรือไม่?


15
  1. บทแทรกของ Neyman-Pearson สามารถนำไปใช้กับกรณีที่เป็นโมฆะง่ายและทางเลือกง่าย ๆ ไม่ได้เป็นของครอบครัวเดียวกันของการแจกแจง? จากการพิสูจน์ฉันไม่เห็นว่าทำไมถึงทำไม่ได้

    ตัวอย่างเช่นเมื่อ Simple Null เป็นการแจกแจงแบบปกติและทางเลือกง่าย ๆ คือการแจกแจงแบบเอ็กซ์โพเนนเชียล

  2. คือการทดสอบอัตราส่วนความน่าจะเป็นวิธีที่ดีในการทดสอบ null คอมโพสิตกับทางเลือกคอมโพสิตเมื่อทั้งสองเป็นของครอบครัวที่แตกต่างกันของการกระจาย?

ขอบคุณและขอแสดงความนับถือ!


ตอนนี้เป็นคำถามที่ดี
Glen_b

1
อย่างที่คุณพูดในคำถามหลักฐานไม่ได้ตั้งสมมติฐานเกี่ยวกับรูปแบบของการแจกแจงสองแบบ เชื่อถือคณิตศาสตร์
สีฟ้า

@Cyan: การทดสอบอัตราส่วนความน่าจะเป็นวิธีที่ดีสำหรับคอมโพสิตโมฆะและทางเลือกคอมโพสิตซึ่งเป็นของครอบครัวที่มีการแจกแจงที่แตกต่างกันหรือไม่?
ทิม

เพื่อชี้แจงความคิดเห็นก่อนหน้าของฉัน: ฉันมักจะเห็นคนพูดว่า "ไม่" - แน่นอนมันดูเหมือนว่าแม้ในเอกสาร : - "[การทดสอบอัตราส่วนความน่าจะเป็น]] ... อาจไม่สามารถใช้เพื่ออ้างถึงรูปแบบการทำงานของการกระจายข้อมูล " มันคงจะดีถ้าการยืนยันแบบนั้นไม่ได้ออกมาบ่อยนัก
Glen_b -Reinstate Monica

3
นี่คือไม่ใช่เพราะคำถามใด ๆทั้งสองแตกต่างกันกระจายและGเป็นส่วนหนึ่งของครอบครัวหนึ่งอย่างต่อเนื่องพารามิเตอร์{ P F + ( 1 - P ) G } , , 0 P 1 FG{pF+(1p)G},0p1
whuber

คำตอบ:


11

ใช่ Neyman Pearson Lemma สามารถนำไปใช้กับกรณีที่เป็นโมฆะง่ายและทางเลือกง่าย ๆ ไม่ได้อยู่ในตระกูลเดียวกันของการแจกแจง

ให้เราต้องการสร้างการทดสอบที่ทรงพลังที่สุด (MP) ของเทียบกับH 1 : X Exp ( 1 )ของขนาดH0:XN(0,1)H1:XExp(1)

สำหรับเฉพาะฟังก์ชันที่สำคัญของเราโดย Neyman Pearson บทแทรกคือk

ϕ(x)={1,f1(x)f0(x)>k0,Otherwise

คือการทดสอบ MP ของเทียบกับ H 1ของขนาดH0H1

นี่

r(x)=f1(x)f0(x)=ex12πex2/2=2πe(x22x)

โปรดทราบว่า ตอนนี้ถ้าคุณวาดรูปภาพของr(x)[ฉันไม่ทราบวิธีการสร้างรูปภาพเป็นคำตอบ] จากกราฟมันจะชัดเจนว่าr(x)>k

r(x)=2πe(x22x)(x1){<0,x<1>0,x>1
r(x)r(x)>kx>c

ดังนั้นสำหรับ particualr φ ( x ) = { 1 , x > 0 , มิฉะนั้น คือการทดสอบ MP ของH oกับH 1ขนาดของมันc

ϕ(x)={1,x>c0,Otherwise
HoH1

คุณสามารถทดสอบ

    1. เทียบกับH1:XCauchy(0,1)H0:XN(0,12)H1:XCauchy(0,1)
    2. เทียบกับ H 1 : X Cauchy ( 0 , 1 )H0:XN(0,1)H1:XCauchy(0,1)
    3. เทียบกับ H 1 : X เลขชี้กำลังสองเท่า( 0 , 1 )H0:XN(0,1)H1:XDouble Exponential(0,1)

โดย Neyman Pearson บทแทรก

θ

นั่นคือทั้งหมดจากฉัน


5

ไตรมาสที่ 2 อัตราส่วนความน่าจะเป็นของสถิติการทดสอบที่สมเหตุสมผล แต่ (a) Neyman-Pearson Lemma ไม่ได้ใช้กับสมมติฐานเชิงประกอบดังนั้น LRT จึงไม่จำเป็นต้องมีประสิทธิภาพมากที่สุด & (b) ทฤษฎีบทของวิลค์ใช้กับสมมติฐานที่ซ้อนกันเท่านั้นดังนั้นหากครอบครัวหนึ่งเป็นกรณีพิเศษของอีกคนหนึ่ง (เช่นเอกซ์โปเนนเชียล / ไวเบอร์ล, ปัวซอง / ลบทวินาม) คุณไม่ทราบการกระจายตัวของอัตราส่วนความน่าจะเป็นภายใต้ค่าว่าง แม้ไม่แสดงอาการ


"... คุณไม่ทราบว่าการกระจายตัวของอัตราส่วนความน่าจะเป็นภายใต้ค่า null แม้จะเป็นแบบไม่แสดงอาการ" นั่นไม่ใช่เรื่องที่น่าเป็นห่วงมากในโลกที่คุณสามารถเขียนรหัสการจำลองภายใต้ค่า Null ได้น้อยกว่า 20 บรรทัดของอาร์
สีฟ้า

@Cyan: การเขียน 20 บรรทัดนั้นอาจต้องใช้ความคิดบ้าง โปรดจำไว้ว่ามันเป็นคอมโพสิตโมฆะโดยทั่วไปเราจะไม่มีเดือยและฉันไม่คิดว่า LR จะต้องเป็นเดือยที่ประมาณ ฉันคิดว่าคุณสามารถเรียนรู้ LR ...
Scortchi - Reinstate Monica

4
  1. คุณพูดถูก ภาพทั่วไปคือ: เราต้องการสถิติทดสอบที่ให้พลังงานสูงสุดในระดับความสำคัญที่กำหนดα. กล่าวอีกนัยหนึ่งคือวิธีการคำนวณค่าφ เพื่อให้จุดส่วนหนึ่งของพื้นที่พารามิเตอร์ที่ φ เกินกว่า αเสื้อชั่วโมง ควอนไทล์ภายใต้ H0 มีน้ำหนักน้อยที่สุดภายใต้ H1. บทสรุปของ Neyman-Pearson แสดงให้เห็นว่าสถิตินั้นเป็นอัตราส่วนความน่าจะเป็น

  2. กระดาษต้นฉบับของ Neyman & Pearson ยังกล่าวถึงการตั้งสมมติฐานร่วมกัน ในบางกรณีคำตอบนั้นตรงไปตรงมา - หากมีทางเลือกของการแจกแจงเฉพาะในแต่ละครอบครัวที่มีอัตราส่วนความน่าจะเป็นแบบอนุรักษ์นิยมเมื่อนำมาใช้กับทั้งครอบครัว นี่คือสิ่งที่เกิดขึ้นบ่อยครั้งตัวอย่างเช่นสำหรับสมมติฐานที่ซ้อนกัน มันง่ายสำหรับสิ่งนี้ที่จะไม่เกิดขึ้น บทความนี้โดย Cox กล่าวถึงสิ่งที่ต้องทำเพิ่มเติม ฉันคิดว่าวิธีการที่ทันสมัยกว่านี้คือการเข้าใกล้แบบเบย์โดยให้นักบวชดูแลทั้งสองครอบครัว


มีการอ้างอิงที่ดี - กระดาษ Cox
Scortchi - Reinstate Monica
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.