$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Dati demografici
Un totale di 665 pazienti con fratture dell'anca soddisfacevano i criteri di inclusione e sono stati arruolati in questo studio (Figura 1). Tra questi, 44 pazienti (6,6%) sono stati diagnosticati con MCVT precoce rilevata all'ecografia di ammissione entro 48 ore dall'infortunio. Un confronto dettagliato delle caratteristiche demografiche e cliniche tra i gruppi MCVT e non MCVT è riassunto nella Tabella 2.
Regressione LASSO ed esiti RCS
Le variabili candidate identificate come statisticamente significative (p-value < 0,01) nell'analisi univariata, insieme ai parametri clinicamente rilevanti, sono state ulteriormente selezionate utilizzando la regressione LASSO (Figura 2). Sono state valutate undici caratteristiche candidate, tra cui nove variabili continue e due variabili fittiche di tipo frattura. Il modello finale di previsione ha mantenuto TFIA, colesterolo totale (TC), Tempo di Tromboplastina Parziale Attivata (APTT), NLR, PLR e TG. La PLR ha mostrato significatività statistica borderline nel modello multivariabile (p-valore = 0,068) ma è stata mantenuta sulla base dello screening LASSO e della plausibilità clinica. Per rilassare l'assunzione di linearità, è stata eseguita un'analisi RCS per i sei predittori continui, aggiustando per le altre variabili identificate (Tabella 3). È stata osservata una significativa relazione non lineare dose-risposta tra TG e rischio di MCVT (valore p complessivo < 0,001). Come illustrato nella Figura 3, il rischio MCVT è aumentato drasticamente nei valori TG bassi e medi e poi è aumentato più gradualmente, con un'incertezza più ampia ai livelli TG più alti. L'analisi RCS veniva eseguita utilizzando 3 nodi. Per TG, i nodi furono posizionati a 0,62 mmol/L, 1,03 mmol/L e 1,64 mmol/L, corrispondenti al 10°, 50° e 90° percentile della distribuzione TG. Il punto di riferimento apparente intorno a 1,03-1,04 mmol/L dovrebbe essere considerato esplorativo piuttosto che un cutoff clinico validato.
Analisi di regressione logistica multivariabile
Incorporando le variabili identificate tramite analisi LASSO e RCS, è stato costruito un modello finale di regressione logistica multivariabile. In questo modello, TFIA, TC, APTT, NLR e PLR erano inclusi come termini lineari, mentre TG era incorporato come spline cubica limitata per tenere conto del suo effetto non lineare (Tabella 4). Sei variabili sono state mantenute nel modello finale di previsione: TFIA (OR = 1,026, IC 95%: 1,004–1,049, valore p = 0,021), TC (OR = 0,244, IC 95%: 0,146–0,406, valore p < 0,001), APTT (OR = 0,856, IC 95%: 0,764–0,959, valore p = 0,007), NLR (OR = 1,189, IC 95%: 1,103–1,280, valore p <0,001), PLR (OR = 1,004, IC 95%: 1.000–1.008, p-valore = 0.068), e TG modellata con RCS. L'effetto complessivo di TG nel modello è stato altamente significativo (p-value < 0,001). Sulla base di queste variabili, è stato sviluppato un nomogramma clinico per facilitare la stima individualizzata del rischio (Figura 4). I coefficienti di regressione completi, l'intercetta del modello, i coefficienti di spline TG e le posizioni dei nodi sono forniti nella Tabella 1. Le definizioni dei predittori sono riassunte nella Tabella 5.
Prestazioni del modello
Le prestazioni discriminative del modello sono state valutate utilizzando l'analisi delle caratteristiche operative del ricevitore (ROC) (Figura 5), ottenendo un AUC di 0,894 (IC 95%: 0,829–0,941). La curva di calibrazione ha mostrato un accordo complessivo tra probabilità previste ed esiti osservati (Figura 6), supportata da un punteggio Brier di 0,041 (IC Bootstrap 95%: 0,031–0,053). La DCA ha suggerito che l'uso del modello potrebbe offrire un beneficio netto superiore rispetto sia alle strategie "treat-all" che "treat-all" in un intervallo di probabilità di soglia (Figura 7). Sono state calcolate ulteriori metriche di calibrazione. L'intercettazione apparente di calibrazione era approssimativamente 0, la pendenza apparente di calibrazione era 1,000 e la calibrazione in grande era circa 0. La pendenza di calibrazione corretta con bootstrap era di 0,896, suggerendo un modesto ottimismo dopo la validazione interna. La validazione interna con 1.000 ricampionamenti bootstrap ha mostrato ottimismo limitato, con un indice C corretto per ottimismo di 0,884 (IC bootstrap 95%: 0,837-0,946) rispetto all'apparente AUC di 0,894. La curva di calibrazione corretta per bias approssimava la linea ideale, con un errore assoluto medio di 0,013. Gli intervalli di confidenza del 95% per le stime del beneficio netto DCA sono mostrati nella Figura 7. Il test di bontà di adattamento di Hosmer-Lemeshow ha mostrato χ2 = 21,77, df = 8, p-valore = 0,005, suggerendo una certa discrepanza tra le probabilità previste e osservate tra i decili di rischio. Pertanto, la calibrazione deve essere interpretata con cautela e rivalutata in coorti esterne.
Disponibilità dei dati:
Il dataset individuale dei partecipanti de-identificato, la checklist completa TRIPOD+AI e il codice completo di analisi statistica R utilizzato per l'elaborazione dei dati, lo sviluppo del modello, la validazione e la valutazione delle prestazioni sono forniti come File Supplementari 1–3. Questi materiali contengono tutti i dati e le risorse analitiche necessarie per riprodurre i risultati riportati in questo studio.

