Subjective sleep quality, quantitative sleep features, and their associations dependent on demographic characteristics, habitual sleep–wake patterns, and distinction of weekdays/weekends

Generally applicable qualifications of sleep and principles for achieving better sleep are difficult to design, because sleep quality can depend on individual demographic characteristics and lifestyles. In this study, the static and dynamic features of sleep–wake patterns were analyzed in association with quantitative sleep-related parameters and self-rated sleep quality to serve as a practical selection of sleep–wake patterns fitted to individual conditions. Data obtained over a 2-week period by actigraphy from university students and information technology workers were measured to obtain a daily subjective rating of sleep quality using the Oguri–Shirakawa–Azumi (OSA) sleep inventory. Qualitative sleep quality in terms of OSA score and quantitative sleep-related and chronobiological features were analyzed with regard to their dependency on the demographic characteristics, habitual sleep–wake patterns (HSWP), and distinction of weekdays/weekends. Multi-factor ANOVA was used to further investigate their dependencies regarding multiple ways of interactions between the demographic characteristics, HSWP, and distinction of weekdays/weekends. Subjective sleep quality and quantitative sleep-related parameters depended on the demographic characteristics, and so did their associations. The classification of day-to-day variations in HSWP showed four clusters that were effective factors for understanding their dependencies. Multi-factor analysis revealed demographic characteristics, HSWP, distinction of weekdays/weekends, and their multi-way interactions up to 3rd order as significant effectors of qualitative and quantitative quality of sleep. This study clarified how quantitative sleep-related parameters, subjective sleep quality, and their associations depended on demographic characteristics. Furthermore, their dependency was understood as a combination of multi-way interactions between the demographic characteristics, HSWP, and the distinction of weekdays/weekends. Our findings could provide a basis for the design of individually matched sleep–wake patterns.


Introduction
Sleep plays an important role in the maintenance of human mental and physical conditions. Poor sleep quality can cause wide range health problems, both acute (e.g., mood swings, slower cognitive response, higher stress levels) and chronic (e.g., obesity, increased risk of cardiometabolic diseases) [1][2][3]. How to achieve good sleep quality has been attracting increasing attention, and numerous theories and methods have been proposed [4]. Among these, wearable devices can inform users of their estimated sleep quality based on daily activity and/or physiological signals such as heart rate, blood pressure, and body temperature [5,6]. However, people used to select their sleep-wake patterns suitable for individual conditions according to one's own empirical ideas. Of 1 3 course, good sleep should take place with appropriate timing and length in a proper environment, with no disruptions [7]. However, because demographic characteristics, working style, individual preferences, and physiological responsiveness in terms of sleep regulation vary from person to person, generally applicable sleep qualifications and principles for achieving better sleep are difficult to design.
Chronotype and social jet lag (SJL) have recently become widely recognized as chronobiological parameters that characterize habitual sleep-wake patterns [2]. People with the early chronotype prefer earlier sleep phases, and vice versa [8][9][10][11]. Naturally, those with the late chronotype suffer from insufficient sleep, and show a tendency to accumulate sleep debt over working days that results in prolonged sleep on the weekend or non-working days, i.e., severe SJL [11]. Several studies have suggested that the late chronotype and severe SJL influence physiological phenomena such as later melatonin secretion [9], later circadian acrophase of activity [10], and lower sleep regularity [12,13], which could underlie the deterioration of sleep quality. Therefore, these parameters may be possible candidates of factors that relate habitual sleep-wake patterns to the quality of sleep. As sleep is dynamically regulated by interactions between circadian rhythms and homeostasis under the influence of environmental and social time cues [14], it might be useful to employ more dynamical features of habitual sleep-wake patterns, rather than static parameters such as chronotype and SJL, in designing sleep-wake patterns for maintaining/improving sleep quality. In fact, the dynamical feature of regularity in the daily sequence of resting states has been suggested as a relevant index in the evaluation of health conditions [3,15].
We analyzed the association of quantitative feature parameters of activity during sleep, static chronobiological parameters such as chronotype and SJL, and self-rated sleep quality with demographic characteristics and day-today variations in habitual sleep-wake patterns. The identified associations were dissected into multiple interactions between these parameters and factors. The present results could contribute to the practical selection of sleep-wake patterns fitted to individual conditions, including demographic characteristics, individual preferences, working style, and chronobiological/homeostatic regulation.

