الإحصاء النفسيالقياس النفسيتصميم التجارب وتحليل البيانات

كيفية إجراء تحليل التباين للقياسات المتكررة في إكسيل

دليل أكاديمي شامل يشرح كيفية إجراء تحليل التباين للقياسات المتكررة (Repeated Measures ANOVA) في إكسيل، وفهم مخرجاته وتفسير نتائجه في البحوث النفسية.

تاريخ النشر

يُمثل تحليل التباين للقياسات المتكررة (Repeated Measures Analysis of Variance – RM-ANOVA) أحد أكثر الأساليب الإحصائية المعلمية تقدماً وأهمية في مجالات البحوث التجريبية والنفسية والعلوم السلوكية والطبية المعاصرة. تكمن القوة المنهجية لهذا الاختبار في قدرته الفريدة على تتبع ورصد استجابات نفس العينة من المفحوصين أو المشاركين عبر مستويات متعددة لمتغير تجريبي معين، أو عبر نقاط زمنية متعاقبة، مما يتيح للباحثين استكشاف ديناميات التغير النمائي والمعرفي والسلوكي بدقة فائقة. وفي حين تفترض الاختبارات الإحصائية التقليدية مثل تحليل التباين الأحادي المستقل استقلالية الملاحظات بين المجموعات، فإن نموذج القياسات المتكررة يستفيد جوهرياً من الارتباط الداخلي بين استجابات الفرد الواحد، محولاً ما يُعد “خطأً عشوائياً” في التصاميم الأخرى إلى تباين يمكن ضبطه وعزله بدقة.

تاريخياً، ارتبط تنفيذ هذا النوع من التحليلات الإحصائية المتقدمة بحزم البرمجيات الاحتكارية المعقدة ومكلفة الثمن مثل SPSS وSAS وStata، أو عبر لغات البرمجة الإحصائية مثل R وPython. ورغم قوة هذه البرمجيات، إلا أن استخدام برنامج مايكروسوفت إكسيل (Microsoft Excel) لتنفيذ تحليل التباين للقياسات المتكررة يمنح الباحثين والأكاديميين ميزة تربوية ومنهجية فريدة؛ إذ يُجبر المحلل على فهم آليات تفكيك التباين خطوة بخطوة، والاطلاع المباشر على البنية الهيكلية لجدول تحليل التباين دون الوقوع في فخ “الصندوق الأسود” البرمجي. وبفضل حزمة أدوات تحليل البيانات المدمجة (Data Analysis ToolPak)، يوفر إكسيل بيئة مرنة، شفافة، وسهلة الوصول تتيح استخراج النتائج بدقة متناهية تتطابق تماماً مع أرقى البرمجيات الإحصائية المتخصصة.

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

1. مفهوم تحليل التباين للقياسات المتكررة وأهميته في البحوث النفسية

1.1 التعريف النظري والرياضي للتحليل

يُعرّف تحليل التباين للقياسات المتكررة (Repeated Measures ANOVA)، والذي يُعرف أيضاً بالتصميم داخل المفحوصين (Within-Subjects Design)، بأنه امتداد اختباري معلمي يُستخدم لمقارنة ثلاثة متوسطات أو أكثر لمجموعة واحدة متصلة من المشاركين تم قياسهم تحت ظروف تجريبية مختلفة أو عند فترات زمنية متتابعة. يقوم الأساس المنطقي لهذا الاختبار على قياس نفس الوحدة التجريبية (غالباً الإنسان أو الحيوان في التجارب السلوكية) بصورة متكررة، مما يجعل كل مفحوص بمثابة ضابط لنفسه (Self-Control). هذا التصميم المنهجي يسمح بالتحكم الدقيق في التباين الفردي الثابت بين الأشخاص، مثل القدرات المعرفية الفطرية، والسمات الشخصية المستقرة، والمحددات الجينية والبيولوجية.

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

Yij = μ + πi + τj + εij

حيث تُمثل Yij الدرجة الملاحظة للمفحوص رقم i في مستوى المعالجة التجريبية رقم j، و μ هو المتوسط العام الكلي للمجتمع، و πi هو الأثر الثابت الخاص بالمفحوص الفردي رقم i (ويمثل الفروق الفردية المستقرة عبر الزمن)، و τj هو الأثر الحقيقي لمستوى المعالجة التجريبية رقم j (وهو المتغير محل الاهتمام والبحث)، وأخيراً εij يُمثل الخطأ العشوائي المتبقي المرتبط بالمفحوص i في المعالجة j، والذي يفترض النموذج أنه موزع توزيعاً طبيعياً بمتوسط صفر وتباين ثابت.

تتجلى العبقرية الرياضية لهذا التحليل في قدرته على عزل تباين الأفراد (πi) من حد التباين الخطئي العام. في التصاميم المستقلة بين المجموعات، يندمج هذا التباين الفردي داخل تباين الخطأ غير المفسر (Within-Groups Error)، مما يؤدي إلى تضخيم مقام النسبة الفائية (F-ratio) ويقلل بالتالي من حساسية الاختبار في اكتشاف الفروق الحقيقية. أما في تحليل القياسات المتكررة، فإن هذا التباين الفردي يُحسب ويُستبعد تماماً، مما يقلص حجم الخطأ المتبقي ويرفع دقة تقدير أثر المعالجة التجريبية (τj) إلى أقصى حد ممكن إحصائياً.

1.2 التطبيقات المنهجية في علم النفس والعلوم السلوكية

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

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

علاوة على ذلك، تستفيد أبحاث التعلم، وعلم النفس العصبي، وتوطيد الذاكرة (Memory Consolidation) من هذه المنهجية لتتبع منحنيات النسيان واسترجاع المعلومات عبر فترات زمنية متباعدة (بعد ساعة، بعد يوم، بعد أسبوع). وفي علم النفس الدوائي والتجريبي، يُستخدم هذا التحليل لمقارنة تأثير جرعات مختلفة من عقار عصبي أو منبه ذهني على الأداء الوظيفي لنفس المتطوعين، مع استخدام فترات زمنية فاصلة (Washout Periods) لضمان زوال الأثر المتبقي للدواء قبل الانتقال للمستوى التجريبي التالي، مما يضمن أعلى درجات الصدق الداخلي للتجربة السيكولوجية.

1.3 المزايا الإحصائية والمنهجية للتصميمات المتكررة

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

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

تتمثل الميزة المنهجية الثالثة في الضبط التجريبي المطلق للمتغيرات الدخيلة (Control of Confounding Variables) المرتبطة بخصائص العينة. في التصاميم بين المجموعات، وحتى مع استخدام التعيين العشوائي (Random Assignment)، يظل هناك احتمال دائم لعدم تكافؤ المجموعات في متغيرات مثل مستوى الذكاء، والحالة الاقتصادية، والتركيبة الجينية، والدافعية الفردية. أما في القياسات المتكررة، فتتطابق هذه الخصائص تماماً لأن الشخص يُقارن بذاته، مما يقضي نهائياً على التباين البيني كمتغير مهدد للصدق الداخلي، ويضمن أن التباين المرصود في المتغير التابع يعود بالدرجة الأولى للتعرض التجريبي المنظم.

2. الفروق الجوهرية بين تحليل التباين المستقل وتحليل القياسات المتكررة

2.1 التصميم بين المجموعات مقابل التصميم داخل المجموعات

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

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

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

2.2 تفكيك مصادر التباين في القياسات المتكررة

لفهم الآلية الرياضية العميقة لتحليل التباين للقياسات المتكررة، يجب استيعاب كيفية تفكيك مجموع المربعات الكلي (Total Sum of Squares – SSTotal). في تحليل التباين الأحادي المستقل، ينقسم التباين الكلي إلى مصدرين أساسيين فقط: التباين بين المجموعات (SSBetween) وهو تباين المعالجة، والتباين داخل المجموعات (SSWithin) الذي يُعامل كلياً كخطأ عشوائي غير مفسر. أما في تحليل القياسات المتكررة، فتجري عملية تفكيك إضافية متعددة المراحل تعيد تقسيم التباين داخل المجموعات إلى مكونات دقيقة للغاية.

تتم عملية التقسيم في القياسات المتكررة وفق المنظومة الرياضية الهيكلية التالية:

  • مجموع المربعات الكلي (SSTotal): يمثل المقدار الإجمالي لتشتت كافة الدرجات الفردية حول المتوسط الحسابي العام لجميع المشاهدات في التجربة.
  • التباين بين المفحوصين (SSSubjects / SSBetween-Subjects): يقيس الفروق الإجمالية بين أداء الأفراد بغض النظر عن المعالجات التجريبية، وهو التباين الناتج عن السمات الذاتية المستقرة لكل مشارك.
  • التباين داخل المفحوصين (SSWithin-Subjects): يمثل التباين المتبقي للدرجات الفردية لكل شخص حول متوسطه الخاص، وينقسم بدوره رياضياً إلى:
    • تباين المعالجة التجريبية (SSTreatment / SSConditions): التباين الناتج مباشرة عن التغير في مستويات المتغير المستقل أو فترات القياس الزمني.
    • تباين الخطأ المتبقي (SSError / SSResidual): وهو التفاعل بين المفحوصين والمعالجة التجريبية (Subjects × Treatment Interaction)، ويُمثل التباين العشوائي الحقيقي غير المفسر بعد عزل الفروق الفردية.

