اختبار دونيت

في الإحصاء ، يُعد اختبار دونيت إجراءً للمقارنة المتعددة [ 1 ] طوره الإحصائي الكندي تشارلز دونيت [ 2 ] لمقارنة كل علاج من بين عدد من العلاجات مع عنصر تحكم واحد. [ 3 ] [ 4 ] وتُعرف المقارنات المتعددة مع عنصر تحكم أيضًا باسم مقارنات متعددة إلى واحد.

تاريخ

تم تطوير اختبار دونيت في عام 1955؛ [ 5 ] تم نشر جدول محدث للقيم الحرجة في عام 1964. [ 6 ]

مشكلة المقارنات المتعددة

تحدث مشكلة المقارنات المتعددة، أو تعدد الاختبارات، عند النظر في مجموعة من الاستدلالات الإحصائية في آنٍ واحد، أو عند استنتاج مجموعة فرعية من المعلمات المختارة بناءً على القيم المرصودة. وتتمثل القضية الرئيسية في أي نقاش حول إجراءات المقارنات المتعددة في مسألة احتمالية الخطأ من النوع الأول. وتنشأ معظم الاختلافات بين التقنيات البديلة من اختلاف المناهج المتبعة في كيفية التحكم في هذه الأخطاء. تُعدّ المشكلة تقنية جزئيًا، ولكنها في جوهرها مسألة ذاتية تتعلق بكيفية تعريف معدل الخطأ، ومدى السماح بحد أقصى لمعدل الخطأ. [ 7 ] يُعدّ اختبار دونيت معروفًا على نطاق واسع ويُستخدم بكثرة في إجراءات المقارنات المتعددة للمقارنة المتزامنة، عن طريق تقدير الفترات أو اختبار الفرضيات، بين جميع العلاجات الفعّالة ومجموعة ضابطة عند أخذ عينات من توزيع يكون فيه افتراض التوزيع الطبيعي معقولًا. صُمّم اختبار دونيت بحيث يكون معدل الخطأ الإجمالي عند أو أقل منα{\displaystyle \alpha }عند إجراء مقارنات متعددة بين مجموعة العلاج والمجموعة الضابطة. [ 7 ]

استخدامات اختبار دونيت

أُنجز العمل الأصلي حول مشكلة المقارنات المتعددة بواسطة توكي وشيفيه . كانت طريقتهما عامة، إذ أخذت في الاعتبار جميع أنواع المقارنات الثنائية. [ 7 ] تسمح طرق توكي وشيفيه بإجراء أي عدد من المقارنات بين مجموعة من متوسطات العينات. من ناحية أخرى، يقارن اختبار دونيت مجموعة واحدة فقط مع المجموعات الأخرى، ويتناول حالة خاصة من مشكلة المقارنات المتعددة، وهي المقارنات الثنائية بين مجموعات علاجية متعددة ومجموعة تحكم واحدة. في الحالة العامة، حيث نقارن كل زوج من الأزواج، نقوم بـك(ك-1)/2{\displaystyle k(k-1){\big /}2}المقارنات (حيث k هو عدد المجموعات)، ولكن في حالة العلاج مقابل الضوابط، سنجري فقط(ك-1){\displaystyle (k-1)}المقارنات. إذا استخدمنا في حالة مجموعات العلاج والمجموعات الضابطة طريقتي توكي وشيفيه الأكثر عمومية، فقد ينتج عنهما فترات ثقة واسعة بشكل غير ضروري. يأخذ اختبار دونيت في الاعتبار البنية الخاصة لمقارنة العلاج بالضابط، مما ينتج عنه فترات ثقة أضيق. [ 5 ] من الشائع جدًا استخدام اختبار دونيت في التجارب الطبية، على سبيل المثال مقارنة قياسات تعداد الدم في ثلاث مجموعات من الحيوانات، إحداها ضابطة بينما عولجت المجموعتان الأخريان بدوائين مختلفين. ومن الاستخدامات الشائعة الأخرى لهذه الطريقة بين المهندسين الزراعيين: فقد يرغبون في دراسة تأثير بعض المواد الكيميائية المضافة إلى التربة على محصول المحاصيل، لذلك يتركون بعض قطع الأرض دون معالجة (قطع ضابطة) ويقارنونها بقطع الأرض التي أضيفت إليها المواد الكيميائية (قطع معالجة).

وصف رسمي لاختبار دونيت

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

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

الافتراضات

