كيفية حساب القيمة الاحتمالية لإحصائية كاي-تربيع في R

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

ومع تطور المنظومات البرمجية مفتوحة المصدر، أصبحت بيئة R الإحصائية المنصة الأكثر مرونة وموثوقية في تنفيذ التحليلات السيكومترية المتقدمة والحسابات الاحتمالية الدقيقة. لا يقتصر التحليل الإحصائي في R على مجرد استخراج قيمة إحصائية الاختبار المحسوبة، بل يمتد إلى استخلاص وتفسير القيمة الاحتمالية (P-Value) التي تمثل حجر الزاوية في اتخاذ القرارات المنهجية وقبول أو رفض الفرضيات الصفرية. إن الفهم الرياضي والمعماري لكيفية حساب هذه القيمة عبر دوال R المخصصة مثل pchisq() ودالة chisq.test() يُكسب الباحث عمقاً إبستيمولوجياً يجنبه الوقوع في المزالق الإحصائية الشائعة ويدعم موثوقية النتائج القابلة لإعادة الإنتاج.

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

1. مقدمة تأسيسية: مفهوم إحصائية كاي-تربيع والقيمة الاحتمالية (P-Value) في التحليل النفسي والإحصائي

1.1 التعريف النظري لاختبار كاي-تربيع واستخداماته المنهجية

يمثل اختبار كاي-تربيع، الذي طوره عالم الإحصاء البريطاني كارل بيرسون في مطلع القرن العشرين، نقلة نوعية في معالجة وتحليل البيانات الفئوية (Categorical Data) والاسمية (Nominal Data). ينتمي هذا الاختبار إلى عائلة الاختبارات اللامعلمية (Non-parametric Tests)، وهي أدوات تحليلية لا تشترط أن تتبع البيانات توزيعاً طبيعياً معلماً أو تجانساً تاماً في التباين المتصل، مما يجعله مثالياً لتحليل الظواهر النفسية والسلوكية المقاسة برتب أو تصنيفات كيفية، مثل الأنماط الشخصية، وفئات التشخيص الإكلينيكي، والتفضيلات السلوكية، والاستجابات الاستبيانية المنفصلة.

في إطار الأبحاث النفسية والتطبيقية، ينقسم استخدام اختبار كاي-تربيع إلى نمطين رئيسيين: الأول هو اختبار جودة المطابقة (Goodness-of-Fit Test)، والذي يُعنى بمقارنة تكرارات عينة واحدة ملاحظة بتوزيع نظري محدد مسبقاً (مثل التأكد مما إذا كانت الإصابة بالاكتئاب تتوزع بالتساوي بين الفصول الأربعة). والنمط الثاني هو اختبار الاستقلالية والتجانس (Test of Independence and Homogeneity)، الذي يبحث في مدى وجود اقتران أو اعتماد متبادل بين متغيرين تصنيفيين في جداول الاقتران (Contingency Tables)، كدراسة العلاقة بين نمط التعلق العاطفي والامتثال لبرامج العلاج النفسي.

تتمحور الغاية المنهجية للاختبار حول تقييم مدى انحراف البيانات الميدانية الملاحظة (Observed Frequencies) عن التكرارات التي كان يُفترض الحصول عليها بالصدفة البحتة (Expected Frequencies) في ظل صحة الفرضية الصفرية ($H_0$). ومن هنا تبرز القيمة الاحتمالية ($p$-value) كأداة كمية حاسمة تتيح للباحث النفسي تقدير مدى معقولية الفرضية الصفرية، مما يمهد الطريق لاتخاذ قرار منهجي صائب يدعم أو يدحض النماذج النظرية المطروحة.

1.2 المفهوم الرياضي للقيمة الاحتمالية (P-Value) في سياق توزيع كاي-تربيع

تُعرّف القيمة الاحتمالية ($p$-value) في الإحصاء الاستدلالي بأنها الاحتمال الرياضي للحصول على قيمة إحصائية مساوية للإحصائية المحسوبة من بيانات العينة أو أكثر تطرفاً منها، بافتراض أن الفرضية الصفرية صحيحة تماماً. وفي إطار فضاء توزيع كاي-تربيع الاحتمالي، تترجم هذه القيمة هندسياً ورياضياً بوصفها المساحة المحصورة تحت منحنى دالة كثافة الاحتمال (Probability Density Function – PDF) الممتدة من قيمة إحصائية الاختبار المحسوبة ($\chi^2$) نحو اللانهاية الموجبة في أقصى الذيل الأيمن للتوزيع.

يرتبط اتخاذ القرار الإحصائي بمقارنة هذه المساحة الاحتمالية بمستوى الدلالة المعياري المحدد مسبقاً والمسمى ألفا ($\alpha$)، والذي يُضبط تقليدياً في العلوم السلوكية والإكلينيكية عند $0.05$ أو $0.01$. إذا كانت القيمة الاحتمالية المحسوبة أصغر من أو تساوي مستوى الدلالة ($p le \alpha$)، فإن ذلك يُعد دليلاً إحصائياً قوياً ضد معقولية الفرضية الصفرية، مما يستوجب رفضها وقبول الفرضية البديلة ($H_1$) التي تشير إلى وجود فروق ذات دلالة إحصائية أو ارتباط حقيقي غير ناتج عن أخطاء المعاينة العشوائية.

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

1.3 مبررات استخدام بيئة R الإحصائية لتحليل البيانات السيكومترية

تُعد لغة R وبيئتها الحسابية الخيار المعياري الذهبي للباحثين والمحللين الإحصائيين في مجالات القياس النفسي (Psychometrics) والعلوم الحيوية والسلوكية. تتميز R عن الحزم الإحصائية التجارية المغلقة بمرونتها البرمجية الفائقة، وشفافيتها المنهجية الكاملة؛ حيث يستطيع الباحث فحص وتعديل الخوارزميات الرياضية الدقيقة الكامنة وراء كل دالة إحصائية، مما يمنع التعامل مع الإحصاء كـ “صندوق أسود”.

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

علاوة على ذلك، تزخر مستودعات R، مثل CRAN، بآلاف الحزم البرمجية المتخصصة في القياس النفسي، مثل حزمة psych وvcd وrcompanion وeffectsize، التي تقدم أدوات متكاملة لحساب حجوم التأثير، وتحليل مصفوفات الاقتران المتعددة، ورسم التوزيعات الاحتمالية بجودة احترافية مخصصة للنشر في الدوريات العلمية المحكمة.

2. البنية الرياضية لتوزيع كاي-تربيع ودور درجات الحرية (Degrees of Freedom)

2.1 خصائص توزيع كاي-تربيع الإحصائي

