كيفية إجراء اختبار ويلكوكسون للرتب الموقعة في R

تُعد معالجة البيانات التجريبية والمقارنة بين القياسات المقترنة من الركائز الأساسية في منهجية البحث العلمي الكمي، لا سيما في مجالات علم النفس، والعلوم الطبية الحيوية، والتربية، والعلوم الاجتماعية. تاريخياً، اعتمد الباحثون بشكل رئيسي على اختبار t للعينات المرتبطة (Paired Samples t-test) لمقارنة المتوسطات الحسابية قبل التدخل وبعده. بيد أن الواقع التطبيقي غالباً ما يفرض تحديات جمّة تنتهك الافتراضات الصارمة للاختبارات المعلمية؛ كعدم تحقق التوزيع الطبيعي للفروق، ووجود قيم شاذة، أو الاعتماد على مقاييس ترتيبية مثل مقياس ليكرت، الأمر الذي يجعل من اللجوء إلى البدائل اللامعلمية ضرورة منهجية ملحة لضمان سلامة الاستنتاجات العلمية.

يمثل اختبار ويلكوكسون للرتب الموقعة (Wilcoxon Signed-Rank Test)، الذي ابتكره عالم الإحصاء الأمريكي فرانك ويلكوكسون عام 1945، الحل الإحصائي الأمثل والأكثر متانة في هذه السيناريوهات. يقوم هذا الاختبار اللامعلمي على تحويل الفروق الكمية بين الأزواج المرتبطة إلى رتب تصاعدية مع مراعاة اتجاه الفرق وإشارته، مما يمنحه قوة استدلالية عالية وحصانة متقدمة ضد الانحرافات التوزيعية دون التضحية بالقدرة على اكتشاف التأثيرات الحقيقية للدراسة.

تُعتبر لغة البرمجة الإحصائية R المنصة الأكاديمية الرائدة لتطبيق هذا الاختبار؛ إذ توفر ترسانة برمجية متطورة تتدرج من الدوال الأساسية المدمجة وصولاً إلى الحزم المتقدمة المتخصصة في الاستدلال الدقيق وحساب أحجام الأثر والتصورات البيانية عالية الجودة. يقدم هذا الدليل الشامل تفكيكاً نظرياً وتطبيقياً متعمقاً لاختبار ويلكوكسون للرتب الموقعة في R، مستعرضاً الأسس الرياضية، والفروض الإحصائية، والشروط المنهجية، والخطوات البرمجية، بالإضافة إلى معالجة الحالات المعقدة وتوثيق النتائج وفق المعايير الأكاديمية العالمية.

1. مقدمة نظرية حول اختبار ويلكوكسون للرتب الموقعة

1.1 المفهوم الإحصائي والتعريف الرياضي للاختبار

يُعرّف اختبار ويلكوكسون للرتب الموقعة (Wilcoxon Signed-Rank Test) بأنه اختبار إحصائي استدلالي لامعلمي يُستخدم لمقارنة مجموعتين من القياسات المقترنة أو المتكررة المأخوذة من نفس أفراد العينة، أو من عينات متطابقة زوجياً (Matched Pairs). يستند المفهوم الإحصائي للاختبار إلى فحص الفروق الزوجية بين درجتي كل مفحوص في القياسين (القبلي والبعدي مثلاً)، ثم ترتيب القيم المطلقة لهذه الفروق تصاعدياً من الأصغر إلى الأكبر، وإعادة إلحاق الإشارة الأصلية (موجبة أو سالبة) برتبة كل فرق.

من الناحية الرياضية، لنفترض وجود أزواج من الملاحظات، حيث يتم احتساب الفروق لكل زوج. يتم استبعاد الفروق التي تساوي صفراً في الإجراء التقليدي، ليتبقى لدينا حجم عينة فعال. تُرتب القيم المطلقة للفروق وتُعطى رتباً تبدأ من 1 لأصغر فرق وصولاً إلى أقصى فرق. بعد ذلك، يُعرّف إحصائي الاختبار الأساسي W بأنه مجموع الرتب الموجبة، أو مجموع الرتب السالبة. في بيئة البرمجيات الإحصائية كـ R، يُحسب إحصائي الاختبار المسمى V غالباً كمجموع رتب الفروق الموجبة فقط.

عند مقارنة البنية الرياضية لاختبار ويلكوكسون مع اختبار الإشارة البسيط (Sign Test)، نجد أن اختبار الإشارة يكتفي برصد اتجاه الفرق (موجب أم سالب) متجاهلاً تماماً الحجم النسبي للفارق بين القيم، مما يؤدي إلى هدر كبير في المعلومات الإحصائية وانخفاض ملحوظ في القوة الاختبارية. في المقابل، يجمع اختبار ويلكوكسون بين اتجاه الفارق وحجمه النسبي عبر منظومة الترتيب بالرتب، مما يجعله أكثر كفاءة ودقة في اتخاذ القرارات الاستدلالية.

1.2 المقارنة بين الاختبارات المعلمية واللامعلمية

تتمحور الفروق الجوهرية بين الاختبارات المعلمية (Parametric Tests) واللامعلمية (Non-parametric Tests) حول طبيعة الافتراضات التوزيعية للبيانات. تتطلب الاختبارات المعلمية، مثل اختبار t للعينات المترابطة، افتراضات صارمة تشمل التوزيع الطبيعي لفروق الدرجات ومستوى قياس فتري أو نسبي. في المقابل، تُصمم الاختبارات اللامعلمية للتعامل مع البيانات التي لا تتبع توزيعاً محدداً مسبقاً، وتعمل بكفاءة تامة مع المستويات الترتيبية والبيانات الملتوية.

فيما يخص القوة الإحصائية (Statistical Power)، يُظهر اختبار ويلكوكسون تفوقاً واضحاً على اختبار t في العينات الصغيرة ذات التوزيعات غير الطبيعية؛ حيث يفقد اختبار t قدرته على اكتشاف الفروق الحقيقية بسبب تضخم خطأ التباين الناتج عن الانحرافات والالتواءات. أما في حال تحقق شرط التوزيع الطبيعي بدقة، فإن اختبار ويلكوكسون يمتلك ما يُعرف بـ الكفاءة المقاربة النسبية (Asymptotic Relative Efficiency – ARE) التي تبلغ حوالي 95.5% مقارنة باختبار t، مما يعني أنه حتى في أفضل الظروف الملائمة للاختبار المعلمي، لا يخسر اختبار ويلكوكسون سوى نسبة ضئيلة جداً من كفاءته.

علاوة على ذلك، تتميز الاختبارات اللامعلمية بمقاومتها الفائقة لتأثير القيم المتطرفة (Outliers). فالقيم الشاذة التي تؤدي إلى تشويه المتوسط الحسابي والانحراف المعياري في اختبار t يتم تحجيم أثرها جذرياً في اختبار ويلكوكسون بمجرد تحويلها إلى رتبة قصوى ثابتة، مما يحمي الباحث من الوقوع في القرارات الإحصائية المضللة الناتجة عن تباينات فردية متطرفة داخل العينة.

1.3 متى نلجأ إلى اختبار ويلكوكسون كبديل لاختبار t للعينات المرتبطة

يمثل اللجوء إلى اختبار ويلكوكسون للرتب الموقعة خياراً منهجياً حتمياً في العديد من الحالات البحثية الواقعية. الحالة الأولى والأكثر شيوعاً هي انتهاك شرط التوزيع الطبيعي لفروق الأزواج المرتبطة، كما يتضح من خلال الاختبارات الاستدلالية مثل اختبار شابيرو-ويلك أو الفحص البصري للمدرجات التكرارية والرسوم الاحتمالية، حيث يصبح استخدام اختبار t مضللاً إحصائياً.

الحالة الثانية تتجلى عند التعامل مع متغيرات تابعة مقاسة على مقاييس ترتيبية (Ordinal Scales)، مثل استبيانات مقياس ليكرت المكونة من نقاط متعددة (غير موافق بشدة إلى موافق بشدة)، أو مقاييس التقييم السريري، حيث تكون المسافات بين الفئات غير متساوية رياضياً، مما يمنع إجراء العمليات الحسابية للمتوسطات والانحرافات المعيارية المعتمدة في الاختبارات المعلمية.

