เหตุใดการควบคุม FDR จึงเข้มงวดน้อยกว่าการควบคุม FWER


16

ฉันได้อ่านแล้วว่าการควบคุม FDR นั้นเข้มงวดน้อยกว่าการควบคุม FWER เช่นในWikipedia :

ขั้นตอนการควบคุม FDR ออกแรงควบคุมที่เข้มงวดน้อยกว่าการค้นพบที่ผิดพลาดเมื่อเปรียบเทียบกับขั้นตอนอัตราข้อผิดพลาดในระดับครอบครัว (FWER) (เช่นการแก้ไข Bonferroni) สิ่งนี้จะเพิ่มพลังงานที่ค่าใช้จ่ายในการเพิ่มอัตราข้อผิดพลาดประเภทที่ 1 เช่นการปฏิเสธสมมติฐานว่างไม่มีผลเมื่อมันควรได้รับการยอมรับ

แต่ฉันสงสัยว่ามันแสดงให้เห็นว่าเป็นจริงทางคณิตศาสตร์ได้อย่างไร

มีความสัมพันธ์ระหว่าง FDR และ FWER บ้างไหม


4
คุณอ่านเอกสารต้นฉบับหรือไม่ เป็นทุกสิ่งที่เราคาดหวังได้จากบทความสถิติ: แนวคิดพื้นฐานเดียวเรื่องราวที่ชัดเจนและรัดกุมในการบอกตัวอย่างที่เป็นประโยชน์และการพิสูจน์ที่แม่นยำ (สั้น ๆ )
พระคาร์ดินัล

คำตอบ:


10

แน่นอน @ cardinal ค่อนข้างถูกต้องว่ากระดาษมีความชัดเจนเท่าที่จะได้รับ ดังนั้นสำหรับสิ่งที่คุ้มค่าในกรณีที่คุณไม่สามารถเข้าถึงเอกสารได้นี่เป็นรุ่นที่อธิบายรายละเอียดเล็กน้อยเกี่ยวกับวิธีที่ Benjamini – Hochberg โต้แย้ง:

FDR คือมูลค่าที่คาดว่าจะมีสัดส่วนของการปฏิเสธเท็จvที่จะปฏิเสธทุกr ตอนนี้Rคือเห็นได้ชัดว่าผลรวมของการปฏิเสธเท็จและถูกต้อง; โทรหลังsQevrrs

โดยสรุป (ใช้ตัวอักษรตัวใหญ่สำหรับตัวแปรสุ่มและตัวอักษรตัวเล็กสำหรับค่าที่รับรู้)

Qe=E(VR)=E(VV+S)=:E(Q).

หนึ่งจะใช้เวลาถ้าR = 0Q=0R=0

ทีนี้มันมีความเป็นไปได้สองอย่าง: โมฆะทั้งหมดเป็นจริงหรือแค่m 0 < mของพวกมันนั้นเป็นจริง ในกรณีแรกที่มีไม่สามารถปฏิเสธที่ถูกต้องเพื่อให้R = V ดังนั้นหากมีการปฏิเสธใด ๆ ( R 1 ) Q = 1มิฉะนั้นQ = 0 ดังนั้นmm0<mr=vr1q=1q=0

FDR=E(Q)=1P(Q=1)+0P(Q=0)=P(Q=1)=P(V1)=FWER

ดังนั้นในกรณีนี้เช่นว่าขั้นตอนใด ๆ ที่ควบคุมF D Rนั้นจะควบคุมF W E R เป็นอย่างมากและในทางกลับกันFDR=FWERFDRFWER

ในกรณีที่สองซึ่งถ้าวี> 0 (ดังนั้นหากมีอย่างน้อยหนึ่งปฏิเสธเท็จ) เราเห็นได้ชัดมี (ซึ่งเป็นส่วนที่มียังมีวีในตัวหาร) ที่วี/ R 1 นี่ก็หมายความว่าฟังก์ชั่นตัวบ่งชี้ที่ใช้ค่า 1 ถ้ามีอย่างน้อยหนึ่งปฏิเสธเท็จ1 V 1จะไม่น้อยกว่าQ , 1 V 1 Q ทีนี้ลองคิดถึงความไม่เท่าเทียมกันทั้งสองข้างซึ่งก็คือความน่าเบื่อของEm0<mv>0vv/r11V1Q1V1QE ทิ้งความไม่เท่าเทียมเอาไว้

E(1V1)E(Q)=FDR

คาดว่าค่าตัวของฟังก์ชั่นตัวบ่งชี้เป็นความน่าจะเป็นของเหตุการณ์ในตัวบ่งชี้ที่เรามีอีกครั้งซึ่งเป็นF W E RE(1V1)=P(V1)FWER

ดังนั้นเมื่อเรามีขั้นตอนที่ควบคุมในแง่ที่ว่าF W E Rอัลฟ่าเราจะต้องมีที่F D Rอัลฟ่าFWERFWERαFDRα

ตรงกันข้ามมีควบคุมในบางαอาจมาพร้อมกับขนาดใหญ่อย่างมีนัยสำคัญF W E R การยอมรับส่วนที่ไม่เป็นศูนย์ของการปฏิเสธที่ผิดพลาด ( F D R ) จากสมมติฐานที่ทดสอบทั้งหมดอาจมีความเป็นไปได้สูงมากที่จะมีการปฏิเสธที่ผิดพลาดอย่างน้อยหนึ่งครั้ง ( F W E R )FDRαFWERFDRFWER

ดังนั้นขั้นตอนจะต้องเข้มงวดน้อยลงเมื่อต้องการเพียงการควบคุมซึ่งเป็นสิ่งที่ดีสำหรับพลังงาน นี่เป็นแนวคิดเดียวกันกับในการทดสอบสมมติฐานขั้นพื้นฐาน: เมื่อคุณทดสอบที่ระดับ 5% คุณจะปฏิเสธบ่อยขึ้น (ทั้งโมฆะที่ถูกต้องและเท็จ) มากกว่าเมื่อทดสอบที่ระดับ 1% เพียงเพราะคุณมีค่าวิกฤตน้อยกว่าFDR


2
(+1) การแสดงออกที่ดี เห็นได้ชัดว่าในกรณีแรกเราสามารถพูดได้ว่าการควบคุม FWER หมายถึงการควบคุม FDR (ซึ่งเป็นเรื่องที่สงสัย) นอกจากนี้มันอาจจะคุ้มค่าที่ชี้ให้เห็นว่าคุณสมบัตินี้มาพร้อมกับสมมติฐานการกระจาย (เช่นความเป็นอิสระ) ในสถิติการทดสอบซึ่งแตกต่างจากขั้นตอนที่ให้ไว้ในกระดาษต้นฉบับสำหรับการควบคุมของ FDR
พระคาร์ดินัล
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.