لماذا تُرفض الأبحاث المتعلقة بالتوزيع العشوائي المندلي رفضًا قاطعًا؟: حسابات القدرة الإحصائية التي لا يُبلغ عنها أحد تقريبًا

تم النشر 19 سبتمبر 2026

يُعدّ حساب القدرة الإحصائية للتوزيع العشوائي المندلي الخطوة الوحيدة في مسار البحث التي لا تُكلّف شيئًا، ومع ذلك يتمّ تجاهلها. قام روزاس-تشافيز وميريمان بمراجعة 86 دراسة توزيع عشوائي مندلي ثنائية العينة حول فرط حمض يوريك الدم والنقرس، نُشرت بين عامي 2016 و2024، ووجدا أن 59% منها تناولت مسألة القدرة الإحصائية، ما يعني أن 41% منها قدّمت تقديرًا سببيًا، إيجابيًا كان أم سلبيًا، دون توضيح ما إذا كانت الدراسة قادرة على رصد التأثير المطلوب (Gout, Urate, and Crystal Depose Disease 2025;3(2):8, doi:10.3390/gucdd3020008 ). تبدو النتيجة السلبية الناتجة عن تصميم ذي قدرة إحصائية منخفضة، والنتيجة السلبية الناتجة عن غياب التأثير فعليًا، متطابقتين ظاهريًا. ولا يُبيّن للقارئ أيّهما يقرأ إلا من خلال حساب القدرة الإحصائية.

هذه المحاكاة هي الثالثة في سلسلة محاكاة التوزيع العشوائي المندلي، وهي تُكمل ما بدأناه في المقالتين السابقتين. أوضحت مقالتنا حول دلالة إحصائية F أن ضعف المتغير المؤثر يُضخّم عرض فاصل الثقة بدلاً من تقدير النقطة. كما أوضحت مقالتنا حول ما يحميك منه كل مُقدِّر من المُقدِّرات الأربعة أن تعدد التأثيرات الجينية قد يجعل أربع طرق مختلفة تتفق على نفس الإجابة الخاطئة. لم تُجب أيٌّ من المقالتين على سؤالٍ أبسط يطرحه المُراجع قبل أن يكون لأيٍّ من هذين التشخيصين أهمية. بالنظر إلى المتغير المؤثر وعينة النتائج المُتاحة، هل كان من الممكن أن تجد دراستك التأثير الذي تدّعي اختباره؟ لقد أنشأنا نموذج قوة التوزيع العشوائي المندلي لعينتين، وتحققنا من صحته بشكل مستقل باستخدام لغتي بايثون وR، وقمنا بتشغيله على نطاق أحجام الدراسات المُستخدمة في الأدبيات المنشورة.

ما الذي يتحقق منه المراجعون بالفعل الآن؟

انتقلت عملية التوزيع العشوائي المندلي من كونها أسلوبًا متخصصًا إلى بند في قائمة مراجعة المراجعين على مدى عقد تقريبًا، وتتألف هذه القائمة من ثلاثة أجزاء. وقد صاغتها سكريفانكوفا وزملاؤها رسميًا في بيان STROBE-MR، وهو امتدادٌ لبيان STROBE يتألف من 20 بندًا، مصمم خصيصًا لدراسات المتغيرات الآلية الجينية، ونُشر بالاشتراك في مجلتي BMJ وJAMA (JAMA 2021؛ 326(16): 1614-1621، doi:10.1001/jama.2021.18236 ). يطلب البند 7 من المؤلفين وصف كيفية تقييمهم لافتراضاتهم، وتُشير وثيقة الشرح صراحةً إلى إحصائية F كأداة لتقييم قوة المتغير الآلي، وهو موضوع الجزء الأول من بحثنا. أما البند 4ب فيطرح سؤالًا منفصلًا. يتطلب ذلك من المؤلفين الإبلاغ عن حجم العينة وتحديد ما إذا تم إجراء حساب القوة أو حجم العينة قبل التحليل الرئيسي، مع ملاحظة أن القوة تُقرأ بشكل أفضل من عرض فاصل الثقة بنسبة 95 بالمائة بمجرد اكتمال الدراسة (BMJ 2021؛375:n2233، doi:10.1136/bmj.n2233 ).