أما الحالة الثالثة فترتبط بصغر حجم العينة (مثلاً أقل من 15 أو 20 زوجاً من الملاحظات)، حيث تعجز البيانات المحدودة عن تحقيق مبرهنة النهاية المركزية (Central Limit Theorem) التي تسمح بالاعتماد على اعتدالية التوزيع للعينات الكبيرة، لا سيما في ظل وجود التواء حاد (Skewness) أو تفلطح شاذ (Kurtosis) في الفروق المسجلة.

2. الفروض الإحصائية لاختبار ويلكوكسون للرتب الموقعة

2.1 فرضية العدم والفرضية البديلة

تنطلق عملية الاستدلال في اختبار ويلكوكسون للرتب الموقعة من صياغة دقيقة للفرضيات الإحصائية. تنص فرضية العدم (Null Hypothesis – H0) على أن وسيط الفروق بين القياسين المرتبطين في المجتمع الإحصائي يساوي الصفر. بصياغة أخرى متقدمة، تفترض H0 أن التوزيع الاحتمالي للفروق بين الملاحظات المقترنة متماثل تماماً ومتمركز حول نقطة الصفر، مما يعني عدم وجود أي أثر ذي دلالة للمتغير المستقل أو التدخل التجريبي.

في المقابل، تنص الفرضية البديلة (Alternative Hypothesis – H1) على أن وسيط الفروق في المجتمع يختلف عن الصفر، مما يشير إلى أن التوزيع الاحتمالي للفروق قد انزاح إيجاباً أو سلباً عن نقطة التمركز المحايدة. من المهم هنا التأكيد على الفرق الجوهري بين هذا الاختبار واختبار t؛ فاختبار t يختبر فرضية تطابق المتوسطات الحسابية، بينما يختبر ويلكوكسون موقع التوزيع والوسيط التوزيعي للفروق (Distributional Location).

تستند الدلالة النظرية لتطابق التوزيع حول الصفر إلى افتراض أن احتمال أن يكون الفرق موجباً يساوي تماماً احتمال أن يكون سالباً لأي زوج يتم اختياره عشوائياً، مع تساوي الأوزان النسبية للرتب الموجبة والسالبة في ظل صحة فرضية العدم، وهو الأساس الذي تُبنى عليه جداول التوزيع الاحتمالي الدقيق للاختبار.

2.2 صياغة الفرضيات أحادية الطرف وثنائية الطرف

يحدد الباحث اتجاه الفرضية الإحصائية استناداً إلى الإطار النظري للدراسة والأدبيات السابقة. تُستخدم الفرضية ثنائية الطرف (Two-tailed Test) عندما يكون الهدف فحص وجود أي فرق دال إحصائياً بين القياسين دون التنبؤ المسبق باتجاه هذا الفرق، وتُصاغ الفرضية البديلة هنا بأن وسيط الفروق لا يساوي الصفر. يُعد هذا الخيار هو الضبط الافتراضي في معظم البرمجيات الإحصائية لكونه أكثر تحفظاً وأقل عرضة للأخطاء الاستدلالية غير المبررة.

من ناحية أخرى، تُصاغ الفرضية أحادية الطرف باتجاه الزيادة (Right-tailed / Greater) عندما يفترض الباحث نظرياً أن التدخل التجريبي سيؤدي حتماً إلى رفع الدرجات (وسيط الفروق أكبر من الصفر)، كأن نختبر فعالية برنامج تدريبي في زيادة مستوى التحصيل الأكاديمي. وعلى العكس، تُصاغ الفرضية أحادية الطرف باتجاه النقصان (Left-tailed / Less) عندما يكون الهدف إثبات انخفاض الدرجات (وسيط الفروق أصغر من الصفر)، مثل دراسة تأثير علاج نفسي في خفض مستويات القلق والتوتر.

ينعكس اتجاه الفرضية بشكل مباشر على حساب وتفسير القيمة الاحتمالية (p-value)؛ حيث تتركز منطقة الرفض في ذيل واحد من التوزيع الاحتمالي في الاختبارات أحادية الطرف، مما يرفع من حساسية الاختبار لرصد التأثير المتوقع، شريطة أن يكون التوجيه مدعوماً بأسس نظرية راسخة وليس مجرد محاولة لتحقيق دلالة إحصائية زائفة.

2.3 تفسير الدلالة الإحصائية ومستوى المعنوية

يُمثل مستوى المعنوية ألفا (Alpha Level)، والذي يُحدد عادة عند 0.05 أو 0.01 في البحوث الأكاديمية، العتبة الاحتمالية القصوى لقبول ارتكاب الخطأ من النوع الأول (Type I Error)، وهو رفض فرضية العدم بينما هي صحيحة في الواقع. تُقارن القيمة الاحتمالية المحسوبة الناتجة عن اختبار ويلكوكسون بقيمة ألفا المحددة مسبقاً لاتخاذ القرار الإحصائي المنهجي.

إذا كانت القيمة الاحتمالية المحسوبة أصغر من أو تساوي مستوى الدلالة، يتم رفض فرضية العدم وتأييد الفرضية البديلة، مما يثبت وجود فروق ذات دلالة إحصائية بين القياسين القبلي والبعدي لا يمكن عزوها للصدفة العشوائية. أما إذا كانت القيمة أكبر من ألفا، فيعجز الباحث عن رفض فرضية العدم، مما يعني عدم كفاية الأدلة الإحصائية لإثبات وجود تأثير للتدخل التجريبي.

يرتبط هذا القرار ارتباطاً وثيقاً بالخطأ من النوع الثاني (Type II Error – Beta)، وهو الفشل في رصد الفرق الحقيقي عند وجوده فعلياً. تزداد احتمالية هذا الخطأ عند صغر حجم العينة أو استخدام اختبارات منخفضة الحساسية، وهنا تتجلى قوة اختبار ويلكوكسون في قدرته على الحفاظ على معدل متوازن للقوة الإحصائية (1 – Beta) مقارنة بالاختبارات اللامعلمية الأضعف كاختبار الإشارة.

3. الشروط والمتطلبات الأساسية لتطبيق الاختبار

3.1 طبيعة البيانات ومستوى القياس

يتطلب تطبيق اختبار ويلكوكسون للرتب الموقعة استيفاء معايير محددة تتعلق بمستوى قياس البيانات. يشترط الاختبار أن يكون المتغير التابع مقاساً على الأقل على مستوى ترتيبي (Ordinal Scale)، أو مستوى كمي فتري (Interval) أو نسبي (Ratio). تكمن أهمية هذا الشرط في ضرورة قدرة الباحث على المفاضلة الترتيبية بين الفروق وحساب اتجاه التغير ومقداره النسبي بين كل زوج من الملاحظات.

تُعد استبيانات مقياس ليكرت (Likert Scales) والاختبارات السيكومترية المقننة، المقاسة بدرجات مجمعة أو فئوية واضحة، من أنسب البيانات لتطبيق هذا الاختبار. تتيح هذه البيانات حساب فارق الدرجات وترتيب تلك الفروق بدقة منطقية. في المقابل، لا يمكن إطلاقاً تطبيق اختبار ويلكوكسون على البيانات الاسمية البحتة (Nominal Data)، مثل تصنيف الحالات إلى (ناجح/راسب) أو (ذكر/أنثى)، حيث يُلجأ في تلك الحالات إلى اختبارات لامعلمية أخرى متخصصة مثل اختبار ماكنيمار (McNemar’s Test) للبيانات الاسمية المزدوجة.

بالإضافة إلى ذلك، يجب أن تكون درجات القياسين القبلي والبعدي مقاسة بنفس وحدة القياس وعلى نفس المدى الرقمي، لضمان أن الفروق المحسوبة تعكس تغيراً حقيقياً في السمة المقاسة وليست ناتجة عن تباين في معايير القياس وأدواته.

3.2 استقلالية الأزواج واقتران الملاحظات

يرتكز التصميم التجريبي لاختبار ويلكوكسون على مفهوم الاقتران التام (Pairing) بين الملاحظات؛ إذ يجب أن تنتمي كل نقطة بيانات في القياس الأول إلى نفس وحدة الملاحظة (نفس الفرد، أو نفس الحيوان المخبري، أو نفس العينة النسيجية) في القياس الثاني، وهو ما يُعرف بتصميم القياسات المتكررة (Repeated Measures Design / Pre-test Post-test).