يُعرف توزيع كاي-تربيع رياضياً بأنه توزيع المتغير العشوائي الناتج عن جمع مربعات عدد $k$ من المتغيرات العشوائية المعيارية المستقلة التي تتبع التوزيع الطبيعي القياسي $N(0, 1)$. وبناءً على هذا التعريف التأسيسي، فإن قيم إحصائية كاي-تربيع تنحصر حصراً في المجال الموجب $[0, \infty)$، مما يجعل وقوع أي قيمة سالبة مستحيلاً رياضياً وإجرائياً.

تعتمد دالة كثافة الاحتمال (PDF) لتوزيع كاي-تربيع اعتماداً كلياً ومباشراً على معلمة وحيدة هي درجات الحرية (Degrees of Freedom – $df$). تتسم هذه الدالة بالالتواء الإيجابي الشديد نحو اليمين عندما تكون درجات الحرية منخفضة ($df = 1$ أو $df = 2$)؛ حيث يتكدس الاحتمال بالقرب من الصفر وتتلاشى المساحة تدريجياً نحو اليمين. ومع تزايد درجات الحرية، يتناقص الالتواء تدريجياً ويقترب شكل التوزيع من التوزيع الطبيعي المتماثل استناداً إلى نظرية النهاية المركزية (Central Limit Theorem).

تُحسب إحصائية كاي-تربيع لبيرسون ($\chi^2$) عبر المعادلة الرياضية التراكمية التالية:

$$\chi^2 = \sum_{i=1}^{k} \frac{(O_i – E_i)^2}{E_i}$$

حيث يمثل $O_i$ التكرار الملاحظ في الفئة أو الخلية $i$، بينما يمثل $E_i$ التكرار المتوقع نظرياً في الخلية ذاتها تحت مظلة الفرضية الصفرية. تعكس هذه المعادلة جوهر الاختبار: فكلما اقتربت التكرارات الملاحظة من التكرارات المتوقعة، صغرت قيمة $\chi^2$ واقتربت من الصفر، مما يرفع القيمة الاحتمالية ويدعم الفرضية الصفرية. وعلى النقيض، كلما تعاظمت الفروق بين الواقع الميداني والافتراض النظري، تضخمت قيمة $\chi^2$ وتحركت عميقاً في الذيل الأيمن للتوزيع، مما يقلص القيمة الاحتمالية ويوجه القرار نحو رفض الفرضية الصفرية.

2.2 حساب درجات الحرية (Degrees of Freedom) للنماذج المختلفة

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

في حالة اختبار جودة المطابقة الأحادي (Goodness-of-Fit) لمتغير تصنيفي ذي $k$ من الفئات، تُحسب درجات الحرية ببساطة من خلال المعادلة:

$$df = k – 1$$

حيث يُطرح الرقم $1$ كقيد إحصائي نظراً لأن مجموع التكرارات المتوقعة يجب أن يساوي تماماً الحجم الكلي للعينة الملاحظة ($N$). أما إذا تطلب حساب التكرارات المتوقعة تقدير معالم إضافية من العينة ذاتها (مثل المتوسط الحسابي أو الانحراف المعياري للتحقق من التوزيع الطبيعي)، فإن كل معلمة مجهولة مقدرة تخصم درجة حرية إضافية، لتصبح المعادلة: $df = k – 1 – m$، حيث $m$ هو عدد المعالم المقدرة.

أما في سياق جداول الاقتران المزدوجة (Contingency Tables) المستخدمة في اختبار الاستقلالية بين متغيرين، حيث يحتوي الجدول على $r$ من الصفوف و $c$ من الأعمدة، فإن درجات الحرية تُحسب بالعلاقة الرياضية التالية:

$$df = (r – 1) \times (c – 1)$$

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

3. الدالة الأساسية لحساب القيمة الاحتمالية في لغة R: تشريح دالة pchisq()

3.1 البنية التركيبية (Syntax) لدالة pchisq() ومحدداتها

تُعد دالة pchisq() الدالة الرياضية الأساسية في لغة R المخصصة لحساب التوزيع التراكمي والمساحات الاحتمالية تحت منحنى توزيع كاي-تربيع. تأتي هذه الدالة مدمجة في حزمة stats التأسيسية التي تُحمّل افتراضياً مع تشغيل R، وتتميز ببنية تركيبية دقيقة وواضحة المعالم:

pchisq(q, df, ncp = 0, lower.tail = TRUE, log.p = FALSE)

يتطلب الاستخدام الفعال لهذه الدالة فهماً عميقاً لجميع معاملاتها (Arguments) البرمجية والإحصائية:

• المعامل q (Quantile): يمثل قيمة إحصائية كاي-تربيع المحسوبة الناتجة من تحليل البيانات، وهي قيمة عددية حقيقية غير سالبة ($q ge 0$).
• المعامل df (Degrees of Freedom): يمثل درجات الحرية المرتبطة بالاختبار، ويجب أن يكون عدداً حقيقياً موجباً يطابق التصميم الإحصائي للبيانات.
• المعامل ncp (Non-centrality Parameter): معلمة اللامركزية، وتُضبط افتراضياً على الصفر ($ncp = 0$) لاختبارات كاي-تربيع المركزية القياسية المستخدمة في اختبار الفرضيات التقليدية.
• المعامل lower.tail: قيمة منطقية (Boolean). عند ضبطها على TRUE (القيمة الافتراضية)، تحسب الدالة الاحتمال التراكمي للذيل الأيسر $P(X le q)$. أما عند ضبطها على FALSE، فإن الدالة تحسب احتمال الذيل الأيمن $P(X > q)$، وهو المطلوب حصراً لحساب القيمة الاحتمالية لاختبار كاي-تربيع.
• المعامل log.p: قيمة منطقية. عند ضبطها على TRUE، تعيد الدالة اللوغاريتم الطبيعي للاحتمال $ln(P)$، وهو خيار فائق الأهمية لتجنب أخطاء التفريغ الحسابي السفلي (Underflow) في النماذج الاحتمالية المعقدة وحسابات المعلومات البيولوجية والوراثية.

3.2 آلية عمل الدالة التراكمية في حساب المساحات الاحتمالية

تعمل دالة pchisq() كدالة توزيع تراكمي (Cumulative Distribution Function – CDF). تقوم الدالة رياضياً بإجراء عملية تكامل عددي متقدم لدالة كثافة الاحتمال الخاصة بتوزيع كاي-تربيع من الصفر وحتى النقطة $q$ المحددة:

$$F(q; df) = \int_{0}^{q} \frac{x^{(df/2) – 1} e^{-x/2}}{2^{df/2} \Gamma(df/2)} , dx$$

