$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Demografische gegevens
In totaal voldeden 665 patiënten met heupfracturen aan de inclusiecriteria en werden ze in deze studie opgenomen (Figuur 1). Onder hen werden 44 patiënten (6,6%) binnen 48 uur na het letsel gediagnosticeerd met vroege MCVT, vastgesteld bij opname-echo. Een gedetailleerde vergelijking van de demografische en klinische kenmerken tussen de MCVT- en niet-MCVT-groepen wordt samengevat in Tabel 2.
LASSO-regressie en RCS-uitkomsten
Kandidaatvariabelen die in univariate analyse als statistisch significant (p-waarde < 0,01) zijn geïdentificeerd, samen met klinisch relevante parameters, werden verder gescreend met behulp van LASSO-regressie (Figuur 2). Elf kandidaatkenmerken werden geëvalueerd, waaronder negen continue variabelen en twee breuktype-dummyvariabelen. Het definitieve voorspellingsmodel behield TFIA, totaal cholesterol (TC), geactiveerde partiële tromboplastinetijd (APTT), NLR, PLR en TG. PLR toonde een grens-statistische significantie in het multivariabele model (p-waarde = 0,068), maar bleef behouden op basis van LASSO-screening en klinische plausibiliteit. Om de lineariteitsaanname te versoepelen, werd RCS-analyse uitgevoerd voor de zes continue voorspellers, waarbij werd aangepast voor de andere geïdentificeerde variabelen (Tabel 3). Er werd een significante niet-lineaire dosis-responsrelatie waargenomen tussen TG en het risico op MCVT (totale p-waarde < 0,001). Zoals geïllustreerd in Figuur 3, nam het MCVT-risico sterk toe over lagere tot midden TG-waarden en steeg daarna geleidelijker, met bredere onzekerheid bij hogere TG-niveaus. RCS-analyse werd uitgevoerd met 3 knopen. Voor TG werden knopen geplaatst op 0,62 mmol/L, 1,03 mmol/L en 1,64 mmol/L, wat overeenkomt met het 10e, 50e en 90e percentiel van TG-verdeling. Het schijnbare referentiepunt rond 1,03-1,04 mmol/L moet als verkennend worden beschouwd in plaats van als een gevalideerde klinische grens.
Multivariabele logistische regressieanalyse
Met de variabelen die via LASSO- en RCS-analyses zijn geïdentificeerd, werd een definitief multivariabel logistisch regressiemodel opgesteld. In dit model werden TFIA, TC, APTT, NLR en PLR als lineaire termen opgenomen, terwijl TG werd opgenomen als een beperkte kubische spline om rekening te houden met het niet-lineaire effect (Tabel 4). Zes variabelen bleven behouden in het uiteindelijke voorspellingsmodel: TFIA (OR = 1,026, 95% BI: 1,004–1,049, p-waarde = 0,021), TC (OR = 0,244, 95% BI: 0,146–0,406, p-waarde < 0,001), APTT (OR = 0,856, 95% BI: 0,764–0,959, p-waarde = 0,007), NLR (OR = 1,189, 95% BI: 1,103–1,280, p-waarde <0,001), PLR (OR = 1,004, 95% BI: 1,000–1,008, p-waarde = 0,068), en TG gemodelleerd met RCS. Het algehele effect van TG in het model was zeer significant (p-waarde < 0,001). Op basis van deze variabelen werd een klinisch nomogram ontwikkeld om individuele risicoschatting te vergemakkelijken (Figuur 4). De volledige regressiecoëfficiënten, modelintercept, TG-splinecoëfficiënten en knooplocaties worden weergegeven in Tabel 1. Voorspellerdefinities worden samengevat in Tabel 5.
Modelprestaties
De discriminatieprestaties van het model werden geëvalueerd met behulp van receiver operating characteristic (ROC) analyse (Figuur 5), wat een AUC van 0,894 (95% BI: 0,829–0,941) opleverde. De kalibratiecurve toonde een algemene overeenstemming tussen voorspelde kansen en waargenomen uitkomsten (Figuur 6), ondersteund door een Brier-score van 0,041 (bootstrap 95% BI: 0,031–0,053). DCA suggereerde dat het gebruik van het model een hoger nettovoordeel kan opleveren dan zowel de "treat-all" als "treat-none" strategieën over een bereik van drempelkansen (Figuur 7). Aanvullende kalibratiemetrieken werden berekend. De schijnbare kalibratie-intercept was ongeveer 0, de schijnbare kalibratiehelling was 1,000 en de kalibratie-in-de-grote was ongeveer 0. De door bootstrap gecorrigeerde kalibratiehelling was 0,896, wat wijst op bescheiden optimisme na interne validatie. Interne validatie met 1.000 bootstrap-resamples toonde beperkte optimisme, met een optimisme-gecorrigeerde C-index van 0,884 (bootstrap 95% BI: 0,837-0,946) vergeleken met de schijnbare AUC van 0,894. De bias-gecorrigeerde kalibratiecurve benadeerde de ideale lijn, met een gemiddelde absolute fout van 0,013. Bootstrap 95% betrouwbaarheidsintervallen voor DCA-nettobateschattingen zijn weergegeven in Figuur 7. De Hosmer-Lemeshow goodness-of-fit-test toonde χ2 = 21,77, df = 8, p-waarde = 0,005, wat wijst op een discrepantie tussen voorspelde en waargenomen waarschijnlijkheden over risicodecilen. Daarom moet kalibratie zorgvuldig worden geïnterpreteerd en opnieuw beoordeeld in externe cohorten.
Beschikbaarheid van gegevens:
De gedeïdentificeerde individuele deelnemerdataset, de voltooide TRIPOD+AI-checklist en de volledige R-statistische analysecode die wordt gebruikt voor gegevensverwerking, modelontwikkeling, validatie en prestatie-evaluatie worden geleverd als aanvullende bestanden 1–3. Deze materialen bevatten alle gegevens en analytische middelen die nodig zijn om de bevindingen uit deze studie te reproduceren.