كما يغطي الاختبار تصميم الأزواج المتطابقة (Matched-Pairs Design)، وفيه يتم مطابقة مفحوص من المجموعة التجريبية مع مفحوص آخر من المجموعة الضابطة بناءً على متغيرات دخيلة متطابقة (كالسن، ومستوى الذكاء، والحالة الاقتصادية)، ليتم التعامل معهما إحصائياً كزوج واحد مترابط من البيانات.

مع تحقق الاقتران الداخلي للزوج، يشترط الاختبار الاستقلالية التامة (Independence) بين أزواج الملاحظات المختلفة؛ بمعنى أن استجابة الزوج الأول يجب ألا تؤثر بأي شكل من الأشكال على استجابة الزوج الثاني. يؤدي انتهاك هذا الشرط، مثل حدوث تفاعل جماعي بين المفحوصين أثناء التجربة أو وجود ارتباط عنقودي غير مضبوط، إلى تضخيم معدل الخطأ من النوع الأول وإفساد مصداقية النتائج الإحصائية للاختبار.

3.3 شرط تماثل توزيع الفروق حول الوسيط

يُعد شرط تماثل توزيع الفروق (Symmetry of Differences) من أدق الشروط المنهجية التي يغفل عنها كثير من الباحثين عند استخدام اختبار ويلكوكسون للرتب الموقعة. إذا كان الهدف من الاختبار هو الاستدلال المباشر على وسيط الفروق (Median Difference)، فيجب أن يكون التوزيع الاحتمالي لفروق الدرجات متماثلاً تقريباً حول وسيطه، حتى وإن لم يكن التوزيع طبيعياً.

يعني التماثل أن شكل ذيل التوزيع على يمين الوسيط يشابه بصرياً ورياضياً شكل الذيل على يساره. إذا تحقق شرط التماثل، فإن وسيط الفروق يتطابق تماماً مع الوسط الحسابي لمركبات والش (Walsh Averages)، مما يتيح تفسير نتائج الاختبار بدقة كاختبار للوسيط. أما إذا كان توزيع الفروق ملتوياً بشدة وغير متماثل، فإن اختبار ويلكوكسون يتحول تقنياً من اختبار لوسيط الفروق إلى اختبار لموقع التوزيع الكلي وشكله (Pseudomedian Location Test).

يمكن التحقق من شرط التماثل باستخدام الأساليب الرسومية المتقدمة كالمدرجات التكرارية لفروق الأزواج، ومخططات الكثافة الاحتمالية، أو عبر الاختبارات الإحصائية الاستدلالية المخصصة لفحص التماثل مثل اختبار مياو-غيلستاين-ستيفنسون (Miao-Gel-Gastwirth Symmetry Test) المتوفر في مكتبات R المتخصصة.

4. فحص وتجهيز البيانات واختبار التوزيع الطبيعي في R

4.1 اختبار شابيرو-ويلك لفحص اعتدالية الفروق

تبدأ الخطوة التشخيصية الأولى في بيئة R بحساب متجه الفروق الزوجية بين درجات القياس القبلي والبعدي. من الأخطاء الشائعة بين المبتدئين فحص التوزيع الطبيعي لكل قياس على حدة؛ فالشرط المعلمي لاختبار t المقترن، والمسوغ الأساسي للتحول إلى بديل ويلكوكسون، يرتكز حصراً على التوزيع الاحتمالي لمتجه الفروق المباشرة بين الأزواج.

يتم تطبيق اختبار شابيرو-ويلك (Shapiro-Wilk Test) باستخدام الدالة المدمجة الرائدة في R عبر تمرير متجه الفروق كمعامل مباشر. تُعد دالة شابيرو-ويلك الاختبار الإحصائي الأكثر قوة وموثوقية في تقييم اعتدالية التوزيع، وتعمل بكفاءة ممتازة للعينات التي يتراوح حجمها من 3 إلى 5000 ملاحظة. تُشير القيمة الاحتمالية الناتجة إذا كانت أصغر من 0.05 إلى انحراف دال إحصائياً عن التوزيع الطبيعي، مما يؤكد انتهاك الافتراض المعلمي ويوجه الباحث نحو اعتماد اختبار ويلكوكسون للرتب الموقعة.

تتميز مخرجات الاختبار في R بالوضوح؛ حيث تقدم قيمة الإحصائي W وقيمة p المقابلة له. يُسهم هذا الاختبار الرقمي، إلى جانب الفحص البصري، في بناء حجة منهجية متماسكة تُبرز التزام الباحث بالدقة الإحصائية وتبرر الانتقال إلى الأساليب اللامعلمية في التقارير الأكاديمية.

4.2 التمثيل البياني لتوزيع الفروق

لا يكتمل التقييم التشخيصي للبيانات دون تدعيمه بالاستكشاف البصري المتعمق؛ حيث تتيح لغة R أدوات رسومية متقدمة لفحص سلوك التوزيع وشكله بدقة. يُعد مخطط الاحتمال الطبيعي (Q-Q Plot) الأداة الأكثر حساسية لرصد انحراف الفروق عن التوزيع المعياري، ويتم إنشاؤه عبر دوال الرسم الأساسية برسم النقاط ومقارنتها بخط الانحدار المرجعي المفترض للتوزيع السوي.

إذا كانت نقاط البيانات تلتف بإحكام حول الخط المستقيم المائل بزاوية 45 درجة، فإن التوزيع يُعد قريباً من النمط الطبيعي. أما في حال انحراف النقاط عند الأطراف العليا أو الدنيا على شكل حرف S أو انحناء مقعر، فإن ذلك يعكس وجود التواء واضح أو ذيول ثقيلة (Heavy Tails) تؤكد عدم ملاءمة الاختبارات المعلمية.

بالإضافة إلى ذلك، يُوفر رسم المدرج التكراري (Histogram) المدمج مع منحنى الكثافة الاحتمالية (Density Curve) رؤية متكاملة لمدى تماثل الفروق حول نقطة الصفر ووجود الالتواءات الإيجابية أو السلبية، وهو أمر بالغ الأهمية في مجالات البحوث النفسية والتربوية التي تتسم بياناتها عادة بتراكم الاستجابات نحو نهايات المقاييس (Ceiling and Floor Effects).

4.3 الكشف عن القيم الشاذة والمتطرفة والتعامل معها

تمثل القيم المتطرفة (Outliers) أحد أبرز التحديات المنهجية في تحليل البيانات المترابطة؛ إذ يمكن لفارق شاذ في زوج واحد أن يُغير متوسط المجموعة بالكامل ويزيد من خطأ التقدير المعياري. يُستخدم المخطط الصندوقي (Boxplot) في R كوسيلة أولية سريعة لرصد تلك الملاحظات المنحرفة التي تقع خارج المدى الربيعي للأعلى أو للأسفل.

توفر حزمة rstatix الحديثة في R أدوات برمجية متخصصة ومبسطة، مثل الدوال المخصصة لتحديد الملاحظات المتطرفة بدقة وتصنيفها إلى قيم متطرفة معتدلة أو قيم شاذة حادة (Extreme Outliers) بناءً على مضاعفات المدى الربيعي (1.5 أو 3 أضعاف IQR)، مما يمنح الباحث تقريراً عددياً شاملاً يحدد هوية المفحوصين ذوي الاستجابات الاستثنائية.

من المنظور الأكاديمي، يتميز اختبار ويلكوكسون بمتانته في التعامل مع هذه القيم؛ حيث لا يُشترط حذفها ميكانيكياً طالما أنها تمثل أداءً حقيقياً غير ناتج عن أخطاء إدخال البيانات. يعمل نظام تحويل الفروق إلى رتب على كبح جماح القيمة الرقمية المتطرفة وحصر أثرها في رتبة نهائية محددة، مما يتيح استبقاء العينة كاملة دون تشويه للواقع التجريبي.

5. بيئة العمل في R والحزم الإحصائية المطلوبة

5.1 الدوال المدمجة في الحزمة الأساسية Stats

تحتوي لغة R في بيئتها التأسيسية على حزمة stats المدمجة افتراضياً، والتي تضم دالة wilcox.test المحورية. تتميز هذه الدالة التأسيسية بسرعتها الفائقة، واعتماديتها البرمجية العالية، وعدم حاجتها لتثبيت أي مكتبات خارجية إضافية، مما يجعلها الخيار الكلاسيكي الأكثر انتشاراً في الأوساط الأكاديمية والتعليمية لتحليل القياسات المستقلة والمقترنة على حد سواء.

