الوحدة التاسعة: الانحدار المتعدد واللوجستي ومقارنة النماذج

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

إعادة تنظيم: انتقل تحليل العوامل إلى الوحدة العاشرة، وانتقلت أخلاقيات p-hacking وحدود NHST إلى الوحدة السادسة.

الإحصاء الوصفي والتحليلي باستخدام البرامج

البرامج الإحصائية تسرع العمليات الحسابية وتقلل الأخطاء. لكن فهم المفاهيم أساسي — البرنامج لا يختار التحليل المناسب!

البرنامج الاستخدام الأمثل نقاط القوة القيد
SPSS البحث في الرياضة والطب واجهة رسومية، سهل التعلم تكلفة، محدود في النماذج المعقدة
R البحث الأكاديمي المتقدم مجاني، مكتبات ضخمة، قابل للبرمجة منحنى تعليمي حاد
Jamovi / JASP تدريس الإحصاء مجاني، واجهة سهلة، نتائج APA جاهزة أقل مرونة من R
Excel البيانات البسيطة والرسوم متاح للجميع محدود في الاختبارات المتقدمة
خطوات العمل على أي برنامج
  1. تنظيف البيانات (تعامل مع القيم المفقودة والشاذة)
  2. فحص افتراضات التحليل (الطبيعية، التباين)
  3. اختيار الاختبار المناسب بناءً على نوع البيانات والسؤال
  4. تشغيل التحليل وتفسير النتائج
  5. إعداد التقارير (Effect size + CI + p-value)

الانحدار الخطي المتعدد Multiple Linear Regression

امتداد للانحدار البسيط: يتنبأ بمتغير تابع من عدة متغيرات مستقلة في آن واحد.

معادلة الانحدار المتعدد
Ŷ = a + b₁X₁ + b₂X₂ + ... + bₙXₙ
حيث: Ŷ = القيمة المتوقعة | X₁, X₂... = المتغيرات المستقلة | b₁, b₂... = معاملات الانحدار
مثال رياضي: التنبؤ بأداء السباق

تنبؤ بزمن 5000م (Y) من:

  • X₁ = VO2max (أقصى استهلاك للأكسجين)
  • X₂ = الوزن النسبي (كتلة العضلات/الوزن الكلي)
  • X₃ = سنوات التدريب

النتيجة: Ŷ = 25 − 0.3(VO2max) − 15(نسبة العضلات) − 0.5(سنوات التدريب)

R² = 0.75: 75% من تباين الأداء يُفسر بالمتغيرات الثلاثة.

شروط وقيود الانحدار المتعدد
  • الخطية: العلاقة خطية بين المتغيرات
  • الاستقلالية: قليل الارتباط بين المتغيرات المستقلة (عدم وجود تعددية خطية multicollinearity)
  • التوزيع الطبيعي: للأخطاء/البواقي
  • حجم العينة: n ≥ 50 + 8k (حيث k = عدد المتغيرات المستقلة) — قاعدة Green (1991). بديل أكثر تحفظاً: 10–20 مشاهدة لكل متغير مستقل (يُفضَّل 20:1 للتحقق المتقاطع الموثوق).

اختيار المتغيرات

الطريقة الوصف متى تُستخدم؟
Enter (إدخال كامل) إدخال جميع المتغيرات دفعة واحدة نظرية قوية تحدد المتغيرات مسبقاً
Stepwise (خطوة بخطوة) إضافة/إزالة المتغيرات تبعاً للمعايير الإحصائية استكشافي — لكن يُنصح بتجنبه في البحث النظري
Hierarchical (هرمي) إدخال المتغيرات في كتل حسب الأولوية التحكم في المتغيرات الدخيلة أولاً
R² المعدَّل (Adjusted R²)
R²_adj = 1 − [(1 − R²)(n − 1) / (n − k − 1)]
n = حجم العينة | k = عدد المتغيرات المستقلة
لماذا R² المعدَّل وليس R² العادي؟

R² العادي يرتفع تلقائياً مع كل متغير مستقل جديد — حتى لو كان المتغير عشوائياً لا فائدة منه. R² المعدَّل يعاقب على إضافة متغيرات ضعيفة:

  • إذا أضفت متغيراً مفيداً: R²_adj يرتفع
  • إذا أضفت متغيراً غير مفيد: R²_adj ينخفض
  • استخدم دائماً R²_adj عند مقارنة نماذج ذات أعداد مختلفة من المتغيرات