تظهر هذه العملية أن مقام النسبة الفائية (F-ratio) في القياسات المتكررة لا يُحسب بقسمة متوسط مربعات المعالجة (MSTreatment) على التباين الإجمالي داخل المجموعات، بل يُقسم حصراً على متوسط مربعات الخطأ المتبقي (MSError):

F = MSTreatment / MSError = (SSTreatment / dfTreatment) / (SSError / dfError)

وبما أن SSError قد جُرد تماماً من SSSubjects، فإن قيمة المقام تصبح أصغر بصورة دراماتيكية مقارنة بالتصاميم المستقلة، مما يؤدي إلى تعظيم قيمة F المحسوبة ورفع الدلالة الإحصائية للنتائج بصورة منهجية متقنة.

2.3 حدود وإمكانات استخدام إكسيل مقارنة بالبرمجيات الإحصائية المتخصصة

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

ومع ذلك، تفرض بيئة إكسيل حدوداً منهجية يجب على الباحث المتقدم إدراكها بوضوح. لا يحتوي إكسيل على أداة مباشرة ومستقلة تحمل اسم “Repeated Measures ANOVA” في واجهته الرسومية، بل يتم تنفيذه بذكاء منهجي عبر توظيف أداة “تحليل التباين ثنائي الاتجاه بدون تكرار” (Two-Factor Without Replication)، حيث يُعامل المشاركون كعامل أول والمعالجات كعامل ثانٍ. هذا التطابق الرياضي يُنتج نفس جدول ANOVA التام، لكنه يتطلب من الباحث فهماً عميقاً لكيفية إعادة تسمية المخرجات وتفسيرها بما يتطابق مع لغة القياس النفسي والتجريبي.

أما القيد الأبرز في إكسيل فيكمن في غياب الفحوصات الإحصائية المؤتمتة للافتراضات المتقدمة؛ إذ لا يقوم إكسيل بحساب اختبار موكلي للكروية (Mauchly’s Test of Sphericity) تلقائياً، ولا يقدم معاملات التعديل الإحصائي المباشرة مثل تصحيح جرينهاوس-جايسر (Greenhouse-Geisser) أو هوينه-فيلدت (Huynh-Feldt) بنقرة زر واحدة. ومع ذلك، يمكن للباحث المتمرس تجاوز هذه القيود كلياً من خلال بناء صيغ إحصائية مخصصة داخل ورقة العمل لفحص تباين الفروق وحساب مؤشرات حجم الأثر والمقارنات البعدية بدقة رياضية متناهية، وهو ما سنستعرضه تفصيلاً في الأقسام اللاحقة من هذا الدليل.

3. الافتراضات الإحصائية والشروط المسبقة للاختبار

3.1 افتراض التوزيع الطبيعي واستقلالية المشاهدات

يستند تحليل التباين للقياسات المتكررة، بوصفه اختباراً معلمياً بارامترياً، إلى جملة من الشروط والافتراضات المسبقة التي يجب فحصها لضمان صدق النتائج وموثوقية القرارات الاستدلالية المتخذة. أول هذه الافتراضات هو التوزيع الطبيعي (Normality)، وتحديداً اعتدالية توزيع درجات المتغير التابع داخل كل مستوى من مستويات المعالجة التجريبية، أو اعتدالية توزيع الفروق بين الأزواج التجريبية في العينات الصغيرة. في البحوث النفسية والسلوكية، يمكن فحص هذا الافتراض عبر تقييم مؤشرات الالتواء (Skewness) والتفرطح (Kurtosis)؛ حيث تُعد البيانات مقبولة التوزيع إذا وقعت قيم الالتواء والتفرطح المعيارية ضمن المجال الموثوق بين -1.96 و +1.96 عند مستوى ثقة 95%.

من الناحية العملية، يتميز اختبار ANOVA بقدر كافٍ من المتانة (Robustness) ضد الانتهاكات الطفيفة للاعتدالية، لا سيما عندما تكون حجوم العينات متساوية عبر كافة الشروط، واستناداً إلى مبرهنة النهاية المركزية (Central Limit Theorem) عند تجاوز حجم العينة 30 مشاركاً. وفي إكسيل، يمكن استخدام دوال التوزيع التكراري وحساب قيم الالتواء عبر الدالة =SKEW() والتفرطح عبر الدالة =KURT() للتحقق من خلو البيانات السلوكية من التشويه التوزيعي الشديد الذي قد يبطل صحة الاختبار المعلمي.

أما افتراض استقلالية الملاحظات بين المشاركين (Independence of Observations Between Subjects)، فيعني أن أداء وسلوك أي مشارك داخل التجربة لا يجب أن يؤثر بأي حال من الأحوال على أداء أو سلوك مشارك آخر. يجب أن تُجمع البيانات من كل فرد بصورة منعزلة ومستقلة تماماً لمنع حدوث التفاعل الاجتماعي أو التأثير المتبادل أثناء الاختبارات السلوكية والمعرفية. تجدر الإشارة هنا إلى أن الاستقلالية مفروضة حصراً بين الأفراد المختلفين، في حين أن الارتباط مسموح به ومطلوب منهجياً بين قياسات الفرد نفسه عبر الزمن أو الظروف التجريبية.

3.2 افتراض الكروية (Sphericity) وأهميته الحرجة

يُعد افتراض الكروية (Sphericity)، والذي يُشار إليه أحياناً بافتراض الدائرية أو المركبة الدائرية، الشرط الإحصائي الأكثر حساسية وخصوصية في تحليل التباين للقياسات المتكررة ذات المستويات الثلاثة فأكثر. ينص هذا الافتراض على تساوي تباينات الفروق بين جميع الأزواج الممكنة من مستويات المتغير المستقل. بتعبير رياضي، إذا كانت لدينا ثلاثة شروط تجريبية (A, B, C)، فإننا نحسب درجات الفروق لكل مشارك بين (A – B) و (A – C) و (B – C)، ويتطلب افتراض الكروية أن تكون التباينات الإحصائية لهذه الفروق الثلاثة متساوية تقريباً في المجتمع الأصلي:

Var(YA – YB) ≈ Var(YA – YC) ≈ Var(YB – YC)

من الضروري التمييز الجوهري بين افتراض تجانس التباين (Homogeneity of Variance) المطلوب في تحليل التباين المستقل وافتراض الكروية في القياسات المتكررة؛ فتجانس التباين يكتفي بطلب تساوي التباين الداخلي لكل مجموعة على حدة، بينما تتطلب الكروية تساوي تباين الفروق الزوجية المشتركة. وتُعد حالة التماثل المركب (Compound Symmetry) حالة خاصة وأكثر صرامة من الكروية، حيث تتطلب تساوي تباينات كافة المستويات وتساوي معاملات الارتباط بين جميع أزواج الشروط التجريبية في آن واحد.

تكمن الخطورة المنهجية الكبرى لانتهاك افتراض الكروية في كونه يؤدي مباشرة إلى تضخم حاد في معدل الخطأ من النوع الأول (α Inflation)؛ أي رفض الفرضية الصفرية وإعلان وجود فروق دالة إحصائياً بينما هي في الواقع فروق وهمية ناتجة عن عدم تجانس تباين الفروق. في حالات الانتهاك الشديد، يمكن أن يرتفع مستوى الدلالة الفعلي من 0.05 الاسمي إلى أكثر من 0.15، مما يُفقد القرارات العلمية موثوقيتها. وفي حين يتم تقييم هذا الافتراض في البرمجيات الكبرى عبر اختبار موكلي (Mauchly’s W)، فإن التعامل معه في إكسيل يتطلب تطبيق تصحيحات تعديل درجات الحرية يدوياً عبر معاملات إبسيلون (ε) لضمان النزاهة الإحصائية التامة.

3.3 مستويات قياس المتغيرات التجريبية

