Machine learning polysomnographically | Shandong Forgingss Group

Scientific Reports volume 13, Article number: 9120 (2023) Cite this article

Metrics details

Excessive daytime sleepiness (EDS) causes difficulty in concentrating and continuous fatigue during the day. In the clinical setting, the assessment and diagnosis of EDS rely mostly on subjective questionnaires and verbal reports, which compromises the reliability of clinical diagnosis and the ability to robustly discern candidacy for available therapies and track treatment response. In this study, we used a computational pipeline for the automated, rapid, high-throughput, and objective analysis of previously collected encephalography (EEG) data to identify surrogate biomarkers for EDS, thereby defining the quantitative EEG changes in individuals with high Epworth Sleepiness Scale (ESS) (n = 31), compared to a group of individuals with low ESS (n = 41) at the Cleveland Clinic. The epochs of EEG analyzed were extracted from a large overnight polysomnogram registry during the most proximate period of wakefulness. Signal processing of EEG showed significantly different EEG features in the low ESS group compared to high ESS, including enhanced power in the alpha and beta bands and attenuation in the delta and theta bands. Our machine learning (ML) algorithms trained on the binary classification of high vs. low ESS reached an accuracy of 80.2%, precision of 79.2%, recall of 73.8% and specificity of 85.3%. Moreover, we ruled out the effects of confounding clinical variables by evaluating the statistical contribution of these variables on our ML models. These results indicate that EEG data contain information in the form of rhythmic activity that could be leveraged for the quantitative assessment of EDS using ML.

Excessive daytime sleepiness (EDS) occurs when staying awake or alert is a progressive challenge for the individual. This condition is not only inherent to hypersomnia disorders such as narcolepsy, idiopathic hypersomnia and sleep disordered breathing, but also can be associated with a range of clinical factors, including metabolic and neurological diseases, ultimately translating to impairment of voluntary activities during the day or night1. EDS has become a significant public concern when associated with fatigue, costing more than $135 billion annually in health-related lost productivity in the United States2. In addition to the financial cost, the individual perception of the difficulty to concentrate and the deterioration of the brain response to audio, visual, and other stimulation motivates the search for a non-invasive biomarker that can help identify EDS to provide effective treatment. Seeking to find associations between sleepiness and its intertwined dynamics in the central nervous system (CNS), we tested the hypothesis that EEG data contain information in the form of rhythmic activity that could be leveraged for the quantitative assessment of EDS using machine learning (ML).

Day-time sleepiness affects the CNS and leads to changes in brain function and rhythms. Previous studies reported desynchronization between the left and right hemispheres under mental fatigue3, and imaging data suggest altered functional connectivity between thalamus and cortex4. In fact, EEG has shown promising results in identifying EDS biomarkers, especially for fatigue versus alert state classification during activities such as driving using portable EEG devices5 and for predicting driving reaction time6. In the clinic, subjective sleepiness symptoms can be assessed by the Epworth Sleepiness Scale (ESS), which is the state-of-practice self-report to quantify EDS or dozing propensity, and in general, it is highly correlated to the standard of care sleepiness measurements such as the multiple sleep latency test7. However, current diagnostic methods remain essentially subjective, because they rely on questionnaires and verbal reports.

In this study, we recorded resting state EEG from awake human subjects obtained from the Cleveland Clinic registry of overnight polysomnograms and pre-processed the data using an automated artifact detection algorithm that our team has previously developed8. Following a statistically-guided approach for EEG feature selection, we trained a ML algorithm to perform a binary classification of low versus high EDS. We further complement our study with a statistical analysis of the contribution of confounding clinical variables to our ML binary classifier.

We leveraged the Cleveland Clinic Sleep Registry, a collection of multimodal physiologic data, including continuous overnight EEG housing sleep studies. Data from this biophysiological repository was extracted for the purposes of this study with a focus on overnight polysomnograms or split night sleep studies. We abstracted polysomnogram data from a merged initiative of Cleveland Clinic's distributed sleep centers to capture those with severe hypersomnolence and those without symptoms of EDS. To serve as the analytic sample for the work, we identified from a total of 72 patients, 31 patients that had severe EDS, defined by an ESS greater or equal to 20 and 41 patients that had no EDS, defined by an ESS less than 5. This study (#22–135) was approved by Cleveland Clinic Institutional Review Board (IRB) Federal Wide Assurance (FWA 00005367) on 2/17/2022 as Exempt Human Subject Research, which granted a waiver of Informed Consent.

