Se muestra un esquema general del flujo de trabajo de preprocesamiento del EEG, extracción de características y clasificación en la Figura 1. Los materiales y el software utilizados en el estudio se enumeran en la Tabla de Materiales.

Figura 1: Flujo de trabajo del preprocesamiento de EEG, extracción de características y clasificación. Descripción esquemática del flujo de trabajo del estudio, que incluye verificación de los datos de EEG, filtrado paso banda y re-referenciación, extracción de las características de raíz media cuadrática (RMS), densidad espectral de potencia (PSD) y entropía, incorporación de la edad, división entre conjuntos de entrenamiento y prueba, desarrollo del modelo de Bosque Aleatorio, análisis de estabilidad de características, análisis SHAP, análisis del tamaño del efecto d de Cohen, selección integrada de características, ablación de características de entropía, reentrenamiento del modelo y evaluación del rendimiento. Haga clic aquí para ver una versión más grande de esta figura.
1. Adquisición de datos
El conjunto de datos de EEG en estado de reposo de acceso público se obtuvo del repositorio OpenNeuro (conjunto de datos ds004504, versión 1.0.9)33. Se verificaron las etiquetas de los participantes y la información demográfica. El conjunto de datos comprendía 88 participantes, incluyendo 36 con enfermedad de Alzheimer (AD), 23 con demencia frontotemporal (FTD) y 29 controles sanos (HC).
2. Importación de la grabación de EEG
Para cada participante, la grabación de EEG se cargó utilizando la biblioteca MNE-Python. Cada archivo de EEG se verificó para confirmar que era accesible, que podía importarse correctamente y que contenía un identificador válido del participante. Se excluyeron las grabaciones con archivos faltantes o dañados. En esta etapa se mantuvo la frecuencia de muestreo original del EEG de 500 Hz. Las señales se submuestrearon para reducir la carga computacional, preservando toda la información clínicamente relevante del EEG para este estudio.
3. Filtrado de paso de banda
Las grabaciones continuas de EEG se filtraron utilizando un filtro de paso de banda de cuarto orden con una frecuencia de corte baja de 0,5 Hz y una frecuencia de corte alta de 40 Hz.
4. Referenciación promedio común y epoquetización
Se aplicó la referenciación promedio común (CAR) a las grabaciones de EEG filtradas. En cada muestra de tiempo, se calculó la señal promedio a través de todos los canales de EEG disponibles y se restó de cada canal. Las señales de EEG referenciadas se conservaron para análisis posteriores. Cada grabación de EEG referenciada se dividió en épocas de duración fija utilizando la función make_fixed_length_epochs() en MNE-Python. Se utilizó una duración de época de 10 s con superposición de 0 s.
5. Validación de épocas
Para cada participante, solo se conservaron épocas completas de 10 s. Cualquier segmento restante de EEG más corto de 10 s al final de una grabación fue descartado. Cada época conservada se utilizó posteriormente como una muestra individual para la extracción de características.
6. Preparación de características
Aunque los datos del Mini-Mental State Examination (MMSE) estaban disponibles en el conjunto de datos original, se excluyeron del presente análisis. Las características extraídas del EEG se combinaron con la información demográfica de los participantes, incluyendo la edad y las etiquetas diagnósticas, obtenidas del archivo de metadatos participants.csv. El identificador del participante se utilizó como clave común, y solo se conservaron aquellos participantes que tenían registros de características del EEG e información demográfica correspondiente.
La matriz de características multidominio incluyó cinco características de raíz cuadrada media (RMS) (delta_rms, theta_rms, alpha_rms, beta_rms y gamma_rms), cinco características de densidad espectral de potencia (PSD) (delta_psd, theta_psd, alpha_psd, beta_psd y gamma_psd) y cinco características de entropía (delta_entropy, theta_entropy, alpha_entropy, beta_entropy y gamma_entropy). La edad se incluyó como una variable demográfica complementaria. Dado que las distribuciones de edad pueden diferir entre los grupos diagnósticos, no se pudo descartar por completo un sesgo demográfico. El grupo diagnóstico, definido como AD, HC o FTD, se asignó como etiqueta objetivo.
7. Preprocesamiento de datos y partición en conjuntos de entrenamiento y prueba
Durante el preprocesamiento, las grabaciones se submuestrearon de 500 Hz a 250 Hz para reducir los requisitos computacionales, manteniendo al mismo tiempo la información de frecuencia de EEG de interés.
El conjunto de datos se dividió en subconjuntos de entrenamiento (80 %) y prueba (20 %) a nivel de sujeto utilizando una estrategia de división agrupada. Las predicciones a nivel de sujeto se obtuvieron posteriormente mediante votación mayoritaria a través de las épocas predichas pertenecientes a cada participante. La evaluación principal utilizó una división agrupada de entrenamiento-prueba a nivel de sujeto para evitar que épocas del mismo participante aparecieran en ambos subconjuntos. Las predicciones a nivel de sujeto se obtuvieron mediante votación mayoritaria a través de las épocas predichas pertenecientes a cada participante. El conjunto de datos se examinó en busca de valores faltantes, y las observaciones faltantes se eliminaron o imputaron según correspondiera. Las etiquetas diagnósticas se codificaron. Un StandardScaler se ajustó con los datos de entrenamiento y luego se aplicó tanto a los conjuntos de datos de entrenamiento como de prueba.
8. Desarrollo del modelo de bosque aleatorio
Se inicializó un clasificador de bosque aleatorio utilizando 80 árboles, una profundidad máxima de árbol de 10, un máximo de cuatro características, un mínimo de cinco muestras por hoja y un estado aleatorio de 42. El clasificador se entrenó utilizando el conjunto de datos de entrenamiento estandarizado.
9. Evaluación del modelo
Se predijeron las etiquetas de clase para ambos conjuntos de datos, de entrenamiento y de prueba. Se construyó una matriz de confusión, y se calcularon la exactitud, la precisión, la recuperación y la puntuación F1 junto con el informe de clasificación. Se registraron las exactitudes de entrenamiento y de prueba.
Para la evaluación por participante, asegúrese de que todas las épocas de un participante determinado se asignen exclusivamente al subconjunto de entrenamiento o al de prueba. Se entrenó el clasificador de bosque aleatorio utilizando los mismos hiperparámetros que en el análisis principal.
10. Validación cruzada
Como procedimiento adicional de evaluación del modelo, se realizó una validación cruzada estratificada de cinco pliegues utilizando shuffle = True y random_state = 30. Se calcularon la precisión media y la desviación estándar a través de los cinco pliegues. Este análisis se consideró por separado de la evaluación principal por sujeto con conjunto de reserva.
11. Análisis de estabilidad de características
La estabilidad de las características se evaluó repitiendo el entrenamiento del bosque aleatorio 10 veces utilizando semillas aleatorias del 0 al 9. En cada ejecución, se registraron la precisión en el conjunto de prueba y las puntuaciones de importancia de las características. Se calcularon la media y la desviación estándar de la puntuación de importancia de cada característica a lo largo de las 10 ejecuciones, y los predictores se clasificaron según su estabilidad. Este análisis se utilizó para evaluar la consistencia de los rankings de características, y no para sustituir la evaluación principal a nivel de sujeto.
12. Análisis de inteligencia artificial interpretable
Se aplicó el TreeExplainer de SHAP al modelo de bosque aleatorio entrenado. Se calcularon los valores de SHAP para estimar la contribución de cada característica a las predicciones del modelo. Se calculó el valor absoluto medio de SHAP para cada característica, y las características se ordenaron según sus contribuciones SHAP. Se identificaron las características con valores de SHAP consistentemente bajos y se compararon con los resultados de los análisis de estabilidad de características y de d de Cohen. Las características que mostraron contribuciones consistentemente bajas se seleccionaron para la ablación subsiguiente y el reentrenamiento del modelo.
13. Análisis del tamaño del efecto estadístico
Se calculó la d de Cohen para cada biomarcador de EEG en las comparaciones de AD versus HC, AD versus FTD y FTD versus HC. Las magnitudes del tamaño del efecto se interpretaron utilizando umbrales de 0,20 para un efecto pequeño, 0,50 para un efecto medio y 0,80 para un efecto grande.
14. Selección integrada de características
Se compararon los resultados del análisis de importancia de características del bosque aleatorio, el análisis SHAP y el análisis de d de Cohen. Los predictores que mostraron consistentemente baja importancia de característica, baja contribución SHAP y tamaños del efecto pequeños se identificaron como candidatos para su eliminación.
15. Ablación de características
Se evaluaron las características de entropía en un análisis exploratorio de ablation, y se construyó una matriz reducida de características que contenía RMS, PSD y edad. Se entrenó nuevamente el clasificador de bosque aleatorio utilizando los mismos hiperparámetros. Se repitieron el entrenamiento del modelo, la prueba, la validación cruzada y el análisis de la curva característica de operación del receptor (ROC) utilizando el conjunto reducido de características.
16. Análisis de ROC
Las probabilidades de clase se obtuvieron a partir del clasificador de Bosque Aleatorio optimizado. Se generaron curvas ROC multiclase utilizando un enfoque de uno contraversus-estrategia de reposo. Se calcularon los valores del área bajo la curva (AUC) específicos para cada clase y el AUC promedio.
17. Comparación del rendimiento
Se comparó el rendimiento del modelo de características completo con el del modelo de características reducido obtenido tras la ablación de características por entropía.