تُعدّ مسألة التعامل مع البيانات الواقعية واحدةً من أعقد التحديات التي تواجه الباحثين في حقول الإحصاء الرياضي، والعلوم السلوكية، والقياس النفسي، والعلوم الطبية الحيوية؛ فالبيانات المستقاة من الميدان نادراً ما تتطابق مع النماذج النظرية المثالية أو تتبع التوزيع الطبيعي المعياري بدقة متناهية. وفي خضم جمع البيانات وتفسيرها، تبرز ظاهرة القيم الشاذة والمتطرفة (Outliers) والانحرافات التوزيعية كعوامل قادرة على إرباك التحليلات الرياضية الكلاسيكية، وتشويه صورة الواقع، وقيادة متخذي القرار إلى استنتاجات مضللة قد تترتب عليها عواقب بحثية وتطبيقية وخيمة.
من هذا المنطلق، برز مفهوم الإحصائية المقاومة (Resistant Statistic) كركيزة منهجية وجوهرية في الفكر الإحصائي المعاصر. تشير المقاومة الإحصائية إلى الحصانة المنهجية والرياضية التي يتمتع بها مؤشر إحصائي معين ضد التأثيرات الجارفة للقيم المتطرفة أو التغيرات الجذرية في نسب صغيرة من مفردات العينة. إنها تعبير عن قدرة المقياس على تقديم تقدير صادق ومستقر للاتجاه العام أو التشتت البنيوي للمجتمع الإحصائي، دون أن ينحرف نحو الأطراف الشاذة الناتجة عن أخطاء القياس أو الطبيعة غير المتجانسة للظواهر المدروسة.
يتناول هذا المقال التأسيسي الشامل مفهوم المقاومة في الإحصاء الرياضي والوصفي بعمق أكاديمي رصين. سنستعرض الجذور النظرية للإحصاء الحصين، ونحلل الخصائص البنيوية للمقاييس المقاومة وغير المقاومة عبر مقارنات رياضية وتطبيقية دقيقة، مدعومة بمحاكاة رقمية مفصلة ودراسة تطبيقية للأثر السيكومتري والبيولوجي، وصولاً إلى بناء مصفوفة منهجية متكاملة تُرشد الباحثين إلى المعايير الدقيقة لاختيار وتطبيق هذه المقاييس وتوثيقها وفق أعلى المعايير الأكاديمية الدولية.
- 1. مدخل نظري: مفهوم المقاومة في الإحصاء الرياضي والوصفي
- 2. طبيعة القيم المتطرفة والشاذة وتأثيرها على الحسابات الإحصائية
- 3. المقارنة الجوهرية: الإحصاءات المقاومة مقابل الإحصاءات غير المقاومة
- 4. مقاييس النزعة المركزية المقاومة: الوسيط كمعيار قياسي
- 5. مقاييس النزعة المركزية غير المقاومة: هشاشة المتوسط الحسابي
- 6. مقاييس التشتت المقاومة: المدى الربيعي ومقاييس الحصانة المتقدمة
- 7. مقاييس التشتت غير المقاومة: قابلية تأثر الانحراف المعياري والمدى
- 8. دراسة حالة تطبيقية: المحاكاة الرقمية لأثر القيمة المتطرفة
- 9. التطبيقات النفسية والسلوكية: لماذا تهم الإحصاءات المقاومة علماء النفس؟
- 10. طرق العرض والتمثيل البياني للإحصاءات المقاومة
- 11. معايير اتخاذ القرار الإحصائي: متى نختار الإحصاء المقاوم ومتى نتجنبه؟
- 12. الخلاصة والتوجيهات المنهجية للباحثين والمحللين
- المراجع الأكاديمية (References)
1. مدخل نظري: مفهوم المقاومة في الإحصاء الرياضي والوصفي
1.1 التعريف الأكاديمي لمفهوم الإحصائية المقاومة (Resistant Statistic)
تُعرّف الإحصائية المقاومة في الأدبيات الإحصائية المتقدمة بأنها المقياس الكمي أو التابع الحسابي الذي لا يطرأ عليه تغير جوهري أو تشوه حاد عند تعديل أو إفساد نسبة محدودة من قيم العينة المدروسة، مهما بلغت درجة تطرف تلك القيم المضافة أو المعدلة نحو اللانهاية الموجبة أو السالبة. ترتبط هذه الخاصية البنيوية ارتباطاً وثيقاً بقدرة المقياس على عزل الضجيج العشوائي والأخطاء المنهجية التي تصيب أطراف التوزيع التكراري، مما يجعله معبراً بصدق عن الكتلة الغالبة أو المركز الحقيقي للبيانات.
ولقياس هذه الحصانة رياضياً، طوّر علماء الإحصاء الاستدلالي مفهوماً تأسيسياً يُعرف باسم نقطة الانهيار (Breakdown Point). تُمثل نقطة الانهيار الحد الأدنى من النسبة المئوية للبيانات الملوثة أو الفاسدة التي يمكن للعينة أن تتحملها قبل أن يُعطي المقدّر الإحصائي قيماً اعتباطية لا معنى لها، أو ينحرف كلياً نحو اللانهاية ($\pm\infty$). فإذا كانت نقطة انهيار مقياس ما تقارب الصفر ($0%$)، كما هي الحال في المتوسط الحسابي الكلاسيكي، فهذا يعني أن قيمة متطرفة واحدة في عينة مكونة من آلاف المفردات قادرة على سحب المقياس وتدمير صلاحيته التفسيرية بالكامل. في المقابل، تتمتع بعض المقاييس المقاومة بنقطة انهيار تصل إلى ($50%$)، وهي أعلى قيمة نظرية ممكنة، مما يعني أن المقياس يظل مستقراً حتى لو تلوّث نصف حجم العينة بالكامل.
من الضروري التفريق من الناحية الاصطلاحية والدلالية بين مفهوم المقاومة (Resistance) ومفهوم القوة الإحصائية أو الحصانة (Robustness). فبينما يُشير مصطلح المقاومة إلى خاصية إجرائية وصفية تعبر عن عدم حساسية المقياس للقيم المتطرفة داخل عينة معينة بصرف النظر عن النموذج التوزيعي، فإن القوة الإحصائية تمثل مفهوماً استدلالياً أوسع صاغه رواد الإحصاء النظري، يعبر عن قدرة الاختبارات والتقديرات الاستدلالية على العمل بكفاءة وثبات حتى في حال انتهاك الافتراضات التوزيعية الصارمة، مثل افتراض التوزيع الطبيعي التام (Normality) أو تجانس التباين (Homoscedasticity).
يتم تقييم مناعة الإحصائية ضد عينات البيانات الملوثة أو غير الطبيعية عبر فحص دوال التأثير (Influence Functions)، وهي أدوات رياضية تقيس معدل التغير اللحظي في قيمة التقدير الإحصائي عند إضافة نقطة بيانات جديدة ذات قيمة شاذة. تتميز الإحصاءات المقاومة بامتلاكها دالة تأثير محدودة (Bounded Influence Function)، مما يضمن عدم تمكن أي قيمة فردية، مهما تعاظمت، من فرض وزن لا نهائي على النتيجة النهائية للتحليل الإحصائي.
1.2 الجذور التاريخية وتطور الإحصاء الحصين (Robust Statistics)
تعود الجذور المعرفية لظهور الإحصاء الحصين والمقاوم إلى منتصف القرن العشرين، كرد فعل منهجي صارم على الهيمنة المطلقة للإحصاء البارامتري التقليدي القائم على افتراضات غاوس (Gauss) وتوزيع لابلاس-غاوس الطبيعي. لقد أدرك كبار علماء الإحصاء أن البيانات المأخوذة من العالم الواقعي—سواء في العلوم الفيزيائية أو الفلكية أو الإنسانية—تحتوي دوماً على توزيعات ذات ذيول ثقيلة (Heavy-tailed Distributions) ونسب تلوث غير قابلة للتفادي، مما يجعل المقاييس المستندة إلى المربعات الصغرى غير واقعية بل ومضللة.
برزت المساهمة التاريخية الأبرز من خلال أعمال عالم الإحصاء الأمريكي جون توكي (John Tukey)، الذي قاد ثورة منهجية في ستينيات وسبعينيات القرن الماضي عبر تأسيس ما يُعرف بـ التحليل الاستكشافي للبيانات (Exploratory Data Analysis – EDA). أكد توكي على ضرورة “ترك البيانات تتحدث عن نفسها” قبل فرض افتراضات مسبقة صارمة، وطور أدوات ومقاييس مقاومة رائدة تعتمد على الأقسام والتصنيفات الرتبية بدلاً من المجاميع الحسابية البحتة.
تزامن ذلك مع العمل النظري الرصين للعالم الرياضي بيتر هوبر (Peter Huber)، الذي نشر في عام 1964 ورقته البحثية الفارقة حول “التقدير الحصين لموقع المعلمة”، والتي أرست دعائم التقديرات المتقدمة المقاومة (مثل M-estimators). وقد تلتها إسهامات فرانك هامبل (Frank Hampel) في السبعينيات عبر إدخال مفهوم دوال التأثير ونقطة الانهيار، مما نقل دراسة المقاومة الإحصائية من مجرد ملاحظات تطبيقية معزولة إلى فرع رياضي قائم بذاته يتمتع ببناء استدلالي محكم.
شهدت العقود الأخيرة تحولاً تدريجياً ولكنه حاسم في النماذج القياسية للعلوم السلوكية والطبية الحيوية نحو تبني هذه الأساليب؛ حيث كشفت أزمات القياس وتكرار النتائج في علم النفس والأبحاث الجينية أن الاستخدام الأعمى للمقاييس غير المقاومة كان أحد المسببات الرئيسية للنتائج الزائفة والتفسيرات غير الدقيقة للظواهر المعقدة، مما جعل الحصانة والمقاومة ضرورة ابستمولوجية لا غنى عنها في البحث الأكاديمي الرصين.
1.3 الخصائص البنيوية للإحصاءات المقاومة
تستمد الإحصاءات المقاومة حصانتها الرياضية من بنيتها الخوارزمية وطريقة تعاملها مع فضاء العينة. السمة الهيكلية الأولى والأهم لهذه المقاييس هي اعتمادها الجوهري على الرتب والمواقع الترتيبية (Rank-based and Positional Mapping) للأعداد، بدلاً من الاعتماد الحصري على القيم الكمية المطلقة للمفردات. فعندما يتم ترتيب البيانات تصاعدياً، تصبح القيمة المتطرفة مجرد “مفردة تحتل الرتبة الأخيرة”، ويفقد وزنها الحسابي الضخم تأثيره المباشر على الحسابات، إذ تتساوى قيمتها الفعلية سواء كانت تزيد عن سابقتها بوحدة واحدة أو بمليون وحدة.
تتجلى السمة البنيوية الثانية في الاستقرار النسبي للمؤشر عند تكرار التجارب وسحب العينات العشوائية من مجتمعات تحتوي على أطراف بيانات مشوهة أو مشوبة بالأخطاء. إن التباين العيني للمقاييس المقاومة في البيئات الملوثة يكون أقل بكثير من التباين العيني للمقاييس غير المقاومة، مما يمنح الباحثين ثباتاً إحصائياً يحول دون تقلب نتائج الدراسات تقلبات حادة ناتجة عن صدف عشوائية غير منضبطة.
تتمثل السمة البنيوية الثالثة في القدرة الفائقة على تمثيل الاتجاه العام للمجتمع دون انحياز أعمى نحو القيم المتطرفة. فالإحصائية المقاومة تحمي النموذج الرياضي من “طغيان الأقلية المتطرفة” في البيانات، وتضمن أن يكون المعيار الإحصائي الناتج ملخصاً أميناً لواقع الشريحة العظمى من الأفراد أو الظواهر الخاضعة للقياس، مع الاحتفاظ في الوقت نفسه بإمكانية رصد وتحديد تلك الشواذ بشكل مستقل ودون تلويث المؤشر الوصفي الأساسي.
2. طبيعة القيم المتطرفة والشاذة وتأثيرها على الحسابات الإحصائية
2.1 تصنيف القيم الشاذة (Outliers) وأسباب نشوئها في البحوث النفسية والبيولوجية
تُمثل القيمة الشاذة (Outlier) نقطة بيانات تبتعد ابتعاداً ملحوظاً عن النمط التوزيعي العام لبقية المشاهدات في مجموعة البيانات. من المنظور المنهجي، تنقسم القيم الشاذة إلى فئتين رئيسيتين: القيم الشاذة الناتجة عن أخطاء، والقيم الشاذة الحقيقية الناتجة عن التباينات الفردية الطبيعية والمتطرفة. وتتطلب كل فئة منهما فهماً سياقياً دقيقاً قبل اتخاذ أي قرار إحصائي حيالها.
تشمل أخطاء القياس الشائعة أخطاء الإدخال المادي (مثل كتابة 100 بدلاً من 10)، وأعطال أجهزة القياس والتسجيل الإلكتروني في مختبرات علم النفس العصبي والفسيولوجيا، أو سوء فهم المبحوث لتعليمات المقياس. في المقابل، تظهر القيم المتطرفة الطبيعية كأمر واقع في دراسات السمات النفسية والقدرات المعرفية المتطرفة؛ مثل رصد أفراد ذوي مستويات ذكاء فائقة (Genius IQ) أو أوقات استجابة سريعة جداً لا تمثل عطلاً تجريبياً بل تعكس سمة بيولوجية حقيقية للمشارك.
في البحوث النفسية والسلوكية على وجه الخصوص، تُفرز ظواهر مثل التراخي في الإجابة، أو عدم الصدق، أو الإجابة العشوائية في استبيانات التقرير الذاتي (Self-report Questionnaires) قيماً متطرفة حادة. كما تؤدي العينات الصغيرة—وهي سمة شائعة في التجارب السريرية وتصوير الدماغ الوظيفي (fMRI)—إلى تفاقم هذه الأزمة، حيث تفتقر العينة الصغيرة إلى القوة الرياضية الكافية لامتصاص التغير المفاجئ الذي تُحدثه قيمة شاذة واحدة، مما يؤدي إلى تشويه التوزيع بالكامل بصورة تفوق بكثير ما يحدث في العينات الضخمة.
2.2 الآليات الرياضية لتأثير القيمة المتطرفة على المؤشر غير المقاوم
لفهم الهشاشة الرياضية للمقاييس غير المقاومة، يجب تحليل آلية تأثير القيمة المتطرفة على مركز الثقل الحسابي. في مصفوفة البيانات المعرفة رياضياً بالعناصر ${x_1, x_2, dots, x_n}$، يعمل المتوسط الحسابي كـ نقطة اتزان ميكانيكية تخضع لقانون الروافع، حيث يكون مجموع الانحرافات الجبرية حول المركز مساوياً للصفر دائماً:
$$\sum_{i=1}^{n} (x_i – \bar{x}) = 0$$
عندما تتباعد قيمة واحدة $x_k$ نحو قيمة متطرفة ضخمة، فإنها تُحدث عزم دوران هائل يجبر نقطة الاتزان بأكملها ($\bar{x}$) على الانزياح نحوها للحفاظ على التوازن الصفري، مما يسحب المؤشر بعيداً عن مركز تركّز بقية البيانات.
يتضاعف هذا الخلل الرياضي بصورة أسية في مقاييس التشتت غير المقاومة مثل التباين والانحراف المعياري، حيث تعتمد هذه المقاييس على تربيع الفروق عن المتوسط الحسابي:
$$SS = \sum_{i=1}^{n} (x_i – \bar{x})^2$$
إن عملية التربيع هذه تُعامل الفروق الكبيرة بطريقة غير خطية؛ فالفارق الذي مقداره 10 وحدات يُسهم بـ 100 وحدة في مجموع المربعات، بينما الفارق الذي مقداره 100 وحدة يُسهم بـ 10,000 وحدة! هذا التضخيم التربيعي يجعل قيمة متطرفة واحدة كافية للسيطرة على أكثر من 90% من قيمة التباين الكلي للمجموعة.
تختلف التداعيات الرياضية أيضاً بحسب طبيعة التلوث في ذيول التوزيع؛ فالتلوث المتماثل على كلا الجانبين الأيمن والأيسر قد يُبقي المتوسط الحسابي في مكانه ظاهرياً ولكنه يدمر التباين ويفسد فترات الثقة، في حين أن التلوث أحادي الجانب (Unilateral Contamination)—وهو النمط الأكثر شيوعاً في البيانات السلوكية الحقيقية—يُحدث التواءً حاداً يسحب مركز النزعة المركزية ومقاييس التشتت في آن واحد.
2.3 مخاطر الاعتماد على مؤشرات غير مقاومة في وجود شواذ
يقود الاستخدام غير الواعي للمؤشرات غير المقاومة في وجود القيم الشاذة إلى كوارث استدلالية تؤثر مباشرة على صدق البحث العلمي وموثوقيته. من أبرز هذه المخاطر الانزلاق في الخطأ من النوع الأول (Type I Error)—أي رفض الفرضية الصفرية وهي صحيحة في الواقع وادعاء وجود فروق أو علاقات زائفة—أو الوقوع في الخطأ من النوع الثاني (Type II Error) المتمثل في الفشل في رصد فروق حقيقية وقائمة بالفعل نتيجة تضخم الخطأ المعياري الناتج بدوره عن تضخم التباين بفعل القيم الشاذة.
في حقل القياس السريري والنفسي، قد يؤدي الاعتماد على المتوسطات الحسابية المشوهة إلى تشخيصات خاطئة وتوصيفات غير دقيقة للحالات المرضية؛ إذ قد تظهر مجموعة كاملة من المفحوصين وكأنها تعاني من مستويات سريرية مرتفعة من القلق أو الاكتئاب، بينما الحقيقة أن فردين أو ثلاثة فقط من أفراد العينة قد أجابوا بأقصى درجات التطرف الممكنة، مما سحب متوسط المجموعة بالكامل نحو المنطقة الحرجة سريرياً، وحجب الواقع الحقيقي للأغلبية الساحقة ذات الدرجات المعتدلة.
يمتد هذا التآكل أيضاً إلى النماذج الإحصائية المتقدمة؛ حيث تفقد معاملات الارتباط الخطي مثل معامل بيرسون (Pearson’s r) ونماذج الانحدار الخطي العادي (OLS) قوتها التفسيرية بصورة كارثية. فقيمة شاذة واحدة ذات نقطة رافعة عالية (High Leverage Point) قادرة على تحويل معامل ارتباط حقيقي موجب وقوي إلى ارتباط صفري أو حتى سالب، والعكس صحيح، مما يؤدي إلى بناء نظريات علمية باطلة تستند إلى سراب إحصائي.
3. المقارنة الجوهرية: الإحصاءات المقاومة مقابل الإحصاءات غير المقاومة
3.1 الفروق الرياضية والمعيارية الأساسية
ترتكز المقارنة المنهجية بين المقاييس المقاومة وغير المقاومة على مجموعة من المعايير المعمارية التي تحكم كيفية معالجة كل رقم داخل الفضاء المتجهي للبيانات. يكمن المعيار الأول في حساسية المقاييس للقيم الطرفية؛ حيث تتعامل الإحصاءات غير المقاومة مع كل مفردة كقيمة وزنية متكافئة وحرة في التأثير، في حين تُخضع الإحصاءات المقاومة تلك القيم لضوابط الموضعة والرتب، مما يضع سقفاً رياضياً صارماً للحد الأقصى من التأثير الذي يمكن لأي نقطة أن تمارسه.
يتعلق المعيار الثاني بـ الكفاءة الإحصائية التقاربية (Asymptotic Efficiency). في ظل التوزيع الطبيعي المثالي والخالي تماماً من أي شائبة أو تلوث، تتفوق المقاييس غير المقاومة (مثل المتوسط والتباين) في الوصول إلى الحد الأقصى من الكفاءة الرياضية المعروف بـ حد كرامر-راو (Cramér-Rao Bound). ومع ذلك، بمجرد حدوث أدنى تلوث في التوزيع (بنسبة لا تتجاوز 1% إلى 5%)، تنحدر كفاءة المقاييس التقليدية بصورة دراماتيكية، بينما تظل المقاييس المقاومة محتفظة بكفاءة عالية ومستقرة، وهو ما يجعل المقاييس المقاومة خياراً عملياً أكثر واقعية وأماناً في التطبيقات الميدانية.
يعالج المعيار الثالث طبيعة التعامل مع كل مفردة؛ فالإحصاءات غير المقاومة تدمج كل معلومة كمية دقيقة في الحساب، وهو ما يُعتبر ميزة في النظم المغلقة والبيانات النقية، ولكنه يتحول إلى عبء منهجي خطير في النظم المفتوحة والمعقدة التي تتعرض لتشوهات مستمرة.
3.2 جدول تصنيفي شامل للمقاييس الإحصائية
يوضح الجدول المقارن التالي التصنيف الدقيق لأهم المقاييس الإحصائية المستخدمة في البحوث الميدانية، مقسمة وفق معيار المقاومة، مع بيان نقطة الانهيار والمعادلات المرجعية لكل مقياس:
| الوظيفة الإحصائية | المقياس المقاوم (Resistant) | نقطة الانهيار | المقياس غير المقاوم (Non-Resistant) | نقطة الانهيار |
|---|---|---|---|---|
| النزعة المركزية | الوسيط (Median) المتوسط المشذب (Trimmed Mean) |
50% $\alpha$ (تصل إلى 50%) |
المتوسط الحسابي (Mean) المتوسط الهندسي (Geometric Mean) |
0% (مفردة واحدة تكسره) 0% |
| التشتت والتباين | المدى الربيعي (IQR) الانحراف المطلق عن الوسيط (MAD) |
25% 50% |
الانحراف المعياري (SD) التباين (Variance) المدى العام (Range) |
0% 0% 0% |
| العلاقة والارتباط | ارتباط الرتب لسبيرمان ($rho$) ارتباط كيندال تاو ($tau$) |
مقاوم للقيم الشاذة مقاوم للقيم الشاذة |
معامل ارتباط بيرسون الخطّي ($r$) | 0% |
| الشكل والالتواء | معامل يولي للالتواء الرتبي (Bowley Skewness) | مقاوم للذيول البعيدة | معامل الالتواء لعزوم بيرسون (Moment Skewness) | 0% |
3.3 الآثار المترتبة على الاستدلال والتعميم في البحث العلمي
إن الاختيار بين المقاييس المقاومة وغير المقاومة ليس مجرد تفضيل حسابي روتيني، بل هو قرار منهجي يحمل آثاراً عميقة على عملية الاستدلال العلمي والقدرة على تعميم النتائج من العينة إلى المجتمع الأصلي. في المجتمعات الإحصائية التي تعاني من الالتواء الشديد والتشتت غير المتجانس، يفشل المتوسط الحسابي في التعبير عن نقطة التمركز الحقيقية للمجتمع، مما يجعل أي استدلال مبني عليه بمثابة تعميم خاطئ لنموذج غير ممثل للواقع.
ينعكس هذا الأثر مباشرة على فترات الثقة (Confidence Intervals)؛ فالفترات المبنية على المتوسط الحسابي والانحراف المعياري تتسع اتساعاً مفرطاً وغير مبرر في وجود القيم المتطرفة، مما يقلل من دقة التقدير ويجعل فترات الثقة واسعة جداً وغير ذات فائدة علمية. في المقابل، تُنتج المقاييس المقاومة فترات ثقة حصينة تحافظ على اتساع منطقي وموثوق، مما يزيد من دقة التقدير الإحصائي.
تكتسب هذه الآثار أهمية بالغة في سياق القرارات السريرية وتخطيط السياسات العامة؛ حيث يؤدي الاعتماد على المقاييس غير المقاومة إلى سوء تقدير الموارد الاقتصادية الموجهة لدعم الفئات المحتاجة، أو الخطأ في تقدير المعايير الوطنية لمستويات الذكاء أو اضطرابات الصحة النفسية، مما يُبرز أن تبني الأساليب المقاومة هو ضمانة أخلاقية وعلمية لسلامة اتخاذ القرار.
4. مقاييس النزعة المركزية المقاومة: الوسيط كمعيار قياسي
4.1 الأساس الرياضي والمنهجي لعمل الوسيط (Median)
يُمثل الوسيط الإحصائي (Median) النموذج المعياري الأكثر شهرة وشيوعاً للمقاييس المقاومة للنزعة المركزية. يُعرف الوسيط رياضياً بأنه القيمة التي تقسم مصفوفة البيانات المرتبة تصاعدياً إلى نصفين متساويين تماماً، بحيث يقع 50% من المشاهدات تحت هذه القيمة و50% فوقها، وهو ما يجعله مكافئاً للمئين الخمسين ($P_{50}$) أو الربيع الثاني ($Q_2$).
تتمثل المعادلة المرجعية لتحديد موقع الوسيط في مصفوفة بيانات مرتبة $X = (x_{(1)}, x_{(2)}, dots, x_{(n)})$ على النحو التالي:
$$\text{Median}(X) = \begin{\cases} x_{\left(\frac{n+1}{2}\right)} &a\mp; \text{إذا كان } n \text{ فردياً} \ \frac{x_{\left(\frac{n}{2}\right)} + x_{\left(\frac{n}{2} + 1\right)}}{2} &a\mp; \text{إذا كان } n \text{ زوجياً} \end{\cases}$$
تنبثق مقاومة الوسيط الخارقة من حقيقة أن حسابه يعتمد على موقع الرتبة الحسابية وليس على القيم الفعلية للبيانات الواقعة عند الأطراف. فلو تخيلنا عينة تحتوي على القيمة القصوى $x_{(n)} = 100$، وتم استبدال هذه القيمة بمليار ($1,000,000,000$)، فلن تتغير قيمة الوسيط بمقدار ذرة واحدة؛ لأن رتبة القيمة وموقعها النسبي لم يتغيرا. ولهذا السبب تحديداً، يمتلك الوسيط أعلى نقطة انهيار ممكنة في علم الإحصاء وهي $50%$، مما يجعله محصناً ضد انهيار تقديره حتى لو تم التلاعب بنصف قيم العينة تقريباً.
4.2 مقارنة سلوك الوسيط والمتوسط في التوزيعات الملتوية
في التوزيعات المتماثلة تماماً الخالية من الالتواء، يتطابق الوسيط مع المتوسط الحسابي، ولكن تظهر القوة الحقيقية للوسيط في التوزيعات الملتوية إيجابياً (Right-Skewed Distribution) وتلك الملتوية سلبياً (Left-Skewed Distribution). في الالتواء الإيجابي—حيث يمتد ذيل البيانات الطويل نحو اليمين كما هو الحال في توزيع الثروات والدخول وأزمنة الرجع العقلية—ينجذب المتوسط الحسابي بقوة نحو الذيل، مقدماً انطباعاً خادعاً بارتفاع المركز، بينما يظل الوسيط مستقراً في قلب التجمع الفعلي للبيانات.
تتضح هذه الظاهرة بجلاء في التطبيقات السيكومترية والاقتصادية؛ فعند قياس مستويات القلق العام في المجتمع، نجد أن الغالبية العظمى من الأفراد يحصلون على درجات منخفضة إلى معتدلة، بينما تسجل قلة قليلة جداً درجات متطرفة للغاية في أقصى اليمين. في هذه الحالة، يعكس الوسيط المستوى الحقيقي للمواطن العادي، في حين يُعطي المتوسط الحسابي تقديراً مضخماً يوحي بأن المجتمع بأكمله على حافة الانهيار النفسي!
الأمر نفسه ينطبق تماماً على قياس أزمنة الاستجابة الحركية والإدراكية (Reaction Times)؛ حيث يتسبب شرود ذهني طفيف لثوانٍ معدودة في قفزة ضخمة في زمن استجابة المحاولة، مما يجعل المتوسط الحسابي غير معبر عن السرعة العقلية الحقيقية للفرد، في حين يستخلص الوسيط القدرة المعرفية الأصيلة بكفاءة متناهية.
4.3 المتوسط المشذب (Trimmed Mean) والمتوسط الموين (Winsorized Mean)
إلى جانب الوسيط، ابتكر الإحصائيون حلولاً توفيقية متطورة تدمج بين الكفاءة الحسابية للمتوسط والمقاومة الاستثنائية للوسيط، ويأتي في مقدمتها المتوسط المشذب (Trimmed Mean) والمتوسط الموين (Winsorized Mean). تقوم فكرة المتوسط المشذب على اقتطاع نسبة مئوية محددة مسبقاً (ولنقل $\alpha = 10%$ أو $20%$) من كلا طرفي التوزيع المرتب (أعلى وأدنى القيم)، ثم حساب المتوسط الحسابي الكلاسيكي للقيم المتبقية في المنتصف.
يُحسب المتوسط المشذب بنسبة $\alpha$ عبر المعادلة التالية (حيث $k = [n \cdot \alpha]$ هو عدد المفردات المحذوفة من كل طرف):
$$\bar{x}_{t\alpha} = \frac{1}{n – 2k} \sum_{i=k+1}^{n-k} x_{(i)}$$
أما في المتوسط الموين (Winsorized Mean)، فلا يتم حذف القيم المتطرفة نهائياً، بل يتم استبدالها بأقرب قيمة غير مقطوعة من الداخل؛ مما يحافظ على حجم العينة الكلي $n$ مع تجريد القيم الشاذة من وزنها المفرط. يمثل هذان المقياسان أداة لا غنى عنها للباحثين الذين يرغبون في الاحتفاظ بالخصائص الجبرية للمتوسط مع تحصين نتائجهم ضد مخاطر الذيول الثقيلة والقيم الشاذة.
5. مقاييس النزعة المركزية غير المقاومة: هشاشة المتوسط الحسابي
5.1 البنية الرياضية للمتوسط الحسابي (Arithmetic Mean)
يُعد المتوسط الحسابي ($\bar{x}$) المقياس الأكثر شهرة واستخداماً في الإحصاء الوصفي والاستدلالي عبر التاريخ، إلا أن بنيته الرياضية تحمله نقطة ضعف هيكلية تجعله هشاً وغير مقاوم على الإطلاق. يُعرف المتوسط الحسابي بأنه حاصل قسمة المجموع الكلي للقيم على عددها:
$$\bar{x} = \frac{1}{n} \sum_{i=1}^{n} x_i$$
تكمن الهشاشة البنيوية في أن البسط الحسابي ($\sum x_i$) يمنح كل مفردة فردية وزناً مساوياً تماماً لقيمتها المطلقة في تكوين الناتج النهائي. وهذا يعني أن أي اضطراب أو تضخم لا نهائي في قيمة مفردة واحدة سينعكس مباشرة وبدون أي كابح رياضي على البسط، مما يؤدي إلى انزياح فوري للمتوسط الحسابي.
ترتبط هذه البنية بخاصية المجموع الصفري للانحرافات، والتي تفرض على المتوسط الحسابي أن يكون المركز الجبري الذي تتوازن عنده الفروق السالبة والموجبة تماماً. هذه الخاصية الحتمية تجعل المتوسط مجبراً على التحرك اللحظي لمطاردة أي قيمة شاذة تحاول الإخلال بهذا التوازن، حتى لو كان ذلك على حساب تمثيل المركز الحقيقي للأغلبية الساحقة من مفردات العينة.
5.2 تحليل أثر القيمة المتطرفة الواحدة على المتوسط الحسابي
لتوضيح الأثر التدميري لقيمة شاذة واحدة على المتوسط الحسابي، دعنا نتأمل عينة صغيرة من درجات مقياس الاكتئاب السريري تتكون من خمسة مفحوصين في عينة غير سريرية: $[10, 12, 11, 13, 14]$. المتوسط الحسابي لهذه المجموعة هو $12$، وهو يمثل مركز البيانات بدقة بالغة. ولكن، إذا تعرض المفحوص الخامس لحالة ذهانية حادة أو أخطأ الفاحص وسجل درجته $140$ بدلاً من $14$، تصبح العينة: $[10, 12, 11, 13, 140]$.
يقفز المتوسط الحسابي في هذه اللحظة من $12$ إلى $37.2$! إن القيمة الجديدة ($37.2$) أعلى من أربعة أضعاف قيم 80% من المفحوصين، وتقع خارج نطاق معظم البيانات الفعلية. هذا ما يُعرف في الأدبيات المنهجية بـ الوهم الإحصائي (Statistical Illusion)؛ حيث يُعطي المتوسط الحسابي انطباعاً خادعاً بأن مستوى الاكتئاب العام للعينة يقع في المنطقة الحادة، بينما 80% من العينة في الواقع طبيعيون تماماً ومستقرون.
يوضح هذا المثال البسيط نقطة الانهيار الصفرية للمتوسط الحسابي؛ حيث تمكنت مشاهدة واحدة شاذة في عينة من 5 مفردات ($20%$) من تدمير المعنى الإحصائي للمقياس بالكامل، بل وسحبه إلى نطاق لا يمت بصلة لأي من المشاهدات الأساسية للمجموعة.
5.3 متى يكون استخدام المتوسط الحسابي مبرراً علمياً؟
على الرغم من الهشاشة الرياضية للمتوسط الحسابي، فإن استخدامَه يظل مبرراً علمياً ومرغوباً بقوة تحت شروط منهجية محددة. يبرز الشرط الأول في التحقق الصارم من اعتدالية التوزيع (Normal Distribution) وخلو البيانات من الذيول الثقيلة والملوثات؛ ففي ظل التوزيع الطبيعي النقي، يُعد المتوسط الحسابي أفضل مقدر غير متحيز ذي تباين أدنى (UMVUE)، مما يجعله يتمتع بالكفاءة الرياضية القصوى التي تفوق كفاءة الوسيط.
يتمثل المبرر الثاني في الحاجة إلى الحسابات الجبرية اللاحقة والنمذجة الخطية المتقدمة؛ فالمتوسط الحسابي يتميز بخصائص رياضية فائقة تسمح بجمعه، وطرحه، وتقسيمه بسهولة في نماذج تحليل التباين (ANOVA)، ونماذج الانحدار المتعدد، ونمذجة المعادلات البنائية (SEM)، وهي عمليات رياضية معقدة يصعب إجراؤها بذات السلاسة الجبرية باستخدام الوسيط أو الرتب.
لذا، فإن القاعدة المنهجية الرصينة تنص على أن المتوسط الحسابي هو المعيار الذهبي عندما تكون البيانات نقية ومستوفية للشروط المعلمية، في حين يتحول إلى نقطة ضعف خطيرة يجب استبدالها بمقاييس مقاومة بمجرد ظهور مؤشرات التلوث أو الالتواء الحاد.
6. مقاييس التشتت المقاومة: المدى الربيعي ومقاييس الحصانة المتقدمة
6.1 المدى الربيعي (Interquartile Range – IQR) وطريقة عمله
كما تحتاج مقاييس النزعة المركزية إلى الحصانة ضد القيم الشاذة، فإن مقاييس التشتت تحتاج إلى ذات الدرجة من المقاومة لضمان التقدير الدقيق لمدى تباين البيانات وتشتتها. يُعد المدى الربيعي (Interquartile Range – IQR) المقياس المقاوم الأكثر انتشاراً في الإحصاء الوصفي، ويُعرف رياضياً بأنه الفرق بين الربيع الثالث والربيع الأول:
$$IQR = Q_3 – Q_1$$
يمثل المدى الربيعي الاتساع الذي تشغله نسبة الـ 50% المركزية من البيانات المرتبة، متجاهلاً بصورة تامة الطرفين الأدنى والأعلى (الـ 25% الدنيا والـ 25% العليا). وبفضل هذه البنية، يمتلك المدى الربيعي نقطة انهيار تبلغ $25%$، مما يعني أن أي قيمة متطرفة تقع في أقصى اليمين أو أقصى اليسار لن يكون لها أي تأثير على قيمة المقياس، طالما أنها لا تتجاوز حدود الربع الأدنى أو الربع الأعلى للعينة.
ابتكر جون توكي استخداماً عبقرياً للمدى الربيعي في بناء الأسوار الإحصائية (Tukey’s Fences) للكشف عن القيم الشاذة بموضوعية ودون تأثر بالقيم الشاذة ذاتها، وذلك عبر تحديد الحدود التالية:
- السور الأدنى (Lower Inner Fence): $Q_1 – 1.5 \times IQR$
- السور الأعلى (Upper Inner Fence): $Q_3 + 1.5 \times IQR$
- السور الخارجي المتطرف (Outer Fences): $Q_1 – 3 \times IQR \quad \text{و} \quad Q_3 + 3 \times IQR$
تُعد أي قيمة تقع خارج الأسوار الداخلية قيمة شاذة محتملة، في حين تُعد القيم الواقعة خارج الأسوار الخارجية قيماً متطرفة حادة تتطلب فحوصاً دقيقة وعزلاً عن التحليلات غير المقاومة.
6.2 الانحراف المطلق عن الوسيط (Median Absolute Deviation – MAD)
إذا كان المدى الربيعي يتمتع بنقطة انهيار تبلغ 25%، فإن الانحراف المطلق عن الوسيط (MAD) يُمثل القمة المطلقة في مقاومة مقاييس التشتت بنقطة انهيار استثنائية تصل إلى 50%. يُعرف مقياس MAD رياضياً بأنه وسيط الانحرافات المطلقة لجميع المفردات عن وسيط العينة الأصلي:
$$MAD = \text{Median}\Big(big|x_i – \text{Median}(X)big|\Big)$$
تتميز آلية حساب MAD بالحصانة المزدوجة؛ فالخطوة الأولى تحسب وسيط المجموعة كمرجع مقاوم للمركز، والخطوة الثانية تحسب وسيط الفروق المطلقة بدلاً من حساب متوسطها أو تربيعها، مما يقفل الباب تماماً أمام أي محاولة من القيم الشاذة لتضخيم مقياس التشتت.
ولكي يتمكن الباحثون من مقارنة مقياس MAD بالانحراف المعياري الكلاسيكي في ظل التوزيعات شبه الطبيعية، يُضرب مقياس MAD بمعامل تصحيح رياضي مستمد من دالة التوزيع الطبيعي ($b = \frac{1}{\Phi^{-1}(0.75)} \approx 1.4826$):
$$\widehat{\sigma}_{MAD} = 1.4826 \times MAD$$
يُستخدم هذا التقدير المصحح ($1.4826 \times MAD$) على نطاق واسع في خوارزميات الكشف الآلي عن الملوثات في القياسات السيكومترية المتقدمة، والرصد الإحصائي الآلي للبيانات البيولوجية في أبحاث الجينوم والعلوم العصبية.
6.3 مقارنة مقاومة المدى الربيعي بغيره من مقاييس التشتت
يُظهر المدى الربيعي استقراراً وثباتاً رياضياً متفوقاً عند مقارنته بمقاييس التشتت التقليدية في البيئات الملوثة. فحينما يُضاف عدد من القيم الشاذة إلى أقصى الذيول، يتضاعف الانحراف المعياري بصورة مفرطة كما سنرى في دراسة الحالة لاحقاً، في حين يُحافظ المدى الربيعي على قيمته دون أي اهتزاز يُذكر.
تتضمن عملية استخدام المدى الربيعي مفاضلة منهجية واعية بين الحصانة الإحصائية وفقدان بعض المعلومات التفصيلية في العينات الصغيرة جداً؛ إذ إن تجاهل الـ 50% الطرفية يعني التضحية بالمعلومات الدقيقة عن شكل الذيول. ومع ذلك، في دراسات الفروق الفردية والقياس النفسي حيث تمثل الفروق الدقيقة بين غالبية المفحوصين جوهر الظاهرة النفسية، يُعد المدى الربيعي الأداة الأصدق لوصف التباين الحقيقي بين الأفراد دون تضخيم زائف ناجم عن استجابات غير مسؤولة أو أخطاء عشوائية.
7. مقاييس التشتت غير المقاومة: قابلية تأثر الانحراف المعياري والمدى
7.1 المدى العام (Range) ونقطة الانهيار الصفرية
يُعد المدى العام (Range) أبسط مقاييس التشتت في الإحصاء الوصفي، ويُحسب بطرح أدنى قيمة في العينة من أعلى قيمة فيها:
$$\text{Range} = x_{(n)} – x_{(1)}$$
على الرغم من سهولة حسابه، فإن المدى العام يُعد أضعف مقاييس التشتت وهشاشته مطلقة؛ إذ يمتلك نقطة انهيار صفرية ($0%$) لاعتماده الكلي والحصري على قيمتين اثنتين فقط في فضاء العينة: القيمة القصوى والدنيا، متجاهلاً بنية التوزيع وتشتت الـ $n-2$ مفردة الباقية بينهما.
بمجرد دخول قيمة شاذة واحدة ملوثة إلى مصفوفة البيانات في أقصى اليمين أو أقصى اليسار، يقفز المدى العام قفزة وهمية هائلة لا تعبر بأي شكل عن التشتت الحقيقي للمجموعة، بل تعكس فقط المسافة بين القيمة الشاذة وبقية البيانات، مما يجعله مقياساً عاجزاً عن تقديم أي معلومة وصفية ذات قيمة علمية رصينة في العينات الملوثة أو الواقعية.
7.2 الانحراف المعياري والتباين (Standard Deviation and Variance)
يُمثل التباين ($\sigma^2$ أو $s^2$) والانحراف المعياري ($\sigma$ أو $s$) الركن الأساسي في الإحصاء البارامتري التقليدي، ويُحسب الانحراف المعياري للعينة عبر الجذر التربيعي لمتوسط مربعات الانحرافات عن المتوسط الحسابي:
$$s = \sqrt{\frac{1}{n-1} \sum_{i=1}^{n} (x_i – \bar{x})^2}$$
تنبع الهشاشة الرياضية القاتلة للانحراف المعياري من مصدرين متزامنين: اعتماده المباشر على المتوسط الحسابي غير المقاوم كمرجع مركزي، واعتماده على آلية تربيع الفروق ($(x_i – \bar{x})^2$). فعندما تبتعد قيمة واحدة عن المتوسط بمسافة كبيرة، يتم تربيع هذا البعد، مما يمنح تلك القيمة المفردة القدرة على السيطرة التامة على بسط التباين وسحبه إلى مستويات فلكية.
تؤدي هذه الخاصية إلى انهيار كارثي للانحراف المعياري في وجود التوزيعات ذات الذيول السميكة (Heavy Tails) والقيم المتطرفة. وتنعكس هذه المشكلة مباشرة على حساب مؤشرات حجم الأثر (Effect Size) في البحوث التجريبية والنفسية، مثل معامل كوهين دي (Cohen’s d) الذي يعتمد على قسمة الفارق بين المتوسطين على الانحراف المعياري المشترك:
$$d = \frac{\bar{x}_1 – \bar{x}_2}{s_{\text{pooled}}}$$
عندما تتسبب قيمة شاذة في تضخيم الانحراف المعياري ($s_{\text{pooled}}$)، يؤدي ذلك إلى تصغير حجم الأثر ($d$) بشكل مصطنع وخادع، مما يوحي بعدم فاعلية تدخل علاجي نفسي أو دواء بيولوجي أثبت كفاءة حقيقية لدى معظم المشاركين، وهو ما يهدد بوقف برامج علاجية واعدة بسبب خلل في المقاييس الإحصائية غير المقاومة المستخدمة.
7.3 معامل الاختلاف (Coefficient of Variation) ومحدوديته
يُستخدم معامل الاختلاف (CV) للمقارنة بين تشتت مجموعتين مختلفتين في وحدات القياس أو المتوسطات، ويُحسب بقسمة الانحراف المعياري على المتوسط الحسابي ومضروباً في 100:
$$CV = \left(\frac{s}{\bar{x}}\right) \times 100$$
يعاني معامل الاختلاف من حساسية مضاعفة ومركبة تجاه القيم الشاذة؛ نظراً لأنه يتكون من بسط غير مقاوم (الانحراف المعياري) ومقام غير مقاوم أيضاً (المتوسط الحسابي). إن حدوث أي تشوه في أحد الطرفين أو كليهما يؤدي إلى تقلبات دراماتيكية في قيمة المعامل، مما يجعل مقارنة التباين النسبي بين المجموعات السلوكية أو الطبية أمراً غير موثوق به في ظل غياب الاعتدالية التامة.
للتغلب على هذه المعضلة، يلجأ الإحصائيون إلى البدائل المقاومة لمعامل الاختلاف، مثل معامل التشتت الربيعي المقاوم (Quartile Coefficient of Dispersion – QCD)، الذي يُحسب عبر قسمة المدى الربيعي على مجموع الربيعين الأول والثالث:
$$QCD = \frac{Q_3 – Q_1}{Q_3 + Q_1}$$
يقدم هذا البديل المقاوم تقديراً محصناً للتجانس والتشتت النسبي دون تأثر بالقيم الشاذة التي تفسد المعامل التقليدي.