حيث يمثل $\Gamma$ دالة غاما الرياضية (Gamma Function). تنتمي دالة pchisq إلى الرباعية الدوالية القياسية في R للتعامل مع توزيع كاي-تربيع، والتي تضم أيضاً:

• dchisq(x, df): لحساب دالة كثافة الاحتمال (PDF) عند نقطة معينة.
• qchisq(p, df): دالة القواسم (Quantile Function) لحساب القيمة الحرجة المقابلة لاحتمال تراكمي محدد (عكس دالة pchisq).
• rchisq(n, df): لتوليد أرقام عشوائية تتبع توزيع كاي-تربيع بدرجات حرية محددة.

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

4. التحكم في اتجاه التوزيع وحساب الذيل الأيمن باستخدام المعامل lower.tail

4.1 أهمية الذيل الأيمن (Right-Tailed Test) في اختبار كاي-تربيع

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

وبناءً على ذلك، فإن القيمة الاحتمالية المطلوبة لاختبار الفرضية الصفرية هي احتمال الحصول على قيمة كاي-تربيع أكبر من أو تساوي القيمة المحسوبة في العينة، أي $P(X ge q)$. وإذا تم استدعاء الدالة بالشكل الافتراضي pchisq(q, df)، فإن ما يتم حسابه فعلياً هو المساحة اليسرى $P(X le q)$، وهي مكمل الاحتمال وليست القيمة الاحتمالية للاختبار.

وللحصول على القيمة الاحتمالية الصحيحة، يمكن للمحلل نظرياً استخدام إحدى طريقتين برمجيتين:

• الطريقة التقليدية: 1 - pchisq(q, df)
• الطريقة المعيارية المباشرة: pchisq(q, df, lower.tail = FALSE)

تُعد الطريقة الثانية المعتمدة على الضبط الصريح للمعامل lower.tail = FALSE هي الخيار الاحترافي الإلزامي في التحليل الإحصائي الرصين في لغة R لأسباب رياضية وحسابية بالغة الحساسية.

4.2 الدقة الرقمية وتفادي أخطاء الطرح الرياضي في R

تخزن الحواسيب الأرقام الكسرية بنظام الفاصلة العائمة المزدوجة الدقة (Double-precision floating-point format وفق معيار IEEE 754)، والذي يوفر حوالي 16 رقماً معنوياً من الدقة. عندما تكون إحصائية كاي-تربيع كبيرة جداً (على سبيل المثال $\chi^2 = 80$ بدرجات حرية $df = 2$)، فإن القيمة الاحتمالية التراكمية في الذيل الأيسر $P(X le q)$ تقترب بصورة شبه مطلقة من الرقم $1$ (مثل $0.99999999999999994$).

إذا قمنا بإجراء عملية الطرح 1 - pchisq(80, df = 2)، فإن الحساب يتعرض لما يُعرف بـ فقدان الدقة الحسابية (Catastrophic Cancellation / Loss of Significance)؛ حيث يؤدي طرح رقم قريب جداً من الواحد من الواحد الصحيح إلى إرجاع قيمة صفرية 0 بسبب حدود تقريب الآلة (Machine Epsilon $\approx 2.22 \times 10^{-16}$).

في المقابل، عندما نستخدم pchisq(80, df = 2, lower.tail = FALSE)، تتجاوز خوارزمية R الداخلية حساب الذيل الأيسر كلياً، وتتجه مباشرة لحساب دالة الذيل الأيمن التكميلية عبر دوال غاما اللوغاريتمية والكسور المستمرة، مما يتيح إرجاع القيمة الاحتمالية الدقيقة جداً (مثل $3.48 \times 10^{-18}$) دون أدنى تشويه رقمي. يوضح الجدول التالي مقارنة مفاهيمية بين الطريقتين الحسابيتين:

الصيغة البرمجية في R	المساحة المحسوبة رياضياً	الدقة عند القيم الكبيرة للإحصائية ($\chi^2 > 40$)	الاستخدام الموصى به
pchisq(q, df)	الذيل الأيسر $P(X le q)$	دقيقة (تقترب من 1)	غير صالحة لاختبار الفرضيات
1 - pchisq(q, df)	الذيل الأيمن $1 – P(X le q)$	معرضة للخطأ والتحول إلى 0	غير موصى بها مع الإحصائيات المرتفعة
pchisq(q, df, lower.tail = FALSE)	الذيل الأيمن المباشر $P(X > q)$	فائقة الدقة الرياضية حتى أرقام متناهية الصغر	المعيار الذهبي الإلزامي في R

5. تطبيق عملي 1: حساب القيمة الاحتمالية لاختبار جودة المطابقة (Goodness of Fit)

5.1 سيناريو بحثي نفسي: توزيع تفضيلات أساليب العلاج النفسي

لفهم الآلية التطبيقية المتكاملة، نفترض سيناريو دراسة إكلينيكية أُجريت في أحد مراكز الاستشارات النفسية الكبرى لمعرفة ما إذا كان لدى المراجعين الجدد تفضيل متساوٍ بين أربعة أساليب رئيسية للعلاج النفسي: العلاج المعرفي السلوكي (CBT)، والعلاج النفسي التحليلي (Psychodynamic)، والعلاج الإنساني الوجودي (Humanistic)، وعلاج القبول والالتزام (ACT). شملت العينة العشوائية $N = 200$ مريضاً، وكانت التكرارات الملاحظة لتفضيلاتهم على النحو التالي:

• العلاج المعرفي السلوكي (CBT): 70 مريضاً
• العلاج التحليلي (Psychodynamic): 40 مريضاً
• العلاج الإنساني (Humanistic): 55 مريضاً
• علاج القبول والالتزام (ACT): 35 مريضاً

تنص الفرضية الصفرية ($H_0$) على أن تفضيلات المرضى تتوزع بالتساوي بين الأساليب العلاجية الأربعة في المجتمع، أي بنسبة 25% لكل أسلوب ($p_1 = p_2 = p_3 = p_4 = 0.25$). وبناءً على ذلك، فإن التكرار المتوقع ($E$) لكل فئة يُحسب بضرب حجم العينة الإجمالي في النسبة المتوقعة: $E_i = 200 \times 0.25 = 50$ مريضاً لكل أسلوب علاجي.

يتم إدخال البيانات في R وحساب إحصائية كاي-تربيع يدوياً عبر الكود البرمجي التالي لتوضيح البنية الرياضية خطوة بخطوة:

observed <- c(70, 40, 55, 35)
expected <- c(50, 50, 50, 50)
chi_sq_stat <- sum((observed - expected)^2 / expected)

