تتضمن هذه الدراسة سجلات طبية سابقة للأشخاص البشر وبيانات سجلات التخدير والتمريض. أجريت الدراسة وفقا للإرشادات المؤسسية ذات الصلة، وتمت مراجعتها واعتمادها من قبل لجنة الأخلاقيات المؤسسية (مراجعة أخلاقيات مستشفى رن رن شو 2025 رقم 2113)، مع تنازل عن الموافقة المكتوبة المستنيرة؛ تم إلغاء تحديد جميع البيانات قبل التصدير واستخدمت فقط لأغراض هذه الدراسة.
تم الإبلاغ عن هذه الدراسة وفقا لبيان TRIPOD لتطوير نموذج التنبؤ متعدد المتغيرات والتحقق الداخلي؛ قائمة التحقق المكتملة ل TRIPOD متوفرة في الجدول التكميلي S1.
1. تصميم البحث وبيئة الدراسة
- عرف تصميم الدراسة كدراسة رصدية بأثر رجعي بمركز واحد باستخدام تصميم مجموعة تاريخية. حدد الإطار الزمني للدراسة ليصبح من 1 يناير 2021 حتى 31 ديسمبر 2024.
- حدد هدف الدراسة كتطوير أداة تقييم تنبؤية للتمريض للمضاعفات المتعلقة بالوضع في غرفة العمليات الهجينة DSA وإجراء التحقق الداخلي وتقييم الأداء. عرف "الجراحة" كوحدة التحليل وخصص لكل عملية معرفا فريدا (مثل رقم المرضى الداخليين + تاريخ العملية + وقت بدء التخدير).
- عندما يخضع نفس المريض لعدة عمليات جراحية مؤهلة، احتفظ فقط بسجل الجراحة الأول المؤهل. عرف "الأول" بأنها الجراحة المؤهلة التي تبدأ فيها التخدير أو التخدير في وقت مبكر ضمن فترة الدراسة.
- حدد نقطة بدء استخراج البيانات بأنها الوقت الذي يدخل فيه المريض غرفة العمليات ويبدأ التخدير أو التخدير. تثبيت نقطة المتابعة عند 72 ساعة بعد العملية، وبالتالي تغطية نافذة الحدوث الأولى للمضاعفات المتعلقة بالوضعية.
- اقتصر إعداد الدراسة على غرفة العمليات الهجينة DSA ووثق أنها مجهزة بنظام تصوير الأوعية الدموية الرقمية القابلة للطرح على الشاشة المسطحة وإمكانية الجراحة المفتوحة. سجل أن الوضعيات الشائعة هي الاستلقاء، الاستلقاء، والجانبية، وأن تعديل الوضعية أثناء العملية محدود بمعدات التصوير، والمجال المعقم، وتخطيط طرق الوصول، مما يؤدي إلى فرص أقل للدوران ومدة العملية أطول نسبيا.
2. مجموعة الدراسة، تدفق فحص الحالات، واستراتيجية حجم العينة
- تصدير جميع إدخالات جراحة "غرفة العمليات الهجينة DSA" ضمن نافذة الدراسة من سجل غرفة العمليات المسجل على التوالي. احتفظ بحقول مطابقة مثل رمز الإجراء، معرف غرفة العمليات، تاريخ العملية، رقم المرضى الداخليين، ووقت بدء التخدير.
- استخدم رقم المريض الداخلي وتاريخ الجراحة للمطابقة الأولى، وقم بإجراء فحص ثانوي باستخدام وقت بدء التخدير ومعرف غرفة العمليات. إذا كانت هناك عدة سجلات تخدير في نفس اليوم لنفس رقم المرضى الداخليين، يرجى تفضيل مطابقة الرقم الذي يحتوي على معرف ثابت لغرفة العمليات ووقت بدء التخدير الأقرب إلى وقت تسجيل الجراحة.
- تشمل المرضى البالغين الداخليين الذين تبلغ أعمارهم ≥18 عاما وتم إجراء جراحاتهم في غرفة العمليات الهجينة DSA وكانت عمليات تدخلية أو أوعية عصبية أو قلبية أو أوعية دمائية طرفية أو هجينة. تأكد من أن التخدير العام أو التخدير قد تم أثناء الجراحة، وأن الوثائق الطبية والتخديرية والتمريضية مكتملة، وأن المتابعة حتى 72 ساعة بعد العملية متاحة.
- استبعد الحالات التي كان فيها حدث النتيجة الأولية موجودا بالفعل قبل العملية لتجنب خط أساس غير صفري للنتيجة. استبعد الحالات التي تم فيها مقاطعة الإجراء أثناء العملية لأسباب غير متعلقة بوضع الوضعية، بحيث لم يكن من الممكن إجراء تقييم إضافي، وكذلك الحالات التي تم نقلها من المستشفى خلال 72 ساعة بعد العملية، مما أدى إلى نتيجة غير محددة.
- تحديد عتبة اختفاء المتغيرات الرئيسية للتنبؤ من خلال حساب نسبة الغياب لكل حالة عبر المتغيرات الرئيسية المحددة مسبقا وتسجيل أسباب الغياب. إذا تجاوزت نسبة الغياب لحالة عبر هذه المتغيرات الرئيسية 30٪، استبعد الحالة وفقا للقواعد المحددة مسبقا ووثق سبب الاستبعاد في سجل الفحص.
- اطلب من باحثين اثنين إكمال الفرز بشكل مستقل وإنشاء سجلات الفحص (بما في ذلك الإدراج والاستبعاد والأسباب). في الحالات المتعارضة، يجب أن يتتبع باحث ثالث إلى السجلات الأصلية للفصل فيها، ويحتفظ بسجل الفصل لإمكانية التتبع.
- لا تحدد مسبقا حدا أعلى لحجم العينة، بل تعتمد استراتيجية شاملة متتالية لكل الحالات. تحكم في تعقيد النموذج بعدد الأحداث، وتحديد عدد المعلمات التي تدخل في النموذج في النهاية وفقا لمبدأ 15 حدثا على الأقل لكل معلمة.
3. مقاييس الدراسة وتعريفات المتغيرات
- أنشئ قاموس للمتغيرات وحدد لكل متغير نظام المصدر، اسم الحقل، قواعد القيم، والتسلسل الزمني. للمتغيرات المستمرة، استخدم أحدث قيمة مختبر أو تقييم قبل الجراحة؛ بالنسبة لمتغيرات الجراحة والتخدير، استخدم أسوأ قيمة مسجلة خلال كامل العملية أو التعرض التراكمي؛ بالنسبة لمتغيرات التمريض والأجهزة، استخدم التطبيق الفعلي داخل العملية. تحديد العتبات المحددة مسبقا وحدود الفئات مسبقا وفقا للمعايير السريرية الشائعة الاستخدام، والتمريض والتخدير حول العمليات، والأدبيات المنشورة. عندما لا يتوفر حد معين مقبول على نطاق واسع، استخدم الفئات العملية لتسهيل التسجيل والتفسير والتطبيق السريري.
- النتيجة الأساسية: النتيجة المركبة لمضاعفات متعلقة بالوضعية
- حدد النتيجة الأساسية كنتيجة مركبة لمضاعفات متعلقة بالوضعية تحدث من أثناء العملية وحتى 72 ساعة بعد العملية. إذا حدثت عدة أحداث مكونة في نفس الحالة، قم بعد النتيجة المركبة مرة واحدة فقط، وسجل وقت الظهور الأول وجميع أنواع المكونات للتحليلات الثانوية.
- عرف إصابة المرحلة الثانية الجديدة أو الإصابة ذات الضغط العالي أو إصابة الأنسجة العميقة خلال 72 ساعة بعد العملية كواحدة من الأحداث المكونة. استخدم نموذج تقييم جلد الرضاعة، والمرحلة، وسجلات التمريض كدليل، واطلب سجلات يمكن تحديدها في مواقع تشريحية محددة.
- عرف إصابة الأعصاب الطرفية بأنها حدوث فقدان حسي، أو تنميل، أو انخفاض في قوة العضلات يستمر لمدة ≥24 ساعة. اطلب أن يكون التوزيع التشريحي متوافقا مع منطقة التلامس/الجذب في الوظيفة في ذلك الوقت، واستخدام ملاحظات التقدم أو سجلات التقييم المتخصصة كدليل رئيسي.
- عرف إصابة المفاصل أو العضلات الهيكلية بأنها ألم موضعي مصحوبا بنتائج تصوير إيجابية أو نتائج فحص بدني متخصص تشير إلى إصابة شد/ضغط. سجل طريقة التصوير أو نتيجة الفحص، وحدد نقاط تطبيق القوة المتعلقة بالموقع.
- عرف الأحداث المتعلقة بالموضع التنفسي أو عدم الاستقرار الديناميكي الدموي بأنها ضغط مجرى الهواء، أو تقييد التنفس، أو تقييد العودة الوريدية الناتجة عن التموضع الذي يتطلب إعادة تموضع، أو ضغط، أو تصعيد ضغط الأوعية، أو مقاطعة الإجراء. استخدم سجلات أحداث نظام التخدير، وسجلات التخدير، وسجلات التمريض/الطبيب التي تحتوي على أوصاف صريحة "متعلقة بالوضع" ومقاييس الإدارة المقابلة كدليل.
- عرف ورم دموي بقطر ≥5 سم على جانب الثقب أو التسرب الذي يتطلب تثبيت ضغط إضافي، مع السجلات التي تشير إلى تفاقم مرتبط بقيود الموضع، كأحد الأحداث المكونة. استخدم سجلات التمريض وسجلات الأطباء في الجناح خلال 72 ساعة بعد العملية كدليل، وسجل الإجراءات الإدارية الإضافية.
- تحديد الصلة الموضعية وفقا لثلاثة معايير: الذكر الصريح للأهمية الموضعية في السجلات، الموقع التشريحي المتوافق مع نقاط التعرض الموضعية، وعدم وجود سبب بديل أكثر إقناعا. إذا كان هناك تعارض في النسبة بين سجلات التمريض والأطباء، ضع الأولوية للسجل الأصلي أقرب إلى وقت الحدث، واجعل اثنين من الحكمين يراجعهم ويتوصلوا إلى توافق.
- تعريفات المتغيرات المتنبئة
- حدد متغيرات مستوى المريض كالعمر، الجنس، مؤشر كتلة الجسم، السكري، الاعتلال العصبي الطرفي، مرض الشرايين الطرفية، مرض الكلى المزمن في المرحلة ≥ 3، فشل القلب الاحتقاني، حالة التدخين، الألبومين قبل العملية، الهيموغلوبين قبل العملية، مؤشر تشارلسون المصاحب، ودرجة برادن قبل العملية. عرف "ما قبل العملية" بأنه أحدث سجل صالح قبل دخول غرفة العمليات، ضمن نافذة التقييم الروتينية للمستشفى حول العملية.
- حدد متغيرات الجراحة والتخدير كفئة الإجراء، نوع التخدير، نوع الوضعية، فئات زاوية دوران الرأس والرقبة، فئات زاوية استخلاء الطرف العلوي، انتزاع أو ثني الطرف السفلي > 45°، مدة العملية، مدة التنظير الفلوري، حالة غير قابلة للدوران أثناء العملية، تقدير فقدان الدم، أقل MAP، الوقت التراكمي مع MAP< 65 مم زئبق، أقل درجة حرارة لنو، استخدام العوامل الفعالة للأوعية، وأقصى ما يعادل النورإبينفرين. صنف الوضعية كوضعية على الأقل/الاستلقاء/الجانبية؛ يصنف دوران الرأس والرقبة إلى ≤30°، 31–60°، و>60°؛ صنف اختطاف الطرف العلوي ب ≤90° و>90°).
- تصدير سلسلة زمنية MAP مع الطوابع الزمنية من نظام إدارة معلومات التخدير. استخدم فرق الزمن بين السجلات المتجاورة كطول الفترة، وجمع مدة كل الفترات مع MAP < 65 مم زئبق للحصول على "MAP < 65 زمن تراكمي"؛ إذا كانت فترة واحدة مفقودة أو كانت الفجوة غير طبيعية (مثلا >5 دقائق بدون تفسير معقول)، ضع علامة على تلك الفترة كمفقودة ولا تدرجها في الإجمالي التراكمي.
- تعريف متغيرات الرضاعة والجهاز كنوع من حشوة الطاولة، وما إذا تم اكتمال حماية نقاط الضغط الرئيسية، وما إذا كانت الأطراف العلوية مثبتة على ألواح ذراع مبطنة، وما إذا كان الرأس مدعوما في وضع محايد، وفئة موقع الثقب، وحجم الغمد، ومدة بقاء الغمد، وطريقة تأمين/إغلاق موقع الطعن، واكتمال توثيق فحص الوضعية أثناء العملية. عرف "اكتمال حماية نقطة الضغط الرئيسية" بأنها مكتملة عندما توجد سجلات حماية في جميع المواقع الأربعة—القذالي، لوح الكتف، منطقة العصب المغطائي، والكعب؛ وإلا، يعتبر غير مكتمل بالكامل.
- سجل ما إذا تم استخدام أي عوامل فعالة للأوعية، وترميز غير المستخدمين ك "غير مستخدم". بالنسبة للمستخدمين، يسجل الحد الأقصى المكافئ للنورإبينفرين (ميكروغرام/كجم/دقيقة); عند كل نقطة زمنية قياسية للتخدير، احسب NE_eq_t = norepinephrine_t + epinephrine_t + 0.01 * dopamine_t + 0.06 * phenylephrine_t + 2.5 * vasopressin_t، حيث يعبر عن النورإبينفرين، الإبينفرين، الدوبامين، والفينيل إفرين بميكروغرام/كجم/دقيقة والفازوبريسين بوحدة U/min; تم ترميز العوامل غير المستخدمة في نقطة زمنية معينة ك 0، و maximum_norepinephrine_equivalent = max(NE_eq_t, na.rm = TRUE) خلال فترة العملية.
ملاحظة: تم تثبيت قواعد التحويل والتجميع مسبقا وتطبيقها في سكريبت التحليل.
4. جمع البيانات ومراقبة الجودة
- توفير تدريب موحد لباحثين مسؤولين عن استخراج البيانات وتوزيع إجراءات التشغيل القياسية (SOP) وقاموس المتغيرات. بعد التدريب، أكمل استخراجا تجريبيا بعينة صغيرة وراجع تعريفات ومسارات استخراج الحقول الغامضة.
- استخلاص التركيبة السكانية والأمراض المصاحبة من السجل الطبي الإلكتروني، ومؤشرات المختبر قبل العملية من النظام المختبري، والعلامات الحيوية والأدوية من نظام التخدير، وتحديد المواقع وإجراءات التمريض من نظام التمريض، ومعلومات مدة التنظير الفلوري والمعلومات المتعلقة بالتصوير من PACS. قم بمطابقة الأنظمة المشتركة باستخدام رقم المرضى وتاريخ الجراحة، واحتفظ بمفاتيح المطابقة وسجل المطابقة.
- عندما تحدث تعارضات ميدانية، أعد إلى سجلات الورق الأصلية واعتبر السجلات الورقية كمعيار. سجل نوع الصراع، والقيمة الأصلية، والقيمة المصححة، ومصدر الأدلة في سجل معالجة النزاعات.
- قم بإدخال بيانات مزدوج من قبل شخصين لجميع الحقول وإنشاء قائمة تناقضات. اطلب من باحث ثالث أن يتحقق من كل بند مع السجلات الأصلية للفصل ويقفل النسخة النهائية من قاعدة البيانات للتحليل.
- إنشاء جدول للنبات السريرية المعقولة للمتغيرات المستمرة وإجراء فحوصات المدى التلقائية بعد الاستيراد. تتبع القيم خارج النطاق للتحقق؛ إذا لم يكن بالإمكان التحقق منها، قم بضبطها على المفقودة واحتفظ بالقيمة الأصلية للرجوع إليها.
- قم باستخراج تكراري لأول 50 حالة وحساب الاتساق، باستخدام كوهين κ للمتغيرات الفئوية وICC (تأثيرات عشوائية ثنائية الاتجاه، توافق مطلق) للمتغيرات المستمرة. إذا كان κ < 0.80 أو ICC < 0.85، يجري إعادة تدريب للمتغيرات ذات الصلة ويعيد التدقيق حتى يتم تحقيق المعيار11.
- احسب عدد ونسبة الحالات المفقودة لكل متغير متنبئ مرشح وأنشئ تقرير عن الاختفاء. بالنسبة للمتغيرات الرئيسية المتنبئة التي لا تتجاوز نسبة الاختفاء 30٪، ابدأ الحساب المتعدد؛ لن يتم احتساب متغير النتيجة الأساسي.
- تفاصيل تنفيذ الحساب المتعدد بواسطة المعادلات المتسلسلة (MICE)
- حدد عدد مجموعات البيانات المحسوبة عند m = 10 واضبط عدد التكرارات على الأقل 20 لتعزيز التقارب. قم بتعيين بذرة عشوائية لضمان تكرار الحساب وحفظ كائن الفأرة لقابلية التدقيق.
- استخدم مطابقة المتوسط التنبؤي (PMM) للمتغيرات المستمرة للحفاظ على الشكل التوزيعي. استخدم الحساب بالانحدار اللوجستي للمتغيرات الثنائية، والحساب اللوجستي متعدد الحدود للمتغيرات متعددة الفئات غير المرتبة (multi categorys)، والحساب اللوجستي الترتيبي للمتغيرات الفئوية المرتبة (التصنيفات). في R، تم تعريف المتغيرات متعددة الفئات غير المرتبة ك nom_vars <- c("smoking_status","procedure_category","anesthesia_type",
"position_type"، "table_padding_type"، "puncture_site_
الفئة","puncture_site_closure_method"), وتم تعريف المتغيرات الفئوية المرتبة ك ord_vars <- c("head_neck_rotation_cat","upper_limb_abduction_cat"); dat[nom_vars] <- lapply(dat[nom_vars]، عامل); dat[ord_vars] <- lapply(dat[ord_vars], أمر بذلك); الميث < - الفئران::make.method(dat)؛ الميث[sapply(dat, is.numeric)] <- "PMM"; الميث[names(dat)[sapply(dat, function(x) is.factor(x) && nlevels(x) == 2)]] <- "logreg"; الميث[nom_vars] <- "بوليريج"; الميث[ord_vars] <- "بولر"؛ وميث[c("inpatient_number","surgery_id","outcome")] <- "".
- تضمين جميع المتغيرات المرشحة المتنبئة ومتغير النتيجة في مصفوفة التنبؤ بالإعدادات لتحسين جودة الحساب، ولكن إجراء الحساب فقط للمتغيرات المتنبئة التي تحتوي على غياب الحساب. اضبط حقول معرف لا يجب أن تعمل كمؤشرات (مثل عدد المرضى الداخليين، معرف الجراحة الفريد) بحيث لا تشارك في نموذج التحديد لتجنب إدخال معلومات غير سريرية. لاتباع هذا البروتوكول، نفذ مصفوفة التنبؤ ك pred <- mice::make.predictorMatrix(dat)؛ pred[,] <- 1; التشخيص (التقسيس) <-0؛ pred[, c("inpatient_number","surgery_id")] <-0; pred[c("inpatient_number","surgery_id","outcome"), ] <-0; وpred[setdiff(colnames(pred), c("inpatient_number","surgery_id","outcome")), "outcome"] <-1، بحيث تستخدم النتيجة فقط كمتنبئ مساعد وليس كشخص مستنسب للنتيجة.
- رسم مخططات تتبع التكرار وتحقق مما إذا كانت المتوسطات والتباين للمتغيرات المستمرة الرئيسية مستقرة لتقييم تقارب الحسابات. تم إجراء التعسيب ك imp < - الفئران::الفئران (dat, m = 10, maxit = 20, method = meth, PredictorMatrix = pred, seed = 20250101, printFlag = FALSE); تم تقييم التقارب باستخدام مخطط (IMP) ومن خلال التحقق من استقرار المتوسطات والتباينات للمتغيرات المستمرة الرئيسية عبر التكرارات.
5. تطوير أداة التقييم التنبؤي التمريضي
- المعالجة المبدئية وترميز المتغيرات المرشحة
- اضبط مستوى المرجع لفئة الإجراءات إلى "عصبي وعائي"، ومستوى المرجع للوضعية على "الاستلقاء". اضبط مستوى المرجع لدوران الرأس والرقبة على "≤30°"، ومستوى المرجع لانحراف الطرف العلوي إلى "≤90°". حدد مستوى المرجع لفئة موقع الثقب إلى "الشريان الكعبري"، ومستوى المرجع لطريقة إغلاق موقع الثقب إلى "جهاز إغلاق الأوعية الدموية".
- عمر المقياس لكل 10 سنوات، مدة العملية لكل 60 دقيقة، مدة التنظير الفلوري كل 10 دقائق، تقدير فقدان الدم لكل 100 مل، وقت الغلاف لكل 60 دقيقة، وحجم الغمد لكل 1 فر. بالنسبة للمتغيرات المفترضة ك "كلما كانت القيمة أقل، زادت المخاطر" (مثل أقل درجة حرارة في النواة، الألبومين، الهيموغلوبين، درجة برادن، أقل MAP)، تطبيق الترميز العكسي بحيث يمثل OR الانحدار "زيادة الخطر لكل وحدة انخفاض".
ملاحظة: يمكن تنفيذ الترميز العكسي بأخذ سالب المتغير الأصلي وتكبير الوحدة مثلا Temp_dec = -(درجة حرارت/0.5) أو Braden_dec = -(برادن/1)، للحفاظ على اتساق اتجاه تفسير OR.
- التحليل الأحادي المتغير والفحص للمتغيرات المرشحة
- ملاءمة نموذج الانحدار اللوجستي أحادي المتغير لكل متنبئ مرشح وتسجيل قيمة OR، و95٪ CI، وقيمة P. للمتغيرات التصنيفية، استخدم الترميز الوهمي مع تقديرات مستوى المرجع وتأثير الإخراج لكل مستوى بالنسبة للمرجع.
- استخدم المنحنيات المكعبة المقيدة لاختبار الخطية للمتغيرات المستمرة، باستخدام أربع عقد موضوعة عند كميات 5، 35، 65٪، و95٪. إذا أشار الاختبار إلى عدم خطية كبيرة، احتفظ بنفس شكل الخطوط في النموذج متعدد المتغيرات التالي.
- خذ اتحاد "المتغيرات المحددة سريريا مسبقا" و"المتغيرات ذات المتغير الأحادي P<0.20" كمجموعة مرشحة متعددة المتغيرات 12. قارن مجموعة المرشحين مع قيد عدد الأحداث، وعند الضرورة، احتفظ تفضيلا بالمتغيرات ذات الصلة السريرية للتحكم في الإفراط في التوافق.
- اختيار اللاسو وبناء النموذج الرئيسي
- تناسب انحدار لوجستي معاقب (LASSO، ألفا=1) على مجموعة المرشحين، مع العائلة = ثنائي الحدين. توحيد المتغيرات المستمرة (توحيد = TRUE) لضمان العقوبة العادلة عبر المتغيرات ذات المقاييس المختلفة.
- قم بإجراء التحقق المتقاطع بعشرة مرات متكررا 1,000 مرة، مع توليد تقسيمات طيات عشوائية في كل مرة وتسجيل انحراف التحقق المتقاطع المقابل ل lambda.min. عرف قوة العقوبة كوسيط lambda.min عبر 1000 تكرار (أو، بالمثل، λ المقابل لأقل متوسط انحراف)، واستخدم هذا λ لتناسب LASSO النهائي للحصول على مجموعة مستقرة من المتغيرات ذات معاملات غير صفرية.
- باستخدام المتغيرات ذات المعاملات غير الصفرية التي اختارها LASSO، يتناسب مع النموذج الرئيسي للانحدار اللوجستي متعدد المتغيرات غير المعاقب في مجموعة البيانات المكتملة. إذا تم الاحتفاظ بمصطلحات الخطوط، فيتم تركيبها في النموذج الرئيسي بنفس شكل الخطوط وتسجيل المعاملات أو مصطلحات الدوال المقابلة.
- ملاءمة النموذج الرئيسي بشكل منفصل في مجموعات البيانات العشر المحسوبة واستخراج معاملات الانحدار والأخطاء المعيارية. اجمع المعاملات والأخطاء القياسية وفترات الثقة باستخدام قواعد روبين، وقدم تقريرا عن قيمة OR وP المجمعة.
- تقييم المخاطر، الرسم النوموغرامي، وتصنيف المخاطر
- قم بمقياس خطي لمعاملات الانحدار الرئيسية للنموذج إلى نقاط صحيحة وفقا لطريقة سوليفان، مع تثبيت المقياس الأساسي وقواعد التقريب. حساب الدرجة الكلية وربطها بالاحتمال المتوقع عبر دالة لوجيت؛ يجب أن تكون صيغة التحويل متسقة مع تقاطع النموذج الرئيسي.
- استخدم حزمة rms لإنشاء صورة نوموغرامية بناء على النموذج الرئيسي، وتصديرها كصورة عالية الدقة أو PDF للنشر. تطبيق تسمية الملفات المعتمدة على مخرجات النوموجرام لتجنب الالتباس الناتج عن التعديلات اللاحقة.
- عرف الاحتمال المتوقع <5٪ كخطر منخفض، 5٪ إلى 14.9٪ كخطر متوسط، و≥15٪ كخطر عالي13. قم بتثبيت مخرجات الأداة إلى تنسيقين—نموذج ورقي من صفحة واحدة وآلة حاسبة إلكترونية—وحدد عملية التقييم إلى ثلاث خطوات لا تزيد عن ثلاث خطوات لتسهيل الاستخدام السريري.
6. التحقق من صحة النموذج وتقييم الأداء
- التحقق الداخلي وتصحيح التفاؤل
- استخدم نظام التمهيد للتحقق الداخلي، واضبط عدد إعادة العينات إلى 1,000 للحصول على تقديرات مستقرة. لكل إعادة عينة، أعد ضبط النموذج في عينة التمهيد، وحساب مقاييس الأداء في كل من عينة التمهيد والعينة الأصلية لتقدير التفاؤل.
- عرف التفاؤل بأنه الفرق المتوسط بين الأداء في عينة الbootstrap والأداء في العينة الأصلية. اطرح التفاؤل من الأداء الظاهر للحصول على الأداء المصحح، وأبلغ عن AUC المصحح بالتفاؤل، واعتراض المعايرة، وميل المعايرة، ودرجة بريير.
- التمييز
- ارسم منحنيات ROC واحسب AUC كمقياس للتمييز. تقدير AUC وCI بنسبة 95٪ باستخدام طريقة DeLong، وإجراء الحسابات على نفس مخرجات التنبؤ (الاحتمالات) لضمان المقارنة.
- عند العتبات المحددة مسبقا (5٪ و15٪)، احسب الحساسية، والنوعية، وPPV، وNPV. لكل عتبة، أخرج أيضا عداد مصفوفة الالتباس للسماح بالتحقق من الحساب.
- المعايرة والخطأ العام
- ملاءمة نموذج معايرة مع المتنبئ الخطي (لوجيت) كمتغير مستقل؛ يشير التقاطع إلى الانحياز المنهجي، والميل يشير إلى درجة الإفراط في التركيب. عرف المعايرة المثالية على أنها تقاطع = 0 وميل = 1، وأبلغ عن المقاييس قبل وبعد التصحيح.
- احسب درجة بريير كمقياس للخطأ الكلي، ورسم منحنيات المعايرة باستخدام تنعيم اللويس بجانب الخط المثالي 45°. تم حساب درجة بريير كمتوسط((نتيجة - predicted_probability)2)، وتم تطبيق تنعيم اللويس كإحصائيات::loess(outcome ~ predicted_probability، span = 0.75، درجة = 2، العائلة = "غاوسي"، السطح = "مباشر").
- منحنيات الفائدة السريرية الصافية واتخاذ القرار
- استخدم تحليل منحنى القرار لتقييم صافي الفائدة عبر احتمالات عتبة تتراوح بين 2٪ إلى 25٪14، ورسم استراتيجيات Treat-All وTreat-Nothing كمراجع. حدد حجم الخطوة العتبية عند 0.01 وحدد ذلك صراحة في السكريبت لضمان إمكانية التكرار.
- قم بتحليل حساسية الفائدة الصافية حول عتبات المخاطر المنخفضة والعالية، وتسجيل اتجاه تأثير تغيرات العتبة على الفائدة الصافية. قدم نتائج تحليل الحساسية جنبا إلى جنب مع العتبات المحددة مسبقا للتحقق من عقلانيتها.
7. التنفيذ الإحصائي وسير عمل السكريبتات القابلة لإعادة الإنتاج
- بيئة البرمجيات وتهيئة المشاريع
- قم بإجراء تحليلات في R 4.3.2 باستخدام حزم pROC وrms وglmnet والفئران وrmda.
- قم بتعيين seed عشوائي ثابت (20250101) في بداية السكريبت، وسجل معلومات الجلسة لتوثيق بيئة البرمجيات.
- استيراد البيانات، إعدادات النوع، وفحوصات الجودة
- استورد البيانات باستخدام read.csv("data/analysis.csv") أو readr::read_csv(). تحقق مما إذا كان عدد الصفوف يتطابق واحدا لواحدا مع معرف الجراحة الفريد، وتحقق من وجود معرفات مكررة.
- حول المتغيرات التصنيفية إلى عوامل واستخدم relevel() لتعيين مستويات مرجعية لتتناسب مع المخطط المسبق. حافظ على المتغيرات المستمرة رقمية، وطبق تحجيم الوحدات والترميز العكسي حسب الإعدادات المسبقة قبل النمذجة.
- بالنسبة للمتغيرات الفئاتية، احسب الترددات والنسب المئوية؛ بالنسبة للمتغيرات المستمرة، احسب المتوسط ± الانحراف المعياري أو الوسيط [IQR]. استخدم الطرق الرسومية و shapiro.test() لتقييم شكل التوزيع، وتحديد النهج الوصفي وفقا لذلك.
- معالجة البيانات المفقودة والحساب المتعدد (m = 10)
- احسب نسبة النقص لكل متغير وتصديرها كإخراج لكل missing_report.csv. تأكد من أن متغير النتيجة لا يحتوي على غياب، وقفل استراتيجية "فرض المتغيرات المتنبئة فقط، لا تحسب النتيجة."
- استخدم mice() مع m = 10 وmaxit ≥ 20 لتشغيل حساب المعادلات المتسلسلة، وحدد طرق الحساب حسب أنواع المتغيرات. استخدم complete(mids, action=k) لتصدير مجموعة البيانات المكتملة kth وحفظها كبيانات/imp_k.csv.
- رسم مخطط (الوسط) للتحقق من استقرار المسارات، ومقارنة توزيعات ما قبل وبعد النسبة للمتغيرات المستمرة الرئيسية. احفظ مخططات التشخيص في الإخراج أو imputation_diagnostics.pdf للتحقق.
- الانحدار الأحادي المتغير، اختبار الخطوط العصبية، وتوليد مجموعة المرشحين
- ملاءمة الانحدار اللوجستي أحادي المتغير في كل مجموعة بيانات محسوبة واستخدام pool() لدمج المعاملات والأخطاء المعيارية. تصدير قيم OR، وقيم CI بنسبة 95٪، وP كنتائج univariate_pool.csv لسجلات فحص المرشحين.
- استخدم حزمة rms rcs() لتناسب مصطلحات spline للمتغيرات المستمرة واختبار الحدود غير الخطية. قم بتثبيت عدد ومواقع العقد وسجلها في تعليقات السكريبت لتجنب الاختلافات بين الجولات.
- خذ اتحاد المتغيرات المحددة سريريا والمتغيرات مع المتغيرات الأحادية P<0.20 لإنشاء قائمة المرشحين، واحفظها كمخرجات/candidate_list.txt. تحقق مما إذا كان عدد المرشحين وعدد الأحداث يحقق قيد التعقيد، وسجل مجموعة المتغيرات النهائية المستخدمة في LASSO.
- اختيار LASSO وتركيب النموذج الرئيسي
- استخدم model.matrix() لإنشاء مصفوفة المتغيرات المستقلة والتأكد من صحة مستويات المرجع. استخدم cv.glmnet(x, y, family="binomial", alpha=1, standardize=TRUE) لإجراء التحقق المتقاطع عشرة أضعاف.
- استخدم حلقة لتكرار تقسيمات الطي العشوائية 1000 مرة، وسجل lambda.min وخطأ التحقق المتقاطع لكل تشغيل، واكتب الرسائل إلى الإخراج/lasso_lambda_trace.csv. خذ وسيط قيمة lambda.min البالغ 1000 كλ المستقر، وقم بتركيب glmnet على مجموعة البيانات الكاملة باستخدام هذه λ للحصول على المتغيرات ذات المعاملات غير الصفرية.
- في مجموعات البيانات العشر المحسوبة، يتناسب مع النموذج الرئيسي glm(family=binomial) باستخدام المتغيرات المختارة. اجمع بين β وSE وOR و95٪ CI باستخدام قواعد روبين، واحفظ النتائج كإخراج أو multivariable_pool.csv.
- مخرجات تحويل الدرجات، الرسم النوموغرامي، وتقسيم المخاطر
- حول β إلى نقاط صحيحة وفقا لمقياس ثابت، وربط النتيجة الكلية إلى اللوجيت والاحتمال باستخدام نفس المقياس. قم بأخذ عينات عشوائية من عدة حالات لمقارنة "احتمال النموذج الأصلي" مع "احتمال تحويل الدرجة" للتأكد من أن الخطأ ضمن النطاق المقبول المحدد مسبقا.
- استخدم حزمة rms لبناء النوموجرام وتصديره كإخراج أو nomogram.pdf. تصدير ورقة النوتة كصيغ إخراج/score_sheet.xlsx أو ما يعادلها لتسهيل إنتاج نموذج ورقي من صفحة واحدة.
- إجراء الطبقات باستخدام عتبات 5٪ و15٪ وإصدار عدد الحالات والحدوث المرصود في كل طبقة. تصدير Se، Sp، PPV، وNPV عند العتبات كمخرجات/threshold_metrics.csv للعرض اللاحق.
- التحقق الداخلي، والمعايرة، ومنحنيات القرار
- اضبط عدد عينات إعادة التمهيد إلى 1000، وفي كل إعادة عينة قم بإعادة التركيب والتنبؤ، مع حساب AUC، ودرجة برير، وتقاطع المعايرة، وميل المعايرة. اكتب النتائج المصححة بالتفاؤل إلى bootstrap_corrected_metrics.csv المخرج/واحتفظ بملخص للنتائج الوسيطة لكل إعادة عينة للتدقيق.
- استخدم حزمة pROC roc() و auc() لحساب AUC، واستخدم ci.auc(method="delong") لإخراج CI بنسبة 95٪. احفظ منحنى ROC كإخراج/roc_curve.png وثبت الدقة والحجم في السكريبت.
- استخدم حزمة RMDA لحساب صافي الفائدة ورسم منحنيات DCA على نطاق العتبة من 2٪ إلى 25٪. تصدير الشكل كناتج/dca_curve.png وثبت حجم خطوة العتبة وتسميات المنحنى لضمان قابلية التكرار.
- تحليلات الحساسية
- استبدل النتيجة ب "مرحلة إصابة الضغط ≥3" وكرر عملية التقييم، واختيار LASSO، وتركيب النموذج الرئيسي، وسير عمل تقييم الأداء. أخرج مقاييس الأداء إلى جانب التحليل الرئيسي لتقييم المتانة.
- استبدل متغير التعرض لانخفاض ضغط الدم ب "الوقت التراكمي مع MAP<60 مم زئبق" وكرر سير عمل النمذجة والتحقق الكامل. قارن اتجاهات المعاملات والتغيرات في AUC، ودرجة بريير، ومقاييس المعايرة لاختبار الحساسية لإعداد العتبة.