$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Het onderzoek werd goedgekeurd door de Biomedische Ethiek Commissie van de Peking Universiteit (IRB00001052-11015). Alle proefpersonen gaven geïnformeerde toestemming.
Studiepopulatie
De gegevens voor dit onderzoek zijn verkregen uit CHARLS 2020, dat 150 eenheden op countyniveau en 450 op dorpsniveau in het hele land omvat, met ongeveer 10.000 huishoudens en 19.395 personen van middelbare leeftijd en ouderen van 45 jaar enouder als 21 jaar. De enquête maakte gebruik van een meerstaps probability sampling-benadering. Individuen werden opgenomen als zij aan de volgende criteria voldeden: leeftijd van 45 jaar of ouder; geselecteerde woning als antwoord op het vragenlijstitem BA007: Wat is het type woning op het woonadres?; en gaf een duidelijk antwoord op de 10-items Center for Epidemiological Studies Depression Scale (CES-D-10). Personen werden uitgesloten als zij voldeden aan een of meer van de volgende criteria: antwoorden die aangaven niet-woonplaats (inclusief gemeenschappen, ziekenhuizen en andere of vermist); gebrek aan informatie over depressie; Gebrek aan informatie over relevante covariaten. Uiteindelijk werden 14.466 volwassenen die aan de oorspronkelijke criteria voldeden in deze studie opgenomen. Het gedetailleerde proces van deelnemersselectie wordt geïllustreerd in Figuur 1.

Figuur 1: Een stroomdiagram voor de selectie van de populatie. Klik hier om een grotere versie van deze figuur te bekijken.
Definitie van de variabelen
Depressieve symptomen werden beoordeeld met behulp van een vereenvoudigde versie van de CES-D-10, die in eerdere studies een hoge betrouwbaarheid en validiteit heeft22,23. De CES-D-10 heeft 10 items met vier toonladders: 0 = nooit, 1 = soms of zelden, 2 = vaak, en 3 = altijd, en de vijfde en achtste items worden omgekeerdmet 24 beoordeeld. De totale score van de depressieschaal wordt verkregen door de beoordelingen van alle 10 items op te tellen, variërend van 0 tot 30; Hoe hoger de score, hoe ernstiger de depressieve symptomen van het individu zijn. Op basis van eerdere studies onder middelbare leeftijds- en oudere bevolkingsgroepenwerden 25,26 deelnemers met CES-D-10 scores ≥10 geclassificeerd als met depressieve symptomen, terwijl scores <10 als niet-depressief werden geclassificeerd.
De informatie die van alle deelnemers aan het onderzoek is verzameld, omvat demografische gegevens (leeftijd, geslacht, burgerlijke staat, woonplaats en opleidingsniveau), leefgewoonten en fysieke gezondheidstoestand (lichamelijke activiteit, roken, drinken, sociale interacties, zelfbeoordeelde gezondheid, de status van chronische ziekte en dutjestijden), en informatie over leefomstandigheden en kinderen (tevredenheid over de leefomstandigheden, de burgerlijke staat van het nageslacht, en de gezondheidstoestand van de nakomelingen). Chronische ziekten worden door artsen vastgesteld en omvatten hypertensie, dyslipidemie, kwaadaardige tumoren, chronische longziekten, hartziekten, beroerte, artritis en reuma.
Statistische analyse
Binaire logistische regressieanalyse
Categorische variabelen worden gepresenteerd als frequenties en percentages. Chi-kwadraattests werden uitgevoerd om verschillen in depressieincidentie te onderzoeken tussen verschillende demografische kenmerken binnen de studiepopulatie. Binaire logistische regressieanalyse werd uitgevoerd om significante factoren te identificeren die geassocieerd zijn met depressieve symptomen. Deze statistische analyses werden uitgevoerd met SPSS-software (versie 21.0).
LASSO-regressie
In R versie 4.3.2 werd de dataset willekeurig gepartitioneerd in een trainingsset en een testset in een verhouding van 7:3, waarbij een vaste willekeurige seed werd gebruikt om reproduceerbaarheid te waarborgen. Specifiek werd 70% van de data toegewezen aan de trainingsset voor modelontwikkeling, en 30% aan de testset voor onafhankelijke validatie. Ten eerste werden variabelen die statistisch significant waren (p < 0,05) in de binaire logistische regressie opgenomen in de LASSO-regressieanalyse om overfitting door het verkleinen van de coëfficiënten te voorkomen. De optimale lambdawaarde werd bepaald via 10-voudige kruisvalidatie. De lambda.1se werd geselecteerd op basis van de ene standaard foutregel om de optimale variabeleselectie te bereiken. LASSO-krimp werd toegepast om overfitting te voorkomen door de coëfficiënten van minder belangrijke variabelen naar nul te verlagen. Voor meer informatie kunt u het Aanvullende Dossier 1 (item 1) raadplegen.
Nomogramconstructie
Het relatieve belang van de geselecteerde voorspellers werd gerangschikt op basis van de absolute grootte van hun gestandaardiseerde coëfficiënten en gevisualiseerd met behulp van een bosplot. Een voorspellend nomogram werd geconstrueerd met behulp van het rms-pakket, gebaseerd op de uiteindelijke logistische regressievergelijking. Voor meer details, raadpleeg het Aanvullende Dossier 1 (item 2).
Modelvalidatie
Discriminatie werd geëvalueerd door het oppervlak onder de receiver operating characteristic curve (AUC) te berekenen met behulp van het pROC-pakket. De AUC werd afzonderlijk berekend voor de trainings- en testset. Typisch duidt een AU groter dan 0,75 op goede discriminatie. Voor meer details, raadpleeg Aanvullend Dossier 1 (item 3).
Kalibratie: Modelkalibratie werd beoordeeld met behulp van de Hosmer-Lemeshow goodness-of-fit test en gevisualiseerd met kalibratiecurves gegenereerd door het rms-pakket, waarbij voorspelde kansen werden uitgezet tegen waargenomen frequenties. Een niet-significante Hosmer-Lemeshow-test (p > 0,05) wijst op een goede calibratie. Voor meer details, raadpleeg het Aanvullende Dossier 1 (item 4).
Klinisch nut: Het klinische nut werd geëvalueerd met behulp van Decision Curve Analysis (DCA) met het RMDA-pakket, waarbij het netto voordeel over drempelkansen van 0% tot 100% werd berekend. Voor meer details, raadpleeg het Aanvullende Dossier 1 (punt 5).
Validatie van willekeurige bossen: Als aanvullende validatie werd een random forest-model geïmplementeerd met behulp van het randomForest-pakket om de stabiliteit en generaliseerbaarheid van het model te beoordelen. De foutsnelheidsconvergentiecurve werd uitgezet om de relatie tussen het aantal bomen en de voorspellingsnauwkeurigheid in zowel trainings- als testsets te evalueren. Een tweezijdige p-waarde < 0,05 werd beschouwd als een statistisch significant verschil. Voor meer details, raadpleeg het Aanvullende Dossier 1 (item 6).