The polysomnogram studies were conducted in accordance with the American Academy of Sleep Medicine (AASM) guidelines using Polysmith software version 10 (Nihon Kohden). Signals were recorded using a standard 10–20 EEG montage. Only 6 referenced channels were selected for analysis, in some subjects, individual channels are referenced by ‘M1’ or ‘M2’. Odd channel numbers (located on the left side of the brain) are referenced by M1 and even channel numbers (located on the right side of the brain) are referenced by M2 in an ideal setting, rarely on cases of reference malfunction, the other side reference was used. The channels used are ‘F3’, ‘F4’, ‘C3’, ‘C4’, ‘O1’, and ‘O2’. The EEG was scanned and only the cleanest signal of 3 continuous minutes interval between "lights out" and the first stage of sleep "N1" was selected, this is the period we define as "awake" epochs in the context for this study. The cleanest signal is based on the number of artifacts. "Lights out" is a log event added by the technician and it means that troubleshooting and the bio calibration has been done, and that the sleep study is ready to be recorded. The awake epochs align with AASM guidelines, which requires that more than 50% of the epoch consists of alpha frequency activity.

Pre-processing of EEG data, feature extraction, statistics, and ML were performed using MATLAB (MathWorks). Since sleepiness is a feature of wakefulness, only EEGs during resting state wakefulness were analyzed, defined as the EEG recording time when subjects were awake prior to the sleep study.

The "awake" epochs selected were scored and chosen as wake manually and annotated electronically in the EEG file. Sleep studies conducted in the Cleveland Clinic Sleep Disorders Laboratory were manually scored for sleep staging including sleep and wake using American Academy of Sleep Medicine scoring rules. After this step, the algorithm developed was based upon EEG data extracted automatically from the wake epochs of the sleep studies. In this respect, we considered the EEG analysis as being automated. Greater than 50% of the epoch was comprised of alpha frequency activity without microsleeps. EEG was collected at a sampling rate of 200 Hz. A high-pass filter with a passband frequency of 1 Hz and a notch filter with a stop-band of 57.5–62.5 Hz were applied to all recordings. All EEG recordings were first visually inspected to confirm overall signal quality for each channel; channels considered to be of low or irretrievable quality were excluded from the study. Waveforms in each channel were further divided into 1-s epochs, and each epoch was tested for the presence of artifacts using a previously validated support vector machine (SVM)8. Epochs containing artifacts were excluded.

From the remaining artifact-free epochs of each recording, the following features were extracted: band-wise PSD for all channels and band-wise Phase-Amplitude Coupling (PAC). To create the band-wise PSD, a periodogram was gathered from artifact-free epochs for each channel, and then these periodograms were averaged together for each channel within each subject. These averaged periodograms were normalized by dividing each frequency bin by the sum of all bins from 3 to 30 Hz. The normalized PSD was used to calculate the band-wise PSD by taking the average of all bins within each of the following four frequency bands: Delta (1–4 Hz), Theta (5–9 Hz), Alpha (10–13 Hz), Beta (14–32 Hz), and low Gamma (33–52 Hz). This yielded 5 PSD features for every channel included.

The Band by Channel PSD, is similar to the band-Wise PSD. The only difference is that channels are not averaged; each channel is considered separately. This yielded a total of 30 features (5 bands × 6 channels) to be considered.

PAC was calculated using the Modulation Index (MI) method9. The center frequencies used for the phase included all the even numbers from 2 to 20. The center frequencies used for amplitude included all multiples of 3 from 30 to 54. MI was measured for each pair of phase and amplitude frequencies (90 total pairs) for each channel, including only those time points for which there were 5 or more consecutive artifact-free eyes-open epochs. This yielded a 9 × 10 MI matrix, for every channel of every subject, with each row corresponding to one phase center frequency and each column corresponding to one amplitude center frequency. This MI matrix was converted to band-wise PAC for the following 4 pairs of bands: Delta—low gamma, Theta—low Gamma, Alpha—low Gamma, and Beta—low Gamma. Other band-wise PAC were computed for the following 4 pairs of bands: Delta—Medium Gamma, Theta—Medium Gamma, Alpha—Medium Gamma, and Beta—Medium Gamma. This conversion was accomplished by averaging across the appropriate regions of the MI matrix. This yielded 4 PAC features for each channel and a maximum of 24 PAC features per subject (4 band-pairs × 6 channels).