يتطلب التطبيق السليم لتحليل التباين للقياسات المتكررة الالتزام الصارم بطبيعة ومستويات قياس المتغيرات الداخلة في النموذج الإحصائي وفقاً لتصنيفات القياس العلمي المعترف بها في العلوم النفسية والتجريبية. يجب أن يكون المتغير التابع (Dependent Variable) مقاساً على مستوى كمي مستمر، أي مقياس فترة (Interval Scale) أو مقياس نسبة (Ratio Scale). تشمل الأمثلة الشائعة في الأبحاث السلوكية: أزمنة الرجع بالمللي ثانية، الدرجات الخام على مقاييس الشخصية المقننة، مستويات الهرمونات في الدم (مثل الكورتيزول)، أو عدد الاستجابات الصحيحة في مهام الذاكرة العاملة.

في المقابل، يجب أن يكون المتغير المستقل (Independent Variable)، والذي يُطلق عليه أيضاً عامل داخل المفحوصين (Within-Subjects Factor)، متغيراً فئوياً تصنيفياً (Categorical/Nominal) أو ترتيبياً (Ordinal) يتضمن ثلاثة مستويات تجريبية أو أكثر (كأن يتضمن ثلاث فترات زمنية: قبلي، بعدي، تتبعي؛ أو ثلاثة أنواع مختلفة من المنبهات الحسية). وإذا كان المتغير المستقل يقتصر على مستويين اثنين فقط، فإن تحليل التباين للقياسات المتكررة يُصبح مكافئاً رياضياً تماماً لاختبار “ت” للعينات المزدوجة والمترابطة (Paired Samples t-test)، وحينها لا تنطبق مسألة الكروية لأن هناك زوجاً واحداً فقط من الفروق.

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

4. تجهيز وإعداد بيئة العمل في مايكروسوفت إكسيل

4.1 تفعيل حزمة أدوات تحليل البيانات (Analysis ToolPak)

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

لتفعيل الحزمة في بيئة نظام التشغيل ويندوز (Windows)، يتبع الباحث الخطوات المنهجية المتسلسلة التالية:

  • فتح برنامج إكسيل والنقر على علامة تبويب ملف (File) في الزاوية العلوية، ثم النقر على خيارات (Options) في أسفل القائمة الجانبية.
  • من النافذة المنبثقة لخيارات إكسيل، يتم اختيار قسم الوظائف الإضافية (Add-ins) من القائمة الجانبية اليسرى.
  • في أسفل النافذة، وضمن القائمة المنسدلة المسماة إدارة (Manage)، يجب التأكد من اختيار وظائف Excel الإضافية (Excel Add-ins)، ثم النقر على زر انتقال (Go…).
  • في المربع الحواري الصغير الذي يظهر، يتم وضع علامة صح داخل مربع الاختيار المجاور لـ Analysis ToolPak (وكذلك Analysis ToolPak – VBA إذا كان الباحث يخطط لاستخدام الماكرو والأتمتة البرمجية)، ثم الضغط على زر موافق (OK).

بمجرد إتمام هذه الخطوات، ستتم إضافة مجموعة أدوات جديدة تُسمى تحليل (Analysis) في أقصى يمين علامة تبويب بيانات (Data) على الشريط الرئيسي (Ribbon)، وتحتوي على زر تحليل البيانات (Data Analysis). يوفر هذا الزر مدخلاً لأكثر من عشرين خوارزمية إحصائية متقدمة، بما فيها الخوارزمية المركزية التي سنعتمد عليها لإجراء القياسات المتكررة: Anova: Two-Factor Without Replication.

4.2 التوافقية وإعدادات التنسيق عبر إصدارات إكسيل

يعمل مايكروسوفت إكسيل عبر منصات وأنظمة تشغيل متعددة تشمل إصدارات Windows، وmacOS، وتطبيقات السحابة عبر Microsoft 365. من الناحية الحسابية والرياضية، تتطابق المعادلات وخوارزميات الحزمة الإحصائية تماماً عبر جميع هذه الإصدارات دون أدنى تفاوت في دقة الأرقام العشرية ومجموع المربعات. ومع ذلك، توجد بعض الفروق الإجرائية الطفيفة في واجهات المستخدم؛ فعلى نظام Mac، يتم الوصول إلى تفعيل الوظائف الإضافية من خلال قائمة أدوات (Tools) العلوية ثم اختيار Excel Add-ins مباشرة وتفعيل Analysis ToolPak.

من المسائل التقنية البالغة الأهمية التي تواجه الباحثين، مسألة ضبط الفواصل العشرية وتنسيق الأرقام (Decimal Separator Settings). تستخدم بعض الأنظمة الإقليمية العربية أو الأوروبية الفاصلة العادية (Comma ,) كفاصل عشري، بينما تستخدم الأنظمة الأنجلوسكسونية النقطة (Dot .). إذا تم لصق بيانات مفصولة بنقاط في بيئة إكسيل معدة لاستخدام الفواصل، فقد يتعامل البرنامج مع الأرقام كنصوص (Text)، مما يؤدي إلى فشل أداة تحليل البيانات وظهور رسائل خطأ تشير إلى أن “نطاق الإدخال يحتوي على بيانات غير رقمية”.

لتجنب هذا الخطأ الكارثي، يجب فحص إعدادات التنسيق عبر الدخول إلى: ملف ← خيارات ← خيارات متقدمة (Advanced)، والتحقق من قسم “استخدام فواصل النظام” (Use system separators)، وضبط الفاصل العشري بدقة. كما يُنصح دائماً بتحديد خلايا البيانات التجريبية والتأكد من تعيين نوع البيانات صراحة على أنه رقم (Number) مع تحديد عدد الخانات العشرية المناسبة، والتحقق من محاذاة الأرقام إلى اليمين تلقائياً داخل الخلايا، وهو المؤشر البصري الصريح على أن إكسيل يقرأ المدخلات كقيم رياضية قابلة للتحليل.

4.3 بناء ورقة عمل منظمة ومعدة للحسابات الإحصائية

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

قبل البدء بالتحليل، يُنصح بتجهيز مناطق الحسابات الرياضية المساعدة مسبقاً، مثل صفوف حساب المتوسطات والانحرافات المعيارية لكل شرط تجريبي باستخدام الدوال =AVERAGE() و =STDEV.S()، وأعمدة حساب متوسطات أداء الأفراد عبر الشروط. كما يُوصى بشدة باستخدام تقنية التنسيق الشرطي (Conditional Formatting) لتحديد القيم الشاذة والمتطرفة (Outliers) بصرياً؛ كأن يتم تلوين الخلايا التي تزيد أو تقل درجاتها المعيارية (Z-scores) عن ±2.5 أو ±3 انحرافات معيارية بلون مميز لفحصها ومراجعة أصولها التجريبية قبل تشغيل نموذج ANOVA.

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

5. هيكلة وتنظيم مصفوفة البيانات التجريبية في أوراق العمل

5.1 تنسيق البيانات بالشكل العريض (Wide Format)

تتطلب أداة تحليل التباين في مايكروسوفت إكسيل هيكلة البيانات حصراً وفق ما يُعرف في علم البيانات بـ التنسيق العريض (Wide Format)، وهو التنسيق القياسي لتحليل القياسات المتكررة التقليدي، على النقيض من التنسيق الطويل (Long Format) المستخدم في النماذج الخطية المختلطة (Linear Mixed Models) في بعض البرمجيات المتقدمة. في التنسيق العريض، تمثل كل وحدة تجريبية (مشارك/مفحوص) وحدة أفقية مستقلة تماماً، بينما تمثل الشروط التجريبية المتغيرات العمودية.

تخضع مصفوفة البيانات في التنسيق العريض للقواعد الصارمة التالية داخل إكسيل:

  • تخصيص الصفوف (Rows): يُخصص كل صف لمشارك تجريبي واحد فقط (Subject/Participant). الصف يحتوي على المعرف الرقمي للمفحوص (مثل P01, P02) متبوعاً بكافة استجاباته عبر جميع الشروط.
  • تخصيص الأعمدة (Columns): يُخصص كل عمود لمستوى واحد من مستويات المعالجة التجريبية أو لفترة زمنية محددة من فترات القياس (مثل Drug_A, Drug_B, Drug_C أو Pre, Post, Follow_up).
  • الخلية المتقاطعة: تمثل القيمة الفردية التي سجلها مشارك محدد في شرط تجريبي محدد حصراً.
  • تجنب دمج الخلايا (No Merged Cells): يجب الامتناع نهائياً عن دمج الخلايا داخل نطاق البيانات أو في صفوف الترويسة، لأن الدمج يعطل قراءة المصفوفات الرقمية بواسطة خوارزمية إكسيل الإحصائية ويؤدي إلى أخطاء برمجية فورية.

5.2 نمذجة مثال تجريبي في علم النفس السلوكي

