$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Dit onderzoek werd goedgekeurd door de National Center for Health Statistics (NCHS) Ethics Review Board, en er werd geïnformeerde toestemming verkregen van alle deelnemers. Publiek beschikbare NHANES-gegevens werden gebruikt in overeenstemming met de relevante richtlijnen.
Studiepopulatie
De National Health and Nutrition Examination Survey (NHANES), uitgevoerd door de NCHS, levert nationaal representatieve gegevens over de gezondheid en voedingsstatus van de niet-geïnstitutionaliseerde Amerikaanse bevolking. Gegevens uit de cycli 2001–2004 en 2009–2010 werden opgehaald. Deelnemers van ≥18 jaar werden opgenomen, terwijl personen met ontbrekende gegevens over chronische LBP, dieetmicronutriënten, opleidingsniveau of belangrijke covariaten — waaronder armoede-inkomensverhouding (PIR), body mass index (BMI), rookstatus, hypertensie, diabetes en alcoholconsumptie — werden uitgesloten. De uiteindelijke analytische cohort bestond uit 8.710 deelnemers (Figuur 1).
Beoordeling van voedingsstoffen
Gegevens over de voedingsinname werden verkregen met behulp van twee 24-uurs dieetherinneringsinterviews. Het eerste interview vond plaats in een mobiel examencentrum, en het tweede werd binnen enkele dagen telefonisch afgerond. De geautomatiseerde meervoudige doorlaatmethode werd toegepast om de nauwkeurigheid van herinneringen te verbeteren en rapportagevooroordelen te minimaliseren. Gedetailleerde AMPM-procedures worden beschreven in de NHANES-voedingsmethodologiedocumentatie.
Beoordeling van chronische LBP
De status van chronische LBP werd vastgesteld met behulp van vragenlijstgegevens. Deelnemers werden geclassificeerd als LBP als zij rugpijn rapporteerden die ≥3 maanden aanhield.
Covariaten
Demografische kenmerken omvatten leeftijd, geslacht, ras/etniciteit (Mexicaans-Amerikaans, Overige Hispanic, Niet-Spaanse blanke, niet-Hispanic zwart, andere ras), opleiding (<9e klas, 9–11e), diploma (middelbare schooldiploma/GED, enkele universitaire of AA-diploma's en universitaire afstudeer), gezinsinkomen-armoedeverhouding (PIR), body mass index (BMI), rookstatus, evenals drinkstatus en geschiedenis van hypertensie en diabetes. Hypertensie werd gedefinieerd op basis van zelfgerapporteerde diagnose door de arts of het huidige gebruik van antihypertensiva. Diabetes werd gedefinieerd met zelfgerapporteerde diagnose door de arts, verhoogde glucosewaarden na een glucosetolerantietest van 2 uur (≥11,1 mmol/L), of geglyceerde hemoglobine ≥6,5% of nuchtere glucose ≥7,0 mmol/L. Prediabetes werd gedefinieerd als door een arts te horen dat je prediabetes had, 2 uur OGTT glucosetolerantiewaarden van 7,8–11,1 mmol/L, nuchtere glucose van 6,1–6,9 mmol/L, of een HbA1c-waarde van tussen 5,7% en 6,4%. De status van roken werd gedefinieerd als (1) Niet-actief vóór de laatste maand, of gestopt met meer dan een jaar, en (2) Roken momenteel (zelfrapportage recent) wordt gedefinieerd als meer dan 2 sigaretten per dag roken na het herwinnen van de gewoonte of vóór het stoppen. Wat betreft de drinkstatus bestaan er twee definities: (1) levenslange niet-drinkers zonder eerder alcoholgebruik (<12 totaal), en (2) huidige drinkers die aangeven ≥12 keer per jaar te hebben gedronken, of minstens zes keer in de afgelopen 12 maanden. BMI berekend als gewicht ÷ lengte2.
Feature preprocessing en selectie voor machine learning
In totaal werden 52 variabelen opgenomen, waaronder 45 continue en 7 categorische variabelen. Alle preprocessing- en featureselectieprocedures—waaronder multicollinearity removal, SMOTE-toepassing en Boruta-gebaseerde featureselectie—werden exclusief uitgevoerd op de trainingsset om datalekkage te voorkomen. Om multicollineariteit te elimineren, werden eerst de variabelen met correlatiecoëfficiënten groter dan 0,9 verwijderd, en daarna de variabelen met VIFs groter dan 3.
Om de klassenonbalans te verbeteren en de herkenning van minderheidsklassen (lees: moeilijkere) klassen te verbeteren, werd de Synthetic Minority Over-sampling Technique (SMOTE) toegepast (Aanvullende Figuur 1), die afstanden creëert tussen steekproeven van minderheidsklassen, de k-dichtstbijzijnde buren voor die afstanden berekent en synthetische data in willekeurige richtingen genereert om klassen in balans te brengen. Alle gebruikte variabelen werden vervolgens gestandaardiseerd.
De selectie van kenmerken werd voltooid met Boruta op Random Forests, waarbij zogenaamde "schaduwkenmerken" werden gegenereerd en deze werden "getest" tegen de invoerkenmerken over 500 iteraties. Degenen die als "bevestigd" werden geclassificeerd, werden behouden om het model te bouwen. Voor gedetailleerde informatie over hyperparameterafstemming, zie Aanvullende Tabel 1.
Statistische analyse
Alle analyses werden uitgevoerd volgens de analytische richtlijnen van NHANES. Continue variabelen werden gerapporteerd als gemiddelde ± standaardafwijkingen (SD), terwijl categorische variabelen werden uitgedrukt als aantal en percentage. Verschillen tussen groepen werden getest met chi-kwadraat- of Student's t-tests waar passend.
Om overfitting te voorkomen, werden de gegevens willekeurig gesplitst in een trainingsset7 en een testset. Met MLR3 werden zes machine learning-algoritmen gebouwd: Random Forest, dat veel beslissingsbomen construeert en hun voorspellingen samenvoegt, met sterke generalisatiekracht en consistentie; LightGBM gebaseerd op gradiënt-boosted decision trees, die is afgestemd op efficiëntie en schaal; K-KNN sorteert monsters op basis van hoe dicht ze bij hun buren liggen; K-KNN heeft doorgaans betere resultaten op kleine niet-lineaire datasets; Naïeve Bayes, die door gebruik te maken van de stelling van Bayes eenvoud geeft en robuust is; SVM dat ideale hypervlakken bezit voor de classificatietaak, zelfs voor monsters in hoge dimensionaliteit; XGBoost is een vorm van gradient boosting ensemble die hoge prestaties heeft, zelfs met zeer grote onvolledige datasets.
Het model werd beoordeeld met behulp van nauwkeurigheid, F-beta, oppervlakte onder de ROC-curves (AUC-ROC), gevoeligheid, specificiteit en AUC-PR. De AUC-ROC is de belangrijkste meetwaarde die zij gebruiken, terwijl de andere meetwaarden diepere inzichten in prestaties geven. Ze valideren hun modellen met tienvoudige kruisvalidatie om stabiliteit te waarborgen. Verschillen in de prestaties van hun modellen werden geëvalueerd met behulp van ANOVA en de Kruskal–Wallis H-test.
Feature-interpreteerbaarheid werd onderzocht met behulp van SHapley Additive exPlanations (SHAP) en Local Interpretable Model-agnostic Explanations (LIME). SHAP leidt schattingen af van de bijdrage van elk kenmerk aan modelvoorspellingen met behulp van Shapley-waarden, gebaseerd op coöperatieve speltheorie, om zowel lineaire als interactieve effecten weer te geven. LIME haalt kenmerkinvloed door complexe modellen te benaderen met eenvoudigere, interpreteerbare surrogaten (bijv. Lasso-regressies) om individuele kenmerken te helpen verklaren. Alle analyses werden uitgevoerd in IBM SPSS Statistics versie 24.0 en in R versie 4.3.0. Een tweezijdige p-waarde < 0,05 wordt als statistisch significant beschouwd.