ينتج عن هذه المعادلة قيمة إحصائية كاي-تربيع تساوي:

$$\chi^2 = \frac{(70-50)^2}{50} + \frac{(40-50)^2}{50} + \frac{(55-50)^2}{50} + \frac{(35-50)^2}{50} = \frac{400}{50} + \frac{100}{50} + \frac{25}{50} + \frac{225}{50} = 8 + 2 + 0.5 + 4.5 = 15.0$$

وبما أن عدد الفئات $k = 4$، فإن درجات الحرية للاختبار هي: $df = 4 – 1 = 3$.

5.2 تنفيذ كود pchisq واستخراج وتفسير القيمة الاحتمالية

لحساب القيمة الاحتمالية المرتبطة بإحصائية الاختبار المحسوبة ($\chi^2 = 15.0$) ودرجات الحرية ($df = 3$) في بيئة R، نطبق الدالة مع ضبط المعامل المعياري للذيل الأيمن:

p_value <- pchisq(q = 15.0, df = 3, lower.tail = FALSE)
print(p_value)

تُرجع الدالة القيمة الاحتمالية الدقيقة: 0.001816912 (أي ما يقارب $p = .0018$).

التفسير الإحصائي والإكلينيكي: بما أن القيمة الاحتمالية الناتجة ($p = .0018$) أصغر بكثير من مستوى الدلالة المعياري ($\alpha = 0.05$)، فإننا نرفض الفرضية الصفرية بثقة إحصائية عالية ونقبل الفرضية البديلة. يشير هذا القرار الإحصائي إلى أن تفضيلات المرضى لأساليب العلاج النفسي لا تتوزع بالتساوي في المجتمع المستهدف، حيث يُظهر الفحص الاستكشافي للتكرارات وجود تفضيل ملحوظ للعلاج المعرفي السلوكي (70 مريضاً مقارنة بـ 50 متوقعاً)، وانخفاضاً واضحاً في الإقبال على علاج القبول والالتزام (35 مريضاً مقارنة بـ 50 متوقعاً)، مما يقدم مؤشرات تخطيطية بالغة الأهمية لمديري مراكز الدعم النفسي عند توزيع الموارد العلاجية.

6. تطبيق عملي 2: حساب القيمة الاحتمالية لاختبار الاستقلالية (Test of Independence)

6.1 سيناريو دراسة ارتباطية: العلاقة بين نوع الشخصية ونمط الاستجابة للضغوط

في دراسة في علم النفس الصحي وسيكولوجيا الشخصية، رغب فريق بحثي في اختبار مدى استقلالية نوع الشخصية (المصنفة إلى ثلاثة أنماط: النمط أ “Type A”، النمط ب “Type B”، النمط ج “Type C”) عن مستوى الاستجابة الفسيولوجية للضغوط النفسية (مصنف إلى مستويين: استجابة منخفضة/معتدلة، واستجابة مرتفعة/شديدة). تم جمع بيانات عينة ممثلة قوامها $N = 300$ مشارك، وجرى تلخيص البيانات في جدول اقتران مزدوج من أبعاد $3 \times 2$:

نوع الشخصية / مستوى الضغط	استجابة منخفضة/معتدلة	استجابة مرتفعة/شديدة	المجموع الهامشي للصفوف
النمط أ (Type A)	30	70	100
النمط ب (Type B)	65	35	100
النمط ج (Type C)	45	55	100
المجموع الهامشي للأعمدة	140	160	الإجمالي الكلي: 300

لحساب التكرارات المتوقعة ($E_{ij}$) لكل خلية في الجدول تحت فرضية الاستقلالية ($H_0$: لا توجد علاقة بين نوع الشخصية ومستوى الاستجابة للضغط)، نطبق القانون الإحصائي لضرب المجاميع الهامشية:

$$E_{ij} = \frac{\text{Row Total}_i \times \text{Column Total}_j}{N}$$

بتطبيق المعادلة على جميع الخلايا، نجد أن:

• الخلية (النمط أ، منخفض): $E_{11} = (100 \times 140) / 300 = 46.67$
• الخلية (النمط أ، مرتفع): $E_{12} = (100 \times 160) / 300 = 53.33$
• الخلية (النمط ب، منخفض): $E_{21} = (100 \times 140) / 300 = 46.67$
• الخلية (النمط ب، مرتفع): $E_{22} = (100 \times 160) / 300 = 53.33$
• الخلية (النمط ج، منخفض): $E_{31} = (100 \times 140) / 300 = 46.67$
• الخلية (النمط ج، مرتفع): $E_{32} = (100 \times 160) / 300 = 53.33$

تُحسب درجات الحرية للجدول ذي الأبعاد $r = 3$ و $c = 2$ كالتالي:

$$df = (3 – 1) \times (2 – 1) = 2 \times 1 = 2$$

6.2 الحساب البرمجي للقيمة الاحتمالية وتفسير الاستقلالية الإحصائية

نقوم بتمثيل جدول الاقتران في بيئة R وحساب إحصائية كاي-تربيع المحسوبة:

contingency_matrix <- matrix(c(30, 65, 45, 70, 35, 55), nrow = 3, ncol = 2, byrow = FALSE)
# حساب المجاميع الهامشية والمتوقعة برمجياً
row_totals <- rowSums(contingency_matrix)
col_totals <- colSums(contingency_matrix)
total_n <- sum(contingency_matrix)
expected_matrix <- outer(row_totals, col_totals) / total_n
chi_square_calc <- sum((contingency_matrix - expected_matrix)^2 / expected_matrix)

ينتج عن هذه الحسابات قيمة إحصائية الاختبار: $\chi^2 = 25.1786$.

لحساب القيمة الاحتمالية الدقيقة عبر دالة pchisq():

p_val_indep <- pchisq(q = chi_square_calc, df = 2, lower.tail = FALSE)
print(p_val_indep)

تُظهر المخرجات الحسابية قيمة احتمالية بالغة الصغر: 3.408428e-06 (أي $p = 0.00000341$).

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

7. الحساب المباشر مقابل الحساب التلقائي باستخدام دالة chisq.test() في R

7.1 استخدام دالة chisq.test() الشاملة في لغة R

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

عند تنفيذ الأمر التالي على مصفوفة دراسة الشخصية والضغوط:

res <- chisq.test(contingency_matrix, correct = FALSE)

تنشئ R كائناً برمجياً متكاملاً من النوع (List) يحتوي على ثروة من المخرجات الإحصائية المفصلة التي يمكن استخراجها عبر المعامل $:

• res$statistic: إحصائية كاي-تربيع المحسوبة ($\chi^2 = 25.1786$).
• res$parameter: درجات الحرية الدقيقة للاختبار ($df = 2$).
• res$p.value: القيمة الاحتمالية النهائية ($p = 3.408 \times 10^{-6}$).
• res$observed: مصفوفة التكرارات الملاحظة الفعلية.
• res$expected: مصفوفة التكرارات المتوقعة نظرياً تحت الفرضية الصفرية.
• res$residuals: البواقي البسيطة لبيرسون (Pearson Residuals) لحساب مقدار ونمط انحراف كل خلية عن المتوقع.
• res$stdres: البواقي المعيارية المنضبطة (Standardized Residuals)، وهي ذات أهمية كبرى في التشخيص الإحصائي لتحديد الخلايا المسؤولة فعلياً عن رفض الفرضية الصفرية (عندما تتجاوز قيمتها المطلقة $\pm 1.96$ أو $\pm 2.58$).

7.2 المقارنة المنهجية بين استخدام pchisq() ودالة chisq.test()

على الرغم من أن دالة chisq.test() تمثل حلاً سريعاً وشاملاً لتحليل البيانات الخام المباشرة، إلا أن دالة pchisq() تظل أداة حيوية ولا غنى عنها لأي باحث أو محلل إحصائي في حالات منهجية متعددة:

• إعادة تحليل الدراسات المنشورة (Secondary Analysis & Meta-Analysis): في كثير من الأحيان، يطالع الباحث دراسات سيكولوجية سابقة تورد قيمة إحصائية كاي-تربيع المحسوبة ودرجات الحرية فقط دون توفير البيانات الخام أو الجداول التفصيلية. في هذا السياق، تمثل pchisq(q, df, lower.tail = FALSE) الوسيلة الوحيدة للتحقق المستقل من صحة القيمة الاحتمالية المذكورة وتدقيقها بدقة رياضية متناهية.
• بناء الاختبارات المخصصة والنماذج المعقدة: عند تطبيق نماذج إحصائية متقدمة مثل نمذجة المعادلة الهيكلية (Structural Equation Modeling – SEM) أو نماذج الانحدار اللوجستي ومقارنة الفروق في نسبة الإمكان (Likelihood Ratio Chi-Square Difference Tests)، تُستخرج إحصائية الفرق في كاي-تربيع ($\Delta \chi^2$) والفرق في درجات الحرية ($\Delta df$)، ويتم تمريرها حصراً إلى دالة pchisq() لتقييم تحسن جودة مطابقة النموذج.
• الأتمتة البرمجية والمحاكاة: عند كتابة خوارزميات مخصصة لاختبار فرضيات معقدة تتضمن توليد آلاف العينات العشوائية، يكون استدعاء دالة pchisq() الخفيفة أسرع حسابياً بمراحل مقارنة بدالة chisq.test() الشاملة التي تستهلك وقتاً إضافياً في بناء الكائنات والمصفوفات المعقدة.

تتطابق مخرجات res$p.value المستخرجة من chisq.test() تطابقاً تاماً مع النتيجة المستخلصة عبر دالة pchisq(q = res$statistic, df = res$parameter, lower.tail = FALSE)، مما يثبت التناسق الرياضي والبرمجي الداخلي لبيئة R الإحصائية.

8. التعامل مع جداول الاقتران 2×2 وتصحيح ييتس للاستمرارية (Yates’ Continuity Correction)

8.1 الأساس الرياضي لتصحيح ييتس في العينات النفسية المتوسطة والصغيرة

ينشأ تحدٍ منهجي جوهري عند تطبيق اختبار كاي-تربيع على جداول الاقتران المزدوجة من أبعاد $2 \times 2$ (درجة حرية $df = 1$)؛ حيث تُعد التكرارات الملاحظة في الخلايا متغيراً عشوائياً منفصلاً (Discrete Random Variable) يتبع توزيع بواسون أو التوزيع فوق الهندسي، في حين أن توزيع كاي-تربيع النظري هو توزيع احتمالي متصل (Continuous Probability Distribution). يؤدي تقريب التوزيع المنفصل بتوزيع متصل في جداول $2 \times 2$ إلى المبالغة في تقدير قيمة إحصائية الاختبار المحسوبة، مما يقلص القيمة الاحتمالية اصطناعياً ويزيد من احتمالية الوقوع في خطأ النوع الأول (Type I Error)؛ أي رفض الفرضية الصفرية وهي صحيحة في الواقع.

لمعالجة هذا التحيز الإحصائي، اقترح عالم الإحصاء فرانك ييتس (Frank Yates) عام 1934 تعديلاً رياضياً يُعرف بـ تصحيح ييتس للاستمرارية (Yates’ Continuity Correction). يقوم هذا التعديل على طرح المقدار الثابت $0.5$ من القيمة المطلقة للفرق بين التكرار الملاحظ والتكرار المتوقع في كل خلية من خلايا الجدول قبل إجراء عملية التربيع، كما توضح المعادلة التالية:

$$\chi^2_{\text{Yates}} = \sum_{i=1}^{2}\sum_{j=1}^{2} \frac{(|O_{ij} – E_{ij}| – 0.5)^2}{E_{ij}}$$

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

8.2 حساب القيمة الاحتمالية مع وبدون تصحيح ييتس في R

تطبق دالة chisq.test() في R تصحيح ييتس تلقائياً وافتراضياً على أي مصفوفة أو جدول من أبعاد $2 \times 2$ من خلال المعامل الافتراضي correct = TRUE. يتيح R للباحث التحكم الكامل في تفعيل أو تعطيل هذا التصحيح لمقارنة النتائج وتقييم الحساسية المنهجية:

# إنشاء جدول اقتران 2x2: التدخل العلاجي (تجريبي/ضابط) مقابل التعافي (نعم/لا)
trial_data <- matrix(c(28, 12, 18, 22), nrow = 2, byrow = TRUE)
# التحليل بتطبيق تصحيح ييتس (الخيار الافتراضي التحفظي)
res_yates <- chisq.test(trial_data, correct = TRUE)
# التحليل بدون تصحيح ييتس (إحصائية بيرسون القياسية)
res_pearson <- chisq.test(trial_data, correct = FALSE)

عند فحص المخرجات نلاحظ الفروق التالية:

• مع تطبيق تصحيح ييتس: $\chi^2_{\text{Yates}} = 4.05$, $df = 1$, $p = .0442$
• بدون تطبيق تصحيح ييتس: $\chi^2_{\text{Pearson}} = 4.88$, $df = 1$, $p = .0272$

إذا رغب الباحث في حساب القيمة الاحتمالية يدوياً عبر pchisq() بعد تطبيق تصحيح ييتس، فإنه يمرر القيمة المصححة مباشرة:

p_val_manual_yates <- pchisq(q = 4.05, df = 1, lower.tail = FALSE)

الجدل الأكاديمي المنهجي: دار نقاش إحصائي مستفيض بين كبار المنهجيين السيكومتريين حول جدوى تصحيح ييتس؛ حيث يرى فريق من الإحصائيين أن تصحيح ييتس قد يكون مفرط التحفظ (Overly Conservative) ويقلل من القوة الإحصائية (Statistical Power) للاختبار، مما يرفع احتمالية الوقوع في خطأ النوع الثاني (Type II Error – الفشل في رصد أثر حقيقي). ولذلك، يُنصح في الممارسات النفسية الحديثة بتعطيل التصحيح (correct = FALSE) إذا كان الحجم الكلي للعينة كبيراً ($N > 100$) وجميع التكرارات المتوقعة تتجاوز 10، بينما يُلجأ إلى اختبار فيشر الدقيق (Fisher’s Exact Test) بدلاً من تصحيح ييتس إذا كانت العينة صغيرة والتكرارات منخفضة.

9. تفسير الدلالة الإحصائية وحجم التأثير (Effect Size) المرتبط بإحصائية كاي-تربيع

9.1 الفصل بين الدلالة الإحصائية (P-Value) والدلالة العملية (Practical Significance)

من أهم الأخطاء الإبستيمولوجية السائدة في الأبحاث السلوكية والنفسية هو الخلط بين الدلالة الإحصائية (Statistical Significance) المعبر عنها بالقيمة الاحتمالية ($p$-value) والدلالة العملية أو الإكلينيكية (Clinical/Practical Significance) المعبر عنها بحجم التأثير وقوة الرابطة. تُعبر القيمة الاحتمالية حصراً عن احتمالية وقوع البيانات تحت فرضية العدم، وتتأثر تأثراً بالغاً وشديداً بحجم العينة الإجمالي ($N$).

عندما تكون العينة ضخمة جداً (مثلاً $N = 50000$ مشارك)، فإن أدنى تباين طفيف أو فرق تافه لا يحمل أي معنى تطبيقي سيكولوجي سينتج عنه قيمة إحصائية كاي-تربيع مرتفعة وقيمة احتمالية بالغة الصغر ($p < .00001$). وعلى العكس من ذلك، في الدراسات الإكلينيكية النادرة ذات العينات المحدودة ($N = 25$)، قد يتواجد أثر علاجي قوي وعلاقة وثيقة بين المتغيرات ولكن القيمة الاحتمالية تفشل في بلوغ مستوى الدلالة ($p > .05$) نتيجة ضعف القوة الإحصائية.

لذلك، تشدد توصيات جمعية علم النفس الأمريكية (APA) على عدم الاكتفاء بتقديم القيمة الاحتمالية كمعيار وحيد لجودة النتيجة، بل يجب إقرانها دوماً بحساب وتفسير مقاييس حجم التأثير (Effect Size Indices) المعيارية التي تعكس قوة الاقتران التجريبي بمعزل عن حجم العينة.

9.2 حساب مقاييس حجم التأثير في R: فاي (Phi) ومعامل كرامر (Cramer’s V)

تعتمد مقاييس حجم التأثير المرتبطة باختبار كاي-تربيع على طبيعة وأبعاد جدول الاقتران المستخدم:

1. معامل فاي (Phi Coefficient – $phi$): يُستخدم حصراً في جداول الاقتران الثنائية $2 \times 2$. ويُحسب رياضياً بجذر حاصل قسمة إحصائية كاي-تربيع على حجم العينة الكلي:

$$\phi = \sqrt{\frac{\chi^2}{N}}$$

2. معامل كرامر (Cramer’s $V$): يمثل التعميم الرياضي لمعامل فاي لجداول الاقتران الأكبر من $2 \times 2$ (مثل الجداول ذات الأبعاد $r \times c$). وتأتي صياغته الرياضية على النحو التالي:

$$V = \sqrt{\frac{\chi^2}{N \times \min(r – 1, c – 1)}}$$

حيث يمثل $min(r-1, c-1)$ القيمة الصغرى بين درجات حرية الصفوف ودرجات حرية الأعمدة. تتراوح قيمة معامل كرامر دائماً بين $0$ (انعدام تام للارتباط) و $1$ (ارتباط تام بين المتغيرين التصنيفيين).

في بيئة R، يمكن حساب هذه المقاييس بسهولة ودقة عبر الأكواد الرياضية المباشرة أو عبر حزم متخصصة مثل rcompanion أو effectsize:

# حساب معامل كرامر لجدول الشخصية والضغط (3 صفوف × 2 أعمدة، N = 300، كاي-تربيع = 25.18)
n_total <- 300
chi2_val <- 25.1786
k_min <- min(3 - 1, 2 - 1) # k_min = 1
cramers_v <- sqrt(chi2_val / (n_total * k_min))
print(cramers_v) # ينتج: V = 0.2897

معايير التفسير النفسي وفق جاكوب كوهين (Cohen’s Benchmarks): يُفسر حجم التأثير لمعامل كرامر في جدول ذي درجة حرية صغرى مساوية للواحد ($min(r-1, c-1) = 1$) كالتالي:

• أثر صغير (Small Effect): $V \approx 0.10$ (يوضح ارتباطاً ضعيفاً ولكنه قابل للرصد).
• أثر متوسط (Medium Effect): $V \approx 0.30$ (يوضح ارتباطاً عملياً واضحاً في الظواهر السلوكية، وهو ما يطابق نتيجتنا $V = 0.29$).
• أثر كبير (Large Effect): $V ge 0.50$ (يوضح علاقة اقترانية قوية جداً وحاسمة).

10. الأخطاء الشائعة واستكشاف المشكلات وإصلاحها عند حساب P-Value في R

10.1 الأخطاء المفاهيمية والبرمجية الشائعة

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

