FDR / FWER
المقارنات المتعددة: FDR مقابل FWER
يؤدي اختبار فرضيات كثيرة إلى زيادة فرص النتائج الإيجابية الخاطئة. تغيّر طرق التصحيح عتبة اعتبار النتيجة ذات دلالة. ينبغي مقارنة ضبط الأخطاء والقدرة على اكتشاف التأثيرات الحقيقية معًا.
لماذا تحتاج الاختبارات المتعددة إلى تصحيح؟
مع 20 فرضية عدم صحيحة ومستقلة وα = 0.05، تبلغ احتمالية ظهور نتيجة إيجابية خاطئة واحدة على الأقل 1 − (1 − 0.05)²⁰ ≈ 64%. والعدد المتوقع للنتائج الإيجابية الخاطئة هو 20 × 0.05 = 1. هاتان كميتان مختلفتان: احتمال وعدد.
ما معدل الاكتشافات الخاطئة؟
FDR = E[V / max(R, 1)]، حيث V عدد الاكتشافات الخاطئة وR إجمالي الاكتشافات. يضبط هدف FDR البالغ 5% هذه النسبة في المتوسط عبر التجارب المتكررة، وفق افتراضات الطريقة. ولا يعني أن احتمال خطأ كل نتيجة ذات دلالة يبلغ 5%.
كيف يختلف FWER؟
FWER = P(V ≥ 1): احتمال حدوث اكتشاف خاطئ واحد على الأقل ضمن مجموعة اختبارات. تستهدف طريقتا بونفيروني وهولم FWER، بينما تستهدف بنجاميني–هوخبرغ FDR. يضع ضبط FWER حدًا لـFDR أيضًا، لكنه غالبًا ما يكتشف عددًا أقل من التأثيرات الحقيقية.
مثال مشروح: بونفيروني مقابل BH
خذ خمس قيم p مرتبة: 0.001، 0.008، 0.039، 0.041، 0.300، مع α = q = 0.05. تستخدم بونفيروني 0.05 / 5 = 0.01 وترفض الفرضيتين الأوليين. تقارن BH الرتب بالعتبات 0.01، 0.02، 0.03، 0.04، 0.05. أكبر رتبة تجتاز العتبة هي 2، لذا ترفض BH الفرضيتين الأوليين أيضًا. لو كانت قيمة p الأخيرة 0.045، لرفضت BH الفرضيات الخمس كلها: اعتمد أكبر رتبة تجتاز العتبة، لا أول رتبة تخفق.
قارن طرق التصحيح الـ11
دون تصحيح
ارفض كل فرضية عندما p < α. لا يوجد ضبط للخطأ على مستوى مجموعة الاختبارات. هذه نقطة مرجعية مفيدة لرؤية تراكم النتائج الإيجابية الخاطئة.
Bonferroni
ترفض بونفيروني عندما p < α/m. تضبط FWER للقيم p الصالحة مهما كانت بنية الاعتماد، لكنها قد تفوّت تأثيرات حقيقية عند إجراء اختبارات كثيرة.
Holm
رتّب قيم p تصاعديًا. قارن الرتبة r مع α/(m−r+1)، وتوقف عند أول إخفاق. تضبط طريقة هولم FWER تحت أي اعتماد، وترفض على الأقل كل ما ترفضه بونفيروني.
Šidák
استخدم العتبة الواحدة 1−(1−α)^(1/m). تضبط طريقة شيداك FWER للاختبارات المستقلة، وهي أقل تحفظًا قليلًا من بونفيروني.
Holm–Šidák
طبّق منطق هولم التدريجي التنازلي باستخدام العتبات 1−(1−α)^(1/(m−r+1)). تستخدم هذه المحاكاة اختبارات مستقلة؛ لا تفترض أن هذه الطريقة تعالج أي اعتماد.
Hochberg
اعثر على أكبر رتبة r تحقق p(r) ≤ α/(m−r+1)، وارفض جميع الرتب حتى r. تضبط طريقة هوخبرغ FWER في حالة الاستقلال أو وفق شروط مناسبة للاعتماد الإيجابي.
Hommel
تستخدم طريقة هومل إجراء الاختبار المغلق المبني على اختبارات سايمز. تضبط FWER في حالة الاستقلال أو في ظروف تسري فيها متباينة سايمز، ولا تقل قوتها عن طريقة هوخبرغ. لكنها أكثر كلفة حسابيًا.
Benjamini–Hochberg
اعثر على أكبر رتبة r تحقق p(r) ≤ (r/m)q، ثم ارفض جميع الرتب حتى r. تضبط طريقة بنجاميني–هوخبرغ FDR مع الاختبارات المستقلة وبعض حالات الاعتماد الإيجابي (PRDS).
Benjamini–Yekutieli
استخدم BH مع q مقسومًا على Hm = Σ(1/k)، حيث k = 1…m. تضبط طريقة بنجاميني–يكوتييلي FDR عند أي اعتماد، عادةً على حساب القوة.
على نهج ستوري (استكشافية)
يقدّر هذا البديل الاستكشافي المبسّط π₀ من نسبة قيم p التي تتجاوز 0.5، ثم يطبّق BH عند q/π̂₀. ليس تنفيذًا كاملًا لقيم q ولا يوفر ضمانًا عامًا لـFDR في العينات محدودة الحجم.
BH التكيفية (استكشافية)
يطبّق هذا البديل الاستكشافي ذو المرحلتين BH عند q أولًا، ثم يقدّر m₀ على أنه max(1, m−R₁)، وبعدها يطبّق BH عند q·m/m̂₀. وهو ليس إجراء BKY ولا نقدّمه على أنه يتمتع بضمانه للعينات محدودة الحجم.
أسئلة شائعة
هل هذه حاسبة FDR لقيم p الخاصة بي؟
هذا محاكي تعليمي، وليس حاسبة لرفع قيم p. يولّد تجارب نعرف حقيقتها كي ترى الاكتشافات الخاطئة والتأثيرات الفائتة، وهي أمور لا يمكن عادةً معرفتها من البيانات الحقيقية وحدها.
هل تضمن قيمة q البالغة 0.05 أن تكون النتائج الإيجابية الخاطئة 5% فقط؟
لا. يتعلق الهدف q = 0.05 بنسبة متوقعة بين الاكتشافات وفق افتراضات الإجراء. وقد تكون نسبة الاكتشافات الخاطئة في تجربة واحدة أعلى أو أدنى.
أي طريقة تصحيح أختار؟
حدّد الأخطاء التي يجب ضبطها وفق تصميم بحثك. تحد طرق FWER من خطر أي اكتشاف خاطئ ضمن مجموعة الاختبارات. وتضبط طرق FDR متوسط نسبة الاكتشافات الخاطئة بين الاكتشافات. تحقّق من افتراضات الاعتماد وحدّد مجموعة الفرضيات قبل فحص النتائج.