Materials and methods
The data were acquired in 2014-2017, partly from a previous study conducted under a different context [16]. The 356 participants with no severe sleep complaints were recruited in the study, including information technology workers (N = 180) or university students (N = 176) aged 20-63 years. All participants were asked to record their physical activity (sampling rate: 30 Hz) using a wrist-mounted accelerometer (Actigraph® GT3X + ; Actigraph LLC, FL, USA) and complete the Oguri-Shirakawa-Azumi (OSA) sleep questionnaire [17] upon waking, for 15 days. The OSA sleep questionnaire quantifies subjective sleep quality into five categories: sleepiness at awakening (OSA-1), sleep maintenance (OSA-2), dreaming (OSA-3), recovery from fatigue (OSA-4), and self-reported sleep duration (OSA-5). Higher OSA scores reflect better sleep quality. All measurement procedures were conducted under the approval of the Ethical Committee of the Graduate School of Information Sciences, Tohoku University, Japan.

Estimation of quantitative sleep-related parameters
Sleep-wake timing and parameters were estimated using ActiLife® (Actigraph LLC, FL, USA) based on the Cole-Kripke sleep-wake identification algorithm [18]. A sleep report using the sensor algorithm comprised wake after sleep onset (WASO), sleep-onset latency (SOL), number of awakenings (NAWK), and sleep efficiency (SE) [19]. Midsleep time (MS), chronotype (MSF sc ; mid-sleep on free days corrected for sleep debt on workdays), and SJL were also calculated as described in [2,8,20]. A simplified version of the Sleep Regularity Index (SRI) was used to estimate a mean coincidence percentage of sleep periods between 2 consecutive days [12].
To extract the onset of posture change during sleep, the data were multiband-filtered (finite impulse response [FIR] filter of order 20; pass band, 0-0.05 and 2-5 Hz) and then segmented into 2-s bins with 50% overlap. The mean and variance of the filtered signals were calculated for each bin. Movement in each bin at time t was detected if the statistical quantities for at least one of three-dimensional signals satisfied either of the following conditions: (1) the absolute difference between the successive mean levels at t and t + 1 was larger than the threshold (= 0.01), or (2) the variance at t was larger than the threshold (= 0.007). Movements lasting for less than 4 s and intermissions of less than 3 s were ignored, thus excluding jitter-like movements and subjecting only major posture changes for analysis. Figure S-1a in Online Resource 1 shows the method used to estimate inter-movement-onset intervals (IMOIs). IMOIs were subjected to logarithmic transformation because of their longtailed distribution, and then, their mean, T mean , and standard deviation (SD), T std , were calculated for each sleep episode. In addition, the number of movements, N m , and movement density, D m (ratio of number of movements to sleep length), were obtained for each sleep episode.
The circadian amplitude and acrophase of daily activity were estimated based on one-component static cosine fitting with 24-h period. Daily activity sequences were band-pass filtered (FIR filter of order 30; pass band: 2-5 Hz), and the locomotor inactivity during sleep (LIDS) ultradian rhythm 1 3 was then extracted following a previously reported method [21]. The parameters obtained were ultradian period, LIDS slope, range of oscillation (RoO), and number of ultradian cycles (N c ). Figure S-1b shows an example of an ultradian rhythm fitted to the LIDS.