التعددية الخطية (Multicollinearity) — VIF والتسامح

التعددية الخطية تحدث حين تترابط المتغيرات المستقلة فيما بينها ارتباطاً عالياً، مما يُضخّم الأخطاء المعيارية ويُعيق تفسير المعاملات. التشخيص يعتمد على:

عامل تضخم التباين (VIF)

VIF = 1 / Tolerance

VIF التفسير
< 5 مقبول
5 – 10 تحذير — راجع النموذج
> 10 إشكالي — تعددية خطية حادة

التسامح (Tolerance)

Tolerance = 1 − R²_j

R²_j = R² من انحدار المتغير j على بقية المتغيرات

Tolerance التفسير
> 0.20 مقبول
0.10 – 0.20 تحذير
< 0.10 إشكالي

مثال رياضي: نموذج للتنبؤ بزمن سباق 5000م يضم VO₂max ونبض القلب الأقصى والعمر. VIF(VO₂max) = 8.3 → تحذير: بعض المتغيرات تقيس جانباً فسيولوجياً مشتركاً → احذف أحد المتغيرات المتداخلة، ادمجها في مؤشر واحد، أو استخدم اختياراً نظرياً أوضح.

مقارنة النماذج AIC / BIC

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

معيار أكايكي للمعلومات (AIC)

AIC = 2k − 2 ln(L)
  • k = عدد المعاملات (بما فيها ثابت النموذج)
  • L = احتمالية النموذج (Likelihood)
  • القيمة الأصغر أفضل
  • مناسب لبيانات بأحجام عينة متوسطة

معيار بيز للمعلومات (BIC)

BIC = k × ln(n) − 2 ln(L)
  • n = حجم العينة
  • يعاقب التعقيد أكثر من AIC
  • القيمة الأصغر أفضل
  • يُفضَّل مع أحجام عينة كبيرة
قراءة فروق AIC/BIC
ΔAIC (أو ΔBIC) قوة الدليل لصالح النموذج الأفضل
0 – 2 لا يوجد فرق جوهري — كلا النموذجَين مقبولان
4 – 7 دليل متوسط لصالح النموذج الأقل قيمة
> 10 دليل قوي لصالح النموذج الأقل قيمة؛ وإذا تقاربت الملاءمة فاختر الأبسط
مثال: مقارنة نماذج التنبؤ بزمن 5000م
النموذج المتغيرات R²_adj AIC BIC
A VO₂max فقط 0.62 145.2 148.7
B VO₂max + نسبة العضلات 0.74 138.5 143.8
C ✓ VO₂max + نسبة العضلات + سنوات التدريب 0.76 136.1 143.2
D النموذج C + 3 متغيرات إضافية غير مهمة 0.74 141.3 152.4

القرار: النموذج C هو الأمثل — أعلى R²_adj، وأدنى AIC وBIC. النموذج D رغم زيادة المتغيرات يُظهر ارتفاعاً في AIC/BIC مما يدل على إفراط في المعاملات (overfitting).

تحليل العوامل Factor Analysis

تحليل العوامل يكشف عن المتغيرات الكامنة (Latent Variables) التي تفسر الارتباطات بين متغيرات مقاسة. يُستخدم كثيراً في بناء الاستبيانات وتحليل بنية القدرات الرياضية.

التحليل الاستكشافي (EFA)

  • متى: عند استكشاف بنية مجهولة سابقاً
  • الهدف: اكتشاف عدد العوامل وبنيتها
  • القرار بعدد العوامل: القيمة الذاتية > 1 (معيار كايزر) أو مخطط الانعطاف (Scree Plot)
  • حمل العامل المقبول: |loading| ≥ 0.40

التحليل التوكيدي (CFA)

  • متى: عند اختبار بنية محددة مسبقاً
  • الهدف: التحقق من ملاءمة النموذج النظري للبيانات
  • مؤشرات الملاءمة:
    • CFI ≥ 0.95 (جيد)
    • RMSEA ≤ 0.06 (جيد)
    • SRMR ≤ 0.08 (جيد)
