$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Esta investigación fue aprobada por la Junta de Revisión de Ética del Centro Nacional de Estadísticas de la Salud (NCHS), y se obtuvo el consentimiento informado de todos los participantes. Los datos públicos de NHANES se utilizaron en cumplimiento de las directrices pertinentes.
Población estudiada
La Encuesta Nacional de Salud y Nutrición (NHANES), realizada por el NCHS, proporciona datos representativos a nivel nacional sobre el estado de salud y nutrición de la población estadounidense no institucionalizada. Se recuperaron datos de los ciclos 2001–2004 y 2009–2010. Se incluyeron participantes de ≥18 años, mientras que se excluyeron individuos con datos faltantes sobre la lombalgia lumbar crónica, micronutrientes dietéticos, nivel educativo o covariables clave —incluyendo la relación pobreza-ingresos (PIR), índice de masa corporal (IMC), estado de tabaquismo, hipertensión, diabetes y consumo de alcohol. La cohorte analítica final comprendió 8.710 participantes (Figura 1).
Evaluación de nutrientes dietéticos
Los datos de ingesta dietética se obtuvieron mediante dos entrevistas de retiro dietético de 24 horas. La primera entrevista se realizó en un Centro Móvil de Exámenes y la segunda se realizó por teléfono en pocos días. Se aplicó el Método Automatizado de Múltiples Pasadas para mejorar la precisión de la llamada y minimizar el sesgo en la notificación. Los procedimientos detallados de AMPM se describen en la documentación de metodología dietética de NHANES.
Evaluación de la lombalgia crónica
El estado de la lombalgia crónica se determinó mediante datos de cuestionarios. Los participantes fueron clasificados como con lombalesis si reportaban dolor de espalda de duración ≥3 meses.
Covariables
Las características demográficas incluían edad, sexo, raza/etnia (mexicoamericano, otros hispanos, blancos no hispanos, negros no hispanos, otra raza), educación (<9º curso, 9º–11º curso), título universitario (diploma/GED, algunos títulos universitarios/AA y graduados universitarios), ratio ingreso-pobreza familiar (PIR), índice de masa corporal (IMC), estado de tabaquismo, así como estado de consumo y antecedentes de hipertensión y diabetes. La hipertensión se definió en función del diagnóstico autoinformado por el médico o del uso actual de medicación antihipertensiva. La diabetes se definió mediante un diagnóstico autoinformado por el médico, niveles elevados de glucosa tras una prueba de tolerancia a la glucosa de 2 horas (≥11,1 mmol/L), hemoglobina glicada ≥6,5% o glucosa en ayunas ≥7,0 mmol/L. La prediabetes se definió como que un médico le dijera que tenías prediabetes, valores de tolerancia a la glucosa OGTT a 2 horas de 7,8–11,1 mmol/L, glucosa en ayunas de 6,1–6,9 mmol/L, o un valor de HbA1c entre el 5,7% y el 6,4%. El estado de fumar se definía como (1) No actual antes del último mes, o dejado de fumar durante más de un año, y (2) fumar actualmente (autoinforme reciente) se define como fumar más de 2 cigarrillos al día al recuperar el hábito o antes de dejar de fumar. En cuanto al estado de consumo de alcohol, existen dos definiciones: (1) no bebedores de por vida sin consumo previo de alcohol (<12 en total), y (2) bebedores actuales que declaran haber bebido ≥12 veces al año, o haber bebido al menos seis veces en los últimos 12 meses. IMC calculado como peso ÷ altura2.
Preprocesamiento y selección de características para aprendizaje automático
Se incluyeron un total de 52 variables, incluyendo 45 continuas y 7 categóricas. Todos los procedimientos de preprocesamiento y selección de características —incluyendo la eliminación de multicolinealidad, la aplicación SMOTE y la selección de características basada en Boruta— se realizaron exclusivamente en el conjunto de entrenamiento para evitar fugas de datos. Para eliminar la multicolinealidad, primero se eliminaron las variables con coeficientes de correlación superiores a 0,9 y luego las variables con FIF mayores a 3.
Para mejorar el desequilibrio de clases y el reconocimiento de clases minoritarias (léase: más difíciles), se aplicó la Técnica de Sobremuestreo de Minorías Sintéticas (SMOTE) (Figura Suplementaria 1), que crea distancias entre muestras de clases minoritarias, calcula los k vecinos más cercanos para esas distancias y genera datos sintéticos en direcciones aleatorias para equilibrar las clases. A continuación, todas las variables utilizadas se estandarizaron.
La selección de características se completó con Boruta en Random Forests, generando las llamadas "características sombra" y "probándolas" frente a las características de entrada a lo largo de 500 iteraciones. Los clasificados como "confirmados" fueron retenidos para construir el modelo. Para información detallada sobre la sintonización de hiperparámetros, consulte la Tabla Suplementaria 1.
Análisis estadístico
Todos los análisis se realizaron de acuerdo con las directrices analíticas de NHANES. Las variables continuas se informaron como media ± desviaciones estándar (DS), mientras que las variables categóricas se expresaron como recuento y porcentaje. Las diferencias entre grupos se evaluaron utilizando chi-cuadrado o pruebas t de Student cuando era apropiado.
Para evitar el sobreajuste, los datos se particionaban aleatoriamente en un conjunto de entrenamiento7 y un conjunto de prueba. Usando MLR3, se construyeron seis algoritmos de aprendizaje automático: Random Forest, que construye muchos árboles de decisión y une sus predicciones, con un fuerte poder de generalización y consistencia; LightGBM basado en árboles de decisión potenciados por gradiente, que está ajustado para eficiencia y escala; K-KNN ordena las muestras según la proximidad que tienen a sus vecinas; K-KNN suele obtener mejores resultados en conjuntos de datos pequeños y no lineales; Naive Bayes, que utilizando el teorema de Bayes da simplicidad y es robusto; SVM que posee hiperplanos ideales para la tarea de clasificación, incluso para muestras de alta dimensionalidad; XGBoost es una forma de conjunto de aumento de gradiente que posee un alto rendimiento, incluso con conjuntos de datos muy grandes e incompletos.
El modelo se evaluó utilizando precisión, F-beta, área bajo las curvas ROC (AUC-ROC), la sensibilidad, especificidad y AUC-PR. El AUC-ROC es la métrica principal que utilizan, mientras que las otras métricas ofrecen una visión más profunda del rendimiento. Validan sus modelos mediante validación cruzada de diez veces para garantizar la estabilidad. Las diferencias en el rendimiento de sus modelos se evaluaron utilizando ANOVA y la prueba de Kruskal–Wallis H.
La interpretabilidad de características se investigó utilizando SHapley Aditivive ExPlanations (SHAP) y Local Interpretable Model-Agnostic Explanations (LIME). SHAP obtiene estimaciones de la contribución de cada característica a las predicciones del modelo utilizando valores de Shapley, basados en la teoría cooperativa de juegos, para representar tanto efectos lineales como interactivos. LIME obtiene la influencia de características aproximando modelos complejos con sustitutos más simples e interpretables (por ejemplo, regresiones de Lazo) para ayudar a explicar la influencia de características individuales. Todos los análisis se realizaron en IBM SPSS Statistics versión 24.0 y en R versión 4.3.0. Un valor p bilateral < 0,05 se considera estadísticamente significativo.