Statistical analysis and classification of habitual sleep-wake patterns (HSWPs)
All statistical analyses were performed using MATLAB (R2017a , where the female workers were not subjected to analysis for stratification reasons. Individual OSA scores were standardized (z-score) before estimating the Spearman correlations separately for gender distinction in the students and age group distinction in the male workers. The male workers were grouped by age into three groups (20 ≤ age < 35, 35 ≤ age < 45, and 45 ≤ age ≤ 63), considering the age distribution of workers in which there was a gap around 35 ( Fig.S-3 of Online Resource 2). This grouping did not so deviate from commonly used every ~ 10 years grouping [22]. The same age grouping was applied in another analysis ( ) who had completed ≥ 14 consecutive days of activity recordings, out of which 14 consecutive days were subjected to the following clustering. K-means clustering was then performed using 10-dimensional features, including the mean and SD of sleep-and wake-onset times, second and third harmonics of the sleepand wake-onset time daily sequences, average MS, and MSF sc . Then, principal component analysis (PCA) was applied to show the clustering result explicitly in a lower dimensional space. Among the variables used for clustering, the former 8 variables were selected as those for effectively characterizing HSWP. The latter two variables were selected out of chronotype (MSFsc), variation of mid-sleep time, SRI, and variation of sleep duration, because including these two variables yielded the best result in terms of percentage of variance explained in PCA. Appropriate numbers of clusters and cluster features were confirmed by hierarchical clustering and the elbow method [23].
To summarize the characterization of subjective sleep quality and quantitative sleep-related features, their multiple dependency on demographic characteristics and HSWP as between-subject factors, and distinction of weekdays/ weekends as a within-subject factor were investigated by mixed factor ANOVA [24], where the data of individual subjects were averaged separately for weekdays or weekends, and a repeat measure model was then generated. The model included both main and interaction effects among factors. In addition, the post hoc analysis was done by Tukey-Kramer test.  Table 1. Characteristically, demographics dependency of quantitative sleep-related parameters was more pronounced than that of subjective sleep quality. Distinctively significant difference in subjective sleep quality was confined to dreaming frequency (OSA-3) between male students and workers, and sleepiness at awakening between female students and workers. Irrespective of gender, students tended to take a later, longer, more irregular, and less disturbed sleep than the workers. Gender difference was observed in a job-specific way, i.e., male students tended to take a later and more disturbed sleep than female students. In contrast, the sleep of male workers tended to suffer more dense disturbances than that of female workers. Because the ages of male workers were distributed over a wide range, they were separated into three generational groups and age dependency was analyzed (Table S-1 of Online Resource 2). Sleep timings tended to be earlier in the eldest group than in other groups, which was also suggested by the chronobiological features. The flat ultradian slope also characterized the sleep of the eldest group. The sleep of the middle-aged group tended to be less disturbed in terms of NAWK and N m , accompanying subjective sleep ratings of highest OSA-2 and OSA-3, and lowest OSA-1 and OSA-4 in comparison with the other groups. The youngest group tended to take a longer and later sleep, which coincided with later chronobiological features, although not significantly. However, sleep in this group was characterized quantitatively by frequent disturbances in terms of NAWK and N m , accompanying subjective sleep ratings of lower OSA-2, 3, and 4.

Demographics dependency of quantitative sleep-related parameters and subjective sleep quality, and their associations
In correlation analyses performed between daily quantitative sleep-related parameters and subjective sleep quality for each demographic group (Tables 2, 3, 4, 5) (Spearman correlations), the student analysis focused on gender difference Table 1 Comparisons with regard to gender difference (male/female) and job difference (students/workers) Values are represented as mean (SD, Gaussian)/median [semi-interquartile range, non-Gaussian]. As post hoc tests, t test, and Mann-Whitney U (MWU) tests were used for Gaussian and non-Gaussian cases, respectively. t-value and degree of freedom, and Z-value are given for t test and MWU test, respectively *** p < 0.001, **p < 0.01, *p < 0.05. MS: mid-sleep time, MSF SC mid-sleep on free days corrected for sleep debt on workdays, SJL social jet lag, SRI sleep regulatory index, SE sleep efficiency, SOL sleep-onset latency, NAWK number of awakenings, WASO wake after sleep onset, N m number of movements, D m movement density, LIDS locomotor inactivity during sleep, RoO range of oscillation in LIDS, N c number of cycles in LIDS. Circadian amplitude and Acrophase are estimated by cosine fitting to the daily activity. These abbreviations and notations are applicable for the following tables  (Tables 2 and 3) and the worker analysis focused on generational difference within the males because of their widely distributed ages and because of the limited number of female worker subjects (Tables 4 and 5; for the middle generation group of workers (35 ≤ age < 45), see Table S-2 of Online Resource 2). The following tendencies were shared by the results for the students. (1) Sleep duration had non-negligible correlations with most of the OSA scores, especially OSA-1 (sleepiness at awakening), 4 (fatigue recovery), and 5 (subjective sleep duration). (2) Parameters characterizing the total disturbances during sleep, such as NAWK, WASO, and N m , exhibited similar profiles to those of sleep duration. (3) Associations with the parameters characterizing the occurrence density of disturbance such as SE and D m were confined to OSA-2 (sleep maintenance), as did the statistics of IMOIs, T mean and T std in the female students. (4) Chronobiological and ultradian rhythm-related parameters such as amplitude and acrophase of circadian rhythm, period of ultradian rhythm, LIDS slope, and RoO were mostly uncorrelated with the  OSA scores, although N c shared a profile with the total disturbance-related parameters. The correlation profile of female students differed from that of males in terms of nonnegligible significant correlations between sleep duration and OSA-3 (dreaming) and between statistics of IMOIs and OSA-2 (sleep maintenance). In addition, the total disturbance-related parameters showed slight correlation with OSA-4 only in female students, and with OSA-1 only in male students. Compared with the students, the male workers showed non-negligible correlations of the parameters related to the total disturbances to OSA-4 as well as OSA-1. Slight but significant OSA-2-specific correlations of the disturbance density-related parameters, as a tendency shared by the students, were found only for the eldest worker group (> 45 years), as did the statistics of IMOIs. In addition, ultradian period had non-negligible correlations with OSA-1 and OSA-4 for the middle-aged group of workers.
Comparison between the students and youngest worker group revealed that sleep duration, the total disturbancesrelated parameters, and N c correlated with more OSA scores in the workers than in the students. In contrast, slight but significant OSA-2-specific correlations of the disturbance density-related parameters were found only in the students.

Classification of HSWPs
Further analyses were performed to clarify how the subjective and quantitative aspects of sleep depend on HSWPs through clustering. The results of HSWP clustering are shown in Fig. 1, which shows day-to-day variations in sleepand wake-onset times superposed around the individual mean level, together with the typical HSWP, for each cluster. Together with the cluster-dependent statistical features of the parameters used for classification and the demographic  (Tables S-3 and S-4 of Online Resource 3), the HSWP clusters were characterized as follows: (1) Cluster 1 (C1, N = 83): The statistics for sleep-and wake-onset times, MS, and MSF sc showed that C1 shared the tendency of earlier habitual sleep with C3. Nevertheless, in contrast with C3, more severe SJL and greater irregularity of daily activity patterns in terms of SRI characterize this cluster. C1 or C3 included more workers than the other clusters, which was statistically suggested by Chi-square test of bias in demographic proportions across clusters. (2) Cluster 2 (C2, N = 34): C2 shared the tendency of later habitual sleep with C4. Extreme day-to-day and weekly variations in sleep-and wake-onset times distinguished C2 from C4, as did the irregularity of daily activity patterns. In C2 and C4, there were more students than workers. This bias of student proportions in C2 and C4 was suggested statistically. PCA conducted to represent the cluster features with a minimal dimension revealed the distribution of day-to-day variations in HSWPs in the third-order principal component (PC) coordinates (the total contribution ratio up to the third-order PC components was 74%) together with the loading vectors, as shown in Fig. 1c-e. From mappings on the PC1-PC2 and PC1-PC3 planes, it was found that PC1 separated C3 and C4, and PC2 separated C1 and C2. From mapping on the PC2-PC3 plane, PC2 separated C1/C4 and C2/C3. The PC1-PC3 vectors are listed in Table S

Multiple dependency of subjective sleep quality and quantitative sleep-related parameters on demographic and HSWP characteristics
The dependency of subjective sleep quality and quantitative sleep-related parameters on demographic characteristics shown in Table 1 was further analyzed by taking into consideration the HSWP clustering results and weekday/ weekend difference. Their multiple dependency on demographics (Gender and Job), Day (distinction of weekdays/ weekends), and HSWP clusters is shown up to 4-way interaction ( Table 6). The significant effectors for the respective features are summarized graphically in Fig. 2. Subsequently, a post hoc analysis by Tukey-Kramer test was performed (Table S-6 of Online Resource 4). Based on these results, the sleep features significantly affected by the multi-way interactions are listed below, and their dependencies are explained, where the highest order of interactions is given in parenthesis with *p < 0.05, **p < 0.01, or ***p < 0.001, and the abbreviations, G for Gender, J for Job, D for Day, and C for HSWP, are used.
OSA-4 and OSA-5 (D*): Subjective evaluation of fatigue recovery and sleep duration tended to be better on weekends in comparison with weekdays.
Sleep and wake onsets, and mid-sleep time (C:D***, J:D*): Sleep phase-related timings were found to exhibit the inter-cluster dependency similar to Table S-3, irrespective of weekdays or weekends. In addition, these timings tended to be significantly later on weekends than weekdays for almost all HSWP clusters. The students tended to wake up significantly later as compared with the workers, irrespective of weekdays or weekends.
Sleep duration (C:G:D*): Significant prolongation of sleep duration tended to take place on weekends for some combinations of cluster (C2, C3, C4) and gender. For the females, weekend prolongation of sleep duration was pronounced for C2 in comparison with C4.
Circadian amplitude and acrophase (G:D*, C:D*): For the males, circadian amplitude tended to be larger on weekends than weekdays, and this weekend increase in circadian amplitude was shown to be more pronounced for the males. Weekday circadian acrophase tended to be earlier for C3 than C2. SE, SOL, NAWK, and D m (J:D**, J:D*, C:G:D**, D*): The workers tended to experience lower weekend SE and longer weekend SOL, as compared with the students. D m tended to be higher on weekends. The male subjects belonging to C3 and C4 tended to experience more NAWK on weekends, as compared with the female workers.
Ultradian period and N c (C:G:D*, C*, D***): Weekday ultradian period of the females belonging to C2 tended to be significantly longer than those to C3. N c tended to be larger in C3 than C4, and on weekends than weekdays.

Discussion
The ability of demographic characteristics such as gender, job, and age to differentiate quantitative sleep-related features and subjective sleep quality and their associations has been shown previously [25][26][27]. In the present study, their demographics dependency was further analyzed. To quantitatively characterize activity during sleep, IMOIs and Table 6 Multiple dependency of subjective sleep quality and quantitative sleep-related parameters on demographic and HSWP characteristics (F test obtained from mixed-effect ANOVA) C HSWP cluster,G gender, J job, D weekday/weekend * p < 0.05, ** p < 0.01, *** p  ultradian inactivity rhythms were estimated, in addition to the conventional parameters produced by the accelerometer algorithm [19,21]. In our analysis, longer IMOIs responsible for body sway or posture change other than jitter-like events [28] were extracted. Ultradian inactivity rhythms were characterized based on a similar range of parameters reported in a previous study [21]. In part, our results shared a previously reported tendency for students to take a later, longer, and more irregular sleep compared with workers, irrespective of gender [10,13,27]. Despite these features, student's sleep was found to be less disturbed than that of the workers. This paradoxical finding might be due to the worker's obligatory regularity of lifestyle, which is suggested by the pronounced SJL for the workers (Table 1), although there is a possibility of confoundment by age. Regarding the gender difference, the male students tended to take a later and more disturbed sleep compared with the female students. There were more dense disturbances in the sleep of male workers than female workers. Although less distinctive, some scores of subjective sleep quality tended to be lower in the female subjects than in the males. Taken together, female subjects have tended to self-evaluate their sleep quality as worse than the quantitative features of less disturbed sleep in comparison with male subjects [21,27,29,30]. Worker's age dependency of sleep features revealed that the middle-aged group tended to have less disturbed sleep in a subjective as well as objective sense, but less satisfactory sleep in comparison with the other generation groups. In addition, the eldest generation group in other studies has also exhibited earlier HSWP and a flatter ultradian slope compared with younger groups [20,21]. Analysis of associations between quantitative sleeprelated parameters and subjective sleep quality for each demographic group showed non-negligible correlations of student's sleep duration with OSA scores except OSA-2 (sleep maintenance), which was also the case for parameters characterizing the total disturbances during sleep. This coincidence could be understood particularly when the occurrence density of disturbances shows not so much variation Dashed lines indicate factors that combine to form multi-way interactions through an individual's sleep. In contrast, there was an association of OSA-2 with parameters characterizing the occurrence density of disturbance. The present study is the first to reveal a categorical profile of associations between the subjective and quantitative features of sleep. Gender difference analysis has added the female-specific characteristics of non-negligible significant OSA-3 (dreaming)-related correlations to some quantitative feature parameters of sleep [31]. Compared with the students, for the male workers, non-negligible correlations of the parameters related to the total disturbances were found additionally for OSA-4 (fatigue recovery), in all generation groups. The eldest worker group (≥ 45 years) exhibited a slight but significant correlation of the disturbance density-related parameters with OSA-2, which was also shared by the students. That this tendency was shared by all generations of workers might suggest their awareness of fatigue recovery. In the job comparison between the male students and workers aged < 35 years, sleep duration was non-negligibly correlated also with OSA-2 for the workers in contrast to the students whose disturbance density-related parameters were found to be correlated only with OSA-2. The students might be more sensitive than the workers to disturbances partly, because the students' sleep is usually less disturbed.
The above analysis showed the fundamental demographics dependency of qualitative and quantitative features of sleep, and their associations. Nevertheless, the demographics dependency was mostly confined to the quantitative features, and the moderate association profiles were not sufficient to fully explain the behavior of sleep quality. To investigate another aspect of sleep, HSWPs were then characterized by clustering. Our classification of HSWPs revealed four clusters, C1-C4, which enabled differentiation of the relationship between quantitative sleep-related features and subjective sleep quality. Among the various indices of the static and dynamical aspects of HSWPs [2,3,8], chronotype was used directly for our classification; the second harmonics of the sleep-and wake-onset times included information on SJL, and SRI was in part shared by SDs and harmonics of the sleep-and wake-onset times. In this sense, our classification was conducted through a composite version of these conventional measures. Testing of bias in demographic proportions across clusters revealed significant biases in job and age group but not in gender (Table S-4 in Online Resource 3). C1 and C3 shared a similar tendency of earlier habitual sleep. In contrast, C2 and C4 shared a similar tendency of later habitual sleep. This cluster dependency might be at least partly due to the demographic structure of clusters that more workers (elderly subjects) in the former clusters and more students (younger subjects) in the latter clusters were included, respectively, considering the job-and age-dependent sleep features known previously [20,32]. In other words, patterns of day-to-day variations in wake-and sleep onsets and SJL, among others, are thought to further differentiate detailed classification. The PCA further provided a reduced dimensional representation of the clusters formed by the composite measures.
By additionally including the results of HSWP classification, subjective sleep quality and quantitative features of sleep were shown to depend in multiple ways on demographics (Gender and Job), HSWP clustering results, and distinction of weekdays/weekends (Day). Further details of dependencies were analyzed by the post hoc test. Sleep phase-related timings were found to exhibit the inter-cluster dependency similar to Table S-3, irrespective of weekdays or weekends. Common habit of SJL was well recognized from the tendency that most of these timings were more delayed on weekends in comparison with weekdays [22,33] In addition, the workers tended to wake up earlier than the students [32]. Chronobiological features of daily activity revealed the gender-specific tendency that the male's circadian amplitude was larger on weekends than weekdays [10,20]. Circadian acrophase characterized C2 as the cluster with significantly later HSWP than the others. Sleep duration was found to be significantly prolonged on weekends within some HSWP clusters in a gender-specific manner. This could be understood at least partly from the common lifestyle causing SJL [22,33]. Increasing tendency of N c on weekends might be associated with the prolongation of sleep duration because of no significant effect of day on the ultradian period. Disturbances were more pronounced on the weekend sleep of workers than the students. Especially for relatively regular HSWP clusters, the male's weekend sleep tended to be more disturbed than the weekday sleep. Such paradoxical tendency of disturbed weekend sleep for relatively regular HSWPs needs further investigation to be clarified. On the other hand, subjective sleep quality in terms of OSA-4 and OSA-5 tended to be better on weekends than weekdays. This tendency suggests multifaceted relationship between quantitative features of sleep and subjective sleep quality.
This study has some limitations. The actigraphy data analyzed here were collected from university students and information technology workers. However, the demographic balance of subjects was not well organized. The ages of female workers were unevenly distributed, which made age-related comparison within them difficult. Age distribution of female workers needs to be further diversified. In addition, the subjective sleep quality and the quantitative sleep features were not fully explained by the demographics and HSWP. Such a disagreement has commonly been recognized, although the sleep quality questionnaires used previously were different from those in the present study [6,34]. More elaborate analysis and additional observations might be necessary to improve the predictability of sleep quality, and will be a topic of future study. There seems to be inconsistency in usage of statistics that the non-parametric test was used in Table 1, and the multi-way ANOVA in Table 6. However, the multi-way ANOVA has been empirically known to be robust against non-normality at least for lower order interactions [24,25]. In addition, we confirmed that the residuals of the multiple dependency analysis did not so seriously deviate from the normality by a quantile-quantile plot. Highly accurate and elaborate analysis taking non-normality into account will be a future subject [35].

Conclusions
Subjective sleep quality was evaluated in light of various aspects-based OSA scores, each of which was associated in a particular way with quantitative sleep-related parameters. The associations were differentiated according to HSWP clusters. Further quantitative dissection of subjective sleep quality, considering multifaceted factors, such as gender, age, and occupation, could improve the validity of the present results, and enable selection of an HSWP fitted to an individual's working style and preferences.