تتعامل الدالة الأساسية بكفاءة مع متجهات البيانات المباشرة وتتيح تخصيص كافة المعاملات اللازمة مثل نوع الاقتران والاتجاه وتصحيح الاستمرارية. وتُنتج الدالة كائناً إحصائياً مركباً من نوع htest يحتوي على كافة المؤشرات التفصيلية مثل قيمة إحصائي الاختبار، والقيمة الاحتمالية، ومحددات فترات الثقة، مما يسمح للباحث باستخراجها وتضمينها برمجياً في التقارير الإحصائية المتقدمة.

مع ذلك، تفرض الدالة التأسيسية بعض القيود البرمجية في معالجة الحالات الخاصة المتقدمة، مثل عدم قدرتها على حساب القيمة الاحتمالية الدقيقة في حال وجود رتب متكررة (Ties)، واكتفائها بالتحول إلى التقريب المقارب المصحوب برسائل تحذيرية، وهو ما يتم تجاوزه باستخدام الحزم المتخصصة الحديثة.

5.2 استخدام حزمة rstatix للتحليلات المتوافقة مع Tidyverse

تمثل حزمة rstatix نقلة نوعية في منهجية التحليل الإحصائي داخل بيئة R؛ حيث صُممت لتتوافق كلياً مع فلسفة معالجة البيانات الحديثة ومنظومة Tidyverse. تقدم الحزمة دالة بديلة تحمل اسم wilcox_test تتميز بقدرتها على العمل بسلاسة تامة مع المعاملات الأنبوبية (Pipe Operators) والجداول المدمجة (Tibbles).

تُسهل هذه الحزمة عملية هيكلة التحليل، وتتيح إجراء الاختبارات الإحصائية المقسمة حسب مجموعات العوامل المتعددة بأسلوب برمجي مقتضب ومباشر. كما تُنتج جداول مخرجات منظمة ومرتبة كأطر بيانات قياسية تحتوي على مسميات واضحة للمتغيرات، وإحصائي الاختبار، وحجم العينة الفعلي، ومستويات الدلالة المشفرة برمز النجوم لتسهيل إدراجها في العروض الأكاديمية.

بالإضافة إلى ذلك، تتكامل الحزمة بصورة عضوية وتلقائية مع حزم الرسم البياني المتطورة مثل ggpubr، مما يتيح نقل مخرجات الدلالة الإحصائية المحسوبة وإسقاطها مباشرة على الأشكال البيانية دون الحاجة إلى إعادة كتابة الأرقام يدوياً، وهو ما يقلل من احتمالات الخطأ البشري في إعداد المسودات العلمية.

5.3 حزمة coin للحسابات الدقيقة والتقريبية المتقدمة

تُعد حزمة coin (Conditional Inference Procedures in a Permutation Test Framework) المرجع الإحصائي الأكثر رصانة وإحكاماً لتنفيذ الاختبارات اللامعلمية المتقدمة في R. تقدم هذه الحزمة إطاراً رياضياً قائماً على التوزيعات الشرطية واختبارات التباديل الدقيقة (Exact Permutation Tests)، وتوفر دالة متخصصة للقياسات المقترنة تُدعى wilcoxsign_test.

تتفوق حزمة coin تفوقاً كاسحاً على الدوال التقليدية في معالجة البيانات المعقدة؛ إذ تمتلك خوارزميات رياضية قادرة على احتساب القيمة الاحتمالية الدقيقة بنسبة 100% حتى في ظل وجود تكرار كثيف للرتب (Ties) أو وجود فروق صفرية متعددة، متجاوزة بذلك أوجه القصور التحذيرية في الدوال الكلاسيكية.

كما تتيح الحزمة إمكانية تفعيل محاكاة مونت كارلو التقريبية (Monte Carlo Approximations) بآلاف التكرارات العشوائية للعينات الضخمة التي يتعذر حساب تباديلها التوافقية كلياً، وتوفر للمحلل الإحصائي استخراج درجات المعيارية المحولة بدقة متناهية، وهو أمر جوهري لحساب أحجام الأثر القياسية وفق المعايير الإحصائية الحديثة.

6. البناء الهيكلي لدالة wilcox.test في R

6.1 معاملات الدالة الأساسية وتحديد نوع الاختبار

يرتكز الاستخدام السليم لدالة wilcox.test في R على الفهم الدقيق لبنيتها الهيكلية ومعاملاتها التشغيلية. في سياق العينات المقترنة، يتم تمرير القياس القبلي كمتجه x والقياس البعدي كمتجه y، ويجب تفعيل المعامل الجوهري paired وضبط قيمته الإحصائية على القيمة المنطقية TRUE. يُعد هذا الضبط هو الفارق الحاسم الذي يُلزم R بتنفيذ اختبار ويلكوكسون للرتب الموقعة بدلاً من تنفيذ اختبار مان-ويتني للعينات المستقلة.

كما تدعم الدالة صيغة المعادلة الرياضية (Formula Interface) في حال كانت البيانات مهيكلة بالتنسيق الطويل؛ حيث يُكتب التعبير بتمرير اسم متغير الدرجات ثم علامة التوافق متبوعة بمتغير المجموعة الزمنية مع إضافة المعامل المقترن وتحديد كائن إطار البيانات المرجعي. يُعد هذا الأسلوب أكثر تنظيماً ومرونة في المشاريع البحثية التي تتضمن قياسات متعددة ومترابطة.

إن إغفال ضبط معامل الاقتران يمثل خطأ منهجياً جسيماً؛ إذ يعامل البرنامج المجموعتين كعينتين منفصلتين تماماً، مما يهدر الارتباط الداخلي للأزواج ويؤدي إلى احتساب خاطئ لدرجات الحرية والرتب، وبالتالي الحصول على نتائج استدلالية غير دقيقة قد تعصف بسلامة الورقة البحثية بأكملها.

6.2 تحديد اتجاه الفرضية عبر معامل alternative

يُدار اتجاه الاختبار الإحصائي داخل دالة wilcox.test من خلال المعامل البرمجي alternative، والذي يتيح ثلاثة خيارات متباينة تعكس الفرضية البحثية المصاغة مسبقاً. الخيار الافتراضي هو الاختبار ثنائي الطرف والمحدد بالقيمة النصية two.sided، وفيه يتم اختبار فرضية أن وسيط الفروق يختلف عن الصفر في أي من الاتجاهين الإيجابي أو السلبي.

عندما يسعى الباحث لاختبار فرضية اتجاهية صريحة تنص على أن القياس الأول (x) أكبر من القياس الثاني (y)، يتم ضبط المعامل على القيمة النصية greater. في هذا السياق، يقوم الاختبار بجمع رتب الفروق الموجبة الناتجة عن طرح المتجه الثاني من الأول لفحص مدى دلالة تفوقها. أما إذا كانت الفرضية تفترض أن القياس الأول أقل من القياس الثاني، فيتم ضبط المعامل على القيمة النصية less ليركز الاختبار على رتب الفروق السالبة.

من الضروري الانتباه إلى الترتيب الرياضي للمتجهات عند استخدام الفرضيات أحادية الطرف؛ فالتبديل بين موضع المتجه الأول والثاني يقلب اتجاه الإشارات تماماً، مما قد يحول النتيجة من دالة إحصائياً إلى غير دالة، وهو ما يفرض على الباحث التحقق البرمجي الدقيق من مواءمة ترتيب المتغيرات مع المنطق النظري للاختبار.

6.3 ضبط تصحيح الاستمرارية والحساب الدقيق

تتضمن الدالة معاملين تقنيين بالغي الأهمية للتحكم في دقة الحسابات الإحصائية: المعامل الأول هو correct، وهو معامل منطقي افتراضي القيمة (TRUE) يقوم بتطبيق تصحيح ييتس للاستمرارية (Continuity Correction). يهدف هذا التصحيح إلى تقليص الفجوة الرياضية الناتجة عن تقريب التوزيع التكراري المنفصل للرتب بواسطة منحنى التوزيع الطبيعي المتصل للعينات، وذلك بخصم أو إضافة 0.5 إلى قيمة إحصائي الاختبار في التقريب المعياري.