لإضفاء الطابع التطبيقي والعملي الكامل على هذا الدليل، سنعتمد نموذجاً تجريبياً واقعياً من علم النفس السلوكي وعلم النفس الدوائي المعرفي. لنفترض أن باحثاً يرغب في دراسة أثر أربعة عقاقير نفسية مختلفة لتحسين التركيز وسرعة المعالجة العصبية (العقار A: دواء وهمي Placebo، العقار B: منبه منخفض الجرعة، العقار C: منبه عالي الجرعة، العقار D: مركب تجريبي جديد) على زمن استجابة المشاركين (مقاساً بالثواني) أثناء أداء مهمة انتباه معقدة. شارك في هذه التجربة خمسة مرضى (N = 5)، وخضع كل مريض للشروط الأربعة كاملة بترتيب عشوائي متوازن (Counterbalanced Order) مع فترات زمنية فاصلة ملائمة.

Raw data in Excel
Raw data in Excel

تم إدخال البيانات في ورقة عمل إكسيل ضمن النطاق A1:E6 كما هو موضح في البناء الجدولي التالي:

المشارك (Subject) العقار A (Placebo) العقار B (Low Dose) العقار C (High Dose) العقار D (Novel Drug)
المريض 1 (P1) 4.5 3.8 3.1 2.9
المريض 2 (P2) 5.2 4.6 3.9 3.4
المريض 3 (P3) 3.9 3.2 2.8
المريض 4 (P4) 6.1 5.4 4.7 4.1
المريض 5 (P5) 4.8 4.1 3.5 3.0

عند بناء هذا الجدول في إكسيل، تحتوي الخلية A1 على تسمية المشاركين، بينما تحتوي الخلايا من B1 إلى E1 على تسميات مستويات المتغير المستقل (أسماء العقاقير). وتتوزع معرفات المشاركين في العمود A من الصف 2 إلى الصف 6، بينما تشغل البيانات الرقمية الخام النطاق المستطيل الصافي من B2 إلى E6، وهو النطاق الرياضي المتكامل الذي ستطبق عليه خوارزمية التحليل الإحصائي.

5.3 تدقيق وتنظيف البيانات السلوكية قبل التحليل

تُعد مرحلة تدقيق وتنظيف البيانات (Data Cleaning and Screening) خطوة إلزامية تسبق أي تحليل معلمي. في تصميمات القياسات المتكررة، يمثل التعامل مع البيانات المفقودة (Missing Data) تحدياً حرجاً؛ فخوارزمية ANOVA في إكسيل تعتمد على مصفوفة مستطيلة متكاملة تماماً متساوية الأبعاد، ووجود حتى خلية فارغة واحدة داخل نطاق البيانات سيمنع تنفيذ الأداة أو يؤدي إلى حسابات خاطئة تماماً.

إذا فقد أحد المشاركين في إحدى جلسات القياس التجريبية، يكون أمام الباحث إجرائياً خياران:

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

يتضمن التدقيق أيضاً التأكد القاطع من عدم وجود أي رموز نصية أو مسافات خفية (Spaces) داخل خلايا الأرقام. يمكن استخدام الدالة =ISNUMBER() في عمود فحص جانبي للتحقق من أن كافة المدخلات تُقرأ كقيم رقمية صريحة. كما يُنصح بفحص الاتساق الداخلي والمنطقي للبيانات للتأكد من أن جميع أزمنة الاستجابة تقع ضمن النطاقات الفسيولوجية الممكنة عقلياً، واستبعاد أي قياس ناتج عن خلل تقني أو عدم انتباه كامل من المشارك.

6. خطوات تنفيذ تحليل التباين ثنائي الاتجاه بدون تكرار كبديل إجرائي

6.1 الربط المنهجي بين تحليل القياسات المتكررة وأداة Two-Factor Without Replication

يتساءل العديد من الباحثين عن السبب المنهجي والرياضي الذي يدفعنا لاستخدام أداة تحمل اسم تحليل التباين ثنائي الاتجاه بدون تكرار (Anova: Two-Factor Without Replication) داخل إكسيل لإجراء تحليل يُفترض أنه أحادي الاتجاه للقياسات المتكررة (One-Way RM-ANOVA). يكمن السر في البنية الجبرية المتطابقة تماماً بين النموذجين الرياضيين؛ حيث يُعامل النموذج في القياسات المتكررة “المفحوصين” كعامل تصنيفي أول (Factor A)، بينما تُمثل “المعالجات التجريبية” العامل التصنيفي الثاني (Factor B).

تسمية “بدون تكرار” (Without Replication) تعني بدقة أنه لا توجد سوى ملاحظة رقمية واحدة فقط في كل خلية ناتجة عن تقاطع مفحوص معين مع شرط تجريبي معين (أي مريض واحد × عقار واحد = زمن استجابة واحد). في هذه الحالة، تتلاشى إمكانية حساب تباين الخطأ الداخلي للخلية الواحدة، ويتحول التفاعل الإحصائي بين المفحوصين والمعالجات (Subjects × Treatments) مباشرة وبصورة جبرية تامة إلى تباين الخطأ المتبقي (Residual Error) للنموذج، وهو التباين نفسه المستخدم لاختبار دلالة أثر المعالجة في القياسات المتكررة.

بناءً على ذلك، فإن مخرجات أداة Anova: Two-Factor Without Replication في إكسيل تقدم تفكيكاً ثلاثياً متكاملاً يعزل تباين الصفوف (المفحوصين) وتباين الأعمدة (المعالجات) وتباين الخطأ المتبقي، مما يمنحنا نفس القيم العددية الدقيقة لمجموع المربعات (SS)، ودرجات الحرية (df)، ومتوسط المربعات (MS)، وقيمة النسبة الفائية (F)، والقيمة الاحتمالية (P-value) التي تُخرجها البرمجيات المتخصصة مثل SPSS تحت قائمة General Linear Model ← Repeated Measures.

6.2 تحديد النطاقات وضبط معاملات الإدخال في إكسيل

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

  • التوجه إلى الشريط العلوي والنقر على علامة تبويب بيانات (Data).
  • النقر على زر تحليل البيانات (Data Analysis) في أقصى يمين الشريط لفتح نافذة الأدوات الإحصائية.
  • التمرير داخل القائمة واختيار الخيار الأول: Anova: Two-Factor Without Replication، ثم النقر على موافق (OK).
Repeated measures ANOVA in Excel
Repeated measures ANOVA in Excel

ستظهر نافذة المعاملات الخاصة بالاختبار، ويتم ضبط الإعدادات بدقة فائقة كالتالي:

  • نطاق الإدخال (Input Range): نقوم بتحديد مصفوفة البيانات كاملة بما في ذلك صف الترويسة وعمود أسماء المشاركين. في مثالنا التطبيقي، نكتب أو نحدد النطاق $A$1:$E$6.
  • التسميات (Labels): يجب تفعيل هذا الخيار بوضع علامة صح داخل المربع. هذا الإجراء الحيوي يخبر إكسيل بأن الصف الأول يحتوي على أسماء المعالجات والعمود الأول يحتوي على أسماء المشاركين، مما يمنع التعامل معها كبيانات رقمية ويضمن تسمية المخرجات بوضوح.
  • مستوى الدلالة ألفا (Alpha): يترك افتراضياً عند القيمة المعيارية 0.05 (وهو احتمال ارتكاب الخطأ من النوع الأول المقبول في العلوم السلوكية)، أو تغييره إلى 0.01 إذا كانت الدراسة تتطلب صرامة إحصائية أعلى.

6.3 اختيار نطاق الإخراج ومعالجة رسائل الخطأ الشائعة

في الجزء السفلي من نافذة إعدادات ANOVA، يحدد الباحث خيارات الإخراج (Output options) المفضلة لعرض النتائج. يتوفر للباحث ثلاثة خيارات أساسية:

  • نطاق الإخراج (Output Range): تحديد خلية بداية فارغة في نفس ورقة العمل الحالية (مثل الخلية $G$1) لعرض الجداول بجوار البيانات الأصلية لتسهيل المقارنة المباشرة.
  • ورقة عمل جديدة (New Worksheet Ply): إنشاء ورقة جديدة وتسميتها تلقائياً لعرض الجداول فيها (وهو الخيار الأكثر تنظيماً ويُنصح به للأبحاث الكبيرة).
  • مصنف جديد (New Workbook): فتح ملف إكسيل جديد بالكامل وتصدير المخرجات إليه.
Repeated measures ANOVA in Excel
Repeated measures ANOVA in Excel

عند النقر على زر “موافق”، قد تظهر أحياناً بعض رسائل الخطأ البرمجية الشائعة. من أبرز هذه الرسائل: Input range contains non-numeric data (نطاق الإدخال يحتوي على بيانات غير رقمية). يحدث هذا الخطأ عادة بسبب نسيان تفعيل خيار Labels مع تحديد نصوص الترويسة، أو بسبب وجود مسافات نصية غير مرئية في خلايا الأرقام. كما قد تظهر رسالة Anova: Two-Factor Without Replication requires at least two rows and two columns of data إذا تم تحديد نطاق غير صحيح لا يشمل المصفوفة الكاملة. بعد التحقق من هذه النقاط والنقر على “موافق”، يقوم إكسيل فورياً بتوليد جداول المخرجات الإحصائية المتكاملة.