مثال: مقياس الإرهاق الرياضي — Raedeke & Smith (2001)

مقياس "إرهاق الرياضي" (15 بنداً) — هل يقيس مفهوماً واحداً أم عدة عوامل؟

EFA
تحليل استكشافي أظهر 3 عوامل ذات قيم ذاتية > 1:
(1) الإجهاد الانفعالي والجسدي | (2) الإحساس المتناقص بالإنجاز | (3) انخفاض قيمة الرياضة
CFA
تحليل توكيدي على عينة مستقلة للتحقق من البنية الثلاثية:
CFI = 0.96، RMSEA = 0.052، SRMR = 0.065 → ملاءمة جيدة

التطبيق: يُستخدم هذا المقياس لتقييم مستوى الإرهاق لدى الرياضيين وتوجيه برامج التعافي النفسي.

متطلبات تحليل العوامل
  • حجم العينة: n ≥ 200 على الأقل (أو نسبة 10:1 مشارك لكل بند)
  • نوع البيانات: مستمرة أو رتبية (مقياس ليكرت ≥ 5 نقاط)
  • فحص KMO: KMO ≥ 0.70 للمصفوفة مناسبة للتحليل
  • Bartlett's Test: p < 0.05 يُثبت وجود ارتباطات كافية بين البنود

الانحدار اللوجستي Logistic Regression

يُستخدم عندما يكون المتغير التابع فئوياً ثنائياً (نعم/لا، فوز/خسارة، إصابة/لا إصابة).

لماذا لا نستخدم الانحدار الخطي العادي؟

لأن المتغير التابع ثنائي (0 أو 1)، والانحدار الخطي يتوقع قيماً مستمرة. الانحدار اللوجستي يتنبأ باحتمال وقوع الحدث.

مثال: التنبؤ بخطر الإصابة

تنبؤ باحتمال إصابة الرباط الصليبي (نعم/لا) من:

  • عمر اللاعب
  • سنوات اللعب السابقة
  • قوة عضلات الفخذ (نسبة مئوية)
  • التوازن أحادي القائمة (الثواني)

النتيجة: "احتمال الإصابة = 15%" → يستدعي برنامج وقائي مكثف.

Odds Ratio (OR)

في الانحدار اللوجستي، OR يخبرنا بكم مرة تزيد أو تنقص الأرجحية (odds)، لا الاحتمال مباشرةً، مع زيادة المتغير المستقل بوحدة واحدة مع ثبات بقية المتغيرات.

  • OR = 1: لا تأثير
  • OR > 1: عامل خطر
  • OR < 1: عامل وقائي

مثال: OR = 2.0 يعني أن أرجحية الإصابة تضاعفت، لكن الاحتمال لا يتضاعف بالضرورة؛ فالانتقال من احتمال 10% إلى أرجحية مضاعفة يعطي احتمالاً يقارب 18%.

التحليلات التنبؤية في الرياضة Predictive Analytics

استخدام البيانات والإحصاء/التعلم الآلي للتنبؤ بالنتائج المستقبلية. تطبيقاتها في الرياضة:

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

نقطة تحقق: النموذج المناسب للسؤال المناسب

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

محتوى انتقالي: الأقسام التالية باقية مؤقتاً لتجنب كسر المراجع القديمة. المسار الرسمي الجديد يدرس تحليل العوامل في الوحدة العاشرة والنزاهة البحثية في الوحدة السادسة.

أخلاقيات البحث الإحصائي Research Ethics & Statistical Integrity

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

ما هو P-Hacking؟

هو إجراء اختبارات متعددة أو تعديل البيانات بطرق مختلفة حتى يصبح p < 0.05، ثم الإبلاغ فقط عن النتيجة الدالة. يُطلق عليه أيضاً "تعذيب البيانات حتى تعترف".

  • مثال 1: تجربة مع 10 متغيرات تابعة — يُبلَّغ فقط عن المتغير الذي أعطى p = 0.03 وتُهمل الـ 9 الأخرى.
  • مثال 2: إزالة المشاركين الذين أعطوا نتائج "غير مريحة" بحجة أنهم قيم شاذة.
  • مثال 3: تجربة اختبارات متعددة (t-test، Mann-Whitney، Wilcoxon) والإبلاغ عن الذي أعطى p < 0.05.