Bandwise coherence was calculated from each of 15 unique channel pairs by averaging the coherence (MATLAB function mscohere) from each artifact-free epoch during the recordings for both channels in a given pair. This average coherence was divided into five bands in the same manner as PSDs. This yielded 5 coherence features for each channel pair, and a maximum of 75 coherence features per subject (5 bands × 15 channel pairs), coherence values were not computed from channel pairs for which one or both of the channels included artifacts.

We used paired two-tailed t tests to compare the band-wise Power Spectral Density (PSDs) between the high ESS and low ESS groups10. We used two-tailed Wilcoxon rank-sum tests to compare the band-wise phase-amplitude coupling (PAC) from the two groups for each of 4 band pairs and 6 channels. We chose a non-parametric statistical test for PAC because values are constrained between 0 and 1 and are therefore less likely to follow a normal distribution, as required by Student's t test. Statistical significance was established throughout at p value < 0.05. As individual testing was conducted, and the purpose of the statistical testing was primarily the selection of features for subsequent ML (rather than testing of a null hypothesis), there was no adjustment for multiple comparisons11,12. Unlike simultaneous family testing with a joint null hypothesis comprising two or more null hypotheses, individual testing is utilized to make a decision about one null hypothesis. As each test provides only one opportunity to make a Type I error, the alpha level does not require lowering.

Although it is common in the literature to include confounding factors as input in a regression model for further adjustment, this approach is inadequate in ML classifiers that can learn complicated non-linear relationships between the input and output. Therefore, we performed the post-hoc analysis described in13, which proposes controlling for confounds by using traditional regression to compare the extent to which a trained ML prediction itself can explain the target variable, in contrast to the independent performance of confounding variables.

Our study population was overall middle-aged (mean age 54 years) with a relatively even distribution of men and women, race-based diversity (34.7% African American) and a mild degree of sleep-disordered breathing (apnea–hypopnea index = 13.4). Those with a higher degree of EDS were more likely to be slightly younger, female and more obese with longer total sleep time and a lower percentage of N1 sleep stage. Of note, as reflective of our pre-specified design, those with a high degree of EDS had a mean ESS of 21 ± 1, in contrast, those without EDS had a mean ESS of 2 ± 1. Detailed demographic information, sleep characteristics, and medical history for overall patients and their subsequent division into high ESS and low ESS groups are shown in Table 1.

Out of the above listed EEG features in Methods, 3 features were selected to create a feature-set for training binary classification algorithms shown to be significantly different between groups; these included 2 from bandwise coherence generated by channels pairs O2 – C4 for both Delta and Theta bands, and channel O2 for the Beta band in PSD.

Following recent trends in best ML practices14, several traditional ML classification algorithms were considered, using cross-validation and a grid-search strategy to find their optimal hyperparameters. The best results were obtained from k-nearest neighbors algorithm (KNN). To ensure data splits were properly distributed between folds, we validated the classifier using a stratified K cross-validation using k = 5. Within our dataset with n = 31 for the high ESS group and n = 41 for the low ESS group, we computed the following metrics for validation: accuracy, area under the receiver operating characteristic curve (AUC-ROC), precision, specificity and recall. Accuracy was calculated within the k-folds cross validation by counting the number of out-of-sample predicted labels that matched the true label of the sample, and dividing this total by the number of samples.

The spectral density for delta, theta, alpha, and beta power bands showed a statistically significant difference in the mean EEG power in high ESS compared to the low ESS group (Fig. 1), including.an increase in delta (0.294 ± 0.010 in high ESS, 0.239 ± 0.005 in low ESS, p < 0.001), in theta (0.208 ± 0.005 in high ESS, 0.181 ± 0.004 in low ESS, p < 0.001) and a decrease in alpha (0.260 ± 0.011 in high ESS, 0.318 ± 0.007 in low ESS, p < 0.001) and beta (0.087 ± 0.003 in high ESS, 0.112 ± 0.007 in low ESS, p < 0.001). Analysis of power in individual 6 channels further showed that significant changes were not localized to particular brain areas. Significant differences in delta, alpha, and beta power bands in each EEG channel were found between low and high ESS groups (Fig. 2).