• إغفال ضبط اتجاه الذيل lower.tail = FALSE: كما أشرنا سابقاً، يؤدي استدعاء pchisq(q, df) دون المعامل إلى حساب المساحة اليسرى التراكمية، مما يعطي ناتجاً مكملاً ($1 – p$). إذا كانت القيمة الحقيقية $p = .02$، فإن الإغفال سيعيد $p = .98$، مما يدفع المحلل إلى قبول الفرضية الصفرية خطأً ورفض اكتشاف علمي حقيقي.
• حساب خاطئ لدرجات الحرية ($df$): إدخال أبعاد المتغيرات بدلاً من طرح الواحد، مثل كتابة $df = 4$ بدلاً من $df = 3$ لمتغير رباعي الفئات في اختبار جودة المطابقة، مما يغير المعلمة التوزيعية لمنحنى كاي-تربيع ويشوه القيمة الاحتمالية.
• انتهاك افتراضات كوشران للتكرارات المتوقعة (Cochran’s Criteria): يشترط اختبار كاي-تربيع التقريبي ألا تقل التكرارات المتوقعة ($E$) في أي خلية عن $1$، وألا تتجاوز نسبة الخلايا التي تحتوي على تكرارات متوقعة تقل عن $5$ نسبة $20%$ من إجمالي خلايا الجدول. يؤدي خرق هذا الافتراض إلى انهيار التقريب التوزيعي لكاي-تربيع، مما يجعل القيم الاحتمالية المستخرجة مضللة تماماً.
• تمرير النسب المئوية بدلاً من التكرارات الخام (Counts): تفترض دوال R مثل chisq.test() أن المدخلات هي أعداد تكرارية حقيقية ($N$). إذا تم إدخال نسب مئوية (مثل $40%$ و $60%$ بدلاً من $400$ و $600$)، فإن الاختبار سيفقد وزنه الإحصائي الحقيقي وتصدر نتائج مشوهة لدرجات الحرية والقيم الاحتمالية.

10.2 التعامل مع رسائل التحذير (Warnings) في بيئة R

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

Warning message: In chisq.test(x) : Chi-squared approximation may be incorrect

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

1. اللجوء إلى اختبار فيشر الدقيق (Fisher’s Exact Test): وهو الخيار المعياري الذهبي لحساب القيمة الاحتمالية الدقيقة (Exact P-value) القائمة على التوزيع فوق الهندسي، خاصة في الجداول الصغيرة ذات التكرارات المحدودة:

fisher_result <- fisher.test(contingency_matrix)
print(fisher_result$p.value)

2. استخدام محاكاة مونت كارلو للتباديل (Monte Carlo Permutations): كما سنفصل في القسم التالي، يتيح R حساب القيمة الاحتمالية التجريبية عبر التوزيعات العشوائية المتكررة دون الاعتماد على التوزيع المقارب.

3. فحص واستبعاد القيم المفقودة (Handling Missing Data): تسبب القيم المفقودة من النوع NA في جداول البيانات تشوهات حسابية. يجب دائماً تنظيف المتغيرات الفئوية باستخدام الدالة na.omit() أو معالجتها ضمن جداول الاقتران باستخدام table(var1, var2, useNA = "no") لضمان بناء المصفوفات على الحالات المكتملة فقط.

11. أساليب المحاكاة وإعادة العينات (Monte Carlo Simulation) لحساب القيمة الاحتمالية

11.1 الحساب الاحتمالي المعتمد على محاكاة مونت كارلو في R

عندما تكون العينات الإكلينيكية نادرة، أو عندما تكون جداول الاقتران كبيرة ومعقدة (مثل جداول $5 \times 4$) مع وجود خلايا فارغة أو ذات تكرارات متوقعة شديدة الانخفاض، يفقد التوزيع النظري لكاي-تربيع صلاحيته الرياضية وتصبح صيغ التكامل الكلاسيكية في pchisq() غير كافية لتمثيل الواقع الاحتمالي بدقة. في هذه السيناريوهات المتقدمة، يُعد أسلوب محاكاة مونت كارلو (Monte Carlo Simulation) المنهجية المثلى لتقدير القيمة الاحتمالية بدقة فائقة.

تقوم خوارزمية مونت كارلو في R على توليد عدد ضخم من الجداول العشوائية التخيلية (Replications $B$) التي تمتلك نفس المجاميع الهامشية للصفوف والأعمدة للجدول الأصلي للبيانات تحت فرضية الاستقلالية الصفرية التامة. تُحسب إحصائية كاي-تربيع لكل جدول عشوائي مولد، ثم تُستخرج القيمة الاحتمالية التجريبية (Empirical P-Value) كنسبة عدد الجداول العشوائية التي أنتجت إحصائية كاي-تربيع مساوية للإحصائية المحسوبة للبيانات الأصلية أو متجاوزة لها، مقسومة على العدد الإجمالي للمحاكاة ($B$).

توفر دالة chisq.test() خياراً مدمجاً ومباشراً لتنفيذ محاكاة مونت كارلو عبر تفعيل المعامل simulate.p.value = TRUE وتحديد عدد مرات التكرار عبر المعامل B (ويُفضل ضبطه على $B = 10000$ أو $B = 100000$ لتحقيق أقصى استقرار إحصائي):

# ضبط بذرة التوليد العشوائي لضمان قابلية إعادة الإنتاج
set.seed(2026)
# تنفيذ اختبار كاي-تربيع بمحاكاة مونت كارلو على جدول متناثر التكرارات
sparse_table <- matrix(c(1, 15, 0, 12, 2, 18), nrow = 3, byrow = TRUE)
mc_res <- chisq.test(sparse_table, simulate.p.value = TRUE, B = 10000)
print(mc_res)

تُظهر المخرجات القيمة الاحتمالية المبنية على المحاكاة: p-value = 0.0384 دون ظهور أي رسائل تحذيرية حول خطأ التقريب، مما يمنح الباحث ثقة منهجية كاملة في القرار الإحصائي المتخذ.

11.2 مقارنة دقة القيمة الاحتمالية للمحاكاة مع التوزيع المقارب

لتوضيح التفوق النوعي لأسلوب المحاكاة عند انتهاك الافتراضات، نقارن بين القيمة الاحتمالية المستخرجة من دالة pchisq() التقليدية المستندة إلى التوزيع المقارب النظري، والقيمة المستخرجة من محاكاة مونت كارلو لبيانات غير متوازنة:

إذا كانت إحصائية كاي-تربيع المحسوبة لجدول متناثر التكرارات هي $\chi^2 = 6.45$ بدرجات حرية $df = 2$، فإن استدعاء دالة التكامل الكلاسيكية pchisq(6.45, df = 2, lower.tail = FALSE) يُعطي قيمة احتمالية تساوي 0.03975. ولكن نظراً لعدم تحقق شرط كوشران في الخلايا، فإن هذه القيمة النظرية قد تكون متفائلة أو مشوهة. عند تطبيق محاكاة مونت كارلو بـ $100,000$ تكرار، قد تبلغ القيمة الاحتمالية التجريبية الدقيقة 0.0542، وهو فارق جوهري ينقل القرار الإحصائي من منطقة الرفض الدال ($p < .05$) إلى منطقة عدم الدلالة ($p > .05$)، مما يجنب الباحث إعلان اكتشاف زائف.

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