7. تفكيك وفهم جدول مخرجات تحليل التباين (ANOVA Output)

7.1 قراءة وتفسير جدول الملخص الإحصائي (Summary Table)

يقوم إكسيل بإنشاء جدولين متتاليين في تقرير المخرجات. الجدول الأول هو جدول الملخص الإحصائي (SUMMARY)، والذي يقدم وصفاً تفصيلياً شاملاً للبيانات مجمعة حسب الصفوف (المشاركون) وحسب الأعمدة (المعالجات التجريبية). يتضمن هذا الجدول أربعة مؤشرات وصفية أساسية لكل مستوى: التعداد (Count)، المجموع (Sum)، المتوسط الحسابي (Average)، والتباين الإحصائي (Variance).

يتيح فحص هذا الجدول للباحث استخلاص المؤشرات الاستكشافية الأولية التالية:

  • ملخص الصفوف (Rows Summary – المفحوصون): يعرض متوسط وتباين أداء كل مريض عبر الشروط الأربعة مجتمعة. نلاحظ تباين المتوسطات العامة للمرضى (مثلاً المريض P4 متوسطه الإجمالي مرتفع 5.075 ثانية، بينما المريض P3 متوسطه أسرع 3.225 ثانية)، وهو ما يعكس بصرياً وجود فروق فردية واضحة بين قدرات المشاركين الأساسية تبرر استخدام تصميم القياسات المتكررة لعزلها.
  • ملخص الأعمدة (Columns Summary – العقاقير): يعرض متوسط زمن الاستجابة لكل عقار عبر كافة المرضى. نلاحظ تدرج المتوسطات: العقار A (4.90 ثانية)، العقار B (4.22 ثانية)، العقار C (3.52 ثانية)، العقار D (3.04 ثانية). هذا الانخفاض التدريجي المنتظم في أزمنة الاستجابة يمنح مؤشراً وصفياً قوياً على فاعلية العقاقير المتزايدة في تسريع المعالجة المعرفية.

7.2 تحليل عناصر جدول التباين الرئيسي (ANOVA Table)

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

Repeated measures ANOVA output in Excel
Repeated measures ANOVA output in Excel

لتوضيح ذلك بدقة، نستعرض جدول المخرجات الفعلي الذي تم حسابه من مصفوفة بيانات أزمنة الاستجابة:

مصدر التباين (Source of Variation) مجموع المربعات (SS) درجات الحرية (df) متوسط المربعات (MS) قيمة F المحسوبة القيمة الاحتمالية (P-value) قيمة F الحرجة (F crit)
الصفوف – المفحوصون (Rows) 8.288 4 2.072 141.27 3.62E-10 3.259
الأعمدة – المعالجة (Columns) 10.378 3 3.459 235.86 6.48E-12 3.490
الخطأ المتبقي (Error) 0.176 12 0.0147
المجموع الكلي (Total) 18.842 19

تخضع العلاقات الرياضية داخل هذا الجدول لقوانين جبرية صارمة:

  • درجات الحرية (Degrees of Freedom – df): لصفوف المشاركين dfRows = n – 1 = 5 – 1 = 4. لأعمدة المعالجات dfColumns = k – 1 = 4 – 1 = 3. لتباين الخطأ dfError = (n – 1)(k – 1) = 4 × 3 = 12. والمجموع الكلي dfTotal = N – 1 = 20 – 1 = 19.
  • متوسط المربعات (Mean Squares – MS): ينتج دائماً عن قسمة مجموع المربعات (SS) على درجات الحرية المقابلة له (MS = SS / df). فمثلاً للأعمدة: MSColumns = 10.378 / 3 = 3.459.
  • النسبة الفائية (F-ratio): تنتج عن قسمة متوسط مربعات المعالجة على متوسط مربعات الخطأ: F = 3.459 / 0.01467 = 235.86.

7.3 مطابقة مخرجات إكسيل مع مصطلحات القياس النفسي

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

  • صف Rows في إكسيل: يُطابق تماماً التباين بين المفحوصين (Between-Subjects Effect / Subjects Variance). يعكس هذا الصف الفروق الفردية الذاتية بين المشاركين، وقيمته الإحصائية تؤكد ما إذا كان المرضى يختلفون جوهرياً عن بعضهم في الأداء العام.
  • عمود Columns في إكسيل: يُطابق تماماً التباين داخل المفحوصين العائد للمعالجة التجريبية (Within-Subjects Treatment Effect). هذا هو الصف الجوهري والأساسي في البحث، حيث تختبر قيمته الفرضية العلمية للدراسة حول ما إذا كانت أنواع العقاقير تؤثر في زمن الاستجابة.
  • صف Error في إكسيل: يُطابق تماماً تباين الخطأ المتبقي / التفاعل داخل المفحوصين (Residual Error / Subjects × Treatment Variance). يمثل التشتت العشوائي غير المفسر، وهو المقام المعتمد لاختبار دلالة المعالجة.

8. تفسير الدلالة الإحصائية واتخاذ القرار العلمي

8.1 اختبار الفرضيات الإحصائية وحساب القيمة الفائية (F-Statistic)

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

  • الفرضية الصفرية (H0): تنص على عدم وجود أي فروق حقيقية بين متوسطات أزمنة الاستجابة للعقاقير المختلفة في المجتمع الأصلي؛ أي أن (μDrug A = μDrug B = μDrug C = μDrug D)، وأن أي اختلافات ظاهرية تعود للمصادفة العشوائية.
  • الفرضية البديلة (H1): تنص على أن متوسطين على الأقل من متوسطات المعالجات التجريبية يختلفان عن بعضهما اختلافاً حقيقياً دالاً إحصائياً.

لاتخاذ القرار العلمي بقبول أو رفض الفرضية الصفرية، نتبع قاعدتين متكافئتين:

القاعدة الأولى (مقارنة القيم الفائية): نقارن قيمة F المحسوبة للأعمدة بقيمة F الحرجة (F crit) المستخرجة عند مستوى دلالة α = 0.05 ودرجات حرية (3, 12). من جدول المخرجات، نجد أن F المحسوبة = 235.86، وهي أعلى بمراحل شاسعة من القيمة الحرجة البالغة 3.49. بناءً على ذلك، تقع قيمة الاختبار عميقاً في منطقة الرفض الحرجة، مما يوجب رفض الفرضية الصفرية رفضاً قاطعاً.

أما بالنسبة لصف المفحوصين (Rows)، فرغم أن قيمة F المحسوبة له بالغة الارتفاع (F = 141.27, p < 0.001) مما يدل على تباين هائل بين قدرات المرضى، إلا أن هذا الأثر لا يُفسر كأثر تجريبي في دراستنا، بل يؤكد صحة القرار المنهجي بعزل هذا التباين الضخم خارج حدود خطأ التجربة.

8.2 تفسير القيمة الاحتمالية (P-value)

تُعد القيمة الاحتمالية (P-value) المقياس الأكثر استخداماً في التعبير عن الدلالة الإحصائية. في مخرجات إكسيل، قد تظهر القيمة الاحتمالية الخاصة بالأعمدة مكتوبة بالترميز العلمي المعتمد على الحرف E (Scientific E-Notation) عندما تكون القيمة بالغة الصغر، مثل: 6.48E-12. يعني هذا التدوين الرياضي ضرب الرقم في 10 مرفوعة للقوة -12، أي بتحريك الفاصلة العشرية 12 خانة إلى اليسار:

P = 6.48 × 10-12 = 0.00000000000648

تخضع قواعد اتخاذ القرار وفق القيمة الاحتمالية للمعايير التالية:

  • إذا كانت P-value ≤ 0.05: نرفض الفرضية الصفرية ونستنتج وجود أثر دال إحصائياً عند مستوى ثقة 95%.
  • إذا كانت P-value ≤ 0.01: نرفض الفرضية الصفرية بدلالة إحصائية عالية جداً عند مستوى ثقة 99%.
  • إذا كانت P-value ≤ 0.001: نرفض الفرضية الصفرية بدلالة قطعية فائقة (Highly Significant).

في تطبيقنا العملي، بما أن P-value (6.48 × 10-12) أصغر بكثير من عتبة الدلالة 0.001، فإننا نرفض الفرضية الصفرية ونقر بوجود أثر فائق الدلالة الإحصائية لنوع العقار النفسي على سرعة الاستجابة المعرفية.

8.3 صياغة الاستنتاج السيكولوجي بناءً على القرار الإحصائي

