$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
La ricerca è stata approvata dal Comitato Etico Biomedico dell'Università di Pechino (IRB00001052-11015). Tutti i soggetti hanno fornito il consenso informato.
Popolazione dello studio
I dati per questo studio sono stati ottenuti da CHARLS 2020, che copre 150 unità a livello di contea e 450 unità a livello di villaggio in tutto il paese, comprendendo circa 10.000 nuclei familiari e 19.395 individui di mezza età e anziani di 45 anni e oltre21 anni. Il sondaggio utilizzava un approccio di campionamento probabilistico a più fasi. Gli individui sono stati inclusi se soddisfacevano i seguenti criteri: età di 45 anni o più; ha selezionato la residenza domiciliare in risposta all'item del questionario BA007: Qual è il tipo di residenza all'indirizzo abitativo?; e ha fornito una risposta chiara alla scala della depressione del Center for Epidemiological Studies (CES-D-10) composta da 10 punti. Gli individui sono stati esclusi se soddisfacevano uno o più dei seguenti criteri: risposte che indicano la residenza non domiciliare (inclusi comunità, ospedali e altri o scomparsi); mancanza di informazioni sulla depressione; mancanza di informazioni sulle covariabili rilevanti. Alla fine, 14.466 adulti che soddisfacevano i criteri iniziali sono stati inclusi in questo studio. Il processo dettagliato di selezione dei partecipanti è illustrato nella Figura 1.

Figura 1: Un diagramma di flusso per la selezione della popolazione dello studio. Clicca qui per visualizzare una versione più grande di questa figura.
Definizione delle variabili
I sintomi depressivi sono stati valutati utilizzando una versione semplificata del CES-D-10, che ha alta affidabilità e validità negli studiprecedenti 22,23. Il CES-D-10 ha 10 elementi con quattro scale: 0 = mai, 1 = a volte o raramente, 2 = spesso e 3 = sempre, e il quinto e l'ottavo item sono valutati alcontrario di 24. Il punteggio totale della scala della depressione si ottiene sommando le valutazioni di tutti e 10 gli elementi, che vanno da 0 a 30; più alto è il punteggio, più gravi sono i sintomi depressivi dell'individuo. Sulla base di studi precedenti su popolazioni di mezza età eanziane, 25,26, i partecipanti con punteggi CES-D-10 ≥10 sono stati classificati come presenti di sintomi depressivi, mentre i punteggi <10 sono stati classificati come non depressi.
Le informazioni raccolte da tutti i partecipanti allo studio includono dati demografici (età, genere, stato civile, residenza e livello di istruzione), abitudini di vita e stato di salute fisica (attività fisica, fumo, alcol, interazioni sociali, autovalutazione della salute, stato di malattia cronica e tempo di pisolino) e informazioni relative alle condizioni di vita e ai figli (soddisfazione delle condizioni di vita, stato civile della prole, e lo stato di salute della prole). Le malattie croniche vengono diagnosticate dai medici e includono ipertensione, dislipidemia, tumori maligni, malattie polmonari croniche, malattie cardiache, ictus, artrite e reumatismo.
Analisi statistica
Analisi di regressione logistica binaria
Le variabili categoriche sono presentate come frequenze e percentuali. Sono stati condotti test chi-quadrato per esaminare le differenze nell'incidenza della depressione tra varie caratteristiche demografiche all'interno della popolazione studiata. È stata eseguita un'analisi di regressione logistica binaria per identificare fattori significativi associati ai sintomi depressivi. Queste analisi statistiche sono state effettuate utilizzando il software SPSS (versione 21.0).
Regressione LASSO
Nella versione 4.3.2 di R, il dataset veniva suddiviso casualmente in un set di addestramento e un set di test con un rapporto 7:3 utilizzando un seed casuale fisso per garantire la riproducibilità. In particolare, il 70% dei dati è stato allocato al set di addestramento per lo sviluppo del modello e il 30% al set di test per la validazione indipendente. Innanzitutto, le variabili statisticamente significative (p < 0,05) nella regressione logistica binaria sono state incluse nell'analisi di regressione LASSO per prevenire l'overfitting riducendo i coefficienti. Il valore lambda ottimale è stato determinato tramite validazione incrociata a 10 volte. Il lambda.1se è stato selezionato sulla base di una regola di errore standard per ottenere la performance ottimale nella selezione delle variabili. Il ritiro LASSO è stato impiegato per prevenire il sovrafitting riducendo i coefficienti di variabili meno importanti verso zero. Per ulteriori dettagli, si prega di consultare il Fascicolo Supplementare 1 (punto 1).
Costruzione del Nomogramma
L'importanza relativa dei predittori selezionati è stata classificata in base all'entità assoluta dei loro coefficienti standardizzati e visualizzata utilizzando un grafico forestale. Un nomogramma predittivo è stato costruito utilizzando il pacchetto rms, basato sull'equazione finale di regressione logistica. Per ulteriori dettagli, si prega di consultare il Fascicolo Supplementare 1 (punto 2).
Validazione del modello
La discriminazione veniva valutata calcolando l'area sotto la curva caratteristica operativa del ricevitore (AUC) utilizzando il pacchetto pROC. L'AUC veniva calcolato separatamente per il set di addestramento e il set di test. Tipicamente, un AUC superiore a 0,75 indica una buona discriminazione. Per ulteriori dettagli, si prega di consultare il Fascicolo Supplementare 1 (punto 3).
Calibrazione: La calibrazione del modello è stata valutata utilizzando il test di buona adattabilità di Hosmer-Lemeshow e visualizzata con curve di calibrazione generate dal pacchetto rms, tracciando le probabilità previste rispetto alle frequenze osservate. Un test Hosmer-Lemeshow non significativo (p > 0,05) indica una buona calibrazione. Per ulteriori dettagli, si prega di consultare il Fascicolo Supplementare 1 (punto 4).
Utilità clinica: L'utilità clinica è stata valutata utilizzando l'analisi della curva decisionale (DCA) con il pacchetto RMDA, calcolando il beneficio netto su probabilità di soglia da 0% a 100%. Per ulteriori dettagli, si prega di consultare il Fascicolo Supplementare 1 (punto 5).
Validazione della foresta casuale: come validazione supplementare, è stato implementato un modello di foresta casuale utilizzando il pacchetto randomForest per valutare la stabilità e la generalizzazione del modello. La curva di convergenza del tasso di errore è stata tracciata per valutare la relazione tra il numero di alberi e l'accuratezza della previsione sia nei set di addestramento che in quelli di test. Un valore p bilaterale < 0,05 è stato considerato indicativo di una differenza statisticamente significativa. Per ulteriori dettagli, si prega di consultare il Fascicolo Supplementare 1 (punto 6).