أما المعامل الثاني فهو exact، وهو معامل منطقي يتحكم في طريقة احتساب القيمة الاحتمالية. عند ضبطه على TRUE، يقوم البرنامج بحساب الاحتمالية الدقيقة عبر إحصاء كافة التباديل الرياضية الممكنة لرتب الفروق، وهو الخيار الأمثل والمفضل في العينات الصغيرة التي يقل حجمها عن 50 زوجاً، شريطة خلو البيانات من الرتب المتساوية تماماً.

في العينات الكبيرة (أكبر من 50 زوجاً)، أو عند وجود رتب مكررة، تقوم لغة R تلقائياً بتعطيل الحساب الدقيق والتحول إلى التقريب المقارب الطبيعي (Asymptotic Approximation) نظراً للتعقيد الحسابي الشديد للتباديل الرياضية، ويُنصح الباحثون في هذه الحالة بالاعتماد على خيارات حزمة coin المتخصصة لضمان أقصى درجات الدقة الإحصائية الممكنة.

7. تطبيق عملي خطوة بخطوة: سيناريو تجريبي واقعي

7.1 إعداد وهيكلة البيانات داخل R

لترسيخ المفاهيم النظرية في إطار تطبيقي شامل، نفترض سيناريو تجريبياً واقعياً لدراسة نفسية-إكلينيكية تهدف إلى تقييم فاعلية برنامج علاجي مستند إلى اليقظة الذهنية في خفض مستويات القلق والتوتر لدى عينة مكونة من 15 مريضاً. تم قياس مستوى القلق لكل مريض قبل بدء البرنامج التدريبي (القياس القبلي)، ثم أُعيد القياس لنفس المرضى بعد إتمام البرنامج العلاجي الذي استمر ثمانية أسابيع (القياس البعدي) باستخدام مقياس تقييم نفسي مكون من درجات من 0 إلى 50.

يتم إدخال البيانات في بيئة R وإنشاء إطار بيانات منظم يحتوي على أرقام تعريفية للمفحوصين، ودرجات القياسين القبلي والبعدي، ثم حساب عمود إضافي للفروق الزوجية الصريحة. بعد بناء البيانات بالتنسيق العريض (Wide Format)، يتم تحويلها أيضاً إلى التنسيق الطويل (Long Format) لتسهيل عمليات النمذجة والرسم البياني المتوافق مع معايير R الحديثة.

تُظهر المعاينة الاستكشافية الأولى للبيانات الوصفية أن وسيط درجات القلق قبل البرنامج كان مرتفعاً، مع تسجيل مدى ربيعي متسع، في حين يظهر القياس البعدي انخفاضاً ملحوظاً في وسيط الدرجات مع تقارب في استجابات المفحوصين، مما يعطي مؤشراً أولياً إيجابياً على فاعلية التدخل التجريبي يستلزم التأكيد عبر الاختبار الاستدلالي المباشر.

7.2 تنفيذ كود الاختبار الإحصائي

يتم تنفيذ الاختبار الإحصائي في R باستخدام الدالة التأسيسية عبر استدعاء المتجهين وضبط الاقتران، وتفعيل استخراج فترات الثقة وحساب الوسيط الزائف لتقدير مقدار التأثير الصافي للبرنامج العلاجي. كما يتم كتابة الشيفرة البرمجية الموازية باستخدام حزمة rstatix لتنفيذ الاختبار بصيغة الصيغة الرياضية الحديثة وحفظ النتائج في كائن إحصائي مهيكل.

تتضمن البرمجة استدعاء دالة wilcox_test مع تحديد فرضية ثنائية الطرف أولاً للتأكد من وجود تغير معنوي في أي اتجاه، ثم إعادة تطبيقها بفرضية أحادية الطرف باتجاه الانخفاض (alternative = ‘less’) لاختبار الفرضية البحثية المباشرة المتمثلة في نجاح البرنامج في خفض حدة القلق لدى أفراد العينة.

يتم تخزين مخرجات الاختبار داخل بيئة العمل لاستخلاص القيم الإحصائية الدقيقة برمجياً، مثل إحصائي الاختبار V ومستوى الدلالة، مما يُسهم في أتمتة عمليات استخراج التقارير وضمان تطابق النتائج الحسابية مع الجداول والرسومات البيانية الملحقة بالدراسة.

7.3 استقراء وفك شفرة المخرجات الإحصائية

عند فحص مخرجات دالة wilcox.test في R، يظهر إحصائي الاختبار تحت مسمى V، وتُشير هذه القيمة العددية إلى مجموع الرتب الإيجابية للفروق بعد استبعاد الفروق الصفرية وترتيب القيم المطلقة. تتوافق القيمة المرتفعة لإحصائي V مع رجحان كفة الفروق في اتجاه محدد، وهو ما يعكس الفارق الجوهري بين الأداء القبلي والبعدي.

تُظهر القيمة الاحتمالية المستخرجة مستوى دلالة أصغر بكثير من مستوى ألفا المعياري (0.05)، مما يقودنا منهجياً إلى رفض فرضية العدم بثقة إحصائية تامة، وقبول الفرضية البديلة التي تؤكد وجود انخفاض دال إحصائياً في مستويات القلق والتوتر لدى المرضى بعد خضوعهم للبرنامج العلاجي المقترح.

بالإضافة إلى ذلك، يُقدم تفعيل معامل فترات الثقة في R تقديراً لما يُعرف بـ الوسيط الزائف (Pseudomedian)، وهو وسيط كافة المتوسطات الممكنة بين كل زوجين من الفروق (Walsh Averages). توفر فترة الثقة المقابلة لهذا التقدير مدى رقمياً دقيقاً لحجم التغير المتوقع في المجتمع، مما يمنح الاستنتاج المبدئي عمقاً إحصائياً يتجاوز مجرد إثبات الدلالة الشكلية إلى تحديد مقدار التحسن الإكلينيكي القابل للقياس.

8. معالجة الحالات الخاصة: الرتب المتساوية والفروق الصفرية

8.1 التعامل مع الفروق الصفرية وتأثيرها على حجم العينة الفعال

تنشأ حالة الفروق الصفرية (Zero Differences) عندما تتطابق درجة المفحوص تماماً في القياسين القبلي والبعدي، مما يجعل قيمة الفرق بينهما مساوية للصفر. في المنهجية الأصلية التي وضعها فرانك ويلكوكسون (Wilcoxon’s Standard Method)، يتم استبعاد كافة الأزواج ذات الفروق الصفرية من التحليل تماماً قبل البدء في حساب وترتيب الرتب، وهو ما تتبعه الدالة الافتراضية في R.

يؤدي استبعاد الفروق الصفرية إلى انخفاض حجم العينة الفعال، مما قد يتسبب في فقدان جزء من القوة الإحصائية للاختبار إذا كان عدد الحالات الثابتة كبيراً داخل العينة. علاوة على ذلك، يمثل استبعاد هذه الحالات تجاهلا لمعلومة تجريبية مفادها أن التدخل لم يُحدث أي أثر لدى هؤلاء المفحوصين، وهو ما قد يؤدي إلى تضخيم غير واقعي لحجم التأثير المحسوب.

كبديل منهجي متقدم، اقترح الإحصائي برات (Pratt’s Method) طريقة يتم فيها تضمين الفروق الصفرية في عملية إعطاء الرتب الأولية مع كافة البيانات، ثم تُحذف الرتب المخصصة للأصفار فقط عند حساب مجموع الرتب النهائي، مما يحافظ على الأوزان النسبية الحقيقية للفروق الأخرى ويقدم تقديراً أكثر متانة، وهي الطريقة التي تتيح حزمة coin في R تطبيقها برمجياً بكل سهولة.

8.2 إدارة الرتب المتساوية والتحذيرات الناتجة

تحدث مشكلة الرتب المقيدة أو المتساوية (Tied Ranks) عندما تتساوى القيم المطلقة للفروق بين زوجين أو أكثر من الملاحظات المختلفة داخل العينة. لمعالجة هذا التساوي، يطبق البرنامج خوارزمية متوسط الرتب (Average Rank Method)، حيث يتم جمع الرتب التي كانت ستُعطى لهذه القيم المتساوية لو كانت متباينة، ثم تُقسم بالتساوي على عدد الحالات المقيدة ليحصل كل فرق على متوسط الرتبة المشتركة.