Figura 1: Diagramma di flusso dei partecipanti. Questo diagramma di flusso riassume i criteri di screening dei pazienti, inclusione ed esclusione, e la coorte finale di studi utilizzata per lo sviluppo del modello. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 2: regressione LASSO per la selezione del predittore. Il grafico mostra la devianza binomiale incrociata come funzione di log(lambda) nella validazione incrociata a 10 volte. Sono state valutate undici caratteristiche candidate, comprendenti nove variabili continue e due variabili fittiche di tipo frattura. I numeri superiori indicano il numero di coefficienti non nulli a ciascun valore lambda. Le barre di errore indicano errori standard. Le linee tratteggiate verticali rappresentano lambda.min (0,00526) e lambda.1se (0,02558). Abbreviazioni: CCI = Indice di Comorbidità di Charlson; TFIA = Tempo dall'infortunio al ricovero; WBC = Globuli bianchi; N = Neutrofilo; L = Linfociti; M = monocito; HB = Emoglobina; HCT = Ematocrito; PLT = Piastrine; NLR = rapporto neutrofili/linfociti; PLR = Rapporto piastrine/linfociti; TP = Proteine totali; Alb = Albumina; Tbil = bilirubina totale; Crea = Creatinina; TC = colesterolo totale; GFR = Tasso di filtrazione glomerulare; IP = Fosfato inorganico; INR = Rapporto internazionale di normalizzazione; Fbg = Fibrinogeno; APTT = Tempo attivo di tromboplastina parziale; TT = Tempo di trombona; PT = Tempo della protrombina. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 3: Analisi RCS di TG e MCVT iniziali. La curva mostra le OR aggiustate per MCVT precoce in base al livello TG (mmol/L). Il modello fu regolato per TC, NLR, TFIA, APTT e PLR e utilizzava 3 nodi situati a 0,62 mmol/L, 1,03 mmol/L e 1,64 mmol/L. Il valore di riferimento era il livello mediano TG (1,03 mmol/L). L'area ombreggiata rappresenta il 95% di IC. TG è stata modellata come un predittore continuo non lineare, e la curva non deve essere interpretata come definendo un cutoff clinico validato. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 4: Nomogramma per prevedere la MCVT precoce rilevata entro 48 ore dopo l'infortunio. Il nomogramma stima la probabilità prevista di MCVT precoce utilizzando TFIA (h), TG (mmol/L), TC (mmol/L), APTT (secondi), NLR e PLR. Le definizioni dei predittori e le unità di misura sono fornite nella Tabella 5. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 5: Curva ROC del modello finale di previsione. La curva ROC mostra la capacità di discriminazione del modello di previsione per la MCVT precoce. L'AUC era di 0,894 (IC 95%: 0,829-0,941). Clicca qui per visualizzare una versione più grande di questa figura.