يُعدّ تدقيق النقرس وفرط حمض اليوريك الفحص الكمّي الوحيد الذي وجدناه لتحديد مدى استيفاء هذا الشرط، ويُظهر النمط في الفحوصات الثلاثة اتساقًا بدلًا من أن يكون مطمئنًا. فقد استوفت 66% من الدراسات الـ 86 شرط الملاءمة المتعلق بقوة الأداة، واستخدمت 51% منها عينات تعرض ونتائج غير متداخلة، وتناولت 59% منها مسألة القدرة الإحصائية، بينما استوفت 17% فقط جميع الافتراضات الأساسية الثلاثة للمتغيرات الآلية في آنٍ واحد. والتزمت 11% من الأبحاث المنشورة بعد عام 2021 بمعايير STROBE-MR التزامًا كاملًا. ولذلك، قد تُشير دراسة ما إلى إحصائية F قوية، وتصميم عينة واحدة دقيق، ومع ذلك لا تُوضّح ما إذا كان لديها أي فرصة لاكتشاف التأثير الذي يناقشه ملخصها.

الصيغة الكامنة وراء حساب قوة التوزيع العشوائي المندلي

قدّم بريون وشاخبازوف وفيشر الإطار الذي لا يزال هذا المجال يستخدمه. انطلاقًا من النظرية التقاربية للتعرض المستمر والمتنبئ الجيني الواحد، بيّنوا أن القوة الإحصائية عند حجم عينة ثابت تختزل إلى كميتين: نسبة التباين في التعرض التي يفسرها المتغير المدروس، والأثر السببي الحقيقي الذي يتم اختباره (المجلة الدولية لعلم الأوبئة 2013؛ 42(5): 1497-1501، doi:10.1093/ije/dyt179 ). وسّع بورغيس هذا النهج ليشمل النتائج الثنائية، وعمل على حجم العينة في الاتجاه المعاكس، من القوة المستهدفة إلى حجم العينة (N) الذي تحتاجه الدراسة (المجلة الدولية لعلم الأوبئة 2014؛ 43(3): 922-929، doi:10.1093/ije/dyu005 ).

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

كيف اختبرناه

قمنا ببناء معادلة حساب قوة التوزيع العشوائي المندلي المغلقة المذكورة أعلاه، وتحققنا منها باستخدام ثلاثة مصادر مستقلة للأدلة، باتباع نفس معايير التحقق التي استخدمناها في الدراستين السابقتين. أولًا، قارنا محاكاة سريعة ومتجهة، تستمد تقديرات التعرض والنتيجة لكل تكرار من توزيعات العينات الخاصة به، مع محاكاة بطيئة تُولّد الأنماط الجينية والتعرضات والنتائج على مستوى الأفراد، وتُجري عمليات الانحدار الفعلية، وذلك على ثلاث مجموعات من المعلمات. وقد اتفقت المعادلتان ضمن هامش خطأ مونت كارلو عند بضع مئات من التكرارات. ثانيًا، قمنا بتشغيل الشبكة الكاملة مرتين، مرة باستخدام بايثون مع NumPy، ومرة ​​أخرى بشكل مستقل باستخدام R، بمعدل 200,000 تكرار لكل خلية، باستخدام مولدات أرقام عشوائية مختلفة ودون استخدام أي شفرة مشتركة. ثالثًا، قارنا كلا المحركين المحاكيين بالمعادلة المغلقة نفسها. يتبع الخطأ المعياري لطريقة دلتا لنسبة والد ذات العينتين، والذي ينشر عدم اليقين من تقديرات التعرض لـ SNP وتقديرات النتائج لـ SNP بدلاً من التعامل مع الأداة على أنها معروفة تمامًا، إطار عمل البيانات الملخصة لبرجس، باتروورث، وتومسون (علم الأوبئة الوراثية 2013؛37(7):658-665، doi:10.1002/gepi.21758 ).

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