عند وجود هذه الرتب المتساوية، تُظهر دالة wilcox.test الأساسية في R رسالة التحذير الشهيرة: cannot compute exact p-value with ties. يرجع سبب هذا التحذير إلى أن التوزيع الاحتمالي الدقيق لتباديل ويلكوكسون يفترض رياضياً عدم وجود أي تكرار في الرتب، وبالتالي فإن وجود القيود يُبطل الصيغة التوافقية الدقيقة المغلقة، مما يدفع الدالة تلقائياً للتحول إلى التقريب الطبيعي المقارب.

يجب على الباحث عدم الخلط بين هذا التحذير البرمجي ورسائل الخطأ القاتلة؛ فالتحذير يُعلم المستخدم فقط بنوع المعالجة التقريبية المطبقة. ومع ذلك، في الأبحاث التي تتطلب صرامة رياضية متناهية، يُستحسن عدم الاكتفاء بالتقريب المعياري العادي والانتقال إلى الدوال المتخصصة في الاستدلال الدقيق المشروط.

8.3 تطبيق اختبار التباديل الدقيق عبر حزمة coin

توفر حزمة coin الحل الرياضي الأمثل والأكثر تقدماً لمعالجة التعقيدات الناتجة عن الرتب المتساوية والفروق الصفرية؛ إذ تتيح دالة wilcoxsign_test حساب التوزيع التبادلي الدقيق المباشر والمشروط بالرتب المقيدة الفعلية الموجودة في العينة، دون اللجوء إلى افتراضات التقريب المعياري غير المبررة.

يتم ضبط خيارات الدالة لتحديد طريقة إدارة الأصفار بدقة، وتحديد خوارزمية حساب الاحتمالية سواء كانت دقيقة كلياً أو تقريبية عبر محاكاة التباديل العشوائية، مما يمنح الباحث قيمة p دقيقة لا تتأثر بالتشوهات التوزيعية. ويُعتبر استخدام هذه الحزمة المعيار الذهبي المفضل لدى المحكمين في المجلات العلمية ذات معامل التأثير المرتفع.

عند مقارنة مخرجات الدالة الأساسية بدالة حزمة coin، غالباً ما نلاحظ تقارباً في النتائج العامة للدلالة في العينات المتوسطة والكبيرة، إلا أن دالة coin تتميز بتقديم قيمة Z المعيارية المصحوبة بضبط كامل للقيود والتكرارات، وهو المتطلب الرياضي الأساسي الذي يُبنى عليه احتساب حجم الأثر اللامعلمي بصورة موثوقة.

9. حساب حجم الأثر للاختبار في R

9.1 أهمية حجم الأثر في التحليلات اللامعلمية

أصبح التوثيق الحصري للقيمة الاحتمالية (p-value) غير كافٍ في المنهجيات البحثية المعاصرة؛ فالقيمة الاحتمالية تتأثر بشدة بحجم العينة، حيث يمكن لفروق تافهة عملياً أن تكتسب دلالة إحصائية في العينات الكبيرة جداً، والعكس صحيح في العينات الصغيرة. من هذا المنطلق، تشترط جمعية علم النفس الأمريكية (APA) توثيق حجم الأثر (Effect Size) كمعيار يوضح القوة العملية والتطبيقية للتدخل التجريبي بصرف النظر عن حجم العينة.

في التحليلات المعلمية، يُستخدم مؤشر كوهين d لقياس الفارق بين المتوسطات بوحدات الانحراف المعياري، إلا أن هذا المؤشر يفقد صلاحيته المنهجية في ظل البيانات اللامعلمية والالتواءات الشديدة. هنا يبرز معامل الارتباط ثنائي المتسلسلة أو معامل الأثر اللامعلمي r كمعيار مثالي يعكس مقدار الارتباط والتغير في الرتب بين شروط القياس المقترنة.

كما يُطرح في الأدبيات الحديثة مفهوم احتمال التفوق (Probability of Superiority / Common Language Effect Size)، والذي يترجم حجم الأثر إلى نسبة مئوية مباشرة تُعبر عن احتمالية أن يُظهر مفحوص تم اختياره عشوائياً تحسناً حقيقياً في القياس البعدي مقارنة بالقياس القبلي، وهو مؤشر بديهي يُسهل تفسير النتائج للممارسين وصناع القرار في الميدان التطبيقي.

9.2 طرق حساب معامل الارتباط r في R

يُحسب معامل حجم الأثر اللامعلمي r رياضياً عن طريق قسمة القيمة المعيارية المطلقة Z الناتجة عن التقريب المعياري لاختبار ويلكوكسون على الجذر التربيعي للعدد الإجمالي للملاحظات المقترنة، وفق المعادلة التالية المعتمدة أكاديمياً:

حجم الأثر r = القيمة المطلقة لـ Z مقسومة على الجذر التربيعي لإجمالي الملاحظات (N).

تجدر الإشارة إلى وجود تدقيق منهجي هام حول تعريف N في العينات المقترنة: يرى بعض الإحصائيين أن N تمثل عدد الأزواج، بينما يرى الاتجاه الأكثر تحفظاً ودقة أن N تمثل إجمالي عدد القياسات الفردية مجتمعة (أي ضعف عدد الأزواج). لحسن الحظ، تتولى دالة wilcox_effsize المدمجة في حزمة rstatix إجراء هذه الحسابات تلقائياً وبأقصى درجات الضبط المنهجي.

تتيح لغة R أيضاً حساب فترات الثقة لمعامل حجم الأثر باستخدام تقنيات إعادة المعاينة غير المعلمية (Bootstrapping)، حيث يتم توليد آلاف العينات العشوائية المتكررة من البيانات الأصلية لاستخراج فترة ثقة بنسبة 95% لمعامل r، مما يوفر تقديراً دقيقاً لمدى استقرار حجم التأثير وقابليته للتعميم في دراسات لاحقة.

9.3 تفسير معايير حجم الأثر اللامعلمي

يتم تفسير القيمة العددية لمعامل حجم الأثر اللامعلمي r بالاعتماد على معايير كوهين (Cohen’s Benchmarks) الكلاسيكية المعدلة للاختبارات الترتيبية والارتباطية. يُصنف حجم الأثر وفق المستويات المعيارية التالية:

• أثر صغير (Small Effect): عندما تتراوح قيمة r بين 0.10 وأقل من 0.30، مما يشير إلى وجود أثر محدود للتدخل التجريبي يتطلب عينات كبيرة للتحقق من استقراره.
• أثر متوسط (Medium Effect): عندما تتراوح قيمة r بين 0.30 وأقل من 0.50، وهو يمثل تأثيراً ملحوظاً وذا قيمة عملية واضحة في معظم التدخلات النفسية والتربوية.
• أثر كبير (Large Effect): عندما تصل قيمة r إلى 0.50 فما فوق، مما يعكس نجاحاً استثنائياً وقوة تأثيرية بالغة للبرنامج التجريبي على سلوك وأداء المفحوصين.

يجب التأكيد على أن هذه العتبات الرقمية تمثل خطوطاً إرشادية عامة وليست قواعد جامدة؛ إذ يجب دائماً ربط تفسير حجم الأثر بالسياق النوعي للدراسة وطبيعة المشكلة المعالجة. ففي بعض السياقات الإكلينيكية الحرجة، كعلاج الاضطرابات النفسية الحادة، قد يُعتبر الأثر الصغير أو المتوسط ذا أهمية علاجية استثنائية إذا ارتبط بتحسن نوعية حياة المريض وانخفاض معدلات الخطر.

10. التصور البياني المتقدم للنتائج في R

10.1 تصميم المخطط الصندوقي للعينات المقترنة باستخدام ggplot2

يمثل التصور البياني الجذاب والمنضبط إحصائياً أحد أهم عناصر الورقة البحثية الرصينة. تُتيح حزمة ggplot2 في R بناء مخططات صندوقية احترافية تُبرز بوضوح التغير في الوسيط والمدى الربيعي بين القياسين القبلي والبعدي. يتم تصميم الصناديق بتدرجات لونية متناسقة مع تخصيص ثيمات جمالية خالية من الشوائب البصرية (Theme Minimal) لتناسب شروط النشر الأكاديمي المرموق.

لتعزيز شفافية العرض وتجنب إخفاء تفاصيل التوزيع الداخلي، يتم دمج نقاط البيانات الفردية للمفحوصين فوق المخطط الصندوقي باستخدام طبقة النقاط المبعثرة (geom_jitter) أو النقاط المباشرة بنصف شفافية ملائمة، مما يسمح للقارئ والمحكم برؤية التوزيع الخام لكافة المفحوصين إلى جانب المؤشرات الإحصائية التلخيصية.