Figura 6: Curva di calibrazione del modello finale di previsione. Il grafico di calibrazione confronta le probabilità previste e osservate della MCVT precoce. La linea ideale rappresenta una calibrazione perfetta; La curva apparente rappresenta le prestazioni del modello sul dataset di sviluppo; e la curva corretta al bias rappresenta la calibrazione corretta dal bootstrap dopo 1.000 ricampionamenti. Il punteggio Brier era 0,041, l'intercettazione apparente di calibrazione era circa 0, la pendenza apparente di calibrazione era 1,000 e la pendenza di calibrazione corretta con bootstrap era 0,896. Clicca qui per visualizzare una versione più grande di questa figura.

Figura 7: DCA del modello finale di previsione. La DCA mostra il beneficio netto del nomogramma attraverso le probabilità soglia. La curva modello rappresenta il beneficio netto dell'uso del nomogramma; la linea tratta-tutto assume che tutti i pazienti sviluppino MCVT; e la linea del trattamento non trattato presuppone che nessun paziente sviluppi MCVT. La curva deve essere interpretata come prova di potenziale utilità clinica solo durante la validazione interna, poiché non è stata stabilita alcuna soglia di trattamento validata. Clicca qui per visualizzare una versione più grande di questa figura.
| Predittore/termine | Trasformazione / definizione | Coefficiente beta |
| Intercettazione | Intercetta del modello | -1.621145541 |
| TFIA | Lineare; Ore dal infortunio al ricovero | 0.025977241 |
| NLR | Lineare; Rapporto neutrofilo/linfociti | 0.172703853 |
| PLR | Lineare; Rapporto piastrine-linfociti | 0.004034446 |
| TG_linear | Termine lineare RCS; TG in mmol/L | 6.24150855 |
| TG_rcs1 | Termine non lineare RCS; nodi a 0,62, 1,03, 1,64 mmol/L | -4.814103739 |
| TC | Lineare; colesterolo totale in mmol/L | -1.411933879 |
| APTT | Lineare; Secondi | -0.155260603 |
| Nodi TG RCS | 0,62, 1,03, 1,64 mmol/L | |
| Probabilità prevista | P = 1 / (1 + exp(-LP)) | |
TABELLA 1: Specifica completa del modello di previsione. Questa tabella presenta il modello previsionale completo, inclusi l'intercetto, i coefficienti di regressione, i coefficienti di spline cubiche limitate per TG e le posizioni dei nodi necessarie per l'implementazione del modello e la validazione indipendente.
| Variabile | Gruppo MCVT | Gruppo non MCVT | Valore p |
| Media ± SD o n (%) | (n = 44) | (n = 621) |
| Genere (maschio/femmina) | 15/29 (34.1/65.9) | 224/397 (36.1/63.9) | 0,872a |
| Età | 79.68 ± 12.03 | 75.66 ± 15.65 | 0.041b |
| BMI | 22.05 ± 2.74 | 21.94 ± 3.12 | 0.796b |
| Ipertensione (sì/no) | 15/29 (34.1/65.9) | 196/425 (31.6/68.4) | 0,739a |
| Diabete (sì/no) | 6/38 (13.6/86.4) | 91/530 (14.7/85.3) | 1.000a |
| Malattia coronarica (sì/no) | 1/43 (2.3/97.3) | 29/592 (4.7/95.3) | 0,713a |
| Tipo di frattura (collo femorale/intertrocanterico/sottotrocanterico) | 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 (ore) | 25.58 ± 18.89 | 16.22 ± 14.60 | 0.002b |
| Conteggio dei globuli bili (x 10⁹/L) | 9.19 ± 2.92 | 8.90 ± 3.18 | 0,529b |
| Conteggio di N (*10⁹/L) | 7.48 ± 2.99 | 4.20 ± 3.11 | <0.001b |
| Conto L (x 10⁹/L) | 0.99 (0.78-1.21) | 1.15 (0.92-1.52) | 0,003d |
| M conteggio (*10⁹/L) | 0.57 (0.44-0.71) | 0.59 (0.43-0.75) | 0.941d |
| Conteggio HB (g/L) | 105.13 ± 22.33 | 110.02 ± 21.16 | 0.165b |
| HCT | 32.07 ± 6.46 | 33.62 ± 6.09 | 0.127b |
| Conteggio PLT (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 |
| Urea (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 (secondo) | 28.16 ± 3.23 | 29.48 ± 3.88 | 0.012b |
| TT (secondo) | 15.95 (15.30-16.95) | 16.10 (15.60-16.80) | 0,533d |
| PT (secondo) | 11.90 (11.40-12.50) | 12.27 (11.70-12.80) | 0,064d |
| D-Dimero (μg/L) | 7186.00 (4,447.25-17,346.75) | 4780.90 (1,980.00-12,140.00) | 0,017d |
| I dati sono presentati in media ± SD, mediana (IQR) o n (%), a seconda di opportunità. aIl test T di Welch; bTest di somma dei rangi di Wilcoxon; cTest del chi-quadrato; dTest esatto di Fiisher (il test esatto di Fisher-Freeman-Halton è stato utilizzato per variabili multi-categoria quando il numero atteso di celle era basso). |
TABELLA 2: Caratteristiche di base dei pazienti con e senza MCVT. Le variabili continue sono presentate come media ± deviazione standard o mediana con intervallo interquartile, a seconda dei conti. Le variabili categoriche sono presentate come numeri e percentuali. La coorte analitica finale non conteneva valori mancanti per i candidati predittori o variabili di risultato; 12 pazienti sottoposti a screening sono stati esclusi a causa della mancanza di dati chiave clinici o di laboratorio.
| Variabile | Nodi (percentile 10, 50, 90) | Valore p non lineare | P-valore complessivo |
| 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 |
TABELLA 3: Associazioni non lineari valutate tramite analisi RCS. La tabella presenta le posizioni dei nodi, i valori p complessivi e i valori p non lineari per i predittori selezionati. I modelli RCS utilizzavano 3 nodi posizionati al 10°, 50° e 90° percentile di ciascuna distribuzione predittiva.
| Variabili | OPPURE | CI-lower | CI-upper | p-valore |
| 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 (modellata con RCS) | - | - | - | <0.001 |
TABELLA 4: Modello di regressione logistica multivariabile per la MCVT precoce. Per termini lineari sono mostrati OR, IC al 95% e valori p. Poiché TG è stata modellata usando RCS, non sono stati riportati un singolo OR e un IC al 95%. La relazione dose-risposta rettificata è mostrata nella Figura 3, e i coefficienti delle spline e le posizioni dei nodi sono forniti nella Tabella 1.
| Predittore | Definizione | Unità | Tempistica |
| TFIA | Tempo dall'infortunio al ricovero | Orari | Esame del sangue di ammissione |
| TG | Trigliceridi | mmol/L | Esame del sangue di ammissione |
| TC | Colesterolo totale | mmol/L | Esame del sangue di ammissione |
| APTT | Tempo di tromboplastica parziale attivato | Secondi | Esame del sangue di ammissione |
| NLR | Conteggio dei neutrofili / conteggio dei linfociti | Rapporto | Esame del sangue di ammissione |
| PLR | Conteggio piastrinico / conteggio dei linfociti | Rapporto | Esame del sangue di ammissione |
TABELLA 5: Definizioni delle variabili del modello di previsione. Questa tabella riassume i predittori inclusi nel modello finale, insieme alle loro definizioni, unità di misura e, dove applicabile, interpretazioni cliniche.
File supplementare 1: Completamento della checklist di segnalazione TRIPOD+AI. Questo file contiene la Segnalazione Trasparente completa di un Modello di Previsione Multivariabile per Prognosi Individuale o Diagnostic Plus Intelligenza Artificiale (TRIPOD+AI), documentando la conformità alle raccomandazioni di reportistica per gli studi di modelli predittivi. Clicca qui per scaricare questo file.
File supplementare 2: Risultati di analisi statistica per lo sviluppo e la validazione dei modelli predittivi. Questo workbook contiene gli output dell'analisi statistica generati durante lo sviluppo e la validazione del modello, inclusi analisi univariate, risultati di regressione LASSO, analisi di spline cubiche limitate, coefficienti di regressione logistica multivariabile, metriche di performance del modello, analisi di calibrazione, risultati del test di buona connessione Hosmer-Lemeshow, analisi della curva decisionale e relative tabelle riassuntive. Clicca qui per scaricare questo file.
File supplementare 3: Script R per l'elaborazione dei dati e l'analisi statistica. Questo file contiene lo script R completo utilizzato per la pre-elaborazione dei dati, la selezione delle variabili, lo sviluppo del modello, la validazione interna, la calibrazione, l'analisi della discriminazione, l'analisi della curva decisionale e la generazione degli output statistici riportati in questo studio. Clicca qui per scaricare questo file.