12. أفضل الممارسات لتوثيق وكتابة نتائج كاي-تربيع وفق دليل جمعية علم النفس الأمريكية (APA)

12.1 صياغة الفقرة الإحصائية لاختبار كاي-تربيع وفق معايير APA 7

يضع دليل النشر العلمي الصادر عن جمعية علم النفس الأمريكية في نسخته السابعة (APA 7th Edition) قواعد وإرشادات صارمة لكيفية كتابة وتوثيق نتائج التحليلات الإحصائية في المتن الأكاديمي. يتطلب التوثيق المعياري لاختبار كاي-تربيع إدراج العناصر التالية بنسق طباعي محدد:

• الرمز الإغريقي للإحصائية بالحرف المائل ($\chi^2$).
• درجات الحرية ($df$) وحجم العينة الإجمالي ($N$) محصورين بين قوسين متجاورين: $\chi^2(df, N = dots)$.
• قيمة إحصائية كاي-تربيع المحسوبة مقربة إلى خانتين عشريتين.
• القيمة الاحتمالية الدقيقة ($p$) مقربة إلى خانتين أو ثلاث خانات عشرية دون وضع صفر قبل الفاصلة (مثال: $p = .015$ وليس $p = 0.015$).
• إذا كانت القيمة الاحتمالية شديدة الصغر (أقل من $0.001$)، تُكتب دائماً بالصيغة ($p < .001$)، ويُحظر تماماً كتابة ($p = .000$) الشائعة خطأً في مخرجات بعض البرمجيات.
• إدراج قيمة حجم التأثير المعياري (معامل فاي $phi$ أو معامل كرامر $V$) مع فترات الثقة (Confidence Intervals) إن أمكن.

نماذج تطبيقية لصياغة النتائج في الأبحاث النفسية:

نموذج 1: اختبار جودة المطابقة (Goodness-of-Fit):
“أظهرت نتائج اختبار كاي-تربيع لجودة المطابقة وجود فروق ذات دلالة إحصائية في تفضيلات المرضى لأساليب العلاج النفسي الأربعة، $\chi^2(3, N = 200) = 15.00, p = .002, V = .27$، مما يشير إلى انحراف واضح عن التوزيع المتساوي المتوقع تحت الفرضية الصفرية.”

نموذج 2: اختبار الاستقلالية (Test of Independence):
“كشف اختبار كاي-تربيع للاستقلالية عن وجود علاقة ارتباطية دالة إحصائياً بين نوع الشخصية ومستوى الاستجابة للضغوط النفسية، $\chi^2(2, N = 300) = 25.18, p < .001, V = .29, 95% \text{ CI } [.18, .40]$. وبناءً على ذلك، تم رفض الفرضية الصفرية التي تفترض استقلالية المتغيرين.”

12.2 إنشاء مخططات بصرية وتوضيحية لتوزيع كاي-تربيع في R باستخدام ggplot2

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

يوضح الكود البرمجي المتكامل التالي كيفية بناء وتصدير المخطط البياني لتوزيع كاي-تربيع بدرجات حرية $df = 3$ مع إحصائية محسوبة $\chi^2 = 15.0$:

# تحميل المكتبات الضرورية
library(ggplot2)
# إنشاء متسلسلة قيم لمحور السينات
x_vals <- seq(0, 25, length.out = 1000)
df_val <- 3
q_stat <- 15.0
# بناء إطار البيانات
plot_data <- data.frame(x = x_vals, y = dchisq(x_vals, df = df_val))
# تصفية البيانات لمنطقة التظليل (الذيل الأيمن: x >= q_stat)
shade_data <- subset(plot_data, x >= q_stat)
# بناء الرسم البياني الاحترافي
p <- ggplot(plot_data, aes(x = x, y = y)) +
geom_line(color = "#1f77b4", size = 1.2) +
geom_ribbon(data = shade_data, aes(ymin = 0, ymax = y), fill = "#d62728", alpha = 0.5) +
geom_vline(xintercept = q_stat, linetype = "dashed", color = "#d62728", size = 1) +
annotate("text", x = q_stat + 2, y = 0.05, label = paste("Chi-Sq =", q_stat, "np = .0018"), color = "#d62728", hjust = 0) +
labs(title = "توزيع كاي-تربيع والمساحة الاحتمالية للذيل الأيمن (df = 3)",
subtitle = "المنطقة المظللة باللون الأحمر تمثل القيمة الاحتمالية (P-Value)",
x = "قيمة إحصائية كاي-تربيع", y = "كثافة الاحتمال") +
theme_minimal(base_size = 14) +
theme(plot.title = element_text(face = "bold", hjust = 0.5),
plot.subtitle = element_text(hjust = 0.5))
# عرض المخطط وتصديره بدقة عالية (300 DPI)
print(p)
ggsave("chi_square_distribution_apa.png", plot = p, width = 8, height = 5, dpi = 300)

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

خاتمة واستنتاجات منهجية

يمثل حساب وتفسير القيمة الاحتمالية لإحصائية كاي-تربيع في بيئة R مهارة تحليلية تأسيسية تجمع بين العمق الرياضي والإتقان البرمجي والوعي السيكومتري التطبيقي. إن إدراك الباحث للاختلافات الجوهرية بين دوال التوزيع التراكمي المباشرة مثل pchisq() والدوال الاختبارية المجمعة مثل chisq.test() يُمكنه من تطويع بيئة R للتعامل مع أكثر التصاميم الإحصائية تعقيداً ودقة.

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

References

• American Psychological Association. (2020). Publication manual of the American Psychological Association (7th ed.). American Psychological Association. https://doi.org/10.1037/0000165-000
• Agresti, A. (2013). Categorical data analysis (3rd ed.). John Wiley & Sons. https://www.wiley.com/en-us/Categorical+Data+Analysis%2C+3rd+Edition-p-9780470463635
• Cohen, J. (1988). Statistical power analysis for the behavioral sciences (2nd ed.). Lawrence Erlbaum Associates. https://doi.org/10.4324/9780203771587
• 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
• Pearson, K. (1900). On the criterion that a given system of deviations from the probable in the case of a correlated system of variables is such that it can be reasonably supposed to have arisen from random sampling. Philosophical Magazine Series 5, 50(302), 157–175. https://doi.org/10.1080/14786440009463897
• 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 (2nd ed.). Springer-Verlag. https://doi.org/10.1007/978-3-319-24277-4
• Yates, F. (1934). Contingency tables involving small numbers and the $\chi^2$ test. Supplement to the Journal of the Royal Statistical Society, 1(2), 217–235. https://doi.org/10.2307/2983604