Figuur 1: Deelnemerstroomdiagram. Dit stroomdiagram vat de screening van patiënten, inclusie- en exclusiecriteria samen, en de uiteindelijke studiecohort die wordt gebruikt voor modelontwikkeling. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 2: LASSO-regressie voor de selectie van voorspellers. De grafiek toont kruisgevalideerde binomiale afwijking als functie van log(lambda) in 10-voudige kruisvalidatie. Elf kandidaatkenmerken werden geëvalueerd, bestaande uit negen continue variabelen en twee breuktype-dummyvariabelen. De bovenste getallen geven het aantal niet-nul coëfficiënten aan bij elke lambda-waarde. Foutbalken geven standaardfouten aan. De verticale stippellijnen stellen lambda.min (0,00526) en lambda,1se (0,02558) voor. Afkortingen: CCI = Charlson Comorbiditeitsindex; TFIA = Tijd van blessure tot opname; WBC = Witte bloedcellen; N = Neutrofiel; L = Lymfocyt; M = Monocyt; HB = Hemoglobine; HCT = Hematocriet; PLT = Bloedplaatje; NLR = Neutrofiel-tot-lymfocytenverhouding; PLR = Bloedplaatje-lymfocytenverhouding; TP = Totaal eiwit; Alb = Albumine; Tbil = Totale bilirubine; Crea = Creatinine; TC = Totaal cholesterol; GFR = Glomerulaire filtratiesnelheid; IP = Anorganisch fosfaat; INR = Internationaal genormaliseerde verhouding; Fbg = Fibrinogen; APTT = Geactiveerde gedeeltelijke tromboplastinetijd; TT = Trombinetijd; PT = Protrombinetijd. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 3: RCS-analyse van TG en vroege MCVT. De curve toont aangepaste OR's voor vroege MCVT volgens TG-niveau (mmol/L). Het model was aangepast voor TC, NLR, TFIA, APTT en PLR en gebruikte 3 knopen bij 0,62 mmol/L, 1,03 mmol/L en 1,64 mmol/L. De referentiewaarde was het mediane TG-niveau (1,03 mmol/L). Het schaduwgebied vertegenwoordigt 95% BI. TG werd gemodelleerd als een continue niet-lineaire voorspeller, en de curve moet niet worden geïnterpreteerd als het definiëren van een gevalideerde klinische cutoff. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 4: Nomogram voor het voorspellen van vroege MCVT, gedetecteerd binnen 48 uur na het letsel. Het nomogram schat de voorspelde kans op vroege MCVT met behulp van TFIA (h), TG (mmol/L), TC (mmol/L), APTT (seconden), NLR en PLR. Voorspellerdefinities en meeteenheden zijn weergegeven in Tabel 5. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 5: ROC-curve van het uiteindelijke voorspellingsmodel. De ROC-curve toont het discriminatievermogen van het voorspellingsmodel voor vroege MCVT. De UK was 0,894 (95% BI: 0,829-0,941). Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 6: Kalibratiecurve van het uiteindelijke voorspellingsmodel. De kalibratiegrafiek vergelijkt voorspelde en waargenomen kansen van vroege MCVT. De ideale lijn staat voor perfecte kalibratie; De schijnbare curve geeft de modelprestaties weer op de ontwikkelingsdataset; en de bias-gecorrigeerde curve vertegenwoordigt bootstrap-gecorrigeerde kalibratie na 1.000 resamples. De Brier-score was 0,041, het schijnbare kalibratie-intercept was ongeveer 0, de schijnbare kalibratie-helling was 1,000 en de bootstrap-gecorrigeerde kalibratie-helling was 0,896. Klik hier om een grotere versie van deze figuur te bekijken.

