Functional forms in joint models for longitudinal and time-to-event data: A practical guide with application and interpretation
تقدم هذه الورقة دليلاً عملياً يوضح أن اختيار الشكل الوظيفي في النماذج المشتركة للبيانات الطولية والبيانات المتعلقة بالوقت حتى وقوع الحدث هو قرار نمذجة حاسم يشكل بشكل أساسي تفسير الارتباطات بين المؤشرات الحيوية والمخاطر، كما هو موضح من خلال تحليل مقارن لمختلف هياكل الارتباط المطبقة على بيانات تجربة الورم الأرومي الدبقي.
في عالم الطب الحديث، غالبًا ما يتتبع الأطباء المرضى على مدار سنوات، مراقبين كيف تتغير أرقام محددة في دمهم أو أنسجتهم يومًا بعد يوم. هذه الأرقام، المعروفة باسم المؤشرات الحيوية، تروي قصة حول كيفية سلوك المرض. وفي الوقت نفسه، يحاول الباحثون فهم متى قد يواجه المريض حدثًا صحيًا خطيرًا، مثل الانتكاسة أو الوفاة. لفترة طويلة، عامل العلماء هاتين القصتين بشكل منفصل: قصة عن الأرقام المتغيرة، وأخرى عن توقيت الحدث. ولكن في الواقع، هاتان القصتان متشابكتان بعمق. فخطر تعرض المريض لنتيجة سيئة لا يتعلق فقط بمكان وصول أرقامه الآن، بل بكيفية وصوله إلى هناك. هل ارتفعت الأرقام ببطء على مدار أشهر، أم أنها قفزت فجأة في الأسبوع الماضي؟ هل تذبذبت صعودًا وهبوطًا بجنون، أم ظلت مستقرة؟ إن فهم الشكل الدقيق لهذه الرحلة أمر بالغ الأهمية للتنبؤ بالمستقبل، ومع ذلك، كان تحديد كيفية ترجمة تاريخ طبي معقد إلى تنبؤ بالمخاطر لغزًا صعبًا.
لقد شرع فريق من الباحثين في حل هذا اللغز من خلال إنشاء دليل عملي لكيفية الربط بين هاتين القصتين. ركزوا على نوع محدد من الأدوات الإحصائية يسمى "النموذج المشترك" (joint model)، وهو مصمم للنظر في القياسات المتكررة وأوقات البقاء معًا. لم يكن جوهر عملهم يتعلق بابتكار رياضيات جديدة، بل بتوضيح كيفية اختيار "العدسة" المناسبة لرؤية البيانات. وجادلوا بأن الطريقة التي يقرر بها الباحث ربط تاريخ المؤشر الحيوي بخطر المريض هي خيار علمي حاسم، وليست مجرد تفصيل تقني. فإذا تم اختيار العدسة الخاطئة، فإن الاستنتاجات الناتجة حول ما يدفع المرض قد تكون مضللة. ولإثبات ذلك، طبقوا دليلهم على بيانات حقيقية من تجربة سريرية كبرى شملت مرضى يعانون من نوع من سرطان الدماغ يسمى "ورم الأرومة الدبقية" (glioblastoma). وقد استخدموا تعداد خلايا الدم البيضاء، وهو مقي de قياسي لوظيفة الجهاز المناعي، ليظهروا كيف يمكن لطرق مختلفة للنظر في نفس البيانات أن تؤدي إلى فهم مختلف، وأحيانًا متضارب، لمخاطر المريض.
حدد الباحثون عدة طرق متميزة لوصف تاريخ المؤشر الحيوي للمريض، حيث تلتقط كل منها جانبًا مختلفًا من رحلته الطبية. النهج الأكثر شيوعًا ينظر ببساطة إلى المستوى الحالي للمؤشر الحيوي في أي لحظة معينة؛ فهو يتساءل: "ما هو مستوى تعداد خلايا الدم البيضاء الآن؟" ويفترض أن هذا الرقم الواحد هو ما يحدد الخطر. ومع ذلك، أظهر الفريق أن هذه الرؤية تغفل الصورة الأكبر. فقد يمكن لمريضين أن يكون لهما نفس تعداد خلايا الدم البيضاء بالضبط في يوم معين، ولكن أحدهما كان في حالة ارتفاع مستمر لأسابيع بينما كان الآخر في حالة انخفاض. إذا كانت سرعة هذا التغيير مهمة، فإن رؤية المستوى الحالي البسيطة ستفشل في التمييز بينهما.
ولمعالجة ذلك، استكشف الباحثون عدسات أكثر ديناميكية. أحد النهج ينظر إلى سرعة التغيير، أو "السرعة المتجهة" (velocity)، متسائلًا عما إذا كان المؤشر الحيوي يتصاعد أو يهبط في هذه اللحظة تحديدًا. والآخر ينظر إلى "التسارع" (acceleration)، أو مدى سرعة تغير تلك السرعة، مما يلتقط ما إذا كان الارتفاع يتسارع أم يتباطأ. كما فحصوا النهج التراكمية، التي تنظر إلى إجمالي التعرض الذي تعرض له المريض للمؤشر الحيوي بمرور الوقت. هذا يشبه حساب المساحة الكلية تحت المنحنى، والتي تمثل عبء المرض على مدار أشهر أو سنوات، بدلاً من مجرد لقطة خاطفة. واختبروا أيضًا طرقًا تركز على التغيرات الأخيرة خلال نافذة زمنية محددة، مثل الثلاثين يومًا الماضية، لمعرفة ما إذا كان التدهور الأخير أكثر خطورة من الاتجاهات طويلة الأمد. وأخيرًا، نظروا في "التباين" (variability)، حيث يقيسون مدى تقلب أرقام المريض حول متوسطه. فقد يكون المريض الذي تتأرجح أرقامه بجنون معرضًا لخطر أكبر من مريض لديه نفس المتوسط ولكن بخط مستقر، لأن عدم الاستقرار هذا قد يشير إلى جسد يكافح لتنظيم نفسه.
عندما طبق الفريق هذه العدسات المختلفة على بيانات تجربة سرطان الدماغ، سلطت النتائج الضوء على أهمية الطريقة المختارة، رغم أن الأدلة الإحصائية كانت غير مؤكدة في كثير من الأحيان. فعندما نظروا فقط إلى تعداد خلايا الدم البيضاء الحالي، أشارت العلاقة بالبقاء على قيد الحياة إلى ارتباط إيجابي محتمل، ولكن مع وجود شك كبير حيث تضمنت فترة الثقة الصفر. وبالمثل، عندما نظروا إلى التعرض التراكمي — أي العبء الإجمالي لخلايا الدم البيضاء بمرور الوقت — كان الارتباط ذا دلالة إحصائية، ومع ذلك تطلب حجم التأثير تفسيرًا دقيقًا. وفي حالات أخرى، قدمت النظر إلى معدل التغيير أو التباين تقديرات عددية مختلفة، لكن فترات الثقة الواسعة لهذه المقاييس أشارت إلى عدم يقين كبير وأدلة ضعيفة على وجود ارتباط واضح في هذه المجموعة من البيانات تحديدًا. ووجدت الدراسة أن هذه الأشكال الوظيفية المختلفة ليست قابلة للتبادل؛ فهي تقيس حقائق بيولوجية مختلفة. فالرقم الكبير في نموذج واحد لا يعني الشيء نفسه في نموذج آخر، لأنها تقيس أشياء مختلفة، مثل المستوى الحالي مقابل التاريخ التراكمي.
أكد الباحثون أنه لا توجد طريقة واحدة "مثلى" لربط نقاط البيانات هذه. فالخيار الصحيح يعتمد تمامًا على السؤال البيولوجي المطروح. فإذا كان المرض مدفوعًا بطفرة مفاجئة، فقد يكون النظر في المستوى الحالي هو الأفضل. وإذا كان مدفوعًا بسنوات من التعرض، فإن المقياس التراكمي يصبح ضروريًا. وإذا كان عدم استقرار النظام هو مصدر الخطر، فإن قياس التباين سيكون هو المفتاح. وتخلص الورقة البحثية إلى أن على الباحثين مواءمة خياراتهم الإحصائية بعناية مع فهمهم العلمي للمرض. ومن خلال معاملة اختيار كيفية ربط البيانات كجزء أساسي من الاستقصاء العلمي، وليس مجرد خطوة تقنية خفية، يمكن للأطباء والعلماء تجنب استخلاص استنتاجات مضللة. تضمن هذه الوضوح أنه عندما يتنبأ نموذج ما بمخاطر المريض، فإن ذلك يستند إلى الميزة المحددة لتاريخه الصحي التي تهم بالفعل، مما يؤدي إلى تنبؤات أكثر دقة ورعاية أفضل.
ملخص تقني: الأشكال الوظيفية في النماذج المشتركة للبيانات الطولية والبيانات المتعلقة بالوقت حتى وقوع الحدث
بيان المشكلة
تُعد النماذج المشتركة للبيانات الطولية والبيانات المتعلقة بالوقت حتى وقوع الحدث أدوات راسخة في الأبحاث السريرية لتحليل العلاقة بين قياسات المؤشرات الحيوية المتكررة ومخاطر وقوع الأحداث (مثل الوفاة أو تطور المرض). ومع ذلك، فإن أحد الجوانب الحرجة التي غالبًا ما يتم إغفالها في هذه النماذج هو اختيار الشكل الوظيفي (أو بنية الارتباط). يحدد هذا الاختيار كيفية ارتباط مسار المؤشر الحيوي بمخاطر الحدث.
في الوقت الحالي، يعتمد الباحثون بشكل متكرر على مواصفات افتراضية (عادةً ما تكون "القيمة الحالية" للمؤشر الحيوي) دون النظر فيما إذا كان هذا الشكل يلتقط الخصائص ذات الصلة سريريًا لتطور المؤشر الحيوي. فالخصائص المختلفة — مثل معدل التغيير، أو التعرض التراكمي، أو التباين — قد تؤثر على التنبؤ بالنتائج بطرق متميزة. إن تجاهل هذه الجوانب يمكن أن يؤدي إلى استنتاجات علمية غير كاملة أو مضللة. علاوة على ذلك، فإن تفسير معاملات الارتباط يعتمد بشكل كبير على الشكل الوظيفي المختار، مما يجعل أحجام التأثير غير قابلة للمقارنة عبر المواصفات المختلفة.
المنهجية
تقدم الورقة نظرة عامة مهيكلة ودليلًا عمليًا لمختلف الأشكال الوظيفية المستخدمة لربط العمليات الطولية وعمليات البقاء. وتتضمن المنهجية ما يلي:
التصنيف النظري: يصنف المؤلفون هياكل الارتباط إلى:
الأشكال اللحظية: ربط المخاطر بميزات يتم تقييمها عند الزمن t (القيمة الحالية، الميل/السرعة المعتمد على الزمن، التسارع/الانحناء).
الأشكال القائمة على التغيير: قياس الانحرافات عن خط الأساس أو النوافذ الزمنية الأخيرة (دلتا، دلتا المقيدة).
التأثيرات العشوائية المشتركة: ربط المخاطر مباشرة بالانحرافات الكامنة الخاصة بكل حالة.
الأشكال القائمة على التباين: ربط المخاطر بعدم الاستقرار أو التباين في مسار المؤشر الحيوي.
التنفيذ والبرمجيات: يركز الدليل على تنفيذ هذه الأشكال باستخدام حزمة JMbayes2 (الإصدار 0.6-0) في لغة R. تفصل الورقة الصيغة (syntax) المطلوبة لتحديد هذه الهياكل، مع ملاحظة أن بعض الأشكال (مثل الميل والتسارع) تتطلب مشتقات عددية (الفروق المحدودة) أو مصفوفات تصميم محددة، بينما تقوم JMbayes2 بمعالجة هذه الحسابات داخليًا.
التوضيح التجريبي: يطبق المؤلفون هذه الأشكال الوظيفية على بيانات من تجربة MIRAGE لسرطان الأرومة الدبقية (EORTC 1709 / Canadian Cancer Trials Group CE.8).
البيانات: 680 مريضًا مصابًا بمرض الأرومة الدبقية حديث التشخيص.
النتيجة الطولية: عدد خلايا الدم البيضاء (WBC)، الذي تم قياسه أسبوعيًا أثناء العلاج الإشعاعي.
النتيجة المتعلقة بالوقت حتى وقوع الحدث: البقاء الإجمالي (OS)، مع التعديل بناءً على المتغيرات الأساسية (العمر، العلاج، مدى استئصال الجراحة، مقياس كارنوفسكي الوظيفي KPS، وحالة MGMT).
النمذجة: تم استخدام نموذج التأثيرات المختلطة الخطية مع استخدام "السبلاينز" التكعيبية الطبيعية (natural cubic splines) للنموذج الفرعي الطولي لالتقاط الاتجاهات غير الخطية لخلايا الدم البيضاء. كما تم استخدام نموذج مخاطر "كوكس" النسبية لنموذج البقاء.
نهج التباين: بالنسبة للشكل القائم على التباين، تم استخدام نهج من خطوتين: أولاً، ملاءمة نموذج مختلط لاستخراج البواقي المطلقة (absolute residuals)؛ ثانيًا، نمذجة هذه البواقي كمخرج طولي جديد مرتبط بالمخاطر.
المساهمات الرئيسية
تصنيف شامل: تصنف الورقة بشكل منهجي وتقارن بين مجموعة واسعة من هياكل الارتباط، بما في ذلك الأشكال الأقل نقاشًا مثل التسارع (الانحناء) والأشكال القائمة على التباين.
دليل تنفيذ عملي: توفر تعليمات صريحة حول صيغة R وسير العمل لـ JMbayes2، مما يوضح كيفية تجاوز الإعدادات الافتراضية لاختبار فرضيات بيولوجية محددة.
إطار التفسير: يوضح المؤلفون أن مقياس والمعنى السريري لمعامل الارتباط (α) جوهريان للشكل الوظيفي المختار. ويؤكدون أن المعاملات من أشكال مختلفة (مثل القيمة الحالية مقابل الميل) ليست قابلة للمقارنة مباشرة من حيث المقدار.
نمذجة التباين: تفصل الورقة نهجًا عمليًا من خطوتين لدمج تباين المؤشر الحيوي (عدم الاستقرار) في برامج النمذجة المشتركة القياسية، مما يعالج نقص الدعم الأصلي للنماذج المشتركة غير متجانسة التباين في العديد من الحزم.
التطبيق التوضيحي: من خلال تطبيق جميع الأشكال على مجموعة بيانات MIRAGE نفسها، توضح الورقة كيف يمكن للأشكال الوظيفية المختلفة أن تؤدي إلى تفسيرات مختلفة للآثار، وفي بعض الحالات، استنتاجات موضوعية مختلفة فيما يتعلق بالعلاقة بين المؤشر الحيوي والمخاطر.
النتائج
أدى تطبيق هذه الأشكال على بيانات تجربة MIRAGE إلى الملاحظات التالية:
القيمة الحالية: ارتبطت زيادة قدرها وحدة واحدة في عدد خلايا الدم البيضاء بزيادة قدرها 1.1% في اللوغاريتم اللحظي لمخاطر الوفاة (α=0.011)، رغم أن فاصل الثقة الائتماني بنسبة 95% شمل الصفر، مما يشير إلى عدم اليقين.
الميل المعتمد على الزمن: أظهر الأشخاص الذين لديهم زيادات سريعة في خلايا الدم البيضاء خطرًا أعلى للوفاة، لكن فاصل الثقة كان واسعًا ([−1.821,4.751])، مما يشير إلى ضعف الأدلة على وجود ارتباط واضح في هذه المجموعة من البيانات.
التسارع: وُجد ارتباط إيجابي بين انحناء مسار خلايا الدم البيضاء ومخاطر البقاء (α=1464.33). ويعكس هذا المقدار الكبير مقياس المشتق الثاني وليس تأثيرًا بيولوجيًا مفرطًا.
التعرض التراكمي: أظهر متوسط التعرض التاريخي لخلايا الدم البيضاء ارتباطًا إيجابيًا ذا دلالة إحصائية مع خطر الوفاة (α=0.051، فاصل ثقة 95% [0.015,0.090]). أما التأثيرات التراكمية المقيدة على مدار 30 يومًا فقد أظهرت ارتباطًا أضعف وغير معتد به إحصائيًا.
التغيير (دلتا): أظهر التغيير عن خط الأساس تقديرًا نقطيًا سالبًا ولكن مع عدم يقين واسع. أما التغيير خلال نافذة زمنية مدتها 30 يومًا فقد أشار إلى ارتباط إيجابي محتمل مع المخاطر، وإن كان فاصل الثقة الائتماني يتضمن الصفر مرة أخرى.
التباين: ارتبط التباين الأكبر داخل الشخص الواحد في خلايا الدم البيضاء إيجابيًا بخطر الوفاة (α=0.236)، رغم أن الفاصل [−0.04 to 0.549] أشار إلى عدم يقين كبير.
تسلط النتـائج الضوء على أن الأشكال الوظيفية المختلفة تستهدف سمات متميزة لمسار المؤشر الحيوي. فعلى سبيل المثال، التقطت الأشكال التراكمية آثار التعرض طويلة المدى التي أغفلتها الأشكال اللحظية، بينما قدمت أشكال التباين إشارة تنبؤية بديلة قائمة على عدم الاستقرار.
الأهمية والادعاءات
تؤكد الورقة أن اختيار الشكل الوظيفي هو قرار نمذجة رئيسي يحدد مباشرة تفسير الارتباط بين المؤشر الحيوي والمخاطر. إنه ليس مجرد تفصيل تقني بل هو قرار علمي يشكل الافتراضات الأساسية حول كيفية تأثير العملية الطولية على المخاطر.
يدعي المؤلفون ما يلي:
المواءمة ضرورية: يجب أن تكون هناك مواءمة دقيقة بين السؤال العلمي (مثل "هل المستوى الحالي مهم؟" مقابل "هل العبء التراكمي مهم؟") والشكل الوظيفي المختار لضمان التفسير الصحيح والتقرير الشفاف.
لا يوجد شكل عالمي: لا يوجد هيكل ارتباط واحد مناسب عالميًا. يجب أن يسترشد الاختيار بالسياق السريري، والمعقولية البيولوجية، والجانب المحدد من مسار المؤشر الحيوي الذي يُعتقد أنه ذو صلة.
القابلية للتفسير: بينما تعد النماذج المرنة (مثل تلك التي تستخدم السبلاينز أو التأثيرات العشوائية) قوية، إلا أن الصيغ القائمة على المسار (التي تربط المخاطر بالقيمة الحالية، والميل، إلخ) غالبًا ما تكون مفضلة لفهم الآليات البيولوجية مقارنة بالتأثيرات العشوائية المشتركة، التي يمكن أن تعيق التفسير.
المنفعة العملية: تسمح مجموعة الأشكال الوظيفية المتاحة للباحثين بالتقاط العلاقات المعقدة الموجودة في البيانات السريرية التي قد تغفلها المواصفات الافتراضية، مما قد يتجنب الاستنتاجات المضللة.
تختتم الورقة بالإقرار بالقيود، مثل الاعتماد على مجموعة بيانات واحدة للتوضيح والطبيعة التقريبية لنهج التباين المكون من خطوتين، لكنها تؤكد أن تقديم دليل عملي لهذه الهياكل المتنوعة يعد خطوة ضرورية لتحسين دقة النمذجة المشتركة في الأبحاث التطبيقية.