يُمثّل المربع الإغريقي اللاتيني أحد أكثر تصاميم القياس والتحكم الإحصائي تطوراً وإحكاماً في مناهج البحث العلمي التجريبي وعلم النفس التجريبي الحديث. يتيح هذا النمط التجريبي المتقدم للباحثين دراسة التأثيرات المستقلة لعامل تجريبي رئيسي، مع ضبط ثلاث مجموعات متمايزة من المتغيرات الدخيلة أو المشوشة في آنٍ واحد بأعلى درجات الكفاءة وأقل استهلاك للموارد. ومن خلال بنيته الرياضية المتوازنة والمتعامدة، يرسخ هذا التصميم جسراً متيناً بين الدقة المنهجية الصارمة وتقليل التكلفة الإجرائية في البحوث الإنسانية والسلوكية المعقدة.
المربع الإغريقي اللاتيني
1. التعريف الموجز والدقيق
المربع الإغريقي اللاتيني هو تصميم تجريبي متعامد متعدد المتغيرات، ينشأ من تراكب مربعين لاتينيين من الرتبة نفسها، بحيث يقترن كل رمز لاتيني بكل رمز إغريقي مرة واحدة فقط عبر كامل الشبكة التجريبية. يُستخدم هذا الترتيب لدراسة متغير تجريبي محدد مع تحييد وعزل ثلاثة متغيرات مشوشة مستقلة (تُمثل بالصفوف، والأعمدة، والحروف الإغريقية، بينما تمثل الحروف اللاتينية المعاملات أو العكس).
في سياق علم النفس القياسي والعلوم السلوكية، يوفر هذا التصميم القدرة على ضبط مصادر التباين غير المرغوب فيها والناجمة عن الفروق الفردية في خصائص العينة، وتأثيرات الترتيب الزمني، وتفاوت السياقات التجريبية أو البيئية، وذلك دون الحاجة إلى تنفيذ تصميم عاملي كامل يتطلب عينات ضخمة قد يستحيل توفيرها ميدانياً أو إكلينيكياً. فهو يمثل حلاً اقتصادياً استراتيجياً يحافظ على القوة الإحصائية وصحة النتائج الداخلية.
2. التأثيل والاشتقاق اللغوي
تعود التسمية الاصطلاحية للمفهوم (Graeco-Latin square) إلى الممارسة الرياضية التاريخية التي دشنها عالم الرياضيات الشهير ليونهارت أويلر في القرن الثامن عشر؛ حيث استخدم الحروف اللاتينية الأبجدية (مثل A، B، C، D) لتمثيل عناصر المجموعة الأولى، بينما استخدم الحروف الإغريقية القديمة (مثل α، β، γ، δ) لتمثيل عناصر المجموعة الثانية. وعند دمج هذين المربعين المتعامدين في شبكة واحدة، تزاوجت الأبجديتان لتشكيل أزواج ثنائية فريدة، ومن هنا نحت المصطلح العلمي المركب في اللغات الأوروبية، وتُرجم بدقة إلى اللغة العربية بـ "المربع الإغريقي اللاتيني" للدلالة على الطبيعة المزدوجة للرموز المستعملة في ترميز المتغيرات.
أما كلمة "مربع" في أصلها اللغوي العربي فمشتقة من الجذر الثلاثي (ر-ب-ع)، وهو ما له أربعة أضلاع متساوية أو مصفوفة تتساوى فيها الأبعاد طولاً وعرضاً، مما يعكس البنية المصفوفية المربعة (k × k) التي تقوم عليها المعالجة الإحصائية الرياضية للتصميم التجريبي، حيث تتطابق درجات مستويات المتغيرات الأربعة المضبوطة في التجربة.
3. النطق والصيغة الصرفية والاستخدام
يُنطق المصطلح باللغة العربية بضم الميم وفتح الراء وتضعيف الباء المفتوحة في كلمة "المُرَبَّع" (اسم مفعول من رَبَّعَ)، وتليها صفة النسبة "الإغريقيّ" بكسر الهمزة وتسكين الغين وكسر الراء وتشديد الياء، ثم "اللاتينيّ" بفتح اللام وكسر التاء والنون المشددة. ويُعامل المصطلح في الكتابات الأكاديمية بوصفه مركباً نعتياً مفرداً مذكراً، وتُستخدم صيغة الجمع السالم والمكسر في حالات التعدد (مثل: المربعات الإغريقية اللاتينية).
شاع استخدام هذا المصطلح في الأدبيات الإحصائية والنفسية كمفهوم تقني متخصص (Jargon) يرتبط ارتباطاً وثيقاً بـ تحليل التباين (ANOVA) وتصميم التجارب المعملية والميدانية المتقدمة، ويُصنف ضمن عائلة "التصاميم كتلية التكافؤ المتقاطعة" (Crossed Block Designs).
4. الشرح المفاهيمي التفصيلي وأبعاده التجريبية
يقوم المربع الإغريقي اللاتيني على مبدأ التعامد الإحصائي التام (Orthogonality)، وهو مبدأ رياضي يضمن أن كل مستوى من مستويات أي متغير من المتغيرات الأربعة يظهر جنباً إلى جنب مع كل مستوى من مستويات المتغيرات الثلاثة الأخرى بالتردد نفسه وبالتساوي الدقيق. هذا التماثل الهندسي الرياضي يلغي تماماً التباين المشترك أو الارتباط المتبادل بين المتغيرات المضبوطة والمتغير المستقل، مما يمنع حدوث ظاهرة الاندماج أو التداخل المربك (Confounding) التي تهدد صدق التجارب النفسية المعملية.
تتمثل الحدود النظرية للتصميم في افتراضه الصارم بعدم وجود تفاعلات ذات دلالة إحصائية (Interaction Effects) بين المتغيرات الخاضعة للدراسة. يُفترض نموذج القياس في المربع الإغريقي اللاتيني أن التأثيرات الملحوظة تجميعية وخطية بحتة (Purely Additive Model)؛ فإذا تفاعل المتغير اللاتيني التجريبي مع أحد متغيرات الحجب (الصفوف أو الأعمدة أو الرموز الإغريقية)، فإن ذلك التفاعل سيندمج مباشرة في حد الخطأ التجريبي أو يشوه التقديرات، مما يؤدي إما إلى تضخيم الخطأ من النوع الثاني وتراجع القوة الإحصائية، أو إلى إصدار استنتاجات خاطئة تماماً حول فاعلية المعالجة.
يشترط التصميم تساوي مستويات جميع المتغيرات الأربعة عند القيمة (k)؛ فإذا كان لدينا أربعة مستويات للمتغير العلاجي النفسي المستقل (k = 4)، فإنه يتحتم بالضرورة توظيف أربعة مستويات للصفوف، وأربعة مستويات للأعمدة، وأربعة مستويات للمتغير الإغريقي. هذا التناظر يفرض قيوداً تصميمية صارمة على الباحث، لكنه بالمقابل يُحقق أعلى كفاءة تقليصية، إذ يتطلب إجراء (k²) من المشاهدات والمحاولات التجريبية فقط، بدلاً من (k⁴) كما يقتضي التصميم العاملي التام، مما يعني تقليصاً استثنائياً في زمن الاختبار، وتكلفة العينة، وإجهاد المفحوصين.
5. التطور التاريخي والمعالم البارزة
يرجع الفضل الأول في ابتكار البنية الرياضية البحتة لهذه المربعات إلى العالم السويسري ليونهارت أويلر في عام 1779، عندما طرح معضلته التاريخية الشهيرة المعروفة بـ "مسألة الضباط الستة والثلاثين"، والتي تساءل فيها عما إذا كان ممكناً ترتيب 36 ضابطاً ينتمون إلى 6 أفواج عسكرية ويحملون 6 رتب مختلفة في مربع أبعاده 6×6، بحيث لا يتكرر أي فوج أو رتبة في أي صف أو عمود. وقد أثبت أويلر استحالة حل هذه المسألة لرتبة المربع 6، وصاغ حدسيته الشهيرة القائلة بعدم وجود مربع إغريقي لاتيني لأي رتبة تُكتب بالصيغة (4n + 2).
ظل هذا المفهوم حبيس النظريات الرياضية البحتة حتى مطلع القرن العشرين، وتحديداً خلال ثلاثينيات القرن، عندما قام عالم الإحصاء والجينات السير رونالد فيشر (Ronald Fisher) بإدخال مفاهيم المربعات اللاتينية وتوسيعاتها الإغريقية اللاتينية إلى حقل التجريب العلمي والزراعي في محطة روثامستد، مؤسساً بذلك الاستخدام التطبيقي الأول لها لضبط مصادر التباين الموضعية.
شهد عام 1959 و1960 تحولاً تاريخياً حاسماً على يد علماء الرياضيات بوز، وشرخاندي، وباركر (Bose, Shrikhande, and Parker)، الذين دحضوا حدسية أويلر بعد قرنين من الزمان، وأثبتوا وجود مربعات إغريقية لاتينية من الرتبة العاشرة وما فوقها، مبرهنين أن الرتبة الوحيدة المستحيلة (بخلاف المربع 2×2) هي الرتبة 6×6 فقط. وخلال النصف الثاني من القرن العشرين، تبنى علماء النفس الرياضي والتجريبي، وفي مقدمتهم وينر (Winer) وإدواردز (Edwards)، هذه التصاميم المتقدمة لتطبيقها على أبحاث الإدراك، والتعلم، والسلوك البشري.
6. الأسس النظرية والمنطلقات المعرفية
يستند المربع الإغريقي اللاتيني في جوهره النظري إلى نموذج التأثيرات الثابتة الخطية (Linear Fixed-Effects Model)، الذي ينص على أن أي ملاحظة س سلوكية أو استجابة نفسية (Y_ijkl) تعود إلى مجموع المكونات التالية: المتوسط العام للأداء، مضافاً إليه التأثير الرئيسي للمتغير التجريبي اللاتيني، وتأثير الصفوف، وتأثير الأعمدة، وتأثير الرموز الإغريقية، وأخيراً مقدار الخطأ العشوائي غير المفسر. هذا الأساس الإبستمولوجي يفترض إمكانية تفكيك السلوك الإنساني بدقة إلى تأثيرات مستقلة قابلة للفرز والعزل.
ينطلق التصميم فلسفياً من مبدأ "التحكم عبر التوازن المنهجي الموضعي" (Local Control) المستمد من نظرية التجريب الكلاسيكية. بدلاً من الاعتماد المطلق على العشوائية الصرفة لمعادلة المتغيرات المشوشة، يقوم المربع الإغريقي اللاتيني بتنظيم التباين تنظيماً قصدياً عبر بنية كتلية متداخلة تلغي الحاجة إلى العينات الهائلة.
يرتبط التصميم بفرع التوافيقيات والرياضيات المنفصلة (Combinatorics)، حيث يتم التعامل مع المتغيرات بوصفها فضاءات احتمالية متعامدة متماثلة تضمن الحصول على أعلى قدر من كفاءة استرجاع درجات الحرية وتخفيض تباين الخطأ التجريبي المتبقي، مما يزيد من حساسية الاختبارات الإحصائية لاكتشاف الفروق الحقيقية الخافتة التي تشيع في أبحاث العلوم النفسية.
7. المكونات الأساسية والأبعاد والتصنيفات
يتشكل تصميم المربع الإغريقي اللاتيني من مجموعة عناصر مصفوفية ورياضية متكاملة تشمل الأبعاد التالية:
- عامل الصفوف (Row Factor): متغير حجب خارجي يمثل مصدراً معروفاً للتباين غير المرغوب، مثل الفئات العمرية للمفحوصين، أو مستويات التحصيل الأكاديمي، أو المجموعات التجريبية المتزامنة.
- عامل الأعمدة (Column Factor): متغير حجب إضافي يمثل غالباً أبعاداً مثل الفترات الزمنية للاختبار (صباحاً، بعد الظهر، مساءً) أو تسلسل ترتيب تقديم المهام النفسية.
- المتغير الإغريقي (α, β, γ, …): متغير حجب ثالث متميز أو متغير مستقل ثانوي متزامن، مثل طبيعة السياق المكاني، أو هوية الفاحص السريري، أو نمط التغذية الراجعة المقدمة للمشاركين.
- المتغير اللاتيني (A, B, C, …): وهو عادة المتغير التجريبي المستقل الخاضع للقياس، مثل أنواع التدخلات العلاجية النفسية، أو جرعات العقاقير السلوكية، أو استراتيجيات التعلم المنظم ذاتياً.
- درجات الحرية (Degrees of Freedom): تُقسم مصفوفة القياس درجات الحرية الإجمالية (k² – 1) بين المتوسط العام، وكل متغير من المتغيرات الأربعة بنصيب (k – 1) لكل واحد منها، بينما تخصص درجات الحرية المتبقية (k – 1)(k – 3) للخطأ العشوائي الصافي.
8. أمثلة وسيناريوهات تطبيقية توضيحية
لتوضيح ذلك تطبيقياً في مختبر أبحاث علم النفس المعرفي، لنفترض رغبة باحث في دراسة أثر أربع استراتيجيات تعليمية معرفية مختلفة (المتغير اللاتيني: A، B، C، D) على سعة الذاكرة العاملة للمفحوصين (متغير الاستجابة). يدرك الباحث وجود مصادر تباين مشوشة خطيرة قد تضلل نتائج البحث إن لم تُضبط بدقة بالغة.
يتم تعيين متغير الحجب الأول (الصفوف) ليمثل أربعة مستويات من الذكاء المقاس (مرتفع جداً، مرتفع، متوسط، منخفض)، وتعيين متغير الحجب الثاني (الأعمدة) ليمثل ترتيب الجلسات المعملية الأربع لتفادي أثر الإجهاد الذهني أو التعلم التراكمي، في حين يمثل المتغير الإغريقي المتعامد نوع المحفز الذهني المستخدم (سمعي، بصري، حركي، حسي مركب: α، β، γ، δ).
تُوزع المعالجات في جدول 4×4 بدقة متناهية؛ بحيث لا يحصل المفحوص في مستوى ذكاء محدد على نفس الاستراتيجية المعرفية ونوع المحفز في نفس الجلسة، ولا يتكرر الاقتران عبر التجربة. يتطلب هذا الإجراء الدقيق 16 خلية اختبارية فقط، في حين كان التصميم العاملي الكلاسيكي يستلزم 256 مجموعة أو خلية مستقلة (4×4×4×4) لبلوغ الضبط الموضعي عينه، مما يعكس توفيراً جباراً في جهد وموارد البحث المعرفي.
9. القياس والأدوات ونماذج التحليل الإحصائي
يخضع جمع البيانات المتولدة عن المربع الإغريقي اللاتيني لاختبار تحليل التباين رباعي الاتجاهات بدون تفاعل (Four-Way ANOVA Without Interaction)، حيث يتم اشتقاق مجموع المربعات الإجمالي (Total Sum of Squares) وتقسيمه حسابياً إلى خمسة مكونات منفصلة تماماً:
- مجموع مربعات الصفوف ($SS_{Row}$) بدرجات حرية ($k – 1$).
- مجموع مربعات الأعمدة ($SS_{Column}$) بدرجات حرية ($k – 1$).
- مجموع مربعات المتغير اللاتيني ($SS_{Latin}$) بدرجات حرية ($k – 1$).
- مجموع مربعات المتغير الإغريقي ($SS_{Greek}$) بدرجات حرية ($k – 1$).
- مجموع مربعات الخطأ التجريبي ($SS_{Error}$) بدرجات حرية ($(k – 1)(k – 3)$).
يُلاحظ إحصائياً أن رتبة المربع الإغريقي اللاتيني يجب ألا تقل عن 4 (أي $k ge 4$)؛ لأن تطبيق التصميم على مربع برتبة $k = 3$ يترك درجات حرية الخطأ التجريبي مساوية للصفر تماماً $((3-1)(3-3) = 0)$، مما يجعل حساب التباين العشوائي واختبار الدلالة بـ (اختبار ف / F-test) أمراً مستحيلاً حسابياً. تُجرى هذه التحليلات بيسر في البرمجيات الإحصائية النفسية المتقدمة مثل SPSS وR من خلال تحديد النموذج الخطي العام (GLM) الخالي من حدود التفاعل الضربِي.
10. التطبيقات والأهمية العملية في الحقول النفسية
يحظى هذا التصميم باهتمام خاص في مجال علم النفس الإكلينيكي والعلاجات السلوكية المعرفية؛ حيث يستخدم لمقارنة فاعلية بروتوكولات تدخل علاجية مختلفة عبر مراحل زمنية متعددة وضمن أقسام استشفائية متباينة وبواسطة معالجين سريريين ذوي خلفيات مختلفة، مع ضمان عدم تداخل كفاءة المعالج أو الترتيب الزمني مع الأثر الحقيقي للتقنية الإكلينيكية المختبرة.
في علم النفس التنظيمي وبيئة العمل، يُطبق المربع الإغريقي اللاتيني لاختبار أثر نظم حوافز وظيفية متنوعة على الرضا والإنتاجية مع حجب عوامل تباين الأقسام الإنتاجية، ونوبات العمل اليومية (صباحي ومسائي وليلي)، والأنماط القيادية للمشرفين، محققاً رؤية دقيقة للمؤسسة بأقل تشويش على بيئة العمل الواقعية.
كما يتألق التصميم في حقل الهندسة الإدراكية وتجربة المستخدم وعلم النفس السيبراني (Cyberpsychology)، حيث يقيس الباحثون سرعة تفاعل وسهولة استخدام واجهات رقمية مختلفة عبر أنواع متعددة من الأجهزة الإلكترونية ومستويات إضاءة مختلفة وفئات مختلفة من الخبرة التقنية لدى المستخدمين، دون إثقال كاهل المشاركين باختبارات مجهدة طويلة.
11. الدلائل والبحوث التجريبية الداعمة
أكدت الأدبيات المنهجية الكلاسيكية، ولا سيما أعمال الباحث فيليب كولز (Cowles) ومستجدات القياس النفسي الرياضي، أن استخدام المربعات المتعامدة يرفع القوة الإحصائية (Statistical Power) للاختبار بمعدل يتراوح بين 30% إلى 60% مقارنة بالتصاميم العشوائية البسيطة التي تتجاهل السيطرة على المتغيرات الحجبية متعددة المصادر.
بينت دراسات علم الأدوية النفسية (Psychopharmacology) المقارنة للتأثيرات المعرفية للمهدئات ومحسنات التركيز، أن استخدام المربع الإغريقي اللاتيني وفّر ضبطاً نموذجياً لتحييد الفروق الفسيولوجية الأساسية للمفحوصين وفترات طرد العقار من الجسم (Washout Periods)، مما مكن الباحثين من استكشاف الفروق الحقيقية الدقيقة بمجموعات اختبار صغيرة جداً مقارنة بما تفرضه التصاميم بين-المفحوصين (Between-Subjects).
أوضحت أبحاث التقويم التربوي والتعلم المنظم المعتمدة على تصميم المربع الإغريقي اللاتيني، أن تحييد متغير صعوبة المادة الدراسية، وخبرة المعلم، وتوقيت الحصة الدراسية أتاح قياساً موضوعياً دقيقاً للأثر المنفرد لتقنيات الذكاء الاصطناعي التوليدي في تنمية التفكير الناقد لدى طلبة المدارس.
12. الاعتبارات الثقافية وعبر الثقافية
يكتسب تطبيق المربع الإغريقي اللاتيني أهمية منهجية مضاعفة في أبحاث علم النفس عبر الثقافي (Cross-Cultural Psychology)؛ حيث تتعدد العوامل الدخيلة المرتبطة باللغة، والطبقة الاجتماعية، ودرجة التحضر، وسياقات التنشئة الاجتماعية للمجتمعات المدروسة. يُمكّن التصميم الباحث من إدراج المجموعات الإثنية أو اللغوية كأحد عوامل الحجب (الصفوف مثلاً)، مما يزيل أثر التباين الجغرافي والثقافي المشوش عند تقييم أثر التدخلات النفسية الأساسية.
ومع ذلك، يجب الحذر المنهجي البالغ من الاختلافات الثقافية في نمط الاستجابة المعرفية أو تأثر المفحوصين بتبدل الباحثين أو بيئات المختبر، إذ قد تتسبب بعض المجتمعات الجمعية في توليد تفاعلات خفية ناجمة عن الرغبة الاجتماعية تتناقض مع فرضية "الجمعية الخطية الخالية من التفاعل" المشترطة في هذا التصميم الرياضي الصارم.
13. الانتقادات والمناقشات والحدود المنهجية
على الرغم من براعة وجمال المربع الإغريقي اللاتيني رياضياً، إلا أنه يواجه انتقادات منهجية متعددة في الأوساط السيكولوجية، أبرزها:
- حتمية غياب التفاعلات (No Interaction Assumption): يُعد هذا الافتراض أكبر كعب أخيل للتصميم؛ فالسلوك البشري بطبيعته مشحون بالتفاعلات المعقدة والمتشابكة. وإذا تداخلت الاستراتيجية المعرفية مثلاً مع الفئة العمرية للمشارك، فإن هذا التصميم يعجز عن رصد التفاعل، بل يؤدي ذلك إلى حشو الخطأ التجريبي بتأثيرات حقيقية مضللة تخفض القوة الإحصائية.
- صعوبة التطبيق العملي والقيود البنيوية: اشتراط تطابق عدد مستويات العوامل الأربعة ($k$) يمثل تحدياً خانقاً؛ فمن النادر جداً في الأبحاث الميدانية الواقعية أن تتساوى مستويات الفئات العمرية مع فترات اليوم مع مستويات التدخل التجريبي بدقة تامة.
- استحالة بعض الرتب الرياضية: كما أثبتت الدراسات، لا يوجد أي مربع إغريقي لاتيني من الرتبة 6×6، مما يعني استحالة استعمال هذا التصميم نهائياً إذا كانت مستويات المتغيرات المراد ضبطها ودراستها تساوي ستة مستويات.
- تضاؤل درجات حرية الخطأ: في المربعات ذات الأبعاد المحدودة مثل $4 imes 4$، تنحصر درجات حرية الخطأ في 3 درجات فقط $((4-1)(4-3) = 3)$، مما يضعف القوة الإحصائية بشكل خطير ما لم يُكرر المربع الإغريقي اللاتيني عدة مرات لمضاعفة حجم العينة ودرجات حرية البواقي.
14. المصطلحات المقارنة والفروق الدقيقة
لتجنب اللبس الشائع في منهجية البحث النفسي، يُقارن المربع الإغريقي اللاتيني بالتصاميم ذات الصلة كالتالي:
- المربع اللاتيني (Latin Square Design): تصميم يضبط متغيرين مشوشين اثنين فقط (صفوف وأعمدة) لدراسة متغير تجريبي واحد، بينما يضبط المربع الإغريقي اللاتيني ثلاثة متغيرات مشوشة مستقلة في آنٍ واحد، لكونه تراكباً متعامداً لمربعين لاتينيين.
- المربع الإغريقي اللاتيني الفائق (Hyper-Graeco-Latin Square): امتداد رياضي أعلى يقوم بتراكب ثلاثة مربعات لاتينية متعامدة أو أكثر لدراسة متغير واحد مع ضبط أربعة متغيرات مشوشة أو أكثر في آن معاً، وتزداد معه الصرامة والتعقيد الحسابي.
- التصميم العاملي الكامل (Full Factorial Design): يدرس كل التوليفات الممكنة بين مستويات المتغيرات ويوفر تحليلاً شاملاً لكافة التفاعلات المتبادلة بينها، بخلاف المربع الإغريقي اللاتيني الذي يضحي بدراسة التفاعلات تماماً في سبيل خفض العينة وضبط التكاليف.
- القطاعات العشوائية الكاملة (Randomized Complete Block Design): تركز على حجب مصدر وحيد فقط للتباين الموضعي غير المرغوب، في حين يتعامل المربع الإغريقي اللاتيني مع ثلاثة أبعاد حجب متقاطعة ومتعامدة هندسياً.
15. خلاصة ومرتكزات أساسية
المربع الإغريقي اللاتيني هو تحفة رياضية ومنهجية تمنح باحثي علم النفس أداة قياس تجريبية فائقة الإحكام لتقليص التحيز وضبط ثلاثة متغيرات مشوشة دفعة واحدة بأقل تكلفة ممكنة من الأفراد والجلسات التجريبية. وتظل نقطة قوته الجوهرية كامنة في التوزيع المتوازن والمتعامد للمستويات التجريبية، بينما تظل نقطة ضعفه ونقده الأساسي تتمثل في افتراضه الصارم الخالي من التفاعلات بين المتغيرات وجمود رتبه المصفوفية. ويظل توظيفه الرشيد مقترناً بالتأكد القبلي من جمعية التأثيرات وامتلاك الباحث لعدد كافٍ من التكرارات لضمان حساسية اختبارات الفروض وصدق القرارات الإحصائية الناتجة.
المراجع
- Cochran, W. G., & Cox, G. M. (1957). Experimental designs (2nd ed.). John Wiley & Sons. https://www.wiley.com/en-us/Experimental+Designs%2C+2nd+Edition-p-9780471162032
- Fisher, R. A. (1935). The design of experiments. Oliver and Boyd. https://digital.library.adelaide.edu.au/dspace/handle/2440/15127
- Kirk, R. E. (2013). Experimental design: Procedures for the behavioral sciences (4th ed.). SAGE Publications. https://us.sagepub.com/en-us/nam/experimental-design/book235437
- Montgomery, D. C. (2017). Design and analysis of experiments (9th ed.). John Wiley & Sons. https://www.wiley.com/en-us/Design+and+Analysis+of+Experiments%2C+9th+Edition-p-9781119113478
- Winer, B. J., Brown, D. R., & Michels, K. M. (1991). Statistical principles in experimental design (3rd ed.). McGraw-Hill. https://archive.org/details/statisticalprinc0000wine