يتناول التحليل الحالة التي تكون فيها نتائج التجربة عددية، وتُجرى التجربة لمقارنة p من المعالجات مع مجموعة ضابطة. ويمكن تلخيص النتائج في مجموعة من(ص+1){\displaystyle (p+1)}المتوسطات المحسوبة لمجموعات المشاهدات،(X0¯،...،Xص¯){\displaystyle ({\bar {X_{0}}},...,{\bar {X_{p}}})}، بينما(X1¯،...،Xص¯){\displaystyle ({\bar {X_{1}}},...,{\bar {X_{p}}})}يشيرون إلى العلاج وX0¯{\displaystyle {\bar {X_{0}}}}يشير ذلك إلى مجموعة الملاحظات الضابطة، وs{\displaystyle s}هو تقدير مستقل للانحراف المعياري المشترك لجميعص+1{\displaystyle p+1}مجموعات من الملاحظات. الكلXأنا¯{\displaystyle {\bar {X_{i}}}}التابعص+1{\displaystyle p+1}يُفترض أن مجموعات الملاحظات موزعة بشكل مستقل وطبيعي مع تباين مشتركσ2{\displaystyle \sigma ^{2}}ويعنيμأنا{\displaystyle \mu _{i}}وهناك أيضاً افتراض بوجود تقدير متاحs2{\displaystyle s^{2}}لσ2{\displaystyle \sigma ^{2}}.

حساب

تعتمد عملية حساب اختبار دونيت على حساب عبارات الثقة حول القيم الحقيقية أو المتوقعة لـص{\displaystyle p}الاختلافاتXأنا¯-X0¯{\displaystyle {\bar {X_{i}}}-{\bar {X_{0}}}}وبالتالي، فإن الفروق بين متوسطات مجموعات العلاج ومتوسطات المجموعة الضابطة. يضمن هذا الإجراء أن احتمالية جميعص{\displaystyle p}البياناتXأنا¯-X0¯{\displaystyle {\bar {X_{i}}}-{\bar {X_{0}}}}أن تكون صحيحاً في آن واحد يساوي قيمة محددة،P{\displaystyle P}عند حساب فاصل الثقة العلوي (أو السفلي) من جانب واحد للقيمة الحقيقية للفرق بين متوسط ​​مجموعة العلاج ومجموعة التحكم ،P{\displaystyle P}يمثل هذا الاحتمال أن تكون هذه القيمة الفعلية أقل من الحد الأعلى (أو أكبر من الحد الأدنى) لتلك الفترة. عند حساب فترة الثقة ذات الجانبين،P{\displaystyle P}يشكل هذا الاحتمال احتمال أن تكون القيمة الحقيقية بين الحدين الأعلى والأدنى.

أولاً، سنرمز إلى المشاهدات المتاحة N بـXأناج{\displaystyle X_{ij}}متىأنا=1،...،ص{\displaystyle i=1,...,p}وج=1،...،شمالأنا{\displaystyle j=1,...,N_{i}}وتقدير التباين المشترك ، على سبيل المثال: s2=أنا=0صج=1شمالأنا(Xأناج-Xأنا¯)2ن{\displaystyle s^{2}={\frac {\sum _{i=0}^{p}\sum _{j=1}^{N_{i}}(X_{ij}-{\bar {X_{i}}})^{2}}{n}}} متىXأنا¯{\displaystyle {\bar {X_{i}}}}هو متوسط ​​المجموعةأنا{\displaystyle i}وشمالأنا{\displaystyle N_{i}}عدد الملاحظات في المجموعةأنا{\displaystyle i}، ون=أنا=0صشمالأنا-(ص+1){\displaystyle n=\sum _{i=0}^{p}N_{i}-(p+1)}درجات الحرية. كما ذكرنا سابقًا، نود الحصول على حدود ثقة منفصلة لكل فرقمأنا-م0،(أنا=1،...،ص){\displaystyle m_{i}-m_{0},(i=1,...,p)}بحيث يكون احتمال أن يكون كلص{\displaystyle p}ستحتوي فترات الثقة على القيم المقابلةمأنا-م0{\displaystyle m_{i}-m_{0}}يساويP{\displaystyle P}.

سننظر في الحالة العامة حيث توجدص{\displaystyle p}مجموعات علاجية ومجموعة ضابطة واحدة. سنكتب:

zأنا=Xأنا¯-X0¯-(مأنا-م0)1شمالأنا+1شمال0{\displaystyle z_{i}={\cfrac {{\bar {X_{i}}}-{\bar {X_{0}}}-(m_{i}-m_{0})}{\sqrt {{\cfrac {1}{N_{i}}}+{\cfrac {1}{N_{0}}}}}}}

دأنا=Xأنا¯-X0¯-(مأنا-م0)s1شمالأنا+1شمال0{\displaystyle D_{i}={\cfrac {{\bar {X_{i}}}-{\bar {X_{0}}}-(m_{i}-m_{0})}{s{\sqrt {{\cfrac {1}{N_{i}}}+{\cfrac {1}{N_{0}}}}}}}}