8. دراسة حالة تطبيقية: المحاكاة الرقمية لأثر القيمة المتطرفة
8.1 تحليل البيانات الأولية (بدون قيمة متطرفة)
لتقديم برهان تجريبي ورياضي ملموس على التباين الجوهري في السلوك الإحصائي بين المقاييس المقاومة وغير المقاومة، سنقوم بتحليل عينة معيارية تمثل درجات 11 مفحوصاً في اختبار لقياس المرونة المعرفية. البيانات مرتبة تصاعدياً على النحو التالي:
$$X_{\text{clean}} = [2,, 5,, 6,, 7,, 8,, 13,, 15,, 18,, 22,, 24,, 29]$$
حجم العينة هنا $n = 11$. سنقوم أولاً بحساب المؤشرات الإحصائية غير المقاومة والمقاومة بدقة متناهية:
- حساب المؤشرات غير المقاومة:
- المتوسط الحسابي ($\bar{x}$):
$$\bar{x} = \frac{2 + 5 + 6 + 7 + 8 + 13 + 15 + 18 + 22 + 24 + 29}{11} = \frac{149}{11} \approx 13.55$$ - التباين العيني ($s^2$):
$$s^2 = \frac{\sum (x_i – \bar{x})^2}{n-1} = \frac{778.73}{10} \approx 77.87$$ - الانحراف المعياري العيني ($s$):
$$s = \sqrt{77.87} \approx 8.82$$ - المدى العام (Range):
$$\text{Range} = 29 – 2 = 27$$
- المتوسط الحسابي ($\bar{x}$):
- حساب المؤشرات المقاومة:
- الوسيط ($\text{Median}$): نظراً لأن $n = 11$ (فردي)، فإن موقع الوسيط هو $\frac{11+1}{2} = 6$. القيمة المقابلة للرتبة السادسة هي:
$$\text{Median} = x_{(6)} = 13$$ - الربيع الأول ($Q_1$): وسيط النصف الأدنى $[2, 5, 6, 7, 8]$ هو:
$$Q_1 = 6$$ - الربيع الثالث ($Q_3$): وسيط النصف الأعلى $[15, 18, 22, 24, 29]$ هو:
$$Q_3 = 22$$ - المدى الربيعي ($IQR$):
$$IQR = Q_3 – Q_1 = 22 – 6 = 16$$ - الانحراف المطلق عن الوسيط ($MAD$): نحسب أولاً الفروق المطلقة عن الوسيط (13):
$$|X – 13| = [11,, 8,, 7,, 6,, 5,, 0,, 2,, 5,, 9,, 11,, 16]$$
برتيب الفروق تصاعدياً: $[0, 2, 5, 5, 6, 7, 8, 9, 11, 11, 16]$، نجد أن وسيطها (الرتبة السادسة) هو:
$$MAD = 7$$
- الوسيط ($\text{Median}$): نظراً لأن $n = 11$ (فردي)، فإن موقع الوسيط هو $\frac{11+1}{2} = 6$. القيمة المقابلة للرتبة السادسة هي:
نلاحظ هنا التوافق والتقارب الشديد بين المقاييس المقاومة وغير المقاومة ($\bar{x} \approx 13.55$ مقابل $\text{Median} = 13$) نظراً للتوزيع شبه المتوازن للبيانات الأصلية.
8.2 إدخال القيمة الشاذة وتحليل الصدمة الإحصائية
سنقوم الآن بمحاكاة حدوث خطأ إدخال أو تسجيل حالة استجابة متطرفة جداً عبر إضافة قيمة شاذة واحدة ($450$) إلى نهاية العينة، لتصبح العينة الجديدة الملوثة مكونة من $n = 12$ مفردة على النحو التالي:
$$X_{\text{contaminated}} = [2,, 5,, 6,, 7,, 8,, 13,, 15,, 18,, 22,, 24,, 29,, 450]$$
دعنا نعيد حساب جميع المقاييس السابقة لنرى التحول الدراماتيكي:
- إعادة حساب المؤشرات غير المقاومة:
- المتوسط الحسابي الجديد ($\bar{x}_{\text{new}}$):
$$\bar{x}_{\text{new}} = \frac{149 + 450}{12} = \frac{599}{12} \approx 49.92$$ - الانحراف المعياري الجديد ($s_{\text{new}}$):
$$s_{\text{new}} \approx 126.27$$ - المدى العام الجديد ($\text{Range}_{\text{new}}$):
$$\text{Range}_{\text{new}} = 450 – 2 = 448$$
- المتوسط الحسابي الجديد ($\bar{x}_{\text{new}}$):
- إعادة حساب المؤشرات المقاومة:
- الوسيط الجديد ($\text{Median}_{\text{new}}$): بما أن $n = 12$ (زوجي)، فالوسيط هو متوسط الرتبتين السادسة والسابعة ($x_{(6)}=13$ و $x_{(7)}=15$):
$$\text{Median}_{\text{new}} = \frac{13 + 15}{2} = 14$$ - الربيع الأول الجديد ($Q_1$): وسيط النصف الأدنى $[2, 5, 6, 7, 8, 13]$ هو متوسط 6 و 7 = $6.5$.
- الربيع الثالث الجديد ($Q_3$): وسيط النصف الأعلى $[15, 18, 22, 24, 29, 450]$ هو متوسط 22 و 24 = $23$.
- المدى الربيعي الجديد ($IQR_{\text{new}}$):
$$IQR_{\text{new}} = 23 – 6.5 = 16.5$$ - الانحراف المطلق عن الوسيط الجديد ($MAD_{\text{new}}$): الفروق المطلقة عن الوسيط (14) مرتبة: $[1, 1, 4, 6, 7, 8, 8, 9, 10, 12, 15, 436]$، ومتوسط الرتبتين السادسة والسابعة هو $\frac{8 + 8}{2} = 8$:
$$MAD_{\text{new}} = 8$$
- الوسيط الجديد ($\text{Median}_{\text{new}}$): بما أن $n = 12$ (زوجي)، فالوسيط هو متوسط الرتبتين السادسة والسابعة ($x_{(6)}=13$ و $x_{(7)}=15$):
8.3 التحليل الاستنتاجي والمقارنة النسبية للتغير
يوضح الجدول الإحصائي التالي المقارنة التحليلية الشاملة لنسب التغير الناتجة عن دخول القيمة المتطرفة الوحيدة:
| المؤشر الإحصائي | القيمة الأصلية النقية | القيمة بعد إدخال الشاذة (450) | الفرق المطلق | نسبة التغير المئوية (%) | طبيعة المقياس |
|---|---|---|---|---|---|
| المتوسط الحسابي (Mean) | 13.55 | 49.92 | +36.37 | +268.4% | غير مقاوم (منهار تماماً) |
| الوسيط (Median) | 13.00 | 14.00 | +1.00 | +7.7% | مقاوم (ثبات استثنائي) |
| الانحراف المعياري (SD) | 8.82 | 126.27 | +117.45 | +1331.6% | غير مقاوم (انفجار تبايني) |
| المدى العام (Range) | 27.00 | 448.00 | +421.00 | +1559.3% | غير مقاوم (مشوه كلياً) |
| المدى الربيعي (IQR) | 16.00 | 16.50 | +0.50 | +3.1% | مقاوم (استقرار تام) |
| الانحراف المطلق (MAD) | 7.00 | 8.00 | +1.00 | +14.3% | مقاوم (محصن رياضياً) |
تكشف هذه النتائج الرقمية بوضوح قاطع البنية المناعية المذهلة للإحصاءات المقاومة؛ فبينما قفز المتوسط الحسابي بنسبة تجاوزت 268%، وانفجر الانحراف المعياري بنسبة صاعقة بلغت 1331%، وتضاعف المدى العام بأكثر من 1559%، حافظ الوسيط على ثباته بنسبة تغير طفيفة جداً (7.7%) وظل المدى الربيعي مستقراً بنسبة تغير لا تتعدى (3.1%).
تثبت هذه التجربة الرقمية أن المفردة الشاذة الواحدة قادرة على تحويل تقرير إحصائي يستند إلى المتوسط والانحراف المعياري إلى وثيقة مضللة علمياً تفقد أي صلة بخصائص العينة الفعلية، في حين يقدم التقرير المستند إلى الوسيط والمدى الربيعي وMAD وصفاً أميناً، ومستقراً، ومعبراً بدقة عن واقع الظاهرة الإنسانية المدروسة.
9. التطبيقات النفسية والسلوكية: لماذا تهم الإحصاءات المقاومة علماء النفس؟
9.1 أبحاث زمن الرجع وسرعة المعالجة العقلية (Reaction Times)
تُمثل تجارب زمن الرجع وسرعة المعالجة الإدراكية (Reaction Times – RT) في علم النفس المعرفي وعلم الأعصاب السلوكي أحد أبرز الميادين التي تفرض حتمية استخدام الإحصاءات المقاومة. تميل بيانات أزمنة الاستجابة بطبيعتها الفسيولوجية إلى اتباع توزيع غاوسي معدل يُعرف بتوزيع إكس-غاوسيان (Ex-Gaussian Distribution) أو التوزيع اللوغاريتمي الطبيعي (Log-normal)، وهو توزيع يتسم بالتواء إيجابي حاد وذيل ممتد نحو اليمين.
في البيئات التجريبية الميدانية، يتعرض المشارك أثناء تنفيذ مئات المحاولات الإدراكية لومضات من السهو اللحظي، أو الرمش بالعين، أو حركات لا إرادية، مما يؤدي إلى توليد محاولات تستغرق 2500 مللي ثانية في سياق مهمة يكون متوسط الاستجابة الطبيعي فيها حوالي 350 مللي ثانية. إن إدراج هذه القيم الشاذة الناتجة عن شرود الانتباه داخل المتوسط الحسابي الكلاسيكي يؤدي إلى تشويه السرعة العصبية الحقيقية للمشارك وتضخيم الفروق بين المجموعات بصورة زائفة.
يتيح استخدام الوسيط الحسابي لكل مشارك والمدى الربيعي (IQR) لعلماء النفس المعرفي استخلاص الكفاءة المعالجاتية الحقيقية للدماغ، وعزل الضجيج العشوائي الناتج عن التشتت اللحظي بدقة متناهية ودون الحاجة إلى الحذف الاعتباطي للبيانات الذي قد ينطوي على تحيز منهجي غير مرغوب فيه.
9.2 المقاييس النفسية السريرية وتشخيص الاضطرابات
في حقل القياس الإكلينيكي وعلم النفس المرضي، تُستخدم مقاييس تقييم اضطرابات المزاج مثل قائمة بيك للاكتئاب (BDI-II) أو مقياس هاملتون للقلق. في الدراسات المسحية غير السريرية التي تستهدف عينات مجتمعية عامة، يحصل معظم الأفراد على درجات منخفضة جداً (صفر إلى 5 درجات)، في حين يُسجل عدد محدود جداً من الأفراد الذين يعانون من نوبات اكتئابية حادة درجات قصوى تتجاوز 45 درجة.
يؤدي الاعتماد على المتوسط الحسابي في هذه الدراسات المسحية إلى رفع المعيار التقديري للمجتمع بالكامل، مما يجعل الباحثين يخلصون إلى أن المجتمع يعاني من “اكتئاب تحت سريري عام”، بينما الواقع يشير إلى وجود عينة سليمة تماماً تحتوي على بؤر مرضية معزولة وشديدة التطرف. هنا تبرز المقاييس المقاومة—وتحديداً المئينيات والوسيط—كأداة تشخيصية لا غنى عنها لبناء المعايير والدرجات المعيارية الوطنية (Psychometric Norms)، وتفادي تشويه الملف النفسي الجماعي بسبب قلة من الحالات الحدية أو الاستجابات المبالغ فيها.
كما تُسهم المقاييس المقاومة في حماية الدراسات السريرية من التضارب في تقييم فاعلية العلاجات النفسية الدوائية والسلوكية؛ فعند مقارنة مجموعة تلقت العلاج المعرفي السلوكي بمجموعة ضابطة، فإن وجود مريض واحد تدهورت حالته لأسباب خارجية غير متعلقة بالتجربة قد يرفع متوسط المجموعة العلاجية ككل، ويوحي بفشل العلاج إذا تم الاعتماد على المتوسط والانحراف المعياري، بينما يكشف الوسيط والمدى الربيعي عن تحسن حقيقي وملموس لدى الأغلبية الساحقة من المرضى.
9.3 الاستبيانات الميدانية ومقاييس ليكرت (Likert Scales)
يدور في الأوساط السيكومترية ومنهجية القياس جدل أكاديمي مستمر حول الطبيعة القياسية للبيانات المستخلصة من مقاييس ليكرت (Likert-type Scales)، وما إذا كانت تمثل بيانات فترية حقيقية (Interval Data) أم أنها بيانات رتبية ترتيبية (Ordinal Data). يحذر علماء القياس الرصين من أن المسافات النفسية بين خيارات مثل “أوافق بشدة” و”أوافق” و”محايد” ليست متساوية رياضياً بالضرورة، مما يجعل تطبيق العمليات الجبرية كالجمع والقسمة لحساب المتوسط الحسابي ممارسة مشوبة بالمخاطر المنهجية.
تتفوق المقاييس المقاومة مثل الوسيط الرتبي والمدى الربيعي في تقديم ملخص صادق لاستجابات الأفراد على مقاييس ليكرت، نظراً لأنها تتعامل مع البيانات وفق طبيعتها الترتيبية الأصيلة دون افتراض تساوي المسافات. بالإضافة إلى ذلك، تُبدي المقاييس المقاومة مناعة عالية ضد ظواهر التحيز السلوكي الشائعة في الاستبيانات، مثل:
- نزعة الاستجابة المركزية (Central Tendency Bias): ميل غالبية المفحوصين لاختيار الإجابات المحايدة لتجنب اتخاذ مواقف حادة.
- نزعة التطرف في الاستجابة (Extreme Response Style): اختيار الأطراف القصوى بصورة منهجية من قبل فئة محدودة من المشاركين نتيجة سمات شخصية تتعلق بالاندفاعية.
يضمن استخدام الوسيط والتحليلات اللامعلمية الحصينة صدق البناء السيكومتري واستقرار النتائج عند مقارنة الاتجاهات النفسية عبر الثقافات والمجموعات الديموغرافية المتعددة.
10. طرق العرض والتمثيل البياني للإحصاءات المقاومة
10.1 مخطط الصندوق والساعدين (Box-and-Whisker Plot)
يُمثل مخطط الصندوق والساعدين (Box-and-Whisker Plot)، الذي ابتكره جون توكي في سبعينيات القرن العشرين، التحفة البصرية والمعمارية الأبرز لتمثيل وتجسيد الإحصاءات المقاومة. يعتمد المخطط في بنيته على ما يُعرف بـ الملخص الخماسي للأعداد (Five-Number Summary)، والذي يتألف حصرياً من مقاييس مقاومة وموضعية:
- القيمة الدنيا الفعلية غير الشاذة (Minimum)، وتحدد نهاية الساعد السفلي.
- الربيع الأول ($Q_1$)، ويمثل الحافة السفلية للصندوق (مستوى 25%).
- الوسيط ($Q_2$)، ويمثل الخط الأفقي البارز داخل الصندوق (مستوى 50%).
- الربيع الثالث ($Q_3$)، ويمثل الحافة العلوية للصندوق (مستوى 75%).
- القيمة القصوى الفعلية غير الشاذة (Maximum)، وتحدد نهاية الساعد العلوي.
يتيح مخطط الصندوق للباحث إجراء كشف بصري فوري ودقيق عن بنية التوزيع؛ فالمسافة الرأسية للصندوق تمثل مباشرة المدى الربيعي ($IQR$)، وموقع خط الوسيط داخل الصندوق يكشف فوراً عن اتجاه الالتواء الداخلي، في حين تُرسم القيم الشاذة والمتطرفة كنقاط أو نجوم معزولة تقع خارج حدود الساعدين، مما يسمح بتشخيصها بصرياً دون أن تشوه الإطار الهيكلي الأساسي للرسم.
عند مقارنة مجموعات تجريبية أو سريرية متعددة، يُتيح رصف مخططات الصندوق جنباً إلى جنب للباحثين والمشرفين السريريين رؤية الفروق في الوسائط ومستويات التشتت وتداخل المديات الربيعية في لمح البصر وبأعلى درجات الموثوقية والمقاومة.
10.2 مقارنة التمثيل البياني بين التوزيعات الطبيعية والملتوية
يكشف التمثيل البياني المتقدم عن عيوب بصرية جوهرية في المخططات البيانية التقليدية المعتمدة على الأعمدة ذات أشرطة الخطأ (Bar graphs with Error Bars) التي تعرض المتوسط الحسابي مضافاً إليه أو مطروحاً منه الانحراف المعياري ($\bar{x} \pm SD$). في التوزيعات الملتوية بشدة، قد تمتد أشرطة الخطأ الكلاسيكية إلى نطاقات رياضية مستحيلة فيزيائياً وسيكومترياً (مثل أن يمتد شريط الخطأ إلى قيم سالبة في قياس زمن الرجع أو عدد الأخطاء المرتكبة!).
في المقابل، يُبرز التمثيل البياني المقاوم المسافة الفاصلة بين المتوسط والوسيط بيانياً كأداة تشخيصية بصرية للالتواء؛ فكلما تباعد الخط الممثل للمتوسط عن خط الوسيط داخل مخطط التوزيع، دل ذلك على قوة السحب التي تمارسها الذيول المتطرفة على النموذج. ووفقاً لأحدث توصيات جمعية علم النفس الأمريكية (APA) ودليل أسلوب النشر الأكاديمي، يُحث الباحثون بشدة على تجنب مخططات الأعمدة المعتمة واستبدالها بمخططات الصندوق المتقدمة، أو مخططات الكمان (Violin Plots) التي تدمج بين تقدير الكثافة الاحتمالية والملخص الخماسي المقاوم.
10.3 التحليل الاستكشافي للبيانات (Exploratory Data Analysis – EDA)
ترتكز الفلسفة المعرفية لـ التحليل الاستكشافي للبيانات (EDA)، كما صاغها جون توكي، على مبدأ توجيه التحقيق الإحصائي لاستكشاف الأنماط غير المتوقعة والبنى الكامنة في البيانات قبل التسرع في فرض نماذج رياضية بارامترية صارمة. يرى توكي أن الإحصاءات المقاومة والمخططات البصرية المرتبطة بها تمثل “خط الدفاع الأول” للباحث ضد التضليل الرياضي.
يتكامل استخدام الأدوات الوصفية المقاومة (كالوسيط، والمدى الربيعي، والرسوم الصندوقية) مع الفحص البصري لمنحنيات الكثافة لإنتاج فهم معمق وشامل لطبيعة الظاهرة النفسية أو الحيوية. إنها خطوة تضمن تطهير الرؤية التحليلية وتمنح الباحث الثقة في اتخاذ القرارات اللاحقة؛ سواء بالمضي قدماً في التحليلات اللامعلمية المقاومة، أو تطبيق تحويلات رياضية للبيانات، أو عزل الأخطاء المادية قبل الدخول في مرحلة النمذجة الاستدلالية المتقدمة.
11. معايير اتخاذ القرار الإحصائي: متى نختار الإحصاء المقاوم ومتى نتجنبه؟
11.1 مصفوفة اتخاذ القرار وفق طبيعة وتوزيع البيانات
إن قرار الانتقال من الإحصاء الكلاسيكي غير المقاوم إلى الإحصاء المقاوم يتطلب استراتيجية منهجية محكمة تستند إلى فحص الخصائص التوزيعية والسياق البحثي. توضح المصفوفة التالية الخطوات المنهجية المتسلسلة لاتخاذ القرار الصائب:
- الخطوة الأولى: اختبار افتراضات الاعتدالية والتماثل:
- إجراء الاختبارات الرسمية للاعتدالية مثل اختبار شابيرو-ويلك (Shapiro-Wilk Test) في العينات الصغيرة والمتوسطة ($n < 50$)، واختبار كولموجوروف-سميرنوف (Kolmogorov-Smirnov) مع تعديل ليلفور للعينات الكبيرة.
- فحص قيم الالتواء (Skewness) والتفرطح (Kurtosis)؛ فإذا تجاوزت الدرجة المعيارية للالتواء النطاق المسموح ($\pm 1.96$ إلى $\pm 2.0$)، يُعد ذلك مؤشراً حاسماً على ضرورة تفضيل المقاييس المقاومة.
- الخطوة الثانية: تقييم حجم العينة ونسبة التلوث:
- في العينات الصغيرة جداً ($n < 30$)، يكون تأثير القيمة الشاذة مدمراً بشكل استثنائي للمقاييس غير المقاومة، مما يجعل الوسيط والمدى الربيعي الخيار الأكثر أماناً واستقراراً.
- في العينات الضخمة جداً ($n > 5000$)، تتفعل نظرية النهاية المركزية (Central Limit Theorem) بالنسبة لتوزيع المتوسطات، ولكن يظل المتوسط الحسابي الوصفي عرضة للتشوه إذا كان التلوث هيكلياً وممتداً في ذيل واحد، مما يستوجب الجمع بين المقاييس المقاومة وغير المقاومة.
- الخطوة الثالثة: تحديد الهدف البحثي والجمهور المستهدف:
- إذا كان الهدف هو الوصف الواقعي الدقيق لخصائص العينة الحالية ومستوياتها المركزية، فإن الإحصاء المقاوم يتفوق بالمطلق.
- إذا كان الهدف هو النمذجة الرياضية متعددة المستويات والتحليلات الجبرية المعقدة في ظل بيانات نقية ومستوفية للشروط، يظل الإحصاء غير المقاوم خياراً ذا كفاءة عالية.
11.2 الموازنة بين المقاومة والكفاءة الإحصائية (Resistance vs Efficiency)
تواجه الباحثين في الإحصاء الرياضي معضلة مفاضلة كلاسيكية تُعرف بـ المفاضلة بين المقاومة والكفاءة (Resistance vs. Efficiency Trade-off). في ظل التوزيع الطبيعي المثالي والخالي تماماً من التلوث، تبلغ الكفاءة التقاربية النسبية للوسيط مقارنة بالمتوسط الحسابي حوالي $63.7%$ ($\frac{2}{\pi} \approx 0.637$). وهذا يعني أن الباحث الذي يستخدم الوسيط في مجتمع اعتدالي تماماً يحتاج إلى عينة أكبر بنحو 36% للحصول على نفس دقة التقدير التي يوفرها المتوسط الحسابي.
تتمثل تكلفة المقاومة هنا في التضحية بجزء من المعلومات التفصيلية التي تقدمها الأرقام الفردية في البيانات النقية. ومع ذلك، تشير الدراسات الإحصائية المتقدمة إلى أن هذا الفقد النظري في الكفاءة يتلاشى بسرعة هائلة بمجرد حدوث أدنى تلوث في الذيول (أقل من 2%)، حيث تنقلب المعادلة وتصبح المقاييس المقاومة أعلى كفاءة بكثير من المقاييس التقليدية المنهارة.
لتجاوز هذه المفاضلة، يُوصى منهجياً بتبني استراتيجية الجمع المنهجي (Dual Reporting Strategy)؛ حيث يقوم الباحث بتقديم المؤشرات المقاومة (الوسيط وIQR) جنباً إلى جنب مع المؤشرات الكلاسيكية (المتوسط وSD) في التقارير المنشورة، مما يمنح القارئ والباحثين الآخرين رؤية بانورامية كاملة ودقيقة لواقع التوزيع.
11.3 المعالجات البديلة للقيم الشاذة قبل اللجوء للمقاييس المقاومة
قبل حسم خيار الاقتصار على المقاييس المقاومة، تتاح أمام الباحثين مجموعة من البدائل المنهجية لمعالجة اضطرابات التوزيع والقيم الشاذة، ومنها:
- التحويلات الرياضية غير الخطية (Data Transformations): تطبيق التحويل اللوغاريتمي ($log(x)$)، أو تحويل الجذر التربيعي ($\sqrt{x}$)، أو تحويلات عائلة بوكس-كوكس (Box-Cox Transformations) لتقليص الذيول الطويلة وإعادة التوزيع إلى الشكل الاعتدالي، مما يسمح باستخدام المقاييس البارامترية بأمان.
- استبعاد القيم الشاذة وفق معايير مسبقة (Outlier Trimming): حذف المشاهدات التي تتجاوز $\pm 3$ انحرافات معيارية أو تقع خارج أسوار توكي، شريطة أن يتم ذلك بشفافية كاملة وتوثيق منهجي صريح في قسم المنهجية لتجنب الوقوع في ممارسات التلفيق العلمي (p-hacking).
- تحليل الحساسية (Sensitivity Analysis): إجراء التحليل الإحصائي مرتين: الأولى مع تضمين جميع البيانات، والثانية بعد استبعاد أو تشذيب القيم المتطرفة ومقارنة النتائج. إذا ظلت النتائج والاستنتاجات الأساسية ثابتة، فهذا يمنح الباحث ثقة مطلقة في متانة خلاصاته العلمية؛ أما إذا انهارت النتائج عند حذف قيمة واحدة، دل ذلك على هشاشة النموذج وحتمية الاعتماد على الأساليب المقاومة.
12. الخلاصة والتوجيهات المنهجية للباحثين والمحللين
12.1 إعادة صياغة الفهم العام لمفهوم المقاومة الإحصائية
إن المقاومة الإحصائية ليست مجرد أداة فنية بديلة أو حيلة حسابية يلجأ إليها المحلل عند تعثر الطرق التقليدية، بل هي فلسفة قياس علمية رصينة تنبع من احترام الطبيعة المعقدة وغير الكاملة للبيانات الواقعية. إن تبني الإحصاءات المقاومة يُحصّن الباحث من الانزلاق في التفسيرات الوهمية، ويضمن تمثيل الحقيقة الإمبيريقية لأغلبية أفراد العينة دون انحياز للتشوهات العشوائية والضجيج المنهجي.
في ظل الأزمة المعاصرة المعروفة بـ أزمة التكرار والموثوقية (Replication Crisis) في العلوم النفسية والطبية والسلوكية، يبرز الإحصاء المقاوم والحصين كأحد الحلول البنيوية الأساسية لاستعادة الصدق التنبؤي للبحوث العلمية؛ إذ يضمن أن النتائج المنشورة تعكس علاقات وفروقاً جوهرية مستقرة وقابلة للتكرار في مختبرات أخرى، وليست مجرد نواتج عرضية لحالات شاذة سيطرت على المتوسطات الكلاسيكية.
12.2 دليل كتابة النتائج الإحصائية وفق دليل أسلوب APA (الإصدار السابع)
عند توثيق النتائج الإحصائية المقاومة في الأوراق البحثية والرسائل الجامعية وفق دليل الجمعية الأمريكية لعلم النفس (APA 7th Edition)، ينبغي على الباحثين اتباع القواعد والمعايير التالية لضمان الدقة الأكاديمية والشفافية المنهجية:
- التوثيق الرمزي الدقيق:
- استخدام الرمز المائل Mdn للدلالة على الوسيط (Median).
- استخدام الرمز المائل IQR للدلالة على المدى الربيعي (Interquartile Range).
- استخدام الرمز MAD للدلالة على الانحراف المطلق عن الوسيط.
- صيغة كتابة التقارير النصية المقاومة:
“A Mann-Whitney U test indicated that cognitive reaction time was significantly faster in the experimental group (Mdn = 345.00 ms, IQR = 42.50 ms) compared to the control group (Mdn = 412.00 ms, IQR = 68.00 ms), U = 142.00, p = .003, r = .45.”
- الترجمة والتطبيق باللغة العربية:
“أظهر اختبار مان-ويتني أن زمن الرجع المعرفي كان أسرع بدرجة دالة إحصائياً لدى المجموعة التجريبية (الوسيط = 345.00 مللي ثانية، المدى الربيعي = 42.50 مللي ثانية) مقارنة بالمجموعة الضابطة (الوسيط = 412.00 مللي ثانية، المدى الربيعي = 68.00 مللي ثانية)، وترافق ذلك مع حجم أثر مرتفع.”
- التعليل المنهجي في قسم المنهجية: يجب على الباحث شرح المبررات العلمية بوضوح؛ مثل الإشارة إلى فشل اختبار شابيرو-ويلك للاعتدالية ($p < .05$) أو ظهور التواء إيجابي حاد، مما استوجب استخدام الوسيط والمدى الربيعي كبدائل مقاومة وأكثر تمثيلاً للنزعة المركزية والتشتت.
12.3 آفاق مستقبلية في القياس النفسي والإحصاء المتقدم
تشهد علوم البيانات الحديثة والقياس النفسي المتقدم طفرة ملحوظة نحو تطوير وتطبيق خوارزميات التقدير المقاوم؛ حيث تتكامل الطرق المقاومة الكلاسيكية اليوم مع خوارزميات إعادة العينات الحصينة (Robust Bootstrapping) التي تتيح بناء فترات ثقة وتوزيعات عينية دقيقة للوسائط والمديات دون فرض أي شروط توزيعية مسبقة على المجتمع الأصلي.
علاوة على ذلك، أصبحت خوارزميات تعلم الآلة والذكاء الاصطناعي الحصين (Robust Machine Learning) تعتمد كلياً على دوال خسارة مقاومة للقيم الشاذة (مثل Huber Loss بدلاً من المربعات الصغرى Mean Squared Error)، لتجنب تسميم النماذج التنبؤية بالبيانات المشوهة أو الملوثة عمداً. نوصي الباحثين وطلاب الدراسات العليا في علم النفس، والعلوم الاجتماعية، والإحصاء الحيوي بجعل الإحصاء المقاوم جزءاً أصيلاً من حقيبتهم المنهجية اليومية، إعلاءً لقيم الصدق العلمي والدقة المعرفية التي لا غنى عنها في بناء المعرفة الإنسانية الرصينة.
المراجع الأكاديمية (References)
- American Psychological Association. (2020). Publication manual of the American Psychological Association (7th ed.). American Psychological Association. https://doi.org/10.1037/0000165-000
- Field, A. (2018). Discovering statistics using IBM SPSS statistics (5th ed.). SAGE Publications.
- Hampel, F. R., Ronchetti, E. M., Rousseeuw, P. J., & Stahel, W. A. (1986). Robust statistics: The approach based on influence functions. John Wiley & Sons. https://doi.org/10.1002/9781118186435
- Huber, P. J. (1964). Robust estimation of a location parameter. The Annals of Mathematical Statistics, 35(1), 73–101. https://doi.org/10.1214/aoms/1177703732
- Huber, P. J., & Ronchetti, E. M. (2009). Robust statistics (2nd ed.). John Wiley & Sons. https://doi.org/10.1002/9780470434697
- Leys, C., Ley, C., Klein, O., Bernard, P., & Licata, L. (2013). Detecting outliers: Do not use standard deviation around the mean, use absolute deviation around the median. Journal of Experimental Social Psychology, 49(4), 764–766. https://doi.org/10.1016/j.jesp.2013.03.013
- Rousseeuw, P. J., & Croux, C. (1993). Alternatives to the median absolute deviation. Journal of the American Statistical Association, 88(424), 1273–1283. https://doi.org/10.1080/01621459.1993.10476408
- Tukey, J. W. (1977). Exploratory data analysis. Addison-Wesley Publishing Company.
- Wilcox, R. R. (2017). Introduction to robust estimation and hypothesis testing (4th ed.). Academic Press. https://doi.org/10.1016/C2015-0-02206-8
- Wilcox, R. R. (2022). Modern statistics for the social and behavioral sciences: A practical introduction (3rd ed.). Chapman and Hall/CRC. https://doi.org/10.1201/9781003184645