يؤكد تحليل النتائج الأولية صحة الصيغة، وليس أي شيء مفاجئ بحد ذاته، وهذا هو الهدف من التحقق من صحتها. عند عينة نتائج قوامها 500,000، وهو نطاق الدراسة الذي تستخدمه معظم تحليلات بنوك البيانات الحيوية الكبيرة حاليًا، تتجاوز القدرة الإحصائية 99% حتى مع أدنى قوة للأداة على الشبكة. أما عند عينة نتائج قوامها 5,000، وهو حجم لا يزال شائعًا في دراسات الحالات والشواهد لنتائج أقل شيوعًا، فلا تصل القدرة الإحصائية إلى 80% في أي مكان على محور R² الذي اختبرناه، بل تبلغ ذروتها عند 60% حتى عندما تفسر أداة واحدة نسبة 10% من تباين التعرض. بين هذين الحدين، يتبادل المحوران بشكل مباشر. تحتاج عينة نتائج قوامها 20,000 إلى أداة تفسر ما يقارب 4 إلى 5% من تباين التعرض للوصول إلى القدرة الإحصائية التقليدية، وهو حد لا تتجاوزه سوى قلة من أدوات SNP الفردية، ويتطلب عادةً تجميع عدة متغيرات ذات دلالة إحصائية على مستوى الجينوم في درجة متعددة الجينات.

ستة تصاميم، حجم تأثير واحد

يُعدّ تحليل المعايير المجردة بمثابة فحص للنموذج. ما يرغب المراجع في معرفته هو ما إذا كان تصميمٌ يبدو عاديًا في المخطوطة المُقدّمة قادرًا على تحقيق النتائج التي يدّعي التوصل إليها. لقد أنشأنا ستة تصاميم دراسات تغطي النطاق المنشور لتحليل الانحدار المندلي ثنائي العينة، ودراسة حالة-مراقبة صغيرة الحجم مع متغير جيني أحادي النوكليوتيد رئيسي واحد، ودراسة مجموعة متوسطة الحجم مع متغير جيني أحادي النوكليوتيد أو مجموع درجات متغيرات جينية متعددة، ودراسة بنك حيوي كبير الحجم، وذلك عبر ثلاث مستويات قوة للأداة. ثمّ حسبنا القدرة الإحصائية لكل تصميم عند ثلاثة أحجام تأثير حقيقية: انحراف معياري متواضع 0.05، وانحراف معياري متوسط ​​0.10، وانحراف معياري كبير 0.20.

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

لماذا يتم نشر اكتشاف ضعيف لـ Null على أنه اكتشاف حقيقي؟

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

هذه هي الآلية الكامنة وراء عدم تطرق 41% من الدراسات في تدقيق روزاس-تشافيز وميريمان إلى مسألة القدرة الإحصائية. فبمجرد أن تُشير ورقة بحثية إلى تقدير نقطي قريب من الصفر وفترة ثقة تتقاطع معه، تُعتبر الورقة مكتملة دون بيان القدرة الإحصائية، ولا يتساءل المؤلفون، ولا معظم المراجعين على ما يبدو، عما إذا كانت الدراسة قد صُممت أصلاً لرصد تأثير بحجم ذي أهمية بيولوجية أو سريرية. وُضع البند 4ب من معيار STROBE-MR تحديدًا لسد هذه الثغرة، إذ يطلب من المؤلفين توضيح ما إذا تم إجراء حساب للقدرة الإحصائية، وتُشير وثيقة الشرح بوضوح إلى أن عرض فترة الثقة بعد إجراء الحساب هو المؤشر الأكثر دقة لما يمكن أن تدعمه الدراسة.

لا يتطلب الحل العملي سوى فقرة واحدة ورقم واحد. قبل إجراء تحليل التوزيع العشوائي المندلي أو تقديمه، قم بإجراء حساب قوة التوزيع العشوائي المندلي للحد الأدنى للتأثير القابل للكشف عند قوة 80%، مع الأخذ في الاعتبار معامل التحديد (R²) المعروف للأداة وحجم عينة النتائج المعروف، باستخدام نفس منهجية معلمات اللامركزية الموضحة في الشكل 1، ثم اذكر هذا الرقم في قسم الطرق بجانب إحصائية F. النتيجة الصفرية المذكورة بجانب حد أدنى للتأثير القابل للكشف يبلغ 0.03 انحراف معياري هي نتيجة حقيقية. أما النتيجة الصفرية المذكورة بجانب حد أدنى للتأثير القابل للكشف يبلغ 0.25 انحراف معياري، بجانب دراسات تشير إلى أن التأثير الحقيقي، إن وُجد، أقرب إلى 0.10 انحراف معياري، فهي قيد تصميمي مُغلّف باستنتاج.