Figuur 7: DCA van het uiteindelijke voorspellingsmodel. DCA toont het netto voordeel van het nomogram over drempelwaarschijnlijkheden. De modelcurve geeft het netto voordeel van het gebruik van het nomogram weer; de treat-all-lijn gaat ervan uit dat alle patiënten MCVT ontwikkelen; en de treat-non-lijn gaat ervan uit dat geen enkele patiënt MCVT ontwikkelt. De curve moet alleen tijdens interne validatie worden geïnterpreteerd als bewijs van potentiële klinische bruikbaarheid, omdat er geen gevalideerde behandelingsdrempel is vastgesteld. Klik hier om een grotere versie van deze figuur te bekijken.
| Voorspeller/term | Transformatie / definitie | Coëfficiënt beta |
| Intercept | Modelinterceptie | -1.621145541 |
| TFIA | Lineair; uren van blessure tot opname | 0.025977241 |
| NLR | Lineair; Neutrofiel-lymfocytenverhouding | 0.172703853 |
| PLR | Lineair; Bloedplaatje-lymfocytenverhouding | 0.004034446 |
| TG_linear | RCS lineaire term; TG in mmol/L | 6.24150855 |
| TG_rcs1 | RCS niet-lineaire term; knopen bij 0,62, 1,03, 1,64 mmol/L | -4.814103739 |
| TC | Lineair; Totaal cholesterol in mmol/L | -1.411933879 |
| APTT | Lineair; Seconden | -0.155260603 |
| TG RCS-knopen | 0,62, 1,03, 1,64 mmol/L | |
| Voorspelde waarschijnlijkheid | P = 1 / (1 + exp(-LP)) | |
TABEL 1: Volledige specificatie van het voorspellingsmodel. Deze tabel toont het volledige voorspellingsmodel, inclusief het intercept, regressiecoëfficiënten, beperkte kubieke splinecoëfficiënten voor TG, en knooplocaties die nodig zijn voor modelimplementatie en onafhankelijke validatie.
| Variabele | MCVT-groep | Niet-MCVT-groep | p-waarde |
| Gemiddelde ± SD of n (%) | (n = 44) | (n = 621) |
| Geslacht (man/vrouw) | 15/29 (34.1/65.9) | 224/397 (36.1/63.9) | 0,872a |
| Leeftijd | 79,68 ± 12,03 | 75,66 ± 15,65 | 0,041b |
| BMI | 22.05 ± 2.74 | 21.94 ± 3.12 | 0,796b |
| Hypertensie (ja/nee) | 15/29 (34.1/65.9) | 196/425 (31.6/68.4) | 0,739a |
| Diabetes (ja/nee) | 6/38 (13.6/86.4) | 91/530 (14.7/85.3) | 1.000a |
| Coronaire ziekte (ja/nee) | 1/43 (2.3/97.3) | 29/592 (4.7/95.3) | 0,713a |
| Breuktype (femorale nek/intertrochanter/subtrochanterisch) | 22/19/3 (50.0/43.2/6.8) | 462/153/6 (74.4/24.6/1.0) | <0,001c |
| CCI | 2.55 ± 2.12 | 2.31 ± 2.09 | 0,472b |
| TFIA (uren) | 25.58 ± 18.89 | 16.22 ± 14.60 | 0,002b |
| WBC-telling (x 10⁹/L) | 9.19 ± 2.92 | 8.90 ± 3.18 | 0,529b |
| N aantal (*10⁹/L) | 7,48 ± 2,99 | 4.20 ± 3.11 | <0,001b |
| L-telling (x 10⁹/L) | 0.99 (0.78-1.21) | 1.15 (0.92-1.52) | 0,003d |
| M-telling (*10⁹/L) | 0.57 (0.44-0.71) | 0.59 (0.43-0.75) | 0,941d |
| HB-telling (g/L) | 105.13 ± 22.33 | 110.02 ± 21.16 | 0,165b |
| HCT | 32.07 ± 6.46 | 33.62 ± 6.09 | 0,127b |
| PLT-aantal (x 10⁹/L) | 195.06 ± 69.21 | 165,82 ± 65,07 | 0,009b |
| NLR | 6.16 (4.30-12.05) | 3.21 (1.56-5.20) | <0,001d |
| PLR | 181.51 (126.19-260.71) | 131.76 (96.00-176.23) | <0,001d |
| TP (g/L) | 61,79 ± 7,25 | 62.33 ± 7.31 | 0,637b |
| Alb (g/L) | 35,77 ± 5,90 | 36.58 ± 4.91 | 0,381b |
| Tbil (μmol/L) | 11.20 (9.23-14.38) | 12.70 (9.00-16.70) | 0,117d |
| Crea (μmol/L) | 70.50 (58.80-89.50) | 67.00 (57.47-82.00) | 0,280d |
| Ureum (mmol/L) | 6.57 (5.28-9.15) | 6.18 (4.92-8.14) | 0,121d |
| TG (mmol/L) | 1.24 (1.03-1.62) | 1.01 (0.77-1.27) | <0,001d |
| TC (mmol/L) | 3.11 ± 0.76 | 3,84 ± 0,95 | <0,001b |
| GFR (ml/min) | 70.1 ± 25.32 | 78.04 ± 24.45 | 0,050b |
| IP (mmol/L) | 1.04 ± 0.24 | 1.06 ± 0.24 | 0,787b |
| INR | 1.04 (0.99-1.07) | 1.05 (1.01-1.10) | 0,094d |
| Fbg (g/L) | 3.7 ± 1.01 | 4.04 ± 1.06 | 0,034b |
| APTT (tweede) | 28.16 ± 3.23 | 29,48 ± 3,88 | 0,012b |
| TT (tweede) | 15.95 (15.30-16.95) | 16.10 (15.60-16.80) | 0,533d |
| PT (tweede) | 11.90 (11.40-12.50) | 12.27 (11.70-12.80) | 0,064d |
| D-dimeer (μg/L) | 7186.00 (4,447.25-17,346.75) | 4780.90 (1,980.00-12,140.00) | 0,017d |
| De gegevens worden gepresenteerd als gemiddelde ± SD, mediaan (IQR), of n (%), afhankelijk van het gepast. aWelch's t-test; bWilcoxon rangsomtest; cChi-kwadraattest; dFiishers exacte test (Fisher-Freeman-Halton exacte test werd gebruikt voor multi-categorie variabelen wanneer de verwachte celtellingen klein waren). |
TABEL 2: Basislijnkenmerken van patiënten met en zonder MCVT. Continue variabelen worden gepresenteerd als gemiddelde ± standaarddeviatie of mediaan met interkwartielbereik, indien van toepassing. Categorische variabelen worden gepresenteerd als getallen en percentages. De uiteindelijke analytische cohort bevatte geen ontbrekende waarden voor kandidaatvoorspellers of uitkomstvariabelen; 12 gescreende patiënten werden uitgesloten vanwege ontbrekende belangrijke klinische of laboratoriumgegevens.
| Variabele | Knopen (10e, 50e, 90e percentielen) | Niet-lineaire p-waarde | Totale p-waarde |
| TC | 2.61, 3.77, 4.92 | 0.369 | <0,001 |
| NLR | 0.708, 3.35, 8.68 | 0.058 | <0,001 |
| TG | 0.62, 1.03, 1.64 | <0,001 | <0,001 |
| TFIA | 2, 12, 48 | 0.327 | 0.072 |
| APTT | 25.5, 28.9, 33.9 | 0.554 | 0.003 |
| PLR | 74, 134, 250 | 0.923 | 0.163 |
TABEL 3: Niet-lineaire associaties beoordeeld met RCS-analyse. De tabel toont knooplocaties, totale p-waarden en niet-lineaire p-waarden voor geselecteerde voorspellers. RCS-modellen gebruikten 3 knopen geplaatst op het 10e, 50e en 90e percentiel van elke voorspellerverdeling.
| Variabelen | OF | CI-lower | CI-upper | p-waarde |
| TFIA | 1.026 | 1.004 | 1.049 | 0.021 |
| TC | 0.244 | 0.146 | 0.406 | <0,001 |
| APTT | 0.856 | 0.764 | 0.959 | 0.007 |
| NLR | 1.189 | 1.103 | 1.28 | <0,001 |
| PLR | 1.004 | 1 | 1.008 | 0.068 |
| TG (gemodelleerd met RCS) | - | - | - | <0,001 |
TABEL 4: Multivariabele logistische regressiemodel voor vroege MCVT. OR's, 95% BI's en p-waarden worden weergegeven voor lineaire termen. Omdat TG werd gemodelleerd met RCS, werden geen enkele OR en 95% BI gerapporteerd. De aangepaste dosis-responsrelatie wordt weergegeven in Figuur 3, en de splinecoëfficiënten en knooplocaties zijn weergegeven in Tabel 1.
| Voorspeller | Definitie | Eenheid | Timing |
| TFIA | Tijd van blessure tot opname | Openingstijden | Toelatingsbloedtest |
| TG | Triglyceriden | mmol/L | Toelatingsbloedtest |
| TC | Totaal cholesterol | mmol/L | Toelatingsbloedtest |
| APTT | Geactiveerde partiële tromboplastinetijd | Seconden | Toelatingsbloedtest |
| NLR | Neutrofieltelling / lymfocytentelling | Ratio | Toelatingsbloedtest |
| PLR | Aantal bloedplaatjes / lymfocytentelling | Ratio | Toelatingsbloedtest |
TABEL 5: Definities van voorspellingsmodelvariabelen. Deze tabel vat de voorspellers samen die in het uiteindelijke model zijn opgenomen, samen met hun definities, meeteenheden en, waar van toepassing, klinische interpretaties.
Aanvullend Dossier 1: TRIPOD+AI-rapportagechecklist voltooid. Dit bestand bevat de voltooide Transparent Reporting of a Multivariable Prediction Model for Individual Prognosis or Diagnosis Plus Artificial Intelligence (TRIPOD+AI) checklist, waarin de naleving van de rapportageaanbevelingen voor voorspellingsmodelstudies wordt gedocumenteerd. Klik hier om dit bestand te downloaden.
Aanvullend Bestand 2: Statistische analyse-uitkomsten voor de ontwikkeling en validatie van voorspellingsmodellen. Dit werkboek bevat de statistische analyse-uitkomsten die zijn gegenereerd tijdens modelontwikkeling en validatie, waaronder univariate analyses, LASSO-regressieresultaten, restrictie-kubische spline-analyses, multivariabele logistische regressiecoëfficiënten, modelprestatie-metrics, kalibratie-analyses, Hosmer-Lemeshow goodness-of-fit testresultaten, beslissingscurve-analyse en gerelateerde samenvattingstabellen. Klik hier om dit bestand te downloaden.
Aanvullend Bestand 3: R-script voor gegevensverwerking en statistische analyse. Dit bestand bevat het volledige R-script dat wordt gebruikt voor datapreprocessing, variabeleselectie, modelontwikkeling, interne validatie, kalibratie, discriminatieanalyse, beslissingscurveanalyse en het genereren van de statistische output die in deze studie worden gerapporteerd. Klik hier om dit bestand te downloaden.