لا تنتهي مهمة الباحث بمجرد إعلان الدلالة الإحصائية الرقمية، بل يجب ترجمة هذا القرار الرياضي إلى استنتاج سيكولوجي ذي معنى علمي ونظري. يؤكد القرار الإحصائي برفض الفرضية الصفرية أن الفروق الملاحظة في أزمنة الاستجابة بين عقاقير (A, B, C, D) ليست وليدة تذبذبات عشوائية في العينة، بل هي نتاج مباشر للتأثيرات الدوائية والعصبية المتباينة لهذه المركبات على الجهاز العصبي المركزي وسرعة المعالجة الحسية الحركية لدى المرضى.

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

9. حساب حجم الأثر والقوة الإحصائية يدوياً في إكسيل

9.1 حساب مربع إيتا الجزئي (Partial Eta Squared – η²p)

يُعد مربع إيتا الجزئي (Partial Eta Squared – η²p) المقياس الأكثر شيوعاً وتوثيقاً في الأبحاث النفسية لتقدير حجم الأثر في تصميمات القياسات المتكررة، وهو المقياس الافتراضي المعتمد في مجلات الجمعية الأمريكية لعلم النفس (APA). يقيس هذا المؤشر نسبة التباين في المتغير التابع التي تُعزى مباشرة إلى المعالجة التجريبية بعد استبعاد وضبط تباين الفروق الفردية بين المفحوصين.

يُحسب مربع إيتا الجزئي من مخرجات جدول ANOVA باستخدام المعادلة الرياضية القياسية التالية:

η²p = SSTreatment / (SSTreatment + SSError)

لتطبيق هذه المعادلة مباشرة في إكسيل، إذا كانت قيمة SS للأعمدة موجودة في الخلية B14 وقيمة SS للخطأ في الخلية B15، نكتب الصيغة الرياضية التالية في خلية حساب مستقلة:

=B14 / (B14 + B15)

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

η²p = 10.378 / (10.378 + 0.176) = 10.378 / 10.554 = 0.983

وفقاً لمعايير جاكوب كوهين (Jacob Cohen) المرجعية لتفسير مربع إيتا الجزئي في العلوم السلوكية:

  • أثر صغير (Small Effect): η²p = 0.01 (يفسر 1% من التباين المتبقي).
  • أثر متوسط (Medium Effect): η²p = 0.06 (يفسر 6% من التباين المتبقي).
  • أثر كبير (Large Effect): η²p = 0.14 (يفسر 14% أو أكثر من التباين المتبقي).

تشير القيمة البالغة 0.983 في تجربتنا إلى حجم أثر هائل واستثنائي؛ حيث يفسر نوع العقار التجريبي ما نسبته 98.3% من التباين الكلي القابل للتفسير في سرعة الاستجابة بعد تحييد الفروق الفردية بين المرضى.

9.2 حساب معامل أوميغا المربع الجزئي (Partial Omega Squared – ω²p)

على الرغم من الانتشار الواسع لمربع إيتا الجزئي، إلا أنه يُعد إحصائياً مقدراً متحيزاً نحو الأعلى (Biased Estimator)؛ حيث يميل لتضخيم حجم الأثر الحقيقي في المجتمع، خاصة في الدراسات ذات العينات الصغيرة جداً (مثل دراستنا الحالية N = 5). لذلك، يوصي كبار الإحصائيين بحساب معامل أوميغا المربع الجزئي (Partial Omega Squared – ω²p) بوصفه تقديراً غير متحيز وأكثر تحفظاً ودقة لخصائص المجتمع الأصلي.

تُعطى المعادلة الرياضية لحساب أوميغا المربع الجزئي بالتالي:

ω²p = [dfTreatment × (MSTreatment – MSError)] / [dfTreatment × MSTreatment + (N – dfTreatment) × MSError]

حيث تُمثل N إجمالي عدد المشاهدات الكلية في المصفوفة (في مثالنا: 5 مفحوصين × 4 شروط = 20 مشاهدة).

يمكن بناء هذه الصيغة المركبة في خلية إكسيل بالربط المباشر مع خلايا الجدول كالتالي:

=(C14*(D14-D15)) / (C14*D14 + (20-C14)*D15)

بالتعويض العددي المباشر من مخرجاتنا:

ω²p = [3 × (3.459 – 0.01467)] / [3 × 3.459 + (20 – 3) × 0.01467] = [3 × 3.44433] / [10.377 + 0.24939] = 10.333 / 10.626 = 0.972

تؤكد القيمة الناتجة (0.972) أن الأثر التجريبي للعقاقير يظل بالغ القوة حتى بعد تطبيق أكثر التصحيحات الإحصائية تحفظاً وتجريداً للتحيز العيني.

9.3 تقييم القوة الإحصائية للتصميم التجريبي

ترتبط القوة الإحصائية (Statistical Power: 1 – β) ارتباطاً وثيقاً بحجم العينة وحجم الأثر الفعلي ومستوى ألفا المحدد. في تصميمات القياسات المتكررة، وبفضل عزل تباين الأفراد، تقترب القوة الإحصائية المحققة في مثل هذه التجارب ذات الأثر الكبير من القيمة القصوى (1.00 أو 100%). هذا يعني أن احتمالية ارتكاب خطأ من النوع الثاني (الفشل في اكتشاف الأثر الحقيقي للعقار) تكاد تنعدم تماماً.

يمكن تقدير القوة الإحصائية يدوياً في إكسيل من خلال حساب معلمة اللامركزية (Non-centrality Parameter – λ) لتوزيع F غير المركزي باستخدام المعادلة: λ = df_treatment × F_observed، ثم استخدام دالة التوزيع التراكمي العكسي لتقدير المساحة الاحتمالية تحت المنحنى. يُعد التقرير المتزامن لكل من قيمة F، ومستوى الدلالة P، وحجم الأثر η²p، والقوة الإحصائية الركيزة المعيارية التي تطلبها كبرى الدوريات العلمية المحكمة لضمان استيفاء شروط النزاهة والشفافية في التقارير البحثية.

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

10.1 مبررات إجراء المقارنات البعدية (Post-Hoc Tests)

يُطلق على اختبار تحليل التباين الكلي اسم الاختبار الإجمالي أو الشامل (Omnibus Test)؛ والسبب في ذلك أنه يخبرنا فقط بوجود فرق دال إحصائياً في مكان ما بين المتوسطات الأربعة، لكنه يعجز بنيوياً عن تحديد أي الأزواج المحددة هي المسؤولة بدقة عن هذا الفرق. هل يختلف العقار A عن B فقط؟ أم أن C يختلف عن D؟ أم أن كافة العقاقير تختلف عن بعضها البعض جوهرياً؟

للإجابة عن هذه التساؤلات التفصيلية، لا بد من إجراء المقارنات البعدية (Post-Hoc Pairwise Comparisons). غير أن إجراء مقارنات ثنائية متعددة باستخدام اختبار “ت” التقليدي دون ضبط إحصائي يؤدي إلى معضلة منهجية خطيرة تُعرف بـ تضخم معدل الخطأ العائلي (Familywise Error Rate Inflation). إذا أجرينا 6 مقارنات ثنائية مستقلة عند مستوى α = 0.05 لكل منها، فإن الاحتمال التراكمي لارتكاب خطأ من النوع الأول لمرة واحدة على الأقل يُحسب بالمعادلة:

αFamilywise = 1 – (1 – α)c = 1 – (1 – 0.05)6 = 1 – 0.735 = 0.265

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

10.2 تطبيق تصحيح بونفيروني (Bonferroni Correction) باستخدام اختبارات t المزدوجة

يُعد تصحيح بونفيروني (Bonferroni Correction) أحد أكثر الأساليب البعدية موثوقية وشهرة في القياسات النفسية. يقوم المنطق الإحصائي لبونفيروني على تقسيم مستوى الدلالة الاسمي (α = 0.05) على إجمالي عدد المقارنات الثنائية الممكنة (c). عدد المقارنات الممكنة لأربعة مستويات يُحسب عبر التوافيق:

c = k(k – 1) / 2 = 4(3) / 2 = 6 مقارنات ثنائية

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

αAltered = 0.05 / 6 = 0.00833

لتطبيق ذلك عملياً في إكسيل، نقوم بإجراء اختبار t للعينات المزدوجة (t-Test: Paired Two Sample for Means) لكل زوج من الأزواج الستة عبر حزمة Data Analysis، أو باستخدام دالة إكسيل المباشرة =T.TEST(array1, array2, 2, 1) (حيث يمثل الوسيط الأخير “1” اختبار العينات المزدوجة Paired، والوسيط قبل الأخير “2” اختبار الطرفين Two-tailed).

نستعرض نتائج المقارنات البعدية الست المطبقة على بياناتنا في الجدول المنظم التالي:

المقارنة الثنائية (Pair) فرق المتوسطات الحسابية القيمة الاحتمالية لاختبار t عتبة بونفيروني المعدلة القرار الإحصائي والتفسير
العقار A مقابل العقار B 0.68 ثانية 0.000045 0.00833 دال إحصائياً (A أبطأ من B)
العقار A مقابل العقار C 1.38 ثانية 0.000008 0.00833 دال إحصائياً (A أبطأ من C)
العقار A مقابل العقار D 1.86 ثانية 0.000002 0.00833 دال إحصائياً (A أبطأ من D)
العقار B مقابل العقار C 0.70 ثانية 0.000031 0.00833 دال إحصائياً (B أبطأ من C)
العقار B مقابل العقار D 1.18 ثانية 0.000012 0.00833 دال إحصائياً (B أبطأ من D)
العقار C مقابل العقار D 0.48 ثانية 0.000185 0.00833 دال إحصائياً (C أبطأ من D)

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

10.3 حساب أقل فرق معنوي (Fisher’s LSD) عبر صيغ إكسيل

يُمثل اختبار أقل فرق معنوي لفيشر (Fisher’s Least Significant Difference – LSD) خياراً بعدياً بديلاً أكثر حساسية، ويُسمح باستخدامه فقط عندما يكون اختبار ANOVA الإجمالي دالاً إحصائياً (Protected LSD). يعتمد هذا الاختبار على حساب قيمة حدية ثابتة للمسافة بين أي متوسطين؛ فإذا تجاوز الفرق المطلق بين متوسطي أي شرطين هذه القيمة الحرجة، يُعتبر الفرق دالاً إحصائياً.

تُعطى المعادلة الرياضية لقيمة LSD لتصميم القياسات المتكررة كالتالي:

LSD = tcritical × √[ (2 × MSError) / n ]

حيث tcritical هي القيمة التائية الحرجة لطرفين عند مستوى α = 0.05 ودرجات حرية الخطأ (dfError = 12)، و n هو عدد المفحوصين (5).

يمكن استخراج قيمة t الحرجة في إكسيل بالدالة: =T.INV.2T(0.05, 12) والتي تعطي 2.1788. بعد ذلك، نطبق صيغة LSD الكاملة في خلية مستقلة كالتالي:

=T.INV.2T(0.05, 12) * SQRT((2 * 0.01467) / 5)

LSD = 2.1788 × √(0.02934 / 5) = 2.1788 × √(0.005868) = 2.1788 × 0.0766 = 0.1669 ثانية

تُظهر المقارنة المباشرة أن جميع الفروق الفعلية بين متوسطات المعالجات (والتي تراوحت بين 0.48 و 1.86 ثانية) أكبر بكثير من قيمة LSD الحرجة البالغة 0.167 ثانية، وهو ما يعزز استنتاجات بونفيروني ويؤكد تفوق كل مستوى علاجي على المستويات السابقة بفارق حقيقي راسخ.

11. فحص انتهاك افتراض الكروية والبدائل التصحيحية

11.1 التشخيص اليدوي لتباين الفروق بين الأزواج

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

لتطبيق ذلك على بياناتنا ذات الشروط الأربعة، نقوم بإنشاء 6 أعمدة حسابية جديدة تمثل الفروق الزوجية لكل مريض:

  • العمود 1: الفروق بين العقار A والعقار B (A – B) ← تباين الفروق =VAR.S() = 0.0070
  • العمود 2: الفروق بين العقار A والعقار C (A – C) ← تباين الفروق = 0.0070
  • العمود 3: الفروق بين العقار A والعقار D (A – D) ← تباين الفروق = 0.0230
  • العمود 4: الفروق بين العقار B والعقار C (B – C) ← تباين الفروق = 0.0150
  • العمود 5: الفروق بين العقار B والعقار D (B – D) ← تباين الفروق = 0.0330
  • العمود 6: الفروق بين العقار C والعقار D (C – D) ← تباين الفروق = 0.0070

يُظهر الفحص البصري لهذه التباينات أنها متقاربة جداً وتتراوح في نطاق ضيق للغاية (بين 0.007 و 0.033)، مما يعطي مؤشراً أولياً مطمئناً على عدم وجود انتهاك جسيم لافتراض الكروية، وأن مصفوفة التباين والتباين المشترك تتسم بقدر عالٍ من التجانس الدائري.

11.2 مفهوم تصحيح جرينهاوس-جايسر (Greenhouse-Geisser Epsilon)

عندما تشير التباينات إلى وجود تفاوت حاد يهدد الكروية، يلجأ الإحصائيون إلى تعديل درجات الحرية للنسبة الفائية عبر ضربها في معامل تصحيحي يُعرف باسم إبسيلون جرينهاوس-جايسر (Greenhouse-Geisser ε). يتراوح هذا المعامل رياضياً بين حد أدنى نظري يمثل أسوأ سيناريو لانتهاك الكروية ويُعطى بالصيغة 1 / (k - 1) (وفي حالتنا مع 4 شروط: 1 / 3 = 0.333)، وحد أقصى يساوي 1.00 ويمثل الكروية التامة والكاملة.

تتم عملية التعديل الإحصائي لدرجات الحرية بضرب درجات حرية البسط والمقام في قيمة ε:

dfTreatment(adj) = ε × dfTreatment = ε × (k – 1)

dfError(adj) = ε × dfError = ε × (k – 1)(n – 1)

إذا افترضنا سيناريو متحفظاً للغاية واعتبرنا أن قيمة إبسيلون التقديرية لدراستنا هي ε = 0.75، فإن درجات الحرية المعدلة تصبح:

  • درجات حرية البسط المعدلة: 3 × 0.75 = 2.25
  • درجات حرية المقام المعدلة: 12 × 0.75 = 9.00

لإعادة حساب القيمة الاحتمالية الدقيقة بعد التصحيح داخل إكسيل، نستخدم دالة التوزيع الفائي لليمين =F.DIST.RT() بتمرير قيمة F الأصلية مع درجات الحرية الكسرية المعدلة:

=F.DIST.RT(235.86, 2.25, 9.00)

تُنتج هذه الدالة قيمة احتمالية بالغة الصغر تساوي 1.43E-09 (أي 0.00000000143). هذه النتيجة تؤكد بشكل قاطع أن أثر المعالجة التجريبية يظل فائق الدلالة الإحصائية حتى تحت أقسى افتراضات انتهاك الكروية وتعديلات درجات الحرية الكسرية.

11.3 تصحيح هوينه-فيلدت (Huynh-Feldt) كبديل أقل تحفظاً

يُعرف تصحيح جرينهاوس-جايسر بكونه شديد التحفظ (Overly Conservative)، مما قد يقلل من القوة الإحصائية للاختبار ويزيد من احتمالية ارتكاب الخطأ من النوع الثاني في العينات المحدودة عندما تكون قيمة إبسيلون قريبة من 0.75 أو أعلى. كبديل منهجي متوازن، طور الباحثان Huynh و Feldt معامل تصحيح معدل يُعرف باسم تصحيح هوينه-فيلدت (Huynh-Feldt Epsilon – ˜ε).

يقوم تصحيح هوينه-فيلدت بتعديل معامل جرينهاوس-جايسر بناءً على حجم العينة ودرجات حرية النموذج، مما يجعله أقل تحفظاً وأكثر ملاءمة للاستخدام عندما يكون الانتهاك طفيفاً إلى متوسط. في التطبيقات البحثية الرصينة، يُنصح بالاعتماد على تصحيح Greenhouse-Geisser إذا كانت قيمة إبسيلون التقديرية أقل من 0.75، والاعتماد على تصحيح Huynh-Feldt إذا كانت إبسيلون أكبر من 0.75، بينما يُعتمد الاختبار غير المصحح كلياً إذا كانت الكروية مستوفاة تماماً (ε = 1.00).

12. صياغة وتوثيق نتائج التحليل وفق معايير الجمعية الأمريكية لعلم النفس (APA 7th)

12.1 القواعد النمطية لتدوين نتائج ANOVA المتكررة

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

F(dfTreatment, dfError) = [F-Value], p = [P-Value], η²p = [Eta-Squared-Value]

تتضمن القواعد التحريرية الصارمة لـ APA النقاط التالية:

  • كتابة الحروف الدالة على الإحصاءات بخط مائل مائل (Italic): مثل F, p, M, SD, t.
  • تقريب قيم F ومؤشرات حجم الأثر إلى منزلتين أو ثلاث منازل عشرية (مثال: F = 235.86, η²p = .98).
  • عدم وضع صفر قبل الفاصلة العشرية للقيمة الاحتمالية (p) وقيم حجم الأثر التي لا يمكن نظرياً أن تتجاوز الواحد الصحيح (تدوين p < .001 أو p = .042 وليس 0.042).
  • الإبلاغ الإلزامي المتزامن عن المتوسطات الحسابية (M) والانحرافات المعيارية (SD) لكل مستوى تجريبي لمساعدة القارئ على استيعاب اتجاه الفروق.

12.2 تصميم الرسوم البيانية التوضيحية لتمثيل التغيرات المتكررة في إكسيل

يُعد التمثيل البصري الجيد ركيزة أساسية في التقارير الأكاديمية. لتمثيل بيانات القياسات المتكررة في إكسيل وفق شروط النشر العلمي، يُفضل إنشاء مخطط خطي بنقاط علام (Line Chart with Markers) يوضح مسار تغير المتوسطات عبر الشروط التجريبية المختلفة.

لإنشاء هذا المخطط الاحترافي في إكسيل، نتبع الخطوات المنهجية التالية:

  • تحديد صف أسماء العقاقير وصف المتوسطات الحسابية المقابلة لها في جدول الملخص.
  • الذهاب إلى قائمة إدراج (Insert) واختيار مخطط خطي (Line with Markers).
  • تنسيق المخطط وفق قواعد APA: إزالة خطوط الشبكة الخلفية (Gridlines)، تحويل خلفية المخطط إلى اللون الأبيض الشفاف، وتعيين لون الخط ونقاط العلام باللون الأسود أو الرمادي الداكن.
  • إضافة أشرطة الخطأ المعياري (Standard Error Bars): النقر على الخط البياني، ثم النقر على زر (+) لإضافة عناصر المخطط واختيار أشرطة الخطأ (Error Bars)، ثم اختيار مخصص (Custom) وتحديد نطاق قيم الخطأ المعياري المحسوبة مسبقاً (SE = SD / √n) لكل شرط لتمثيل عدم اليقين الإحصائي بدقة متناهية.

12.3 نماذج عملية لصياغة فقرة النتائج الأكاديمية

النموذج التطبيقي لتوثيق نتائج دراستنا الحالية (الدالة إحصائياً):

“أُجري تحليل التباين للقياسات المتكررة أحادي الاتجاه (One-Way Repeated Measures ANOVA) لتقييم أثر نوع العقار النفسي على زمن الاستجابة في مهمة الانتباه المعرفي عبر أربعة شروط تجريبية (العقار A: الدواء الوهمي، العقار B: جرعة منخفضة، العقار C: جرعة عالية، العقار D: المركب التجريبي الجديد). أظهرت النتائج وجود أثر رئيسي فائق الدلالة الإحصائية لنوع العقار على زمن الاستجابة، F(3, 12) = 235.86, p < .001, η²p = .983, ω²p = .972. وتراوحت المتوسطات الحسابية لأزمنة الاستجابة بين 4.90 ثانية (SD = 0.81) في شرط الدواء الوهمي (العقار A)، و 4.22 ثانية (SD = 0.81) في العقار B، و 3.52 ثانية (SD = 0.74) في العقار C، و 3.04 ثانية (SD = 0.65) في العقار D. وأكدت المقارنات البعدية باستخدام تصحيح بونفيروني (αadjusted = .0083) وجود فروق دالة إحصائياً بين جميع أزواج العقاقير (p < .001 لكافة المقارنات)، مما يثبت أن العقار D يؤدي إلى أعلى تسريع ذي دلالة في المعالجة المعرفية مقارنة بكافة الشروط الأخرى.”

النموذج التطبيقي في حالة النتائج غير الدالة إحصائياً:

“أشارت نتائج تحليل التباين للقياسات المتكررة إلى عدم وجود أثر رئيسي دال إحصائياً للبرنامج التدريبي عبر فترات القياس الثلاث (القبلي: M = 45.20, SD = 5.12؛ البعدي: M = 46.10, SD = 4.88؛ التتبعي: M = 45.80, SD = 5.05)، حيث بلغت قيمة F(2, 18) = 1.14, p = .342, η²p = .112. وتشير هذه النتائج إلى أن التغيرات الطفيفة الملاحظة في درجات الأداء عبر الزمن لا تتجاوز حدود التباين العشوائي المتوقع داخل العينة.”

يوضح الجدول الأكاديمي النهائي التالي التلخيص الشامل الموصى به للنشر المكتمل:

الشرط التجريبي (العقار) المتوسط (M) الانحراف المعياري (SD) الخطأ المعياري (SE) مجال الثقة 95% [CI]
العقار A (Placebo) 4.90 0.812 0.363 [3.89, 5.91]
العقار B (Low Dose) 4.22 0.814 0.364 [3.21, 5.23]
العقار C (High Dose) 3.52 0.743 0.332 [2.60, 4.44]
العقار D (Novel Drug) 3.04 0.650 0.291 [2.23, 3.85]

خاتمة

يُشكل تحليل التباين للقياسات المتكررة (Repeated Measures ANOVA) ركيزة منهجية متينة في ترسانة الباحث العلمي، لما يوفره من قدرة فائقة على عزل تباينات الفروق الفردية وتعظيم القوة الإحصائية للتجارب السلوكية والنفسية والبيولوجية. وقد أثبت هذا الدليل المرجعي أن مايكروسوفت إكسيل، عبر توظيف أداة Anova: Two-Factor Without Replication وبناء المعادلات التكميلية المخصصة، يمثل بيئة إحصائية كاملة، دقيقة، وشفافة قادرة على مجاراة أعتى البرمجيات الإحصائية المتخصصة، مع منح الباحث فهماً عميقاً وتفاعلياً لحركة البيانات ومصادر التباين.

إن التطبيق المنهجي الناجح لهذا الاختبار يتجاوز مجرد استخراج قيمة F والقيمة الاحتمالية P-value؛ إذ يقتضي الالتزام الصارم بفحص الافتراضات الإحصائية كالكروية والاعتدالية، وتقدير حجوم الأثر الحقيقية عبر مؤشرات η²p و ω²p، وتطبيق ضوابط بونفيروني الصارمة في المقارنات البعدية، وتوثيق النتائج وفق نسق APA 7th المعياري. يتيح هذا التكامل بين الصرامة الرياضية والوضوح الإجرائي في إكسيل للباحثين والأكاديميين إنتاج أبحاث رصينة تلبي أعلى معايير الجودة والشفافية في المجتمع العلمي الدولي.

المراجع (References)

  • American Psychological Association. (2020). Publication manual of the American Psychological Association (7th ed.). American Psychological Association. https://doi.org/10.1037/0000165-000
  • Cohen, J. (1988). Statistical power analysis for the behavioral sciences (2nd ed.). Lawrence Erlbaum Associates.
  • Field, A. (2018). Discovering statistics using IBM SPSS statistics (5th ed.). SAGE Publications.
  • Greenhouse, S. W., & Geisser, S. (1959). On methods in the analysis of profile data. Psychometrika, 24(2), 95–112. https://doi.org/10.1007/BF02289823
  • Huynh, H., & Feldt, L. S. (1976). Estimation of the Box correction for degrees of freedom from sample data in randomized block and split-plot designs. Journal of Educational Statistics, 1(1), 69–82. https://doi.org/10.2307/1164736
  • Keppel, G., & Wickens, T. D. (2004). Design and analysis: A researcher’s handbook (4th ed.). Pearson Prentice Hall.
  • Maxwell, S. E., & Delaney, H. D. (2004). Designing experiments and analyzing data: A model comparison perspective (2nd ed.). Lawrence Erlbaum Associates.
  • Microsoft Corporation. (2023). Use the Analysis ToolPak to perform complex data analysis in Excel. Microsoft Support. https://support.microsoft.com/
  • Olejnik, S., & Algina, J. (2003). Generalized eta and omega squared statistics: Measures of effect size for some common research designs. Psychological Methods, 8(4), 434–447. https://doi.org/10.1037/1082-989X.8.4.434
  • Tabachnick, B. G., & Fidell, L. S. (2019). Using multivariate statistics (7th ed.). Pearson.

اقتباس هذا المقال

looti, M. (2026, أغسطس 26). كيفية إجراء تحليل التباين للقياسات المتكررة في إكسيل. عرب سايكلوجي. https://arabpsychology.com/statistics/how-to-perform-repeated-measures-anova-excel/
looti, Mohammed. “كيفية إجراء تحليل التباين للقياسات المتكررة في إكسيل.” عرب سايكلوجي, 26 أغسطس 2026, https://arabpsychology.com/statistics/how-to-perform-repeated-measures-anova-excel/.
looti, Mohammed. “كيفية إجراء تحليل التباين للقياسات المتكررة في إكسيل.” عرب سايكلوجي. أغسطس 26, 2026. https://arabpsychology.com/statistics/how-to-perform-repeated-measures-anova-excel/.