التوتر مع الطرفين الآخرين

يُعدّ رفع القدرة الإحصائية بإضافة المزيد من المتغيرات الجينية إلى الأداة استجابةً بديهيةً لانخفاض معامل التحديد (R²)، ويتفاعل هذا مع نمطي الفشل الآخرين في هذه السلسلة بدلاً من تجنّبهما. فكل متغير جيني إضافي يُدمج في النتيجة يرفع معامل التحديد (R²) وبالتالي القدرة الإحصائية، وهو أمرٌ مفيدٌ للمشكلة المطروحة في هذا المقال، ولكن كل متغير إضافي يُمثّل أيضاً فرصةً إضافيةً للتأثير الجيني المتعدد الأفقي، وهو نمط الفشل الذي حدّده الجزء الثاني من المقال، حيث يمكن لمسارٍ جيني مشترك أن يُحرّك جميع المُقدِّرات الأربعة المشتركة نحو نفس الإجابة الخاطئة ظاهرياً. فالأداة الكبيرة والقوية ليست بالضرورة أداةً صالحة، ولا تقلّ أهمية تشخيصات التأثير الجيني المتعدد المذكورة في ذلك المقال مع ازدياد حجم الأداة.

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

الدعم المهني لدراسة التوزيع العشوائي المندلي

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

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

خاتمة

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

لا يصف أيٌّ من العناصر الثلاثة في هذه السلسلة منهجًا معيبًا. فإحصائية F، والمُقدِّرات الأربعة المقاومة لتعدد التأثيرات، وحساب القدرة الإحصائية، جميعها تؤدي وظيفتها كما تنص عليه نظريتها الإحصائية. لكن ما يجمعها هو أن أيًّا منها لا يُغني عن العنصرين الآخرين، وأنّ ورقة بحثية تُشير إلى عنصر واحد فقط، كإحصائية F أعلى من 10 دون بيان القدرة الإحصائية، أو فرضية العدم الواضحة دون أدنى تأثير قابل للكشف، تكون قد استوفت الشروط الشكلية لقائمة STROBE-MR دون الإجابة على السؤال الذي يطرحه المُراجع المُدقِّق. إنّ ذكر العناصر الثلاثة جميعها لا يتطلب سوى فقرة قصيرة. وإغفال أحدها هو ما يجعل قصور التصميم استنتاجًا منشورًا.

مراجع حسابات

1. بريون إم جيه إيه، شاخبازوف كيه، فيشر بي إم. حساب القوة الإحصائية في دراسات التوزيع العشوائي المندلي. المجلة الدولية لعلم الأوبئة. 2013؛42(5):1497-1501. دوى: 10.1093 / عيي / dyt179
2. بورغيس إس. حسابات حجم العينة وقوة الاختبار في التوزيع العشوائي المندلي مع متغير آلي واحد ونتيجة ثنائية. المجلة الدولية لعلم الأوبئة. 2014؛43(3):922-929. doi:10.1093/ije/dyu005
3. بورغيس إس، باتروورث إيه، طومسون إس جي. تحليل التوزيع العشوائي المندلي مع متغيرات جينية متعددة باستخدام بيانات ملخصة. علم الأوبئة الوراثية. 2013؛37(7):658-665. doi:10.1002/gepi.21758
4. سكريفانكوفا VW، ريتشموند RC، وولف BAR، وآخرون. تعزيز الإبلاغ عن الدراسات الرصدية في علم الأوبئة باستخدام التوزيع العشوائي المندلي: بيان STROBE-MR. JAMA. 2021؛326(16):1614-1621. doi:10.1001/jama.2021.18236
5. سكريفانكوفا VW، ريتشموند RC، وولف BAR، وآخرون. تعزيز الإبلاغ عن الدراسات الرصدية في علم الأوبئة باستخدام التوزيع العشوائي المندلي (STROBE-MR): شرح وتفصيل. المجلة الطبية البريطانية. 2021؛375:n2233. doi:10.1136/bmj.n2233
6. روزاس-تشافيز إف، ميريمان تي آر. دراسات التوزيع العشوائي المندلي: مقياس لتقييم الجودة. مرض ترسب بلورات اليورات النقرسي. 2025؛3(2):8. doi:10.3390/gucdd3020008