يتيح هذا النمط من الرسوم إبراز التباين والانكماش في درجات العينة بوضوح؛ حيث يُظهر هبوط الصندوق البعدي وانخفاض وسيطه بالمقارنة مع الصندوق القبلي دليلاً بصرياً قاطعاً يدعم الاستنتاج الاستدلالي لاختبار ويلكوكسون.

10.2 رسم مخطط المسارات المتطابقة Paired Spaghetti Plot

على الرغم من جودة المخطط الصندوقي، إلا أنه يعامل المجموعتين بصرياً كعينات مستقلة ويفشل في إبراز التغير الحاصل على المستوى الفردي لكل مفحوص على حدة. هنا تبرز الأهمية الفائقة لمخطط المسارات المتطابقة (Paired Spaghetti Plot) كأفضل أداة بصرية مخصصة للبيانات المترابطة.

يقوم هذا المخطط بربط نقطة القياس القبلي بنقطة القياس البعدي لكل مريض عبر خط اتجاهي فردي مستقل باستخدام دالة geom_line(aes(group = id)). يتيح هذا التصور الرسومي رصد التباينات الفردية في الاستجابة للتدخل العلاجي؛ حيث تظهر الخطوط المنحدرة لأسفل استجابة المرضى الإيجابية بالعلاج، بينما تكشف الخطوط المستقيمة أو المرتفعة عن الحالات التي لم تستجب للبرنامج أو تراجعت حالتها.

يمكن تطوير هذا المخطط برمجياً في R بتلوين الخطوط الاتجاهية وفق طبيعة التغير (أخضر للتحسن، وأحمر للتراجع، ورمادي للثبات)، مما يحول الشكل البياني إلى لوحة تشخيصية ثرية بالمعلومات السريرية الدقيقة التي تضفي مصداقية استثنائية على نتائج البحث.

10.3 دمج القيم الإحصائية ومؤشرات الدلالة على الرسم

تصل قوة التصور البياني في R إلى ذروتها عند دمج نتائج الاختبار الإحصائي ومؤشرات الدلالة وحجم الأثر مباشرة فوق الرسوم البيانية. توفر حزمة ggpubr دوال متقدمة تتيح إضافة أقواس المقارنة التلقائية (Comparison Brackets) وكتابة قيم إحصائي V وقيمة p المشفرة برمجياً بدقة فائقة.

تتم محاذاة هذه المؤشرات الرياضية على ارتفاعات مدروسة أعلى الصناديق لتجنب التداخل البصري مع نقاط البيانات، مع إمكانية إضافة شروح نصية مخصصة توثق قيمة حجم الأثر r وحجم العينة الفعال، مما يجعل الشكل البياني مستقلاً بذاته وقادراً على إيصال الرسالة العلمية الكاملة دون حاجة القارئ للرجوع المطول للمتن النصي.

في المرحلة النهائية، يتم تصدير هذه المخططات من R باستخدام دالة ggsave بجودة نقطية فائقة (300 إلى 600 نقطة في البوصة DPI) وبصيغ احترافية مثل PDF للرسوم المتجهية أو TIFF عالية الدقة، وهي الصيغ الصارمة التي تشترطها دور النشر العالمية للأوراق البحثية المحكمة.

11. كتابة التقرير وتوثيق النتائج وفق معايير APA

11.1 المعايير الأكاديمية لتوثيق الاختبارات اللامعلمية

تفرض الإصدارة السابعة من دليل النشر العلمي لجمعية علم النفس الأمريكية (APA 7th Edition) قواعد محددة لتوثيق نتائج الاختبارات اللامعلمية؛ إذ لا يجوز الاكتفاء بذكر المتوسطات الحسابية والانحرافات المعيارية عند الاعتماد على اختبار ويلكوكسون، بل يجب تقديم الإحصاءات الوصفية اللامعلمية المناسبة والمتمثلة في الوسيط (Median) والمدى الربيعي (Interquartile Range – IQR) لكل من شروط القياس المقترنة.

يجب أن يتضمن التوثيق الإحصائي الدقيق العناصر التالية مدمجة في سياق النص الأكاديمي: اسم الاختبار المطبق بالكامل، وحجم العينة الفعال بعد ضبط الأصفار، وقيمة إحصائي الاختبار (V أو W) بحروف مائلة، والقيمة الاحتمالية الدقيقة p محسوبة حتى ثلاث خانات عشرية (مع كتابة p < .001 إذا كانت متناهية الصغر بدلاً من كتابة أصفار مطلقة)، بالإضافة إلى توثيق قيمة حجم الأثر r وفترة الثقة الخاصة به إن وُجدت.

كما يُلزم الدليل الباحثين بالإشارة الصريحة إلى المسوغ المنهجي الذي استدعى اللجوء لاختبار ويلكوكسون بدلاً من اختبار t (مثل عدم تحقق التوزيع الطبيعي للفروق أو طبيعة القياس الترتيبي)، وتوضيح الطريقة المتبعة في إدارة الرتب المتساوية والفروق الصفرية لضمان الشفافية وقابلية التكرار العلمي.

11.2 جداول النتائج والمخرجات الأكاديمية

يمثل الجدول الإحصائي المنظم وسيلة مثالية لتلخيص المقارنات الزوجية، لا سيما في الدراسات التي تتضمن متغيرات تابعة متعددة مقاسة قبل وبعد التدخل. يجب أن يلتزم الجدول بتنسيق APA الصارم؛ حيث يخلو من الخطوط العمودية كلياً ويقتصر على ثلاثة خطوط أفقية رئيسية تفصل العناوين، ومحتوى البيانات، وتذييل الجدول.

يشتمل الجدول النموذجي على أعمدة تحدد أسماء المتغيرات، ووسيط ومدى كل قياس (القبلي والبعدي)، وقيمة إحصائي ويلكوكسون V، ومستوى الدلالة الإحصائية المقابل، وقيمة معامل حجم الأثر اللامعلمي r، وتصنيفه التفسيري (صغير، متوسط، كبير). يتم تصدير هذه الجداول برمجياً من R مباشرة بصيغة جداول مايكروسوفت وورد محررة بالكامل عبر حزم متطورة مثل flextable أو sjPlot.

يُختتم الجدول بهامش توضيحي سفلي (Table Note) يشرح الرموز والاختصارات الإحصائية المستخدمة، ويشير إلى الفرضية المطبقة (أحادية أم ثنائية الطرف)، مع توضيح مستويات الدلالة المشفرة برمز النجوم، مما يمنح الجدول استقلالية وقابلية كاملة للقراءة السريعة.

11.3 نماذج نصية جاهزة للتقارير والأبحاث النفسية

لتسهيل عملية الصياغة الأكاديمية على الباحثين، نورد فيما يلي نماذج نصية قياسية باللغة العربية مطابقة لمتطلبات APA في مختلف سيناريوهات القرارات الإحصائية:

نموذج التقرير في حال وجود فروق دالة إحصائياً (فرضية موجهة):
“أُجري اختبار ويلكوكسون للرتب الموقعة لتقييم أثر البرنامج العلاجي المستند إلى اليقظة الذهنية على خفض مستويات القلق لدى أفراد العينة (N = 15). كشفت النتائج عن انخفاض دال إحصائياً في مستويات القلق بعد التدخل العلاجي (الوسيط = 18.00، المدى الربيعي = 6.50) مقارنة بمستويات القلق المسجلة قبل التدخل (الوسيط = 32.00، المدى الربيعي = 8.00)، حيث بلغت قيمة إحصائي الاختبار V = 120.00، والقيمة الاحتمالية p < .001 (أحادي الطرف). كما أظهر التحليل حجم أثر لامعلمي كبير جداً (r = 0.88)، مما يؤكد الفاعلية العملية والتطبيقية العالية للبرنامج العلاجي في خفض حدة الأعراض.”