Power Spectral Density (mean of 6 EEG channels) in high ESS (n=31) and low ESS (n=41) subjects in the 0-50 Hz frequency range (upper row). Mean power in the frequency bands delta (1-4 Hz), theta (5-9 Hz), alpha (10-13 Hz), beta (14-32 Hz) and low gamma (33-52 Hz) (lower row).

Histograms show Power Spectral Density in 6 individual EEG channels in high ESS (n=31) and low ESS (n=41) subjects in the frequency bands delta (1-4 Hz), theta (5-9 Hz), alpha (10-13 Hz), beta (14-32 Hz) and low gamma (33-52 Hz). Heat map shows corresponding t-test p-values for individual channels in each band (red hue indicates p < 0.05, highlighted cell indicates a feature selected for the training of a ML algorithm).

For Phase-amplitude coupling (PAC) in 6 individual EEG channels, there was no significant difference between low gamma and delta, theta, alpha, beta respectively, as well as between medium gamma and delta, theta, alpha, beta respectively (Fig. 3).

Phase-amplitude coupling (PAC) in 6 individual EEG channels in high ESS (n=31) and low ESS (n=41) subjects between low gamma and delta, theta, alpha, beta respectively (upper row), as well as between medium gamma and delta, theta, alpha, beta respectively (lower row). No statistically significant difference was noted between groups in any individual channel.

In Fig. 4, we compared the coherence in 6 individual EEG channels across each frequency band. Significant difference was found for the high ESS and low ESS groups and between coherence in some of the individual channels (C4 and O2).

Coherence in 6 individual EEG channels in high ESS (n=31) and low ESS (n=41) subjects. The values are symmetrical across one diagonal (the 2 highlighted cells indicate features that were selected for the training of a ML algorithm).

We computed PSD by Hz (Fig. 5), which is similar to Band by Channel (Fig. 2 bottom right), except that instead of bands, we are using bins from 1 to 16 Hz. This yielded to 96 features (16 Hz × 6 channels) to be considered. In Fig. 2 (Heat map (bottom right)), we compared average PSD in 6 individual EEG channels across each frequency bin. Significant difference was found between the high ESS and low ESS groups in channel O2 in Beta band.

Heat map shows t-test p-values for individual channels, in high ESS (n=31) and low ESS (n=41) subjects, in each frequency bin from 1 Hz to 16 Hz. (red hue indicates p<0.05).

The 3 features generated from the PSD and Coherence analysis were used to train k-NN binary classifier of high versus low ESS. Our model reached an accuracy of 80.2%, an AUC-ROC of 79.5%, a precision of 79.2% a recall of 73.8%, and specificity of 85.3%.

In our analysis of confounding variables that could influence our EEG findings, we compared the predictions of our ML classifier with the potential confounding variables: age, BMI, and gender. After training our model, in each fold of our cross-validation, we used a generalized linear model to predict EDS on the test set in 3 scenarios: Using only our ML output, using only the confounding variables, using both the confounding variables and the ML output. To aggregate the results for the 5 folds, we computed the mean Pseudo-R2 and the mean p-values. In Fig. 6, we show the mean Pseudo-R2, whereas our ML only reached 41% and the confounding only reached 8%. When combining both ML and confounding variables reached 56%. Table 2 shows the mean p values from the generalized linear regression for the same 3 scenarios for each variable. We observed that we have a statistically significant predictor when only the ML output is used (p value < 0.05).

Mean Pseudo R² computed by a generalized linear model on the test sets of cross-validation for 3 scenarios: using only the machine learning prediction, using only the confounding variables age, BMI and gender and using the combination of both machine learning and confounding.

Clinically, EDS overlaps with other common sleep and mood disorders that are considered causal to EDS, such as obstructive sleep apnea (OSA), narcolepsy, depression, and post acute sequelae of COVID-1915, with an estimated prevalence of 20% of adults in the United States16. These causes are usually misdiagnosed or undiagnosed in the population, even when a patient performs a sleep study7,17.