التحيز في أبحاث الرياضة
  • تحيز النشر (Publication Bias): الدراسات ذات النتائج الدالة تُنشر أكثر من الدراسات غير الدالة — مما يُشوّه صورة الأدلة المتاحة.
  • تحيز الكيان الراعي: الدراسات الممولة من شركات معدات أو مكملات غذائية تُظهر نتائج مفيدة لها بنسب أعلى من الدراسات المستقلة.
  • تحيز الاختيار: اختيار الرياضيين الأفضل أداءً لإثبات فاعلية برنامج، ثم تعميم النتيجة على الجميع.
  • انتقائية البيانات: الإبلاغ عن نقطة زمنية معينة بعد التدخل لأن النتائج أفضل فيها.
أشكال P-hacking
  • اختبارات متعددة: اختبار عشرات الفرضيات حتى نجد واحدة دالة
  • إيقاف مبكر: إيقاف جمع البيانات عندما نحصل على p < 0.05
  • تغيير المجموعات: تعديل التصنيفات بعد رؤية البيانات
  • إزالة القيم الشاذة: إزالة بيانات "لأنها لا تناسب"
  • HARKing: تقديم فرضية بعد رؤية النتائج (Hypothesizing After Results are Known)
مخاطر التلاعب بالبيانات
  • قرارات تدريبية خاطئة تضر الرياضيين (زيادة الإجهاد، التغذية غير المناسبة)
  • إهدار ميزانيات الأندية على تدخلات غير فعالة
  • ضرر بالسمعة الأكاديمية والمهنية عند الكشف عن التلاعب
  • تراجع الثقة في أبحاث علوم الرياضة عموماً
كيف تحمي نزاهة بحثك؟
1
التسجيل المسبق (Pre-Registration): سجّل فرضياتك ومتغيراتك في منصة مثل OSF قبل جمع البيانات.
2
مستوى الدلالة المحدد مسبقاً: حدد α = 0.05 قبل بدء الدراسة، ولا تغيّره بعد رؤية النتائج.
3
تصحيح المقارنات المتعددة: إذا أجريت اختبارات متعددة، طبّق تصحيح Bonferroni (α المعدَّل = α / عدد الاختبارات).
4
الشفافية الكاملة: أبلّغ عن جميع الاختبارات التي أجريتها، وليس فقط الدالة منها.
5
تحليل حجم العينة مسبقاً (Power Analysis): احسب العينة المطلوبة لقوة 0.80 قبل بدء الجمع.
مبادرات تحسين البحث
  • Open Science: مشاركة البيانات والكود لإعادة الإنتاج
  • Registered Reports: قبول البحث للنشر بناءً على المنهجية لا النتائج
  • Confidence Index: تقييم جودة الأدلة بدلاً من الثنائية (دال/غير دال)
أسئلة تقويمية — الفصل السابع
مفاهيمي ما الفرق بين الانحدار الخطي واللوجستي؟ متى تستخدم كل منهما؟
تطبيقي في نموذج تنبؤ بالإصابات، حصلت على OR = 2.5 لـ "عدم التوازن". كيف تفسر هذا؟
بحثي ناقش: لماذا يُعتبر Stepwise regression مُشكلاً في البحث العلمي الرصين؟
نقدي ما هي "أزمة إعادة الإنتاج" في العلوم؟ وما علاقتها بـ p-hacking؟
حسابي نموذج انحدار للتنبؤ بالحمل التدريبي اليومي يضم 4 متغيرات: R² = 0.81، n = 45. احسب R²_adj وفسّر الفرق.
تطبيقي في نموذج انحدار متعدد لتنبؤ بالأداء: VIF(VO₂max) = 12.4. ماذا يعني ذلك؟ وما الحل المقترح؟
مقارنة نموذجان: الأول AIC = 152، الثاني AIC = 141. ما النموذج الأفضل؟ وما قوة الدليل بناءً على الفرق؟
بحثي مقياس "الثقة بالنفس الرياضية" مكوّن من 20 بنداً. اشرح كيف تستخدم EFA ثم CFA للتحقق من بنيته.