نموذج التقرير في حال عدم وجود فروق دالة إحصائياً (فرضية ثنائية الطرف):
“أظهرت نتائج اختبار ويلكوكسون للرتب الموقعة عدم وجود فروق ذات دلالة إحصائية بين درجات الرضا الوظيفي المقاسة قبل البرنامج التدريبي (الوسيط = 45.00، المدى الربيعي = 10.00) والدرجات المقاسة بعده (الوسيط = 47.00، المدى الربيعي = 9.00)، V = 48.50، p = .421 (ثنائي الطرف)، بحجم أثر صغير بلغ r = 0.15. تشير هذه النتيجة إلى أن التدريب الإداري لم يُحدث تغيراً جوهرياً في مستوى الرضا العام لدى الموظفين.”

12. الأخطاء الشائعة واستراتيجيات استكشاف الأخطاء وإصلاحها

12.1 الأخطاء المنهجية في اختيار وتطبيق الاختبار

يقع بعض الباحثين في منزلقات منهجية تؤثر سلباً على سلامة التحليل الإحصائي؛ ومن أبرز هذه الأخطاء الخلط بين اختبار ويلكوكسون للرتب الموقعة للعينات المقترنة واختبار مان-ويتني (المعروف أحياناً باسم اختبار ويلكوكسون لمجموع الرتب Wilcoxon Rank-Sum) المخصص للعينات المستقلة. يؤدي تطبيق الاختبار المستقل على بيانات مقترنة إلى تجاهل الاقتران الداخلي للأفراد وفقدان جزء جوهري من حساسية التحليل.

الخطأ المنهجي الثاني هو إغفال فحص شرط تماثل توزيع الفروق حول الوسيط والاستمرار في تقديم النتائج على أنها فحص مباشر للوسيط الحسابي، في حين أن الالتواء الحاد يجعل الاختبار معبراً عن موقع التوزيع العام (Pseudomedian) وليس وسيط الفروق البسيط، وهو ما يتطلب تدقيقاً في المصطلحات المستخدمة في قسم مناقشة النتائج.

كما يُعد تطبيق الاختبار على عينات ذات أحجام ضخمة جداً (آلاف الملاحظات) دون فحص لحجم الأثر من الأخطاء التفسيرية الشائعة؛ فمع العينات الهائلة، قد تُظهر الفروق البسيطة جداً التي لا تتجاوز كسوراً عشرية غير ذات قيمة سريرية دلالة إحصائية فائقة (p < .001)، مما يفرض على الباحث الحذر التام وعدم الانجراف وراء الدلالة الرقمية المجردة دون تقييم الأثر الميداني العملي.

12.2 الأخطاء البرمجية الشائعة في لغة R وحلولها

أثناء كتابة وتنفيذ الشيفرات البرمجية في R، قد تظهر رسائل خطأ متكررة تتطلب فهماً تقنياً لمعالجتها. من أكثر الأخطاء شيوعاً رسالة: ‘x’ and ‘y’ must have the same length. ينشأ هذا الخطأ عند وجود تباين في عدد عناصر المتجهين القبلي والبعدي، وغالباً ما ينتج عن وجود قيم مفقودة (Missing Values – NA) في أحد القياسات دون الآخر. يتم حل هذه الإشكالية برمجياً بحذف الأزواج غير المكتملة مسبقاً عبر دوال استبعاد القيم الناقصة مثل na.omit لضمان التوافق التام لأطوال المتجهات المقترنة.

خطأ برمجي آخر يرتبط باستيراد البيانات من ملفات خارجية كـ Excel أو SPSS؛ حيث يتم أحياناً استيراد الأعمدة الرقمية كمتغيرات نصية (Characters) أو عوامل اسمية (Factors). يؤدي تمرير هذه المتغيرات غير الرقمية إلى دالة الاختبار لظهور رسائل خطأ تفيد بعدم ملاءمة نوع البيانات، ويكمن الحل في التحويل البرمجي الصريح للأعمدة إلى الصيغة الرقمية عبر الدالة as.numeric قبل تنفيذ الاختبار.

كذلك، يقع المبرمجون في خطأ إغفال ضبط ترتيب العوامل (Factor Levels) عند استخدام صيغة المعادلة الطويلة (Formula Interface)، مما يجعل R يطرح القياس البعدي من القبلي بشكل معكوس يقلب إشارات الفروق، وهو ما يستوجب إعادة ترتيب مستويات المتغير الفئوي بضبط المستويات مسبقاً لضمان توجيه رياضي منضبط لعمليات الطرح التلقائية داخل الدالة.

12.3 البدائل الإحصائية المتقدمة عند فشل شروط ويلكوكسون

في الحالات التي تفشل فيها البيانات في تلبية الشروط التأسيسية لاختبار ويلكوكسون ذاته، تتوفر ترسانة من البدائل الإحصائية المتقدمة. إذا كان توزيع الفروق غير متماثل على الإطلاق ويعاني من التواء شديد وتراكم شاذ لا يمكن علاجه، فإن اختبار الإشارة (Sign Test) يُعد البديل اللامعلمي الأنسب، حيث يتخلى كلياً عن شرط التماثل ويكتفي بفحص اتجاهات الفروق دون أوزانها.

أما في التصاميم التجريبية المعقدة التي تتضمن قياسات متكررة عبر أكثر من نقطتين زمنيتين مع وجود بيانات مفقودة ومتغيرات وسيطة ودخيلة، فإن النماذج الخطية المعممة المختلطة (Generalized Linear Mixed Models – GLMM) أو النماذج الترتيبية المختلطة (Cumulative Link Mixed Models) تمثل الحل الإحصائي الحديث الأكثر دقة ومرونة في نمذجة البنى المعقدة.

أخيراً، يمكن اللجوء إلى اختبارات إعادة التوزيع الحر والتباديل العشوائية النقية (Permutation Paired Tests)، أو استخدام أسلوب التحويلات الرياضية المتقدمة للبيانات (مثل تحويلات بوكس-كوكس Box-Cox Transformations) لإعادة البيانات إلى حظيرة التوزيع الطبيعي بما يتيح العودة الآمنة لتطبيق اختبار t المعلمي الكلاسيكي بكامل قوته الاستدلالية.

خاتمة شاملة وتوصيات منهجية

يُمثل اختبار ويلكوكسون للرتب الموقعة أداة استدلالية لامعلمية بالغة الأهمية والقوة في الترسانة المنهجية للباحث المعاصر؛ إذ يجمع بين المرونة التوزيعية والحصانة المتقدمة ضد القيم الشاذة، مع الحفاظ على كفاءة وقوة إحصائية تقارب الاختبارات المعلمية الكلاسيكية. توفر لغة البرمجة R بيئة متكاملة تتيح تنفيذ هذا الاختبار بمرونة فائقة، بدءاً من المعالجة التشخيصية الأولية وحتى إنتاج التقارير البصرية والأكاديمية المتقدمة وفق معايير النشر العالمية.

نوصي الباحثين في العلوم السلوكية والطبية بضرورة تبني الفحص التشخيصي المزدوج (الرقمي والبصري) لفروق القياسات المقترنة قبل الاستقرار على نوع الاختبار، والحرص الدائم على توثيق أحجام الأثر اللامعلمية وفترات الثقة إلى جانب القيم الاحتمالية. كما نشدد على ضرورة استخدام الحزم المتخصصة في الاستدلال الدقيق عند التعامل مع الرتب المقيدة والفروق الصفرية، بما يضمن سلامة القرارات العلمية ويعزز من موثوقية وتراكمية المعرفة الأكاديمية.

References

• Conover, W. J. (1999). Practical nonparametric statistics (3rd ed.). John Wiley & Sons. https://www.wiley.com/en-us/Practical+Nonparametric+Statistics%2C+3rd+Edition-p-9780471160687
• Field, A., Miles, J., & Field, Z. (2012). Discovering statistics using R. SAGE Publications. https://www.sagepub.com/en-us/nam/discovering-statistics-using-r/book236067
• Hothorn, T., Hornik, K., van de Wiel, M. A., & Zeileis, A. (2006). A Lego system for conditional inference. The American Statistician, 60(3), 257–263. https://doi.org/10.1198/000313006X118430
• Kassambara, A. (2020). rstatix: Pipe-Friendly Framework for Basic Statistical Tests (R package version 0.7.0). https://CRAN.R-project.org/package=rstatix
• R Core Team. (2023). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/
• Wickham, H. (2016). ggplot2: Elegant graphics for data analysis. Springer-Verlag New York. https://ggplot2.tidyverse.org/
• Wilcoxon, F. (1945). Individual comparisons by ranking methods. Biometrics Bulletin, 1(6), 80–83. https://doi.org/10.2307/3001968