Several approaches have been investigated to develop an objective neurophysiologic biomarker capable of capturing symptoms of EDS. For example, the least absolute shrinkage and selection operator (LASSO) was used to predict ESS from EEG signals collected from train drivers, but with varying degrees of success and requiring more complex computational techniques compared to our study, which was guided by a statistical approach for the selection of ML features18. Another proposed sleepiness biomarker is the odds ratio product, computed from the delta, theta, alpha-sigma, and beta frequency bands from EEG signals, and its association to ESS19. Despite previous work, however, the direct classification of high or low ESS using EEG lacks transparency in the ML approach.

In this study, we demonstrate the feasibility of an automated analytical pipeline, using resting state EEG during wakefulness and ML, of accurately classifying EDS as low versus high. Our results showed that the average power spectrum across EEG channels in low ESS patients is significantly enhanced in the alpha and beta bands and attenuated in the delta and theta bands compared to high ESS subjects.

After analyzing sleep EEG in two groups of patients, we identified significant differences in the delta, theta, alpha, and beta frequency bands between those who reported high versus low sleepiness. Aiming to estimate daytime sleepiness at the individual level, we trained a k-NN classifier reaching 80.2% accuracy, 79.2% precision, 73.8% recall, and 85.3% specificity. Thus, our work demonstrates a potential for EEG analysis to generate biomarkers for excessive daytime sleepiness and the use of signal processing and ML to classify sleepiness at the individual level. Moreover, we emphasized the novelty of these biomarkers when we showed that even after controlling for confounding variables, ML predictors alone significantly explain EDS.

These results suggest that short duration, resting state, wake EEG contains information that could be leveraged to reliably assess EDS, thereby enhancing clinical care. These findings should be contextualized with the historical standard of utilizing EEG patterns to qualitatively and to a lesser extent quantitatively define sleep stages as per the American Academy of Sleep Medicine and assessment of symptoms of excessive daytime sleepiness with time-intensive objective tests such as the multiple sleep latency test20.

Our results show that low ESS subjects have higher EEG activity in the beta band in the occipital region than high ESS subjects. Previous research has also highlighted the significant differences between these biomarkers in groups of EDS patients21. Others have suggested that neural mechanisms related to visual tasks and attention may be linked to these biomarkers. For instance, an increase in response time for visual tasks among students in a real classroom is associated with a decrease in the mean baseline beta band power over the occipital region22. Also, in23, beta band activity in the occipital region was linked to attention deficits in elderly subjects, with increased activity preceding correct responses for visual tasks.

Some limitations of our work deserve attention. Since the sample size could be considered relatively small (n = 72) in the field of EEG/ML14, our model would benefit from validation against a prospective dataset as per ML best practices, as well as across different geographical sites. Our approach still depends on manual sleep staging annotation since the automated pipeline does not determine "wake" epochs. Still, the "wake" period is determined after we scan the data looking for the manual annotations of "lights out," "wake," and sleep stage "N1". Moreover, according to AASM, the awake epochs can range from full alertness through early stage of drowsiness. We also used signal processing techniques to extract EEG features, for example power in predefined delta, theta, alpha, beta, and gamma bands. This may have limited our feature space for biomarkers in contrast to other techniques that are based on training deep learning models directly from EEG raw signals24, although such techniques are non-transparent and more computationally demanding. We also acknowledge that some demographic variables might have an impact on EEG, such as age25,26,27,28 and gender29,30, and we did not specifically control for such parameters, although mean demographic values were overall comparable.

We chose extremes of subjective perception of dozing propensity defined by ESS as an exploratory use case for algorithm development. The sleep studies were not selected based upon obstructive sleep apnea status, but rather patients were most often referred for the sleep study based upon suspected OSA with varying levels of self-reported sleepiness. Therefore, there are opportunities to build upon this work in the future by examining EEG-based biomarkers which reflect the spectrum of sleepiness.

Finally, we envision that our automated and quantitative method for assessing EDS can be operationalized by adding our method's output into patients' Electronic Health Record after performing a polysomnogram test. Future investigation should focus on whether a point of care measure such as this could be used to risk stratify patients to identify those for example who may need drowsy driving education or to identify sleep phenotypes more responsive to specific interventions such as pharmacotherapeutics to treat hypersomnia disorders.

In this study, we investigated the potential of using EEG signals recorded during the awake period in polysomnograms as biomarkers for excessive daytime sleepiness. We identified significant differences in the delta, theta, alpha, and beta frequency bands between those who reported high versus low sleepiness measured by the widely used ESS. Aiming to estimate daytime sleepiness at the individual level, we trained a k-NN ML classifier that reached 80.2% accuracy, 79.2% precision, 73.8% recall, and 85.3% specificity in a retrospective cohort. After controlling for potential confounding variables, we show that this study innovates in building a direct association between EEG and ESS. Ultimately, we provide powerful techniques for those areas that can leverage advanced EEG analysis, such as underdiagnosed EDS studies.

The data that support the findings of this study and the datasets generated and analyzed are not publicly available due to human subject concerns, but are available from the corresponding author on reasonable request.

Jaussent, I. et al. Incidence, worsening and risk factors of daytime sleepiness in a population-based 5-year longitudinal study. Sci. Rep. 7(1), 1372 (2017).

Article ADS CAS PubMed PubMed Central Google Scholar

Ricci, J. A. et al. Fatigue in the US workforce: Prevalence and implications for lost productive work time. J. Occup. Environ. Med. 49(1), 1–10 (2007).

Article MathSciNet PubMed Google Scholar

Jarchi, D., B. Makkiabadi, & S. Sanei. Mental fatigue analysis by measuring synchronization of brain rhythms incorporating enhanced empirical mode decomposition. In 2010 2nd International Workshop on Cognitive Information Processing, pp. 423–427 (2010).

Killgore, W. D. S. et al. Daytime sleepiness is associated with altered resting thalamocortical connectivity. NeuroReport 26(13), 779–784 (2015).

Article PubMed Google Scholar

Ren, Z. et al. EEG-based driving fatigue detection using a two-level learning hierarchy radial basis function. Front. Neurorobot. 15, 3–3 (2021).

Article Google Scholar

Ko, L.-W. et al. Eyeblink recognition improves fatigue prediction from single-channel forehead EEG in a realistic sustained attention task. J. Neural Eng. 17(3), 036015–036015 (2020).

Article ADS PubMed Google Scholar

Aurora, R. N. et al. Correlating subjective and objective sleepiness: Revisiting the association using survival analysis. Sleep 34(12), 1707–1714 (2011).

Article PubMed PubMed Central Google Scholar

Levitt, J. et al. Automated detection of electroencephalography artifacts in human, rodent and canine subjects using machine learning. J. Neurosci. Methods 307, 53–59 (2018).

Article PubMed Google Scholar

Tort, A. B. L. et al. Measuring phase-amplitude coupling between neuronal oscillations of different frequencies. J. Neurophysiol. 104(2), 1195–1210 (2010).

Article PubMed PubMed Central Google Scholar

Niven, D. J. et al. Matched case-control studies: A review of reported statistical methodology. Clin. Epidemiol. 4, 99 (2012).

Article PubMed PubMed Central Google Scholar

Greenland, S. Analysis goals, error-cost sensitivity, and analysis hacking: Essential considerations in hypothesis testing and multiple comparisons. Paediatr. Perinat. Epidemiol. 35(1), 8–23 (2021).

Article PubMed Google Scholar

Rubin, M. When to adjust alpha during multiple testing: A consideration of disjunction, conjunction, and individual testing. Synthese 199, 10969–11000 (2021).

Article Google Scholar

Dinga, R., Schmaal, L., Penninx, B. W., Veltman, D. J., & Marquand, A. F. (2020). Controlling for effects of confounding variables on machine learning predictions. BioRxiv (2020).

Steyerberg, E. W. Validation in prediction research: The waste by data splitting. J. Clin. Epidemiol. 103, 131–133 (2018).

Article PubMed Google Scholar

Pérez-Carbonell, L. et al. Impact of the novel coronavirus (COVID-19) pandemic on sleep. J. Thorac. Dis. 12, 163–175 (2020).

Article Google Scholar

Pagel, J. F. Excessive daytime sleepiness. Am. Fam. Phys. 79(5), 391–396 (2009).

CAS Google Scholar

Sahni, A. S. et al. Management of excessive sleepiness in patients with narcolepsy and OSA: Current challenges and future prospects. Nat. Sci. Sleep 11, 241–252 (2019).

Article ADS PubMed PubMed Central Google Scholar

Lees, T. et al. Electroencephalography as a predictor of self-report fatigue/sleepiness during monotonous driving in train drivers. Physiol. Meas. 39(10), 105012 (2018).

Article PubMed Google Scholar

Lechat, B. et al. A novel EEG marker predicts perceived sleepiness and poor sleep quality. Sleep 45(5), 051 (2022).

Article Google Scholar

Academy of sleep medicine A. The AASM manual for the scoring of sleep and associated events summary of updates in version 2.5. J. Clin. Sleep Med. 16(4), 605–607 (2020).

Melia, U. et al. Mutual information measures applied to EEG signals for sleepiness characterization. Med Eng Phys 37(3), 297–308 (2015).

Article PubMed Google Scholar

Ko, L. W. et al. Sustained attention in real classroom settings: An EEG study. Front Hum Neurosci 11, 388 (2017).

Article PubMed PubMed Central Google Scholar

Gola, M. et al. EEG beta band activity is related to attention and attentional deficits in the visual performance of elderly subjects. Int. J. Psychophysiol. 89(3), 334–341 (2013).

Article PubMed Google Scholar

Nygate, Y. et al. 543 EEG-based deep neural network model for brain age prediction and its association with patient health conditions. Sleep 44(Suppl 2), 214 (2021).

Article Google Scholar

Hughes, J. R. & Cayaffa, J. J. The EEG in patients at different ages without organic cerebral disease. Electroencephalogr. Clin. Neurophysiol. 42(6), 776–774 (1977).

Article CAS PubMed Google Scholar

Nguyen, P. et al. (EDS) Age and gender classification using EEG paralinguistic features. In International IEEE/EMBS Conf. on Neural Engineering, NER (2013).

Vlahou, E. L. et al. Resting-state slow wave power, healthy aging and cognitive performance. Sci. Rep. 4, 5101 (2014).

Article CAS PubMed PubMed Central Google Scholar

Zoubi, O. A. et al. Predicting age from brain EEG signals-a machine learning approach. Front. Aging Neurosci. 10, 184 (2018).

Article PubMed PubMed Central Google Scholar

Hu, J. An approach to EEG-based gender recognition using entropy measurement methods. Knowl. Syst. 140, 134–141 (2018).

Article Google Scholar

Van Putten, M. J. A. M., Olbrich, S. & Arns, M. Predicting sex from brain rhythms with deep learning. Sci. Rep. 8(1), 3069 (2018).

Article ADS PubMed PubMed Central Google Scholar

Download references

The pre-processing automation stage was performed on Oracle Cloud platform.

These authors contributed equally: Araujo Matheus and Ghosn Samer.

Sleep Disorders Center, Neurological Institute, Cleveland Clinic Foundation, Cleveland, OH, USA

Araujo Matheus, Wells Samantha & Mehra Reena

Department of Biomedical Engineering, Cleveland Clinic Foundation, Cleveland, OH, USA

Ghosn Samer & Saab Y. Carl

Quantitative Health Sciences, Lerner Research Institute, Cleveland Clinic Foundation, Cleveland, OH, USA

Wang Lu

Metro Health Medical Center, Cleveland, OH, USA

Hariadi Nengah

Department of Biomedical Engineering, Brown University, Providence, RI, USA

Saab Y. Carl

Respiratory Institute, Cleveland Clinic Foundation, Cleveland, OH, USA

Mehra Reena

Cardiovascular and Metabolic Sciences, Lerner Research Institute, Cleveland Clinic Foundation, Cleveland, OH, USA

Mehra Reena

Heart and Vascular Institute, Cleveland Clinic Foundation, Cleveland, OH, USA

Mehra Reena

You can also search for this author in PubMed Google Scholar

M.A., S.G., R.M. and C.S. wrote the main manuscript text. S.W. and N.H. identified subjects and collected data, L.W. extracted subjects data and computed Statistics in Table. 1. M.A. prepared the Machine Learning confounding study, added Table 2 and Fig. 6. S.G. developed the EEG data pipeline, prepared data analysis, machine learning and Figs. 1, 2, 3, 4 and 5.

Correspondence to Mehra Reena.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Matheus, A., Samer, G., Lu, W. et al. Machine learning polysomnographically-derived electroencephalography biomarkers predictive of epworth sleepiness scale. Sci Rep 13, 9120 (2023). https://doi.org/10.1038/s41598-023-34716-5

Download citation

Received: 05 October 2022

Accepted: 05 May 2023

Published: 05 June 2023

DOI: https://doi.org/10.1038/s41598-023-34716-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.