سنكتب أيضًا:دأنا=zأناs{\displaystyle D_{i}={\frac {z_{i}}{s}}}، والذي يتبع توزيع إحصائية t للطالب بدرجة حرية n . حدود الثقة الدنيا مع معامل الثقة المشتركP{\displaystyle P}لـص{\displaystyle p}آثار العلاجمأنا-م0،(أنا=1،...،ص){\displaystyle m_{i}-m_{0},(i=1,...,p)}سيتم تقديمها من قبل:

Xأنا¯-X0¯-دأناs1شمالأنا+1شمال0،أنا=1،...،ص{\displaystyle {\bar {X_{i}}}-{\bar {X_{0}}}-d_{i}'s{\sqrt {{\frac {1}{N_{i}}}+{\frac {1}{N_{0}}}}},i=1,...,p}

وص{\displaystyle p}الثوابتدأنا{\displaystyle d_{i}'}يتم اختيارهم بحيثPرoب(ت1<د1،...،تص<دص)=P{\displaystyle Prob(t_{1}<d_{1}',...,t_{p}<d_{p}')=P}وبالمثل، ستُعطى الحدود العليا بما يلي:

Xأنا¯-X0¯+دأناs1شمالأنا+1شمال0،أنا=1،...،ص{\displaystyle {\bar {X_{i}}}-{\bar {X_{0}}}+d_{i}'s{\sqrt {{\frac {1}{N_{i}}}+{\frac {1}{N_{0}}}}},i=1,...,p}

لتحديد الحدودمأنا-م0{\displaystyle m_{i}-m_{0}}في كلا الاتجاهين، يمكن أخذ الفترة التالية:

Xأنا¯-X0¯±دأنا"s1شمالأنا+1شمال0،أنا=1،...،ص{\displaystyle {\bar {X_{i}}}-{\bar {X_{0}}}\pm d_{i}''s{\sqrt {{\frac {1}{N_{i}}}+{\frac {1}{N_{0}}}}},i=1,...,p}

متىدأنا"{\displaystyle d_{i}''}يتم اختيارها لإرضاءPرoب(|ت1|<د1،...،|تص|<دص)=P{\displaystyle Prob(|t_{1}|<d_{1}',...,|t_{p}|<d_{p}')=P}الحل لتلك القيم المحددة لـدأنا"{\displaystyle d_{i}''}للاختبار ذي الجانبين ودأنا{\displaystyle d_{i}'}تُعطى القيم لاختبار من جانب واحد في الجداول. [ 5 ] نُشر جدول مُحدّث للقيم الحرجة في عام 1964. [ 6 ]

مثال: قوة تحمل القماش

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

قوة الكسر (بالرطل)
معيارالعملية 1العملية 2العملية 3
155555550
247644944
348645241
وسائل50615245
التباين1927921

يمكن حساب اختبار دونيت باتباع الخطوات التالية:

1. إدخال البيانات مع المتوسطات والتباينات:

  • اجمع القياسات لكل مجموعة (العمليات القياسية وعمليات المعالجة). راجع البيانات في الجدول أعلاه لمعرفة الأرقام الخام والمتوسطات والتباينات لكل مجموعة.

2. حساب التباين المجمعs2{\displaystyle s^{2}}:

552+472+482+552+...+412-3(502+612+522+452)8=1528=19{\displaystyle {\frac {55^{2}+47^{2}+48^{2}+55^{2}+\ldots +41^{2}-3(50^{2}+61^{2}+52^{2}+45^{2})}{8}}={\frac {152}{8}}=19}.

3. حساب الانحراف المعياريs{\displaystyle s}:

  • احسب الجذر التربيعي لمتوسط ​​التباين. على سبيل المثال،
s=19=4.36{\displaystyle s={\sqrt {19}}=4.36}.

4. حساب الخطأ المعياري:

  • الصيغة التالية تعطي الخطأ المعياري لفرق المتوسطين. على سبيل المثال،
s2شمال=4.3623=3.56{\displaystyle s{\sqrt {\frac {2}{N}}}=4.36{\sqrt {\frac {2}{3}}}=3.56}.

5. تحديد القيمة الحرجةت{\displaystyle t}:

  • استخدم جداول دونيت للعثور علىت{\displaystyle t}بالنسبة لدرجات الحرية ومستوى الثقة المحددين. على سبيل المثال،
لص=3{\displaystyle p=3}وdf=8{\displaystyle {\text{d.f.}}=8}:
من جانب واحد:ت=2.42{\displaystyle t=2.42}
ذو وجهين:ت=2.88{\displaystyle t=2.88}

6. يُرمز إلى الكمية التي يجب إضافتها و/أو طرحها من الفروق الملحوظة بين المتوسطات لإعطاء حدود ثقتها بـأ{\displaystyle A}(وقد أطلق توكي على هذا "البدل" )، ويمكن حسابه على النحو التالي:

  • ضاعفت{\displaystyle t}باستخدام الخطأ المعياري لفرق المتوسطين. مثال:
من جانب واحد:
أ=2.42×3.56=8.61{\displaystyle A=2.42\times 3.56=8.61}
ذو وجهين:
أ=2.88×3.56=10.25{\displaystyle A=2.88\times 3.56=10.25}

7. حساب حدود الثقة:

  • احسب حدود الثقة لكل عملية مقارنةً بالمعيار. على سبيل المثال،
حدود أحادية الجانب:
العملية 1:
61-50-8.61=2.39رطل{\displaystyle 61-50-8.61=2.39\,{\text{lbs}}}
العملية الثانية:
52-50-8.61=-6.61رطل{\displaystyle 52-50-8.61=-6.61\,{\text{lbs}}}
العملية 3:
45-50-8.61=-13.61رطل{\displaystyle 45-50-8.61=-13.61\,{\text{lbs}}}
الحدود ذات الجانبين:
العملية 1:
61-50±10.25=[0.75،21.25]رطل{\displaystyle 61-50\pm 10.25=[0.75,21.25]\,{\text{lbs}}}
العملية الثانية:
52-50±10.25=[-8.25،12.25]رطل{\displaystyle 52-50\pm 10.25=[-8.25,12.25]\,{\text{lbs}}}
العملية 3:
45-50±10.25=[-15.25،5.25]رطل{\displaystyle 45-50\pm 10.25=[-15.25,5.25]\,{\text{lbs}}}

8. استخلاص النتائج:

  • استنادًا إلى حدود الثقة المحسوبة، استخلص استنتاجات حول كل عملية مقارنةً بالمعيار. على سبيل المثال،
من جانب واحد:
العملية 1: تتجاوز قوة الكسر المعيار بما لا يقل عن 2.39 رطل.
العملية 2: قوة الكسر لا تتجاوز المعيار (قيمة سالبة).
العملية 3: قوة الكسر لا تتجاوز المعيار (قيمة سالبة).
ذو وجهين:
العملية 1: تتجاوز قوة الكسر المعيار بما يتراوح بين 0.75 رطل و 21.25 رطل.
العملية 2: قوة الكسر تتراوح بين -8.25 رطل و 12.25 رطل (قد تتجاوز المعيار أو لا).
العملية 3: قوة الكسر تتراوح بين -15.25 رطل و 5.25 رطل (قد تتجاوز المعيار أو لا).

مراجع

  1. Upton G. & Cook I. (2006.) قاموس الإحصاء ، الطبعة الثانية، مطبعة جامعة أكسفورد، أكسفورد، المملكة المتحدة.
  2. رومسي، ديبورا (19 أغسطس 2009). الإحصاء 2 للمبتدئين . وايلي. ص 186. تم الاطلاع عليه بتاريخ 22 أغسطس 2012. تم تطوير اختبار دونيت بواسطة. 
  3. إيفريت بي إس وشروندال أ. (2010). قاموس كامبريدج للإحصاء ، الطبعة الرابعة، مطبعة جامعة كامبريدج، كامبريدج، المملكة المتحدة.
  4. "البرمجيات الإحصائية | قسم تكنولوجيا المعلومات بجامعة كنتاكي" . Uky.edu. مؤرشف من الأصل بتاريخ 31 يوليو 2012. تم الاطلاع عليه بتاريخ 22 أغسطس 2012 .
  5. 1 2 3 4 دونيت، سي دبليو (1955). "إجراء مقارنة متعددة لمقارنة عدة علاجات مع مجموعة ضابطة" . مجلة الجمعية الإحصائية الأمريكية . 50 : 1096-1121. doi : 10.1080/01621459.1955.10501294 .
  6. 1 2 Dunnett CW (1964.) "جداول جديدة للمقارنات المتعددة مع عنصر تحكم"، Biometrics ، 20 : 482-491 .
  7. 1 2 3 هاول، ديفيد سي. الأساليب الإحصائية لعلم النفس (الطبعة الثامنة ). 
  8. اختبار دونيت ، هايبرستات أونلاين: كتاب تمهيدي في الإحصاء وبرنامج تعليمي عبر الإنترنت للمساعدة في دورات الإحصاء
  9. آليات الاختبارات المختلفة - الإحصاء الحيوي BI 345 ، مؤرشف بتاريخ 1 يونيو 2010 في أرشيف الإنترنت ، كلية سانت أنسيلم
  10. فيلارز، دونالد ستاتلر (1951). التصميم الإحصائي وتحليل التجارب لبحوث التنمية . دوبوك، أيوا: شركة ويليام سي براون.