ORIGINAL RESEARCH article

Front. Sleep, 17 April 2025

Sec. Insomnia

Volume 4 - 2025 | https://doi.org/10.3389/frsle.2025.1459854

Gratitude and sleep disturbance in primary care patients: the mediating roles of health self-efficacy, health behaviors, and psychological distress

  • 1. Department of Psychiatry and Behavioral Sciences, Johns Hopkins University School of Medicine, Baltimore, MD, United States

  • 2. Department of Psychiatry and Behavioral Sciences, Quillen College of Medicine, East Tennessee State University, Johnson City, TN, United States

  • 3. Ärzte der Welt e.V. [Doctors of the World], Munich, Germany

  • 4. Division of Health Promotion, Faculty for Applied Natural Sciences and Health, Coburg University of Applied Sciences and Arts, Coburg, Germany

  • 5. Institute of General Practice and Family Medicine, Faculty of Medicine, Ludwig Maximilian University of Munich, Munich, Germany

  • 6. Department of Psychology, Luther College, Decorah, IA, United States

  • 7. Department of Psychology, Durham University, Durham, United Kingdom

  • 8. Gasteiner Heilstollen, Böckstein, Austria

  • 9. Institute of Ecomedicine, Paracelsus Medical University, Salzburg, Austria

Abstract

Introduction:

Sleep disturbances are prevalent among primary care patients, and psychological dysfunction, including stress, anxiety, and depression, are robust contributors to poor sleep health. Yet, the presence of potential protective characteristics, such as health self-efficacy and engaging in adaptive health behaviors, may mitigate such outcomes. Gratitude (i.e., recognition and appreciation of experiences, relationships, and surroundings), a positive psychological cognitive-emotional characteristic, may serve as a catalyst of these beneficial downstream effects, given its association with improved health functioning and sleep.

Methods:

In a sample of primary care patients (N = 869, M age = 53; 60.7% female) from 50 urban and 30 rural practices in Germany, health self-efficacy (i.e., belief in ability to perform necessary actions to manage health) and constructive health behaviors (i.e., actions taken to modify health positively), separately and together as parallel first-order mediators, and stress, anxiety, and depression, as parallel second-order mediators, were investigated as potential serial mediators of the association between gratitude and sleep disturbances. Participants completed self-report measures in person and online.

Results:

Significant serial mediation was observed across models, although effects varied. In the first model, gratitude was associated with greater health self-efficacy and, in turn, to less stress (a1d21b4 = −0.019, 95% CI [−0.039, −0.002], SE = 0.010), anxiety (a1d31b5 = −0.026, 95% CI [−0.045, −0.008], SE = 0.009), and depression (a1d41b6 = −0.020, 95% CI [−0.040, −0.003], SE = 0.009), and fewer consequent sleep disturbances. In the second model, health behaviors, and anxiety (a1d31b5 = −0.009, 95% CI [−0.019, −0.002], SE = 0.004) and depression (a1d41b6 = −0.007, 95% CI [−0.016, −0.001], SE = 0.004), were serial mediators, but health behaviors and stress were not. In a final combined model, serial mediation occurred on two pathways, health self-efficacy and anxiety (a1d41b6 = −0.026, 95% CI [−0.046, −0.009], SE = 0.009), and health self-efficacy and depression (a1d51b7 = −0.019, 95% CI [−0.037, −0.003], SE = 0.009), and a specific indirect effect was found for health behaviors (a2b4= −0.086, 95% CI [−0.140, −0.030], SE = 0.026), but not self-efficacy.

Discussion:

Overall, gratitude was associated with reduced sleep disturbances through positive health behavior engagement, and via the serial mediation effects of greater health self-efficacy and lower psychological distress. Clinical interventions that enhance gratitude (e.g., gratitude listing or diaries), self-efficacy (e.g., disease self-management programs), or health behavior engagement (e.g., weight management programs) may promote favorable downstream effects on psychological distress and sleep disturbances among primary care patients.

1 Introduction

Sleep disturbances are prevalent worldwide, with 44% of adults reporting difficulty falling asleep and 35% indicating difficulty maintaining sleep (Aernout et al., 2021). Forty percent of primary care visits address sleep-related complaints (Arroll et al., 2012). In the United States, the number of outpatient ambulatory visits with sleep disturbance listed as the chief complaint increased 29% in a decade (Ford et al., 2014). Problems initiating and maintaining sleep are the defining criteria of the most prevalent sleep disorder, insomnia, which afflicts 22% of the global population and is more common (1.6:1) in females than males (Zeng et al., 2020). Insomnia symptoms occur in approximately one-third of primary care patients and difficulty maintaining sleep is the most prevalent characteristic, present in 80% of primary care patients with insomnia (Léger et al., 2010). Both oversleeping and undersleeping are associated with musculoskeletal, endocrine, respiratory, and neurological disorders, inflammation, higher cardiometabolic risk, and mortality (Åkerstedt et al., 2017; Irwin et al., 2016; Kanagasabai and Chaput, 2017; Ohayon et al., 2013; Shorofsky et al., 2019).

Given the deleterious impact of sleep disturbance on health, detection and intervention are imperative (Chattu et al., 2018). Individual-level cognitive-emotional factors are often more amenable to clinical intervention than fixed biological contributors (Hale et al., 2020). Gratitude is one such factor, and has been linked to improvements in health functioning, including sleep (Jackowska et al., 2016; Ng and Wong, 2013). Gratitude is an affective trait or emotion that encompasses appreciation for benefits received, including thankfulness for the present moment, nature, beauty, and life circumstances (McCullough et al., 2002; Wood et al., 2008). Several studies have demonstrated that gratitude can improve sleep. A systematic review of randomized controlled trials revealed that gratitude interventions (e.g., gratitude diaries, daily gratitude lists) improved subjective sleep quality in five samples, including individuals diagnosed with anxiety and depression, patients with neuromuscular disease, community adults, and college students (Boggiss et al., 2020). In these and other studies, the influence of gratitude on sleep disturbances is rarely explored without evaluating cognitive-emotional mechanisms or psychopathology. For example, in an RCT among German community members, a gratitude intervention improved gratitude and reduced worry, anxiety, depression, and insomnia symptoms post-intervention and at 3- and 6- month follow-ups (Heckendorf et al., 2019). In another RCT, individuals endorsing anxiety and depression diagnoses reported reduced stress, anxiety, depression, and sleep difficulties after completing a 3-week gratitude diary and, at 3-week follow-up, stress scores were maintained and anxiety scores had significantly improved (Southwell and Gould, 2017). Finally, among young adults, depressive symptoms explained the effect of gratitude on longer sleep duration and higher daytime energy (Alkozei et al., 2019), and among Austrian patients with rheumatic and musculoskeletal disease, gratitude was associated with less stress, anxiety, and depressive symptoms and, in turn, to better sleep quality and less functional impairment (Hirsch et al., 2021).

Extant research provides support for a beneficial effect of gratitude on amelioration of stress, anxiety, and depression, and sleep disturbances, both directly and indirectly. Yet, the mechanisms responsible for this protective influence have remained largely unexplored. In our theoretical model (see Figure 1), we propose that health self-efficacy, or belief in ability to perform health-supporting behaviors (Sirois, 2007), and engagement in wellness activities, are two potential explanatory linkages (Cousin et al., 2021). Drawing from the broaden-and-build theory (Fredrickson, 2004), gratitude may promote wellbeing by broadening one's awareness, thoughts, and behaviors over time, thereby promoting the development of resources (e.g., health, relationships, skills) that can be utilized to cope with stress (Fredrickson, 2004; Wood et al., 2010). Additionally, the coping hypothesis suggests that gratitude encourages adaptive coping strategies, such as planning and social support, which protect against stress and other adverse mental health outcomes (Wood et al., 2007). Finally, according to the positive affect hypothesis, gratitude, as a pleasant emotion, may produce additional constructive emotions that promote wellbeing and life satisfaction (Wood et al., 2010).

Figure 1

Consistent with these theories, prior studies independently link gratitude to health-promoting beliefs (e.g., health self-efficacy), behaviors (e.g., exercise, diet), and reduced psychological distress (DuBois et al., 2016; Grant and Gino, 2010). For example, among patients with heart failure, gratitude was associated with improved self-efficacy to maintain physical functioning and, in turn, to better medication adherence (Cousin et al., 2020), and with greater self-efficacy to preserve heart function, less depression, better sleep, and lower levels of inflammatory markers (Mills et al., 2015). Further, gratitude for one's health predicted heightened physical activity, as assessed via accelerometer, and gratitude for one's life predicted medical adherence after 6 months, among acute coronary syndrome patients (Legler et al., 2019). In primary care samples, the roles of health self-efficacy and health behavior engagement in promoting psychological and sleep health is well-established (Alfaris et al., 2015; Bluestein et al., 2010). However, no studies have integrated these pathways into a unified model, nor clarified the distinct and shared roles of cognitive-emotional and behavioral mediators in the gratitude-sleep relationship.

Given that nearly 70% of primary care visits involve psychological complaints, including depression and anxiety (Hunter and American Psychological Association, 2009), significant comorbidity exists between psychological and physical conditions (Gili et al., 2010), and emerging conditions are most often identified and treated in primary care settings (Schoen et al., 2008), it is important to identify risk and protective factors that can be targeted for clinical intervention in primary care. Gratitude is one such clinically malleable factor (Heckendorf et al., 2019; Yuan et al., 2021), with robust linkages to adaptive health outcomes.

In the current study, we examined the potential serial mediation effects of health self-efficacy and health behaviors, and stress, anxiety, and depression, on the relation between gratitude and sleep disturbances. As perceived stress, anxiety, and depression frequently co-occur and interact and, in conjunction, exacerbate sleep disturbance (Demyttenaere and Heirman, 2020; Kalmbach et al., 2018; Karing, 2021) due to shared cognitive-emotional and biological mechanisms (Boggero et al., 2017; Packard et al., 2016; Thorsteinsson et al., 2019), we examined them as parallel, rather than serial or separate, mediators in our models. We hypothesized that health self-efficacy and health behaviors, separately and together as parallel first-order mediators, and stress, anxiety, and depression, as parallel second-order mediators, would serially mediate the association between gratitude and sleep disturbances, such that greater gratitude would be associated with higher levels of health self-efficacy and/or health behaviors and, in turn, to less stress, anxiety and depressive symptoms, and sleep disturbances.

2 Materials and methods

Over a period of 2 months, patients from 50 urban primary care practices in the Munich area and 30 rural primary care practices in the Upper Franconia region of Germany were recruited in person by clinic staff and providers, and via informational posters, to complete a questionnaire battery. Paper surveys could be returned to the practice in person or to either study center at Coburg University of Applied Sciences or the University of Munich via postal mail, and questionnaires could also be completed online. The original study was approved by the Ethics Commission of the Faculty of Medicine at the Ludwig Maximilian University of Munich, and participation was voluntary. To ensure privacy, the Institutional Review Board waived consent requirements and no written informed consent was obtained. To participate, respondents were required to be age 18 or older, to have sufficient knowledge of German, and to possess a willingness to participate. Participants under age 18 and with severe psychiatric disorders were excluded.

2.1 Participants

Participants (N = 869) were 39.1% male (n = 343) and 60.7% female (n = 533), ranging from 20 to 92 years of age (M = 53.00, SD = 11.96). Race and ethnicity data were not collected, as legal restrictions and sociopolitical pressures enjoin researchers from collecting ethnoracial data in Germany and much of Western Europe (Roig, 2017; Simon, 2012). Regarding education, most participants reported having an intermediate school certificate (n = 287; 32.7%), followed by a high school diploma (n = 252; 28.7%), some high school completed (n = 196; 22.3%), university entrance qualification (n = 75; 8.5%), and other certificate (n = 48; 5.5%). Most participants reported one chronic medical condition (n = 305; 34.7%), followed by zero (n = 206; 23.5%), two (n = 182; 20.7%), three (n = 97; 11.0%), four (n = 46; 5.2%), and five or more (n = 41; 4.6%).

2.2 Measures

In addition to measures assessing this study's variables, participants responded to demographic questions (e.g., age, family composition, zip code, marital status, height/weight), meant for characterizing the sample and to serve as covariates.

2.2.1 Sleep disturbances

Frequency of sleep disturbances was assessed using Item 3 from the Patient Health Questionnaire-9, which was developed to detect depressive symptoms among primary care patients (Kroenke et al., 2001). Item 3 was extracted from the German translation of the PHQ-9, developed by Gräfe et al. (2004) using internationally accepted translation methods and available on Pfizer's Patient Health Questionnaire website. Participants indicated, on a scale from 0 (not at all) to 3 (nearly every day), the frequency with which they experience “Trouble falling or staying asleep, or sleeping too much”, and a higher score indicates greater sleep difficulty. This item is considered a valid measure of sleep disturbances that can be used as a substitute for longer instruments. For example, in a study aiming to evaluate the item's utility for screening for sleep disturbances in primary care patients, the PHQ item 3 demonstrated convergent validity (r = 0.75, p < 0.001) with the Insomnia Severity Index (ISI; Morin, 1993), and a cutoff score of 1, which indicates sleep problems were present several days over the past 2 weeks, yielded the optimal balance of sensitivity (82.5%) and specificity (84.5%; MacGregor et al., 2012). Further, among German cancer patients, this item was strongly correlated with the ISI (r = 0.72, p < 0.001) and performed similarly to the ISI; the ISI and Item 3 exhibited similar correlations with the Patient Health Questionnaire-2 (PHQ-2; Kroenke et al., 2003; r = 0.45, 0.42), Generalized Anxiety Disorder 2-item scale (GAD-2; Kroenke et al., 2007; r = 0.45, 0.39), and had comparable magnitudes and patterns of association with sex, age, and tumor type (Schulte et al., 2021).

2.2.2 Perceived stress

The Perceived Stress Scale–4 (PSS-4; Cohen and Williamson, 1988), German version (Stächele and Volz, 2013), was used to evaluate stress. On a Likert scale, ranging from 0 (never) to 4 (very often), respondents indicated the degree of appraised stress over the past month via statements such as, “In the last month, how often have you felt that you were unable to control the important things in your life?” and “In the last month, how often have you felt that things were going your way?” After recoding two items to indicate greater perceived stress, items are summed to generate a total score between 0 and 16. The four items comprising the scale were extracted from the original scale, the PSS-14 (Cohen et al., 1983). In the current sample, the PSS-4 demonstrated acceptable reliability (Cronbach's α = 0.73; McDonald's ω = 0.74), and acceptable (α = 0.74) to good (α = 0.84) reliability in prior health samples (Vallejo et al., 2018; Wu and Amtmann, 2013).

2.2.3 Anxiety

The Generalized Anxiety Disorder 2-item scale (GAD-2; Kroenke et al., 2007) measures anxiety on a Likert scale ranging from 0 (not at all) to 3 (nearly every day), via items representing core anxiety symptoms from the Generalized Anxiety Disorder 7-item scale (GAD-7; Spitzer et al., 2006). The two items utilized were from the German translation of the GAD-7, developed by Löwe et al. (2008) using internationally accepted translation methods and downloaded from Pfizer's Patient Health Questionnaire website. Participants indicate how often over the last 2 weeks they were bothered by “Feeling nervous, anxious or on edge” and “Not being able to stop or control worrying”. The total score, which ranges from 0 to 6, is obtained by summation of item scores. Higher scores denote greater anxiety symptoms, with a cutoff of 3 providing good sensitivity and specificity to detect clinically significant anxiety symptoms (Kroenke et al., 2007), which was supported in a systematic review (N = 5,223) of validation studies (Plummer et al., 2016). Internal consistency reliability was good in the current sample (α = 0.80; ω = 0.82). The GAD-2 exhibited acceptable reliability and strong correlations with the GAD-7 in samples of community-dwelling adults in Germany (α = 0.72; r = 0.87; Hinz et al., 2017) and multiple sclerosis patients (α = 0.77; r = 0.94; Hughes et al., 2018).

2.2.4 Depression

The Patient Health Questionnaire-2 (PHQ-2; Kroenke et al., 2003), comprised of the first two items of the original measure, the Patient Health Questionnaire-9 (Kroenke et al., 2001), was used to assess depressive symptoms over the past 2 weeks. The two items were taken from Gräfe et al.'s (2004) German translation of the PHQ-9. Respondents indicated how often they were bothered by “Little interest or pleasure in doing things” and “Feeling down, depressed, or hopeless”, on a scale ranging from 0 (not at all) to 3 (nearly every day). Responses were summed for a total score ranging from 0 to 6, with higher scores indicating greater depressive symptoms. A score of 3 is considered the cutoff for detecting clinically significant depression (Staples et al., 2019), although some research suggests a more modest cutoff of 2 (Manea et al., 2016). In the current study, reliability of the PHQ-2 was good (α = 0.79; ω = 0.83), and in a sample of German COPD inpatients, composite reliability was good (CR = 0.89; Schuler et al., 2018). Among 1,619 German primary care outpatients, the PHQ-2 was strongly correlated with the PHQ-9 (r = 0.87; Löwe et al., 2005).

2.2.5 Gratitude

Gratitude was measured using the Gratitude Questionnaire-Six Item Form (GQ-6; McCullough et al., 2002), German translation (Personality and Assessment Group, Department of Psychology, University of Zurich, n.d.), which includes six statements (e.g., “I am grateful to a wide variety of people”; “If I had to list everything that I felt grateful for, it would be a very long list”) assessing appreciation. Responses are rated on a scale ranging from 1 (strongly disagree) to 7 (strongly agree) and are summed to compute a total score between 6 and 42, after reverse scoring two items so that higher scores denote higher gratitude levels. In the current sample, internal consistency reliability was acceptable (α = 0.77; ω = 0.77). In a meta-analysis, internal consistency reliability was acceptable in European studies (α = 0.75), and good when averaging all 74 studies (α = 0.82; Card, 2019).

2.2.6 Health self-efficacy

Health self-efficacy, or the degree of confidence one possesses in their ability to engage in the necessary actions to control their health, was assessed using six items from the 8-item Health Self-Efficacy/Mastery Beliefs subscale of the Control Beliefs Inventory (HSE-CBI; Sirois, 2003a), a 26-item measure of perceived control. Two items assessing aspects of health self-efficacy incorporated in the remaining items with alternative wording were excluded to reduce participant burden. As one item is reverse scored, the omission is not expected to significantly affect internal consistency. The original measure was translated from English into German, and then back-translated to English by a panel comprising members of the research team who are bilingual native German speakers fluent in English. Inconsistencies were reconciled between the researchers. Respondents indicate on a scale ranging from 1 (strongly disagree) to 6 (strongly agree) the degree to which they agree with statements such as “I am confident that I can successfully look after my health” and “I am able to meet the challenge of following a healthy routine”. After reverse scoring two negatively worded items, the average of the six items is computed to calculate an HSE-CBI total score between 1 and 6, with a higher score indicating a greater degree of health self-efficacy. In the current study, the HSE-CBI has acceptable reliability (α = 0.77; ω = 0.73), and good reliability (α = 0.84) in two samples of community adults (Sirois, 2004). The HSE-CBI demonstrated convergent validity with the General Perceived Self-Efficacy Scale (Schwarzer and Jerusalem, 1995) among community adults (Sirois, 2003b).

2.2.7 Health behaviors

Participants completed the Wellness Behaviors Inventory (WBI), a measure that assesses the frequency of engagement in health-promoting behaviors, such as exercising, relaxation, and healthy eating (Sirois, 2001). The WBI was translated into German using the same procedure as the HSE-CBI. Responses to items such as “I eat breakfast” and “I exercise for 20 continuous minutes or more, to the point of perspiration” are provided on a 5-point Likert scale ranging from 1 (less than once a week or never) to 5 (every day of the week). Although the WBI consists of 12 items, only the 10 items utilized to compute the total score were included in this study. The two excluded items refer to vitamin and supplement use and form a separate optional index of health behaviors (Sirois, 2001). Two items are reverse scored before calculating the total mean score, which ranges from 1 to 5, with higher scores denoting greater engagement in health behaviors. In the current primary care sample, the WBI had questionable internal consistency reliability (α = 0.63; ω = 0.64), consistent with prior research; in a meta-analysis of 15 undergraduate and community adult samples, reliability ranged from questionable to acceptable (α = 0.64–74; Sirois et al., 2015), likely due to its checklist status.

2.3 Statistical analyses

2.3.1 Bivariate analyses

All analyses were conducted in R Version 2022.07.1 (R Core Team, 2022). To establish the presence of bivariate relationships before testing mediation effects, Pearson's product-moment and Spearman's rank-order correlation coefficients were calculated to assess the associations among gratitude, health self-efficacy, health behaviors, stress, anxiety, and depression, and the correlations between each of these variables and sleep disturbances were evaluated with Spearman's rank-order correlation coefficients, using the R package psych (Revelle, 2022). A multicollinearity cutoff of r > 0.80 between study variables was used, as associations exceeding this threshold may indicate lack of independence and would be interpreted with caution (Field, 2013).

2.3.2 Serial multivariate mediation analyses

Prior to conducting the mediation analyses, skewness, kurtosis, and heteroscedasticity of the data were evaluated using base R functions, to ensure the assumptions of linear regression were met (R Core Team, 2022). All predictors were standardized (i.e., converted to z-scores) because they have no natural or common metric, allowing the relative influence of each predictor to be meaningfully compared in models containing both linear and ordinal logistic regression analyses (Agresti, 1996; Menard, 2004). Sleep disturbances, being an ordered polytomous variable, was not standardized.

Three serial mediation path analyses were conducted. In the first model, the relation between gratitude and sleep disturbances, via the potential mediating effects of health self-efficacy as a first-order mediator, and stress, anxiety, and depression, as second-order parallel mediators, was assessed. In the second model, health self-efficacy was replaced with health behaviors as the first-order mediator. Finally, in the third model, health self-efficacy and health behaviors were evaluated as parallel first-order mediators. Age, sex, education level, and number of chronic diseases were included as covariates in each model.

2.3.2.1 Analytic procedure

As the outcome variable, sleep disturbances, is not continuous and, instead, was measured using one item with four levels, the final regression equation of the serial mediation model (i.e., sleep disturbances regressed on gratitude and all mediators and covariates) was calculated using cumulative odds ordinal logistic regression with the VGAM package (Yee, 2015). All other pathways of the serial mediation model were analyzed using ordinary least squares linear regression, via the stats package (R Core Team, 2022).

For the ordinary least squares regression steps of each model, a change of one standard deviation in the predictor is associated with a change of β standard deviations in the endogenous variable. Regarding the ordinal logistic regression step of each model, a one standard deviation difference in the predictor is associated with a β-unit difference in ordered logits of sleep disturbance frequency (Menard, 2004). To simplify interpretation of the ordinal logistic regression results, the standardized coefficients were converted to odds ratios (OR), so that for each one standard deviation increase in the predictor, the odds of being more likely to endorse a higher frequency of sleep disturbances is multiplied OR times, holding constant all other variables.

Regarding the final step of each model, ordinal logistic and multinomial logistic regression models were compared, using the likelihood ratio test, to determine which type of model fit the data better, via the VGAM package (Yee, 2015). Next, the proportional odds assumption of ordinal logistic regression (i.e., stating each predictor has an equivalent effect, or parallel slopes, across all levels of the outcome variable; McCullagh, 1980) was evaluated, using likelihood ratio tests in the ordinal package (Christensen, 2019). If the assumption was violated, a partial proportional odds model was fit, which relaxes the proportional odds assumption and allows multiple regression coefficients only for predictors that have differential effects on the levels of the outcome variable (Peterson and Harrell, 1990).

Subsequently, non-significant covariates were removed from the linear and ordinal logistic regression analyses in a backward stepwise manner to preserve power, unless their inclusion improved overall model fit (Faraway, 2015). Finally, a function was written to calculate all mediation effects for each model, by multiplying coefficients for each path (e.g., the serial mediation effect of gratitude on sleep through health self-efficacy and stress is calculated as a1d31b5). As each model used standardized regression coefficients and continuous mediators, the mediated effect may be obtained by multiplying coefficients (i.e., using the product of coefficients method) from both ordinary least squares and ordinal logistic regression pathways and bootstrapping the confidence intervals (MacKinnon et al., 2007). Bootstrapping using 10,000 simulated samples and 95% confidence intervals, per the percentile method, was utilized in each model, to estimate the sampling distribution of the mediation effects. All p-values for the regression coefficients were adjusted using the p.adjust function in the stats package (R Core Team, 2022), applying the False Discovery Rate (FDR) method (Benjamini and Hochberg, 1995) to control for multiple comparisons across models. Pathways were considered statistically significant if the FDR-adjusted p-values were below the conventional threshold of 0.05. To correct the confidence intervals for the mediation effects, the False Coverage Rate (FCR) method was applied at a threshold of 0.97, which was calculated using the proportion of ratio of null effects to total effects across models. This method controls for false coverage in the estimation of CIs, reducing the risk of obtaining misleading conclusions about significance (Benjamini and Yekutieli, 2005). Mediation effects were considered statistically significant if the corrected 97% CIs did not include zero. Finally, multicollinearity for all predictors in each model was assessed by calculating the variance inflation factor (VIF) with the R package car (Fox and Weisberg, 2019). The VIF of each coefficient should be < 5 (Hair et al., 2019), and ideally around 1 (Akinwande et al., 2015).

Of note, path analysis was determined to be a more appropriate statistical approach than structural equation modeling (SEM) because we were interested in the relationships between variables and the mediation pathways, rather than any latent constructs or the measurement properties of the scales, and our one-item outcome variable precludes the use of comprehensive latent SEM. Further, the robust estimators available in R (maximum likelihood, weighted least squares mean and variance), cannot be used with bootstrapping for models with ordinal outcome variables, limiting methods for confidence interval estimation, and cannot be modeled to allow multiple slopes for a predictor when the proportional odds assumption is violated (Rosseel, 2012). However, as a supplementary analysis, we conducted SEM to obtain overall fit statistics for each model, dichotomizing the outcome variable and specifying psychological distress as a latent variable comprised of stress, anxiety, and depression. More details regarding the procedure and results can be found in the Supplementary material.

3 Results

3.1 Bivariate analyses

Pearson's bivariate and Spearman's rank-order correlation analyses indicated that gratitude was positively correlated with health self-efficacy (r = 0.27, p < 0.001) and health behaviors (r = 0.09, p < 0.01), and negatively related to stress (r = −0.17, p < 0.001), anxiety (r = −0.12, p < 0.001), depression (r = −0.15, p < 0.001), and sleep disturbances (r = −0.10, p < 0.01). Similarly, health self-efficacy and health behaviors were positively related to each other (r = 0.30, p < 0.001) and negatively related to stress, anxiety, depression, and sleep disturbances (all p < 0.001). Bivariate correlations between all continuous and ordinal study variables are provided in Table 1. No correlations approached or exceeded the multicollinearity cutoff of r > 0.80.

Table 1

Mean (SD)HSEHBStressAnxietyDepressionSleepaAgeEducationaConditions
GS33.44(6.53)0.270.09#−0.17−0.12−0.15−0.10#−0.07*0.17−0.16
HSE4.50(0.75)0.30−0.31−0.20−0.23−0.25−0.070.23−0.28
HB3.55(0.65)−0.21−0.13−0.16−0.230.180.06−0.07*
Stress5.81(3.09)0.490.510.29−0.04−0.180.28
Anxiety1.35(1.42)0.620.31−0.05−0.050.26
Depression1.28(1.38)0.29−0.06−0.110.32
Sleepa1.19(0.97)0.06−0.120.28
Age53.00(11.96)−0.210.25
Educationa−0.20
Conditions1.79(1.50)

Means, standard deviations, and bivariate correlations of study variables.

GS, Gratitude: Gratitude Questionnaire–6; HSE, Health self-efficacy: Health Self-Efficacy/Mastery Beliefs subscale, Control Beliefs Inventory; HB, Health behaviors: Wellness Behaviors Inventory; Stress, Perceived stress: Perceived Stress Scale–4; Anxiety, Generalized Anxiety Disorder 2-item scale; Depression, Patient Health Questionnaire–2; Sleep, Sleep disturbances: Patient Health Questionnaire–9, Item 3.

Covariates: Age = participant age; Education = participant education level; Conditions = number of chronic medical conditions.

aSpearman's rank-order correlation coefficients.

*p < 0.05.

#p < 0.01.

p < 0.001.

3.2 Serial multivariate mediation analyses

Examination of skewness, kurtosis, and heteroscedasticity indicated the assumptions of linear regression were met. Table 2 shows the frequencies of the ordinal outcome variable, sleep disturbances, by response category, and normality data for the item.

Table 2

ResponseFrequencySkewnessKurtosis
n%ValueSEValueSE
Not at all21224.6
Several days39645.9
More than half the days13115.2
Nearly every day12314.3
Total8621000.560.08−0.600.17

Frequencies, skewness, and kurtosis for sleep disturbances.

Skewness and kurtosis values indicate data are normally distributed.

3.2.1 Model 1

Regarding model 1, the ordinal logistic regression model fit the data significantly better than the multinomial model ( = 538.92, p ≤ 0.001) and was more parsimonious (i.e., contained the fewest regression coefficients). However, the proportional odds assumption was violated for health self-efficacy ( = 7.312, p = 0.026). Therefore, a partial proportional odds model incorporating three regression coefficients for health self-efficacy was used.

Linear and ordinal logistic regression results and fit statistics for each pathway, including significant covariates identified through backward stepwise selection, are provided in Table 3. Variance inflation factors are also provided in Table 3, and no predictors exceed the multicollinearity cutoff. Regression results excluding covariate effects are illustrated in Figure 2. Gratitude was related to higher health self-efficacy (β = 0.273, 95% CI [0.209, 0.337], SE = 0.032, t = 8.412, padjusted < 0.001) which was, in turn, related to lower stress (β = −0.320, 95% CI [−0.387, −0.253], SE = 0.034, t = −9.413, padjusted < 0.001), anxiety (β = −0.218, 95% CI [−0.288, −0.147], SE = 0.036, t = −6.078, padjusted < 0.001), and depression (β = −0.255, 95% CI [−0.322, −0.189], SE = 0.034, t = −7.530, padjusted < 0.001). A one standard deviation increase in stress was associated with a 24% greater likelihood of endorsing a higher frequency of sleep disturbances (OR = 1.242, 95% CI [1.091, 1.414], padjusted = 0.002), while equivalent increases in anxiety (OR = 1.553, 95% CI [1.357, 1.778], padjusted < 0.001) and depression (OR = 1.337, 95% CI [1.162, 1.540], padjusted < 0.001) were associated with 55% and 34% greater likelihood of sleep problem frequency, respectively.

Table 3

EffectaβSEβ95% CI βtpAdj. pVIFFit statistics
LLUL
HSER2 = 0.201;R2adjusted = 0.197F(4,807) = 50.74,p < 0.001
(Intercept)− 0.3610.179− 0.713− 0.010− 2.0180.0440.055
Gratitude0.2730.0320.2090.3378.412< 0.001< 0.0011.052
Age0.0080.0030.0030.0142.9380.0030.0051.158
Education0.0770.0230.0310.1233.3000.0010.0021.074
Conditions− 0.1960.023−0.242−0.151− 8.488< 0.001< 0.0011.181
StressR2 = 0.256;R2adjusted = 0.252F(4,814) = 69.96,p < 0.001
(Intercept)0.4490.1390.1760.7223.2260.0010.002
HSE− 0.3200.034− 0.387− 0.253− 9.413< 0.001< 0.0011.239
Gratitude− 0.1540.033− 0.218− 0.090− 4.713< 0.001< 0.0011.141
Age− 0.0140.003− 0.019− 0.008− 4.927< 0.001< 0.0011.149
Conditions0.1470.0230.1020.1936.335< 0.001< 0.0011.300
AnxietyR2 = 0.178;R2adjusted = 0.172 F(5,816) = 35.13, p < 0.001
(Intercept)0.1280.159− 0.1840.4390.8050.4210.421
HSE− 0.2180.036− 0.288− 0.147− 6.078< 0.001< 0.0011.242
Gratitude− 0.0680.034− 0.136− 0.001− 1.9870.0470.0571.145
Age− 0.0110.003− 0.017− 0.006− 3.953< 0.001< 0.0011.164
Sex (female)0.2550.0660.1250.3863.844< 0.001< 0.0011.037
Conditions0.1820.0240.1340.2307.455< 0.001< 0.0011.301
DepressionR2 = 0.252;R2adjusted = 0.248F(5,816) = 55.09,p < 0.001
(Intercept)0.3290.1500.0350.6232.1940.0290.038
HSE− 0.2550.034− 0.322− 0.189− 7.530< 0.001< 0.0011.241
Gratitude− 0.1220.032− 0.186− 0.059− 3.766< 0.001< 0.0011.145
Age− 0.0150.003− 0.020− 0.010− 5.490< 0.001< 0.0011.163
Sex (female)0.1480.0630.0250.2712.3590.0190.0271.037
Conditions0.2080.0230.1630.2549.048< 0.001< 0.0011.300
EffectbβSEβ95% CI βzpAdj. pVIFOR95% CI OR
LLULLLUL
SleepLL = −886.06, df = 2,351
0|10.8170.2810.2661.3682.9060.0040.0062.2641.3053.927
1|2− 1.5970.280− 2.146− 1.048− 5.701< 0.001< 0.0010.2020.1170.351
2|3− 2.7100.290− 3.280− 2.141− 9.331< 0.001< 0.0010.0670.0380.118
Gratitude0.0860.053− 0.0170.1901.6310.1030.1111.1791.0900.9831.209
HSE1− 0.1000.090− 0.2760.077− 1.1090.2670.2671.397c0.9050.7591.080
HSE2− 0.3650.081− 0.522− 0.207− 4.527< 0.001< 0.0010.6940.5930.813
HSE3− 0.2060.103− 0.408− 0.003− 1.9920.0460.0550.8140.6650.997
Stress0.2170.0660.0870.3463.2760.0010.0022.0131.2421.0911.414
Anxiety0.4400.0690.3050.5756.392< 0.001< 0.0012.2351.5531.3571.778
Depression0.2910.0720.1500.4324.047< 0.001< 0.0012.3841.3371.1621.540
Age0.0090.0040.0000.0182.0550.0400.0541.2361.0091.0001.018
Education− 0.0740.036− 0.145− 0.003− 2.0390.0410.0541.1100.9290.8650.997
Conditions0.1770.0370.1050.2504.783< 0.001< 0.0011.4161.1941.1101.284

Regression coefficients of serial indirect effects Model 1 for gratitude, health self-efficacy, psychological distress, and sleep disturbances, with covariates.

Standardized coefficients are reported.

CI, confidence interval; LL, lower limit; UL, upper limit. Adj. p, Adjusted p-value (FDR correction); VIF, variance inflation factor; OR, odds ratio; HSE, health self-efficacy.

Bolded p-values indicate those that became non-significant after adjustment. 0|1: sleep disturbances threshold; the likelihood of scoring a 0 on sleep disturbances as opposed to 1−3. 1|2: sleep disturbances threshold; the likelihood of scoring a 0−1 on sleep disturbances as opposed to 2−3. 2|3: sleep disturbances threshold; the likelihood of scoring 0−2 on sleep disturbances as opposed to 3. HSE1: The likelihood of a one standard deviation increase in health self-efficacy being associated with higher sleep disturbances scores of 1−3 vs. 0. HSE2: The likelihood of a one standard deviation increase in health self-efficacy being associated with higher sleep disturbances scores of 2-3 vs. 0−1. HSE3: The likelihood of a one standard deviation increase in health self-efficacy being associated with a higher sleep disturbances score of 3 vs. 0−2.

Covariates: Age = participant age; Sex = participant sex; Education = participant education level; Conditions = number of chronic medical conditions.

aMultiple linear regression analyses.

bOrdinal logistic regression analysis.

cVariance inflation factor for health self-efficacy.

Figure 2

Mediation effect results for model 1 are provided in Table 4. As hypothesized, health self-efficacy as a first-order mediator, and stress (a1d21b4= −0.019, 95% CI [−0.039, −0.002], SE = 0.010), anxiety (a1d31b5 = −0.026, 95% CI [−0.045, −0.008], SE = 0.009), and depression (a1d41b6 = −0.020, 95% CI [−0.040, −0.003], SE = 0.009) as parallel second-order mediators, serially mediated the relation between gratitude and sleep disturbances. Additionally, specific indirect effects linking gratitude and sleep disturbances were observed for one of the three health self-efficacy pathways, HSE2 (a1b2 = −0.099, 95% CI [−0.166, −0.043], SE = 0.031), and stress (a2b4 = −0.033, 95% CI [−0.072, −0.003], SE = 0.018) and depression (a4b6 = −0.036, 95% CI [−0.078, −0.004], SE = 0.019). Overall, the total indirect effect (β = −0.347, 95% CI [−0.530, −0.066], SE = 0.105) and the total effect of gratitude on sleep disturbances (c = −0.261, 95% CI [−0.478, −0.030], SE = 0.112) were significant, while the direct effect was non-significant (c' = 0.086, 95% CI [−0.066, 0.235], SE = 0.078), indicating mediation. After adjusting the confidence intervals to 0.97 to account for false coverage rate correction, the conclusions regarding all indirect effects remained consistent.

Table 4

EffectβSEβ95% CI β
LLUL
Gratitude → self-efficacy1 → sleepa1 b1−0.0270.028−0.0840.027
Gratitude → self-efficacy2 → sleepa1 b2−0.0990.031−0.166−0.043
Gratitude → self-efficacy3 → sleepa1 b3−0.0560.046−0.1370.064
Gratitude → stress → sleepa2 b4−0.0330.018−0.072−0.003
Gratitude → anxiety → sleepa3 b5−0.0300.019−0.0720.003
Gratitude → depression → sleepa4 b6−0.0360.019−0.078−0.004
Gratitude → self-efficacy → stress → sleepa1 d21 b4−0.0190.010−0.039−0.002
Gratitude → self-efficacy → anxiety → sleepa1 d31 b5−0.0260.009−0.045−0.008
Gratitude → self-efficacy → depression → sleepa1 d41 b6−0.0200.009−0.040−0.003
Totalc−0.2610.112−0.478−0.030
Total indirect−0.3470.105−0.530−0.066
Directc0.0860.078−0.0660.235

Indirect effects of health self-efficacy and psychological distress on the relation between gratitude and sleep disturbances (Model 1).

Covariates include age, sex, education level, and number of chronic diseases. Standardized coefficients are reported.

CI, confidence interval; LL, lower limit; UL, upper limit. Self-efficacy = health self-efficacy.

Self-efficacy1: Specific indirect effect of gratitude related to sleep disturbances scores of 1–3 vs. 0, via health self-efficacy. Self-efficacy2: Specific indirect effect of gratitude related to sleep disturbances scores of 2–3 vs. 0–1, via health self-efficacy. Self-efficacy3: Specific indirect effect of gratitude related to sleep disturbances scores of 3 vs. 0–2, via health self-efficacy. All effects remained consistent after adjusting the CIs for the FCR at 0.97.

3.2.2 Model 2

Regarding model 2, the ordinal logistic regression model fit the data significantly better than the multinomial model ( = 98.897, p ≤ 0.001) and was more parsimonious. However, the proportional odds assumption was violated for health behaviors ( = 6.558, p = 0.038). Therefore, a partial proportional odds model including three regression coefficients for health behaviors was used.

Linear and ordinal logistic regression results and fit statistics for each pathway, including significant covariates identified through backward stepwise selection, are provided in Table 5. Variance inflation factors are also provided in Table 5, and no predictors exceed the multicollinearity cutoff. Regression results excluding covariate effects are illustrated in Figure 3. Gratitude was related to better health behavior engagement (β = 0.173, 95% CI [0.106, 0.239], SE = 0.034, t = 5.081, padjusted < 0.001) which was, in turn, related to lower stress (β = −0.209, 95% CI [−0.274, −0.144], SE = 0.033, t = −6.290, padjusted < 0.001), anxiety (β = −0.113, 95% CI [−0.180, −0.046], SE = 0.034, t = −3.314, padjusted = 0.002), and depression (β = −0.143, 95% CI [−0.207, −0.079], SE = 0.033, t = −4.393, padjusted < 0.001). A one standard deviation increase in stress was associated with a 17% greater likelihood of endorsing a higher frequency of sleep disturbances (OR = 1.168, 95% CI [1.026, 1.330], padjusted = 0.024), while equivalent increases in anxiety (OR = 1.602, 95% CI [1.394, 1.840], padjusted < 0.001) and depression (OR = 1.344, 95% CI [1.166, 1.549], padjusted < 0.001) were associated with 60% and 34% greater likelihood of sleep problem frequency, respectively.

Table 5

EffectaβSEβ95% CI βtpAdj. pVIFFit statistics
LLUL
HBR2 = 0.126;R2adjusted = 0.122F(4,817) = 29.38,p < 0.001
(Intercept)− 1.4810.187− 1.848− 1.114− 7.922< 0.001< 0.001
Gratitude0.1730.0340.1060.2395.081< 0.001< 0.0011.051
Age0.0260.0030.0210.0328.970< 0.001< 0.0011.158
Education0.0740.0250.0250.1222.9950.0030.0041.075
Conditions− 0.0970.024− 0.144− 0.050− 4.028< 0.001< 0.0011.179
StressR2 = 0.211;R2adjusted = 0.206F(5,820) = 43.8, p < 0.001
(Intercept)0.0820.162− 0.2350.4000.5080.6120.629
HB− 0.2090.033− 0.274− 0.144− 6.290< 0.001< 0.0011.120
Gratitude− 0.2210.033− 0.286− 0.157− 6.729< 0.001< 0.0011.082
Age− 0.0090.003− 0.015− 0.003− 3.0950.0020.0031.269
Sex (female)0.1130.065− 0.0150.2421.7330.0830.1041.043
Conditions0.1870.0230.1420.2328.180< 0.001< 0.0011.205
AnxietyR2 = 0.148;R2adjusted = 0.142F(5,826) = 28.58,p < 0.001
(Intercept)− 0.0070.166− 0.3320.318− 0.0410.9670.967
HB− 0.1130.034− 0.180− 0.046− 3.314< 0.0010.0021.123
Gratitude− 0.1140.034− 0.181− 0.048− 3.385< 0.0010.0011.084
Age− 0.0100.003− 0.016− 0.004− 3.2690.0010.0021.262
Sex (female)0.2570.0670.1250.3893.827< 0.001< 0.0011.041
Conditions0.2070.0230.1610.2538.844< 0.001< 0.0011.202
DepressionR2 = 0.217;R2adjusted = 0.213F(5,826) = 45.88,p < 0.001
(Intercept)0.1600.157− 0.1490.4691.0180.3090.348
HB− 0.1430.033− 0.207− 0.079− 4.393< 0.001< 0.0011.120
Gratitude− 0.1750.032− 0.238− 0.112− 5.433< 0.001< 0.0011.084
Age− 0.0130.003− 0.019− 0.007− 4.562< 0.001< 0.0011.259
Sex (female)0.1590.0640.0330.2842.4800.0130.0171.041
Conditions0.2420.0220.1980.28610.821< 0.001< 0.0011.201
EffectbβSEβ95% CI βzpAdj. pVIFOR95% CIOR
LLULLLUL
SleepLL = −888.96, df = 2,417
0|1− 0.2570.256− 0.7590.244− 1.0050.3150.3150.7730.4681.277
1|2− 2.7480.265− 3.267− 2.229− 10.374< 0.001< 0.0010.0640.0380.108
2|3− 3.9310.280− 4.480− 3.382− 14.024< 0.001< 0.0010.0200.0110.034
Gratitude0.0620.053− 0.0410.1661.1840.2360.2561.1531.0640.9601.180
HB1− 0.3900.093− 0.572− 0.208− 4.201< 0.001< 0.0011.175c0.6770.5640.812
HB2− 0.6540.077− 0.804− 0.503− 8.499< 0.001< 0.0010.5200.4470.605
HB3− 0.5360.097− 0.726− 0.345− 5.514< 0.001< 0.0010.5850.4840.708
Stress0.1560.0660.0260.2852.3500.0190.0241.9751.1681.0261.330
Anxiety0.4710.0710.3320.6106.650< 0.001< 0.0012.2691.6021.3941.840
Depression0.2950.0730.1530.4384.070< 0.001< 0.0012.4171.3441.1661.549
Age0.0210.0050.0120.0304.516< 0.001< 0.0011.3061.0211.0121.031
Sex (female)0.2160.1010.0180.4142.1350.0330.0391.0691.2411.0181.513
Conditions0.2360.0370.1630.3096.356< 0.001< 0.0011.3941.2661.1771.361

Regression coefficients of serial indirect effects Model 2 for gratitude, health behaviors, psychological distress, and sleep disturbances, with covariates.

Standardized coefficients are reported.

CI, confidence interval; LL, lower limit; UL, upper limit; Adj. p, Adjusted p-value (FDR correction); VIF, variance inflation factor; OR, odds ratio; HB, health behaviors.

0|1: sleep disturbances threshold; the likelihood of scoring a 0 on sleep disturbances as opposed to 1−3. 1|2: sleep disturbances threshold; the likelihood of scoring a 0−1 on sleep disturbances as opposed to 2−3. 2|3: sleep disturbances threshold; the likelihood of scoring 0−2 on sleep disturbances as opposed to 3. HB1: The likelihood of a one standard deviation increase in health behaviors being associated with higher sleep disturbances scores of 1−3 vs. 0. HB2: The likelihood of a one standard deviation increase in health behaviors being associated with higher sleep disturbances scores of 2−3 vs. 0−1. HB3: The likelihood of a one standard deviation increase in health behaviors being associated with a higher sleep disturbances score of 3 vs. 0−2. Covariates: Age = participant age; Sex = participant sex; Education = participant education level; Conditions = number of chronic medical conditions.

aMultiple linear regression analyses.

bOrdinal logistic regression analysis.

cVariance inflation factor for health behaviors.

Figure 3

Mediation effect results for model 2 are provided in Table 6. Despite significant regression results, the multivariate hypothesis for model 2 was only partially supported. Serial mediation was observed through health behaviors as a first-order mediator, and anxiety (a1d31b5= −0.009, 95% CI [−0.019, −0.002], SE = 0.004) and depression (a1d41b6= −0.007, 95% CI [−0.016, −0.001], SE = 0.004), but not stress, as parallel second-order mediators. Additionally, specific indirect effects linking gratitude and sleep disturbances were observed for all three health behavior pathways, health behaviors1 (a1b1= −0.067, 95% CI [−0.120, −0.026], SE = 0.024), health behaviors2 (a1b2= −0.113, 95% CI [−0.180, −0.059], SE = 0.030), health behaviors3 (a1b3= −0.092, 95% CI [−0.158, −0.040], SE = 0.032), as well as anxiety (a3b5= −0.054, 95% CI [−0.104, −0.015], SE = 0.022) and depression (a4b6= −0.052, 95% CI [−0.103, −0.011], SE = 0.024). Overall, the total indirect effect (β = −0.435, 95% CI [−0.624, −0.268], SE = 0.093) and the total effect of gratitude on sleep disturbances (c = −0.372, 95% CI [−0.607, −0.146], SE = 0.115) were significant, while the direct effect was non-significant (c' = 0.062, 95% CI [−0.090, 0.209], SE = 0.076), indicating mediation. After adjusting the confidence intervals to 0.97 to account for false coverage rate correction, the conclusions regarding all indirect effects remained consistent.

Table 6

EffectβSEβ95% CI β
LLUL
Gratitude → health behaviors1 → sleepa1 b1−0.0670.024−0.120−0.026
Gratitude → health behaviors2 → sleepa1 b2−0.1130.030−0.180−0.059
Gratitude → health behaviors3 → sleepa1 b3−0.0920.032−0.158−0.040
Gratitude → stress → sleepa2 b4−0.0340.022−0.0810.007
Gratitude → anxiety → sleepa3 b5−0.0540.022−0.104−0.015
Gratitude → depression → sleepa4 b6−0.0520.024−0.103−0.011
Gratitude → health behaviors → stress → sleepa1 d21 b4−0.0060.004−0.0140.001
Gratitude → health behaviors → anxiety → sleepa1 d31 b5−0.0090.004−0.019−0.002
Gratitude → health behaviors → depression → sleepa1 d41 b6−0.0070.004−0.016−0.001
Totalc−0.3720.115−0.607−0.146
Total indirect−0.4350.093−0.624−0.268
Directc0.0620.076−0.0900.209

Indirect effects of health behaviors and psychological distress on the relation between gratitude and sleep disturbances (Model 2).

Covariates include age, sex, education level, and number of chronic diseases. Standardized coefficients are reported.

CI, confidence interval; LL, lower limit; UL, upper limit.

Health behaviors1: Specific indirect effect of gratitude related to sleep disturbances scores of 1–3 vs. 0, via health behaviors. Health behaviors2: Specific indirect effect of gratitude related to sleep disturbances scores of 2–3 vs. 0–1, via health behaviors. Health behaviors3: Specific indirect effect of gratitude related to sleep disturbances scores of 3 vs. 0–2, via health behaviors. All effects remained consistent after adjusting the CIs for the FCR at 0.97.

3.2.3 Model 3

Regarding model 3, the ordinal logistic regression model fit the data significantly better than the multinomial model ( = 549.04, p ≤ 0.001) and was more parsimonious. However, the proportional odds assumption was violated for both health self-efficacy and health behaviors (HSE: = 6.990, p = 0.030; health behaviors: = 6.712, p = 0.035; both: = 10.397, p = 0.034). As all three models fit significantly better than the parallel model, the partial proportional odds model containing multiple slopes for both health self-efficacy and health behaviors was then compared to models containing multiple slopes for either health self-efficacy or health behaviors, to determine the most parsimonious model. The model including multiple slopes for both health self-efficacy and health behaviors did not fit the data significantly better than the health self-efficacy-only ( = 3.407, p = 0.182) or health behaviors-only ( = 3.685, p = 0.158) models. Next, the health self-efficacy-only and health behaviors-only models were compared. The health self-efficacy-only model fit better than the health behaviors-only model ( = 0.278, p < 0.001); thus, the health self-efficacy-only partial proportional odds model was used.

Linear and ordinal logistic regression results and fit statistics for each pathway, including significant covariates identified through backward stepwise selection, are provided in Table 7. Variance inflation factors are also provided in Table 7, and no predictors exceed the multicollinearity cutoff. Regression results excluding covariate effects are illustrated in Figure 4. Gratitude was related to better health self-efficacy (β = 0.273, 95% CI [0.209, 0.337], SE = 0.032, t = 8.412, padjusted < 0.001) which was, in turn, related to lower stress (β = −0.276, 95% CI [−0.350, −0.203], SE = 0.038, t = −7.352, padjusted < 0.001), anxiety (β = −0.202, 95% CI [−0.280, −0.123], SE = 0.040, t = −5.064, padjusted < 0.001), and depression (β = −0.230, 95% CI [−0.303, −0.157], SE = 0.037, t = −6.169, padjusted < 0.001). Gratitude was also related to better health behavior engagement (β = 0.173, 95% CI [0.106, 0.239], SE = 0.034, t = 5.081, padjusted < 0.001) which was, in turn, related to lower stress (β = −0.107, 95% CI [−0.177, −0.037], SE = 0.036, t = −3.001, padjusted = 0.005), but not anxiety and depression. A one standard deviation increase in stress was associated with a 17% greater likelihood of endorsing a higher frequency of sleep disturbances (OR = 1.168, 95% CI [1.023, 1.333], padjusted = 0.038), while equivalent increases in anxiety (OR = 1.602, 95% CI [1.394, 1.840], padjusted < 0.001) and depression (OR = 1.343, 95% CI [1.163, 1.551], padjusted < 0.001) were associated with 60% and 34% greater likelihood of sleep problem frequency, respectively.

Table 7

EffectaβSEβ95% CI βtpAdj. pVIFFit statistics
LLUL
HSER2 = 0.201;R2adjusted = 0.197F(4,807) = 50.74,p < 0.001
(Intercept)− 0.3610.179− 0.713− 0.010− 2.0180.0440.065
Gratitude0.2730.0320.2090.3378.412< 0.001< 0.0011.052
Age0.0080.0030.0030.0142.9380.0030.0061.158
Education0.0770.0230.0310.1233.3000.0010.0021.074
Conditions− 0.1960.023− 0.242− 0.151− 8.488< 0.001< 0.0011.181
HBR2 = 0.126;R2adjusted = 0.122F(4,817) = 29.38,p < 0.001
(Intercept)− 1.4810.187− 1.848− 1.114− 7.922< 0.001< 0.001
Gratitude0.1730.0340.1060.2395.081< 0.001< 0.0011.051
Age0.0260.0030.0210.0328.970< 0.001< 0.0011.158
Education0.0740.0250.0250.1222.9950.0030.0051.075
Conditions− 0.0970.024− 0.144− 0.050− 4.028< 0.001< 0.0011.179
StressR2 = 0.267;R2adjusted = 0.261F(6,808) = 49.01,p < 0.001
(Intercept)0.2120.157 −0.0970.5221.3490.1780.199
HSE− 0.2760.038− 0.350− 0.203− 7.352< 0.001< 0.0011.515
HB− 0.1070.036− 0.177− 0.037− 3.0010.0030.0051.369
Gratitude− 0.1560.033− 0.220− 0.092− 4.779< 0.001< 0.0011.146
Age− 0.0100.003− 0.016− 0.005− 3.611< 0.001< 0.0011.269
Sex (female)0.1180.063− 0.0060.2431.8730.0610.0861.041
Conditions0.1470.0230.1010.1926.329< 0.001< 0.0011.303
AnxietyR2 = 0.171;R2adjusted = 0.164F(7,795) = 23.4, p < 0.001
(Intercept)− 0.1850.202− 0.5810.212 −0.9140.3610.386
HSE− 0.2020.040− 0.280− 0.123− 5.064< 0.001< 0.0011.518
HB− 0.0560.038− 0.1300.018− 1.4810.1390.1601.389
Gratitude− 0.0620.035− 0.1300.006− 1.7800.0750.1021.146
Age− 0.0090.003− 0.015− 0.003− 2.8830.0040.0071.300
Sex (female)0.2780.0670.1470.4104.148< 0.001< 0.0011.042
Education0.0480.0240.0010.0961.9990.0460.0661.096
Conditions0.1710.0250.1230.2206.966< 0.001< 0.0011.292
DepressionR2 = 0.256;R2adjusted = 0.251F(6,813) = 46.72,p < 0.001
(Intercept)0.2410.155− 0.0630.5461.5550.1200.142
HSE− 0.2300.037− 0.303− 0.157− 6.169< 0.001< 0.0011.512
HB− 0.0600.035− 0.1290.009− 1.7200.0860.1131.365
Gratitude− 0.1210.032− 0.185− 0.058− 3.756< 0.001< 0.0011.146
Age− 0.0130.003− 0.019− 0.008− 4.756< 0.001< 0.0011.259
Sex (female)0.1570.0630.0340.2802.5100.0120.0191.038
Conditions0.2060.0230.1610.2519.006< 0.001< 0.0011.300
EffectbβSEβ95% CI βzpAdj. pVIFOR95% CI OR
LLULLLUL
SleepLL = -879.75, df = 2,386
0|1− 0.1820.258− 0.6870.324− 0.7040.4820.5190.8340.5031.383
1|2− 2.6600.266− 3.181− 2.139− 10.006< 0.001< 0.0010.0700.0420.118
2|3− 3.8370.281− 4.388− 3.286− 13.653< 0.001< 0.0010.0220.0120.037
Gratitude0.0760.054− 0.0290.1821.4130.1580.2001.1901.0790.9711.199
HSE10.0740.094− 0.1110.2590.7860.4320.5041.633c1.0770.8951.296
HSE2− 0.1830.086− 0.351− 0.015− 2.1400.0320.0450.8320.7040.985
HSE3− 0.0170.109− 0.2300.196− 0.1570.8750.8750.9830.7941.217
HB− 0.4960.057− 0.608− 0.384− 8.684< 0.001< 0.0011.3770.6090.5450.681
Stress0.1550.0680.0230.2882.2950.0220.0382.0441.1681.0231.333
Anxiety0.4710.0710.3320.6106.638< 0.001< 0.0012.2811.6021.3941.840
Depression0.2950.0740.1510.4394.013< 0.001< 0.0012.4351.3431.1631.551
Age0.0200.0050.0110.0304.344< 0.001< 0.0011.3101.0211.0111.030
Sex (female)0.2260.1010.0270.4252.2270.0260.0401.0651.2531.0271.529
Conditions0.2070.0380.1320.2825.406< 0.001< 0.0011.4481.2301.1411.325

Regression coefficients of serial indirect effects Model 3 for gratitude, health self-efficacy, health behaviors, psychological distress, and sleep disturbances, with covariates.

Standardized coefficients are reported.

CI, confidence interval; LL, lower limit; UL, upper limit. Adj. p, Adjusted p-value (FDR correction); VIF, variance inflation factor; OR, odds ratio; HSE, health self-efficacy; HB, health behaviors.

Bolded p-values indicate those that became non-significant after adjustment. 0|1: sleep disturbances threshold; the likelihood of scoring a 0 on sleep disturbances as opposed to 1−3. 1|2: sleep disturbances threshold; the likelihood of scoring a 0−1 on sleep disturbances as opposed to 2−3. 2|3: sleep disturbances threshold; the likelihood of scoring 0−2 on sleep disturbances as opposed to 3. HSE1: The likelihood of a one standard deviation increase in health self-efficacy being associated with higher sleep disturbances scores of 1−3 vs. 0. HSE2: The likelihood of a one standard deviation increase in health self-efficacy being associated with higher sleep disturbances scores of 2−3 vs. 0−1. HSE3: The likelihood of a one standard deviation increase in health self-efficacy being associated with a higher sleep disturbances score of 3 vs. 0−2.

Covariates: Age = participant age; Sex = participant sex; Education = participant education level; Conditions = number of chronic medical conditions.

aMultiple linear regression analyses.

bOrdinal logistic regression analysis.

cVariance inflation factor for health self-efficacy.

Figure 4

Mediation effects results for model 3 are provided in Table 8. The multivariate hypothesis for model 3 was only partially supported. No serial mediation was found when health behaviors was included as a first-order mediator and, as a result, health self-efficacy and health behaviors did not function as parallel first-order mediators. Only health self-efficacy and anxiety (a1d41b6 = −0.026, 95% CI [−0.046, −0.009], SE = 0.009), and health self-efficacy and depression (a1d51b7 = −0.019, 95% CI [−0.037, −0.003], SE = 0.009), serially mediated the relation between gratitude and sleep disturbances. Specific indirect effects linking gratitude and sleep disturbances were observed for health behaviors (a2b4 = −0.086, 95% CI [−0.140, −0.030], SE = 0.026) and depression (a5b7= −0.036, 95% CI [−0.079, −0.005], SE = 0.019). Despite serial mediation occurring in only two of six potential paths, the sum of all indirect effects of gratitude on sleep disturbances via health self-efficacy, health behaviors, and psychological distress was significant (β = −0.276, 95% CI [−0.451, −0.053], SE = 0.093), as was the total effect (c = −0.200, 95% CI [−0.412, −0.001], SE = 0.105). The direct effect was non-significant (c' = 0.076, 95% CI [−0.078, 0.227], SE = 0.079), indicating mediation. After adjusting the confidence intervals to 0.97 to account for false coverage rate correction, the conclusions regarding all indirect effects remained consistent. However, the total effect of gratitude on sleep became non-significant after adjustment, indicating a suppression effect in which the positive direct and negative indirect effects cancel each other out (MacKinnon, 2000), and likely reflecting the lack of serial indirect effects for health behaviors.

Table 8

EffectβSEβ95% CI β
LLUL
Gratitude → self-efficacy1 → sleepa1 b10.0200.030−0.0370.079
Gratitude → self-efficacy2 → sleepa1 b2−0.0500.030−0.1130.009
Gratitude → self-efficacy3 → sleepa1 b3−0.0050.041−0.0820.081
Gratitude → health behaviors → sleepa2 b4−0.0860.026−0.140−0.030
Gratitude → stress → sleepa3 b5−0.0240.017−0.0610.006
Gratitude → anxiety → sleepa4 b6−0.0290.020−0.0730.006
Gratitude → depression → sleepa5 b7−0.0360.019−0.079−0.005
Gratitude → self-efficacy → stress → sleepa1 d31 b5−0.0120.008−0.0290.003
Gratitude → self-efficacy → anxiety → sleepa1 d41 b6−0.0260.009−0.046−0.009
Gratitude → self-efficacy → depression → sleepa1 d51 b7−0.0190.009−0.037−0.003
Gratitude → health behaviors → stress → sleepa2 d32 b5−0.0030.002−0.0080.001
Gratitude → health behaviors → anxiety → sleepa2 d42 b6−0.0050.004−0.0130.002
Gratitude → health behaviors → depression → sleepa2 d52 b7−0.0030.003−0.0090.001
Totalc−0.2000.105−0.4120.001
Total indirect−0.2760.093−0.451−0.053
Directc0.0760.079−0.0780.227

Indirect effects of health self-efficacy, health behaviors, and psychological distress on the relation between gratitude and sleep disturbances (Model 3).

Covariates include age, sex, education level, and number of chronic diseases. Standardized coefficients are reported.

CI, confidence interval; LL, lower limit; UL, upper limit; Self-efficacy, health self-efficacy.

Self-efficacy1: Specific indirect effect of gratitude related to sleep disturbances scores of 1–3 vs. 0, via health self-efficacy. Self-efficacy2: Specific indirect effect of gratitude related to sleep disturbances scores of 2–3 vs. 0-1, via health self-efficacy. Self-efficacy3: Specific indirect effect of gratitude related to sleep disturbances scores of 3 vs. 0–2, via health self-efficacy. All effects remained consistent after adjusting the CIs for the FCR at 0.97 except for bolded CI, which became positive, indicating a non-significant total effect.

Finally, the results of our secondary SEM analyses provided additional support for our models, with modified models (i.e., dichotomized outcome, creation of psychological distress latent variable) exhibiting good fit and comparable patterns of direct and indirect effects. Please refer to Supplementary material for detailed results.

4 Discussion

In a sample of German primary care patients, we examined the potential serial mediating effects of health self-efficacy (HSE) and health behaviors as first-order mediators, and stress, anxiety, and depression as second-order mediators, on the relation between gratitude and sleep disturbances. Consistent with our theoretical model, in model 1, the serial mediation hypothesis was supported; gratitude was related to greater health self-efficacy and, in turn, to less stress, anxiety, and depression in parallel, with subsequent diminished sleep disturbances. However, in model 2, although gratitude was related to health behaviors, and health behaviors were related to stress, anxiety, and depression, the hypothesis was only partially supported, in that serial mediation was observed via health behaviors and anxiety and depression, but not via health behaviors and stress. In model 3, health self-efficacy and health behaviors did not function as parallel first-order mediators as hypothesized. Instead, serial mediation emerged in only two of six potential paths, through health self-efficacy and anxiety, and health self-efficacy and depression.

Support for our theoretical model arises from existing research outlining the benefits of gratitude for health self-efficacy and health behaviors. Gratitude is a pleasant emotion that, when evoked, facilitates the development of additional positive emotions (i.e., positive affect hypothesis; Wood et al., 2010), and the building of resources to be called upon during stressful times (i.e., broaden-and-build theory, Fredrickson, 2004; Wood et al., 2010), thereby contributing to health and wellbeing. Grateful individuals employ health-promoting behaviors as a positive coping strategy and are less likely to engage in maladaptive health behaviors, such as substance use, according to the coping hypothesis (Wood et al., 2007). Further, according to the cognitive framework, gratitude stimulates a positive perception of one's abilities, strengths, and circumstances, engagement in adaptive emotion regulation strategies, and a tendency to focus on and recall encouraging experiences which, in turn, enhances resilience against stress and improves psychological and physical wellbeing (Alkozei et al., 2018). Bandura's (1989) social cognitive theory describes self-efficacy as the belief in one's ability to control life circumstances, and health self-efficacy specifically refers to belief in ability to perform necessary behaviors to attain health goals (Sirois, 2007). Self-efficacy influences health behaviors through goals, perceptions of barriers and facilitators, and outcome expectations (Bandura, 2004). Experience with prior successes helps bolster belief in capability to exert control despite setbacks (Bandura, 1977, 2004); thus, gratitude's tendency to promote a positive perspective of capabilities and experiences may improve sense of control over health (Swain et al., 2020).

Gratitude can foster self-efficacy and mitigate psychological distress and sleep problems, consistent with the significant serial indirect effects findings of models 1 and 3. For example, among first-year U.S. undergraduates with high, but not low, gratitude levels, self-efficacy for COVID-19 safety predicted less depression after 2 months (Ang et al., 2022). Further, gratitude interventions increased self-efficacy and reduced perceptions of stress in German young adults (Lorenz et al., 2022), improved pain self-efficacy and pain-related anxiety in individuals with arthritis (Swain et al., 2020), and improved environmental mastery (i.e., perceived ability to effectively change surrounding circumstances and events; Ryff and Keyes, 1995), perceived stress, anxiety, depression, and insomnia in adults (Czyżowska and Gurba, 2022). Finally, in patients with asymptomatic heart failure, gratitude was related to greater self-efficacy to preserve heart function, less depression, better sleep, and lower levels of inflammatory markers (Mills et al., 2015).

Additionally, consistent with model 2's findings, engagement in adaptive health behaviors can positively impact psychological health and sleep (Su et al., 2021), and numerous studies offer support for gratitude's salubrious influence on these associations. In a global ecological momentary assessment study of adults (N = 4,825), greater gratitude was associated with increased exercise, lower blood pressure, heart rate, and stress, more daily positive expectations and reflections, and better sleep quality (Newman et al., 2021). Similarly, after a gratitude diary intervention, women reporting emotional distress and moderate sleep disturbances experienced more positive emotions, reduced blood pressure, and better sleep quality (Jackowska et al., 2016). In a large online survey of 59,985 respondents from ~160 countries, gratitude was associated with increased likelihood of engaging in healthy eating, social activities, self-care, and exercise, and reduced risk for smoking, depression, anxiety, poor physical health, and insufficient sleep (Weziak-Bialowolska et al., 2023). Finally, in a large multinational online survey of pregnant women, participants endorsing greater gratitude reported less detrimental impact of COVID-19 on diet, fitness, and sleep, and did not have anxiety or depression (Choi et al., 2022).

Importantly, it may be necessary to improve self-efficacy to affect positive changes in health behaviors, and health self-efficacy should be investigated as an antecedent to health behaviors in our model, rather than a parallel mechanism. There is some precedent for this assertion, which could explain the lack of parallel and serial indirect effects for health behaviors and substantiation of serial indirect effects for health self-efficacy in model 3. Rosenstock et al. (1988) revised the health belief model to incorporate self-efficacy, along with perceived benefits and barriers, perceived vulnerability to a severe health threat, and sufficient motivation or health concern, as important explanatory factors in predicting health-related behavior among individuals with chronic illnesses (Rosenstock, 1966). As an example, in a sample of Ukrainian primary care patients with both chronic kidney disease and hypertension, self-efficacy accounted for 76% of the variance in quality of life associated with health behaviors (i.e., physical activity, diet, stress management, smoking, medication adherence, blood pressure control; Korzh et al., 2022). Further, in German primary care patients with heart failure, low self-efficacy predicted both poor lifestyle and medication adherence (Eisele et al., 2020).

Given that affective traits are known to influence health behaviors, the self-regulation resource model posits that affective (i.e., high positive and low negative affect) and social-cognitive variables (i.e., self-efficacy) work in conjunction to bolster the self-regulation capacity needed to perform health-promoting behaviors (Sirois, 2015). Gratitude is one possible positive affective trait that can be drawn upon, alongside health self-efficacy, to facilitate behavior change (Klibert et al., 2019; Wood et al., 2010). Self-efficacy mediated the effect of gratitude on medication adherence in heart failure patients (Cousin et al., 2020), and, in a systematic review, gratitude for the medical team and patient's donor, and self-efficacy, emerged as the primary factors contributing to medical self-management among kidney transplant recipients (Jamieson et al., 2016).

Multiple patterns emerged across models that warrant further discussion. In model 2, health behaviors were an independent mediator at all severity levels of sleep disturbances whereas, in model 3, health behaviors were not a significant mediator in serial pathways but emerged as an independent mediator of the relation between gratitude and sleep disturbances. This pattern of findings indicates that there is a strong independent influence of health behaviors on sleep. Engagement in both positive and negative health behaviors can have direct effects on sleep. For example, caffeine use can increase arousal and disrupt sleep onset, duration, and efficiency, and quality of slow-wave sleep (Clark and Landolt, 2017). Similarly, cocaine and nicotine intoxication can impair REM sleep, sleep onset time, and sleep duration (Garcia and Salloum, 2015). On the other hand, physical exercise and weight control can directly improve sleep quality, efficiency, duration, and onset latency (Alfaris et al., 2015; Rogers et al., 2017). Other health behaviors, such as a healthy diet, can help to regulate circadian rhythms through nutrients that affect the release of hormones, such as melatonin, growth hormone, and serotonin (Vernia et al., 2021).

The final notable finding, across models, was that, although gratitude and the first-order mediators (i.e., M1; health self-efficacy and/or health behaviors) predicted greater reductions in stress than anxiety and depression, the effect of stress on sleep disturbances was weaker than the effects of anxiety and depression. As well, serial mediation for the stress pathway occurred only in model 1, which was focused on health self-efficacy. Thus, although the broaden-and-build theory (Fredrickson, 2004) indicates that gratitude, health self-efficacy, and health behavior engagement may lead to stress reduction, these salubrious effects are not transmitted to sleep, as with anxiety and depression, unless health behaviors are excluded from the model and only health self-efficacy is considered. Because biological stress is considered a significant contributor to the development of anxiety and depression (Ross et al., 2017), physiological stress and, as a result, perceived stress, may exacerbate anxiety and depression, rather than escalating concurrently.

Overall, our findings indicate that gratitude has a greater effect on psychological distress through health self-efficacy than through health behaviors, and that health self-efficacy is the primary mechanism through which gratitude predicts reductions in psychological distress and downstream sleep problems.

4.1 Limitations and directions for future research

It is important that our findings are interpreted within the context of several study limitations. To begin, this study's cross-sectional design precludes determinations of causality (O'Laughlin et al., 2018) and, thus, bidirectionality of variables is a possibility. First, sleep disturbances and depression are known to have a bidirectional association in primary care patients (Bouwmans et al., 2017). This may occur because, in the context of a maladaptive stress response, a persistent feedback loop of anxiety and depressive symptoms, and sleep problems, can manifest (Gold, 2015). Additionally, sleep is considered a health behavior that, independently and through its impact on other health behaviors, such as diet and exercise, contributes to the development and progression of depression (Dzierzewski et al., 2014; Lopresti et al., 2013).

In adults, anxiety and depressive disorders predict unfavorable health behaviors, such as smoking, low physical activity, and poor diet (Difrancesco et al., 2022; Gall et al., 2016; Walsh et al., 2023), and self-efficacy and smoking behavior may be bidirectional (Clyde et al., 2019). In sum, given the current study's cross-sectional, survey-based design, future randomized controlled trials in which study variables are manipulated and prospective, longitudinal studies may help to elucidate these linkages.

Further, the use of self-report measures may limit external validity, as responses are affected by participants' understanding of the items, level of personal insight, and potential social desirability and response biases (Demetriou et al., 2015). Future studies should employ objective or physiological assessments, such as medical record reviews, actigraphy, polysomnography, and ecological momentary assessment, to improve validity (Jackson et al., 2018; Prince et al., 2020; Scarlett et al., 2021). Relatedly, the outcome variable of sleep disturbances was measured using only one item, item 3 of the PHQ-9, which precludes estimates of internal consistency and may not adequately capture a construct as multifaceted as sleep disturbances (Allen et al., 2022). However, this item has been validated in primary care patients as a brief screening tool for sleep disturbances (MacGregor et al., 2012). In future studies, utilizing a multidimensional measure, such as the Pittsburgh Sleep Quality Index (Buysse et al., 1989), may enhance reliability and validity and allow for the use of SEM or ordinary least squares regression at each stage of the model, simplifying interpretation of the regression results and potentially reducing measurement error.

Notably, as most participants (58%) reported zero or one medical condition, the results may not fully capture the experiences of individuals with multimorbidity. Nearly 40% of German adults experience multimorbidity (i.e., have two or more chronic medical conditions simultaneously), including half of adults between 50 and 59 years old, and 47% of adults between 40 and 49 years of age with low education levels (Puth et al., 2017). More chronic conditions predict development of stress, anxiety, and depression in representative population-based samples (Hajek and König, 2020; Liu et al., 2021), and sleep disturbances in primary care patients (Ullmann et al., 2022). Number of chronic diseases was covaried in all analyses and emerged as a significant predictor on almost all pathways of every model.

Finally, race and ethnicity statistics are not collected in German studies due to legal prohibitions and political factors, thus precluding inferences of generalizability for these demographic factors. Yet, it is likely that disparities exist, based on immigration research findings indicating that first-generation migrants in Germany have significantly higher rates of depression and generalized anxiety (Beutel et al., 2016). Future research conducted with additional demographic groups, chronic disease populations (e.g., cardiovascular, rheumatological, neurodegenerative), and sleep disorder populations (e.g., insomnia, sleep apnea, narcolepsy) is needed to improve generalizability and substantiate our findings across diverse groups.

4.2 Implications

Despite limitations, this study's novel findings may help to inform the selection of targeted clinical interventions to improve sleep in primary care patients. To begin, self-management interventions are effective for improving gratitude, health behaviors, self-efficacy, psychological distress, and sleep, including in individuals living with cancer (Martin et al., 2020), sedentary European adults with insomnia (Hartescu et al., 2015), and primary care patients from 10 practices in Germany (Zimmermann et al., 2016) and enrolled in a global cardiovascular risk reduction program (Prince et al., 2017). Importantly, according to a systematic review, to promote health self-efficacy, successful self-management programs in primary care should be tailored to meet individual needs, and should include independent symptom monitoring with self-treatment, regular communication with one's provider, stress management and coping strategies for addressing psychological distress, and promotion of responsibility for health behavior choices (Dineen-Griffin et al., 2019).

Further, systematic reviews and meta-analyses indicate beneficial effects of gratitude interventions for self-efficacy (Grant and Gino, 2010), stress, anxiety, and depression (Carr et al., 2021; Cregg and Cheavens, 2021), and subjective sleep quality (Boggiss et al., 2020). In a systematic review of positive psychological interventions (PPIs) in which 16 of 27 studies included gratitude-specific exercises, PPIs were responsible for favorable effects on health behaviors (Feig et al., 2022). Additional evidence-based interventions that may improve self-efficacy, health behaviors, psychological distress, and sleep disturbance include cognitive behavioral therapy (e.g., changing maladaptive thoughts to alter emotional responses) and acceptance and commitment therapy (e.g., fostering present-moment awareness, acceptance, and values-driven behavior), which can be efficaciously implemented in primary care settings (Arnold et al., 2022; Benninghoven et al., 2022; Hopko et al., 2009; Jonkers et al., 2012; Siengsukon et al., 2020). Behavioral activation for depression (i.e., identifying and intentionally scheduling enjoyable activities), a component of cognitive behavioral therapy, is particularly well-suited to cultivating engagement in salubrious behaviors, such as reducing problem drinking, facilitating smoking cessation, and improving physical activity and depression (May et al., 2024). In sum, as psychological distress and sleep disturbances are most often identified in primary care populations, providers have an opportunity to facilitate the initial catalyst of advantageous downstream effects that begins with promotion of gratitude.

4.3 Conclusion

In a sample of German primary care patients, health self-efficacy as a first-order mediator, and stress, anxiety, and depression, as parallel second-order mediators, serially mediated the relation between gratitude and sleep disturbances. In addition, health behaviors, and anxiety and depressive symptoms, functioned as serial mediators of the gratitude-sleep linkage. Finally, we found that health self-efficacy and health behaviors were not parallel first-order mediators, and only health self-efficacy and anxiety, and health self-efficacy and depression, serially mediated the relation between gratitude and sleep disturbances. Future prospective research utilizing objective assessments and across diverse samples is needed to substantiate these findings. However, despite limitations, this study provides preliminary support for the role of gratitude as a catalyst for health self-efficacy and health behaviors and, in turn, for its beneficial influence on psychological distress and sleep disturbances.

Statements

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving humans were approved by Ethics Commission of the Faculty of Medicine at the Ludwig Maximilian University of Munich. The studies were conducted in accordance with the local legislation and institutional requirements. The Ethics Committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants' legal guardians/next of kin because participation was voluntary and, to ensure privacy, no personally identifiable information was obtained.

Author contributions

HA: Conceptualization, Data curation, Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing. JH: Conceptualization, Data curation, Methodology, Supervision, Validation, Writing – review & editing. AW: Conceptualization, Investigation, Project administration, Writing – review & editing, Data curation. NK: Conceptualization, Investigation, Project administration, Resources, Supervision, Writing – review & editing, Methodology. JS: Conceptualization, Investigation, Project administration, Resources, Supervision, Writing – review & editing. LT: Conceptualization, Methodology, Writing – review & editing. FS: Conceptualization, Methodology, Writing – review & editing. MO: Conceptualization, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frsle.2025.1459854/full#supplementary-material

References

  • 1

    ÅkerstedtT.NarusyteJ.AlexandersonK.SvedbergP. (2017). Sleep duration, mortality, and heredity—A prospective twin study. Sleep. 40:zsx135. 10.1093/sleep/zsx135

  • 2

    AernoutE.BenradiaI.HazoJ.-B.SyA.Askevis-LeherpeuxF.SebbaneD.et al. (2021). International study of the prevalence and factors associated with insomnia in the general population. Sleep Med.82, 186192. 10.1016/j.sleep.2021.03.028

  • 3

    AgrestiA. (1996). An Introduction to Categorical Data Analysis. New York, NY: Wiley.

  • 4

    AkinwandeM. O.DikkoH. G.SamsonA. (2015). Variance inflation factor: As a condition for the inclusion of suppressor variable(s) in regression analysis. OJS5, 754767. 10.4236/ojs.2015.57075

  • 5

    AlfarisN.WaddenT. A.SarwerD. B.DiwaldL.VolgerS.HongP.et al. (2015). Effects of a 2-year behavioral weight loss intervention on sleep and mood in obese individuals treated in primary care practice. Obesity23, 558564. 10.1002/oby.20996

  • 6

    AlkozeiA.SmithR.KillgoreW. D. S. (2018). Gratitude and subjective wellbeing: a proposal of two causal frameworks. J. Happiness Stud.19, 15191542. 10.1007/s10902-017-9870-1

  • 7

    AlkozeiA.SmithR.KotzinM. D.WaugamanD. L.KillgoreW. D. S. (2019). The association between trait gratitude and self-reported sleep quality is mediated by depressive mood state. Behav. Sleep Med.17, 4148. 10.1080/15402002.2016.1276017

  • 8

    AllenM. S.IliescuD.GreiffS. (2022). Single item measures in psychological science: a call to action. Eur. J. Psychol. Assess.38, 15. 10.1027/1015-5759/a000699

  • 9

    AngJ. Y.-Z.MonteV.TsaiW. (2022). First-year college students' adjustment during the COVID-19 pandemic: the protective roles of hope and gratitude. Transl. Issues Psychol. Sci.8, 375388. 10.1037/tps0000320

  • 10

    ArnoldT.HaubrickK. K.Klasko-FosterL. B.RogersB. G.BarnettA.Ramirez- SanchezN. A.et al. (2022). Acceptance and commitment therapy informed behavioral health interventions delivered by non-mental health professionals: a systematic review. J. Context. Behav. Sci.24, 185196. 10.1016/j.jcbs.2022.05.005

  • 11

    ArrollB.FernandoA.FalloonK.Goodyear-SmithF.SamaranayakeC.WarmanG. (2012). Prevalence of causes of insomnia in primary care: a cross-sectional study. Br. J. Gen. Pract.62, e99e103. 10.3399/bjgp12X625157

  • 12

    BanduraA. (1977). Self-efficacy: toward a unifying theory of behavioral change. Psychol. Rev.84, 191215. 10.1037//0033-295X.84.2.191

  • 13

    BanduraA. (1989). “Social cognitive theory,” in Annals of Child Development, ed. R. Vasta (Greenwich, CT: JAI Press), 160.

  • 14

    BanduraA. (2004). Health promotion by social cognitive means. Health Educ. Behav.31, 143164. 10.1177/1090198104263660

  • 15

    BenjaminiY.HochbergY. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B Stat. Methodol.57, 289300. 10.1111/j.2517-6161.1995.tb02031.x

  • 16

    BenjaminiY.YekutieliD. (2005). False discovery rate–adjusted multiple confidence intervals for selected parameters. J. Am. Stat. Assoc.100, 7181. 10.1198/016214504000001907

  • 17

    BenninghovenD.MenkeE.ChinaC.SchroederF.BethgeM. (2022). Implementation of a behavioral medicine-oriented concept in cardiological rehabilitation. Psychother. Psychosom. Med. Psychol.72, 429437. 10.1055/a-1749-6379

  • 18

    BeutelM. E.JüngerC.KleinE. M.WildP.LacknerK. J.BlettnerM.et al. (2016). Depression, anxiety and suicidal ideation among 1st and 2nd generation migrants – Results from the Gutenberg health study. BMC Psychiatry16:288. 10.1186/s12888-016-0995-2

  • 19

    BluesteinD.RutledgeC. M.HealeyA. C. (2010). Psychosocial correlates of insomnia severity in primary care. J. Am. Board Fam. Med.23, 204211. 10.3122/jabfm.2010.02.090179

  • 20

    BoggeroI. A.HostinarC. E.HaakE. A.MurphyM. L. M.SegerstromS. C. (2017). Psychosocial functioning and the cortisol awakening response: meta-analysis, P-curve analysis, and evaluation of the evidential value in existing studies. Biol. Psychol.129, 207230. 10.1016/j.biopsycho.2017.08.058

  • 21

    BoggissA. L.ConsedineN. S.Brenton-PetersJ. M.HofmanP. L.SerlachiusA. S. (2020). A systematic review of gratitude interventions: effects on physical health and health behaviors. J. Psychosom. Res.135:110165. 10.1016/j.jpsychores.2020.110165

  • 22

    BouwmansM. E. J.ConradiH. J.BosE. H.OldehinkelA. J.De JongeP. (2017). Bidirectionality between sleep symptoms and core depressive symptoms and their long-term course in major depression. Psychosom. Med.79, 336344. 10.1097/PSY.0000000000000407

  • 23

    BuysseD. J.ReynoldsC. F.MonkT. H.BermanS. R.KupferD. J. (1989). The Pittsburgh sleep quality index: a new instrument for psychiatric practice and research. Psychiatry Res.28, 193213. 10.1016/0165-1781(89)90047-4

  • 24

    CardN. A. (2019). Meta-analyses of the reliabilities of four measures of gratitude. J. Posit. Psychol.14, 576586. 10.1080/17439760.2018.1497690

  • 25

    CarrA.CullenK.KeeneyC.CanningC.MooneyO.ChinseallaighE.et al. (2021). Effectiveness of positive psychology interventions: a systematic review and meta-analysis. J. Posit. Psychol.16, 749769. 10.1080/17439760.2020.1818807

  • 26

    ChattuV. K.Md ManzarD.KumaryS.BurmanD.SpenceD. W.Pandi-PerumalS. R. (2018). The global problem of insufficient sleep and its serious public health implications. Healthcare7:1. 10.3390/healthcare7010001

  • 27

    ChoiK. W.KimH. H.BasuA.KwongA. S. F.Hernandez-DiazS.WyszynskiD. F.et al. (2022). COVID-19 perceived impacts on sleep, fitness, and diet and associations with mental health during pregnancy: a cross-national study. J. Affect. Disord. Rep.7:100288. 10.1016/j.jadr.2021.100288

  • 28

    ChristensenR. H. B. (2019). Ordinal—Regression Models for Ordinal Data, Software.

  • 29

    ClarkI.LandoltH. P. (2017). Coffee, caffeine, and sleep: a systematic review of epidemiological studies and randomized controlled trials. Sleep Med. Rev.31, 7078. 10.1016/j.smrv.2016.01.006

  • 30

    ClydeM.PipeA.ReidR.ElsC.TullochH. (2019). A bidirectional path analysis model of smoking cessation self-efficacy and concurrent smoking status: impact on abstinence outcomes. Addict. Biol.24, 10341043. 10.1111/adb.12647

  • 31

    CohenS.KamarckT.MermelsteinR. (1983). A global measure of perceived stress. J. Health Soc. Behav.24:385. 10.2307/2136404

  • 32

    CohenS.WilliamsonG. (1988). “Perceived stress in a probability sample of the United States,” in The Social Psychology of Health, eds. S. Spacapan and S. Oskamp (Newbury Park, CA: Sage Publications, Inc.), 3167.

  • 33

    CousinL.BuckH.BenitezB.MillsP.RedwineL. (2020). A structural equation model of gratitude, self-efficacy, and medication adherence in patients with stage B heart failure. J. Cardiovasc. Nurs.35, E18E24. 10.1097/JCN.0000000000000721

  • 34

    CousinL.RedwineL.BrickerC.KipK.BuckH. (2021). Effect of gratitude on cardiovascular health outcomes: a state-of-the-science review. J. Posit. Psychol.16, 348355. 10.1080/17439760.2020.1716054

  • 35

    CreggD. R.CheavensJ. S. (2021). Gratitude interventions: effective self-help? A meta-analysis of the impact on symptoms of depression and anxiety. J. Happiness Stud.22, 413445. 10.1007/s10902-020-00236-6

  • 36

    CzyżowskaN.GurbaE. (2022). Enhancing meaning in life and psychological well-being among a European cohort of young adults via a gratitude intervention. Front. Psychol.12:751081. 10.3389/fpsyg.2021.751081

  • 37

    DemetriouC.OzerB. U.EssauC. A. (2015). “Self-report questionnaires,” in The Encyclopedia of Clinical Psychology, eds. R. L. Cautin and S. O. Lilienfeld (Wiley), 16. 10.1002/9781118625392.wbecp507

  • 38

    DemyttenaereK.HeirmanE. (2020). The blurred line between anxiety and depression: hesitations on comorbidity, thresholds and hierarchy. Int. Rev. Psychiatry32, 455465. 10.1080/09540261.2020.1764509

  • 39

    DifrancescoS.PenninxB. W. J. H.MerikangasK. R.Van HemertA. M.RieseH.LamersF. (2022). Within-day bidirectional associations between physical activity and affect: a real-time ambulatory study in persons with and without depressive and anxiety disorders. Depres. Anxiety39, 922931. 10.1002/da.23298

  • 40

    Dineen-GriffinS.Garcia-CardenasV.WilliamsK.BenrimojS. I. (2019). Helping patients help themselves: a systematic review of self-management support strategies in primary health care practice. PLoS ONE14:e0220116. 10.1371/journal.pone.0220116

  • 41

    DuBoisC. M.MillsteinR. A.CelanoC. M.WexlerD. J.HuffmanJ. C. (2016). Feasibility and acceptability of a positive psychological intervention for patients with type 2 diabetes. Prim. Care Companion CNS Disord.18:10.4088/PCC.15m01902. 10.4088/PCC.15m01902

  • 42

    DzierzewskiJ. M.BumanM. P.GiacobbiP. R.RobertsB. L.Aiken-MorganA. T.MarsiskeM.et al. (2014). Exercise and sleep in community-dwelling older adults: evidence for a reciprocal relationship. J. Sleep Res.23, 6168. 10.1111/jsr.12078

  • 43

    EiseleM.HarderM.RakebrandtA.BoczorS.MarxG.BlozikE.et al. (2020). Association of depression and anxiety with adherence in primary care patients with heart failure—cross-sectional results of the observational RECODE-HF cohort study. Fam. Pract.37, 695702. 10.1093/fampra/cmaa042

  • 44

    FarawayJ. J. (2015). Linear Models with R. 2nd Edn. Boca Raton, FL: CRC Press.

  • 45

    FeigE. H.MadvaE. N.MillsteinR. A.ZambranoJ.AmonooH. L.LongleyR. M.et al. (2022). Can positive psychological interventions improve health behaviors? A systematic review of the literature. Prev. Med.163:107214. 10.1016/j.ypmed.2022.107214

  • 46

    FieldA. P. (2013). Discovering Statistics Using SPSS. 4th Edn. Los Angeles, CA: Sage.

  • 47

    FordE. S.WheatonA. G.CunninghamT. J.GilesW. H.ChapmanD. P.CroftJ. B. (2014). Trends in outpatient visits for insomnia, sleep apnea, and prescriptions for sleep medications among US adults: findings from the National Ambulatory Medical Care Survey 1999-2010. Sleep37, 12831293. 10.5665/sleep.3914

  • 48

    FoxJ.WeisbergS. (2019). An R Companion to Applied Regression. 3rd Edn. Thousand Oaks, CA: Sage.

  • 49

    FredricksonB. L. (2004). “Gratitude, like other positive emotions, broadens and builds,” in The Psychology of Gratitude, eds. R. A. Emmons and M. E. McCullough (New York, NY: Oxford University Press), 144166.

  • 50

    GallS. L.SandersonK.SmithK. J.PattonG.DwyerT.VennA. (2016). Bi-directional associations between healthy lifestyles and mood disorders in young adults: The Childhood Determinants of Adult Health Study. Psychol. Med.46, 25352548. 10.1017/S0033291716000738

  • 51

    GarciaA. N.SalloumI. M. (2015). Polysomnographic sleep disturbances in nicotine, caffeine, alcohol, cocaine, opioid, and cannabis use: a focused review. Am. J. Addict.24, 590598. 10.1111/ajad.12291

  • 52

    GiliM.ComasA.García-GarcíaM.MonzónS.AntoniS.-B.RocaM. (2010). Comorbidity between common mental disorders and chronic somatic diseases in primary care patients. Gen. Hosp. Psychiatry32, 240245. 10.1016/j.genhosppsych.2010.01.013

  • 53

    GoldP. W. (2015). The organization of the stress system and its dysregulation in depressive illness. Mol. Psychiatry20, 3247. 10.1038/mp.2014.163

  • 54

    GräfeK.ZipfelS.HerzogW.LöweB. (2004). Screening psychischer Störungen mit dem “Gesundheitsfragebogen für patienten (PHQ-D)”. Diagnostica50, 171181. 10.1026/0012-1924.50.4.171

  • 55

    GrantA. M.GinoF. (2010). A little thanks goes a long way: explaining why gratitude expressions motivate prosocial behavior. J. Pers. Soc. Psychol.98, 946955. 10.1037/a0017935

  • 56

    HairJ. F.RisherJ. J.SarstedtM.RingleC. M. (2019). When to use and how to report the results of PLS-SEM. EBR31, 224. 10.1108/EBR-11-2018-0203

  • 57

    HajekA.KönigH.-H. (2020). Prevalence and correlates of individuals screening positive for depression and anxiety on the PHQ-4 in the German general population: Findings from the nationally representative German Socio-Economic Panel (GSOEP). IJERPH17:7865. 10.3390/ijerph17217865

  • 58

    HaleL.TroxelW.BuysseD. J. (2020). Sleep health: an opportunity for public health to address health equity. Annu. Rev. Public Health41, 8199. 10.1146/annurev-publhealth-040119-094412

  • 59

    HartescuI.MorganK.StevinsonC. D. (2015). Increased physical activity improves sleep and mood outcomes in inactive people with insomnia: a randomized controlled trial. J. Sleep Res.24, 526534. 10.1111/jsr.12297

  • 60

    HeckendorfH.LehrD.EbertD. D.FreundH. (2019). Efficacy of an internet and app-based gratitude intervention in reducing repetitive negative thinking and mechanisms of change in the intervention's effect on anxiety and depression: results from a randomized controlled trial. Behav. Res. Ther.119:103415. 10.1016/j.brat.2019.103415

  • 61

    HinzA.KleinA. M.BrählerE.GlaesmerH.LuckT.Riedel-HellerS. G.et al. (2017). Psychometric evaluation of the Generalized Anxiety Disorder Screener GAD-7, based on a large German general population sample. J. Affect. Disord.210, 338344. 10.1016/j.jad.2016.12.012

  • 62

    HirschJ. K.AltierH. R.OffenbächerM.ToussaintL.KohlsN.SiroisF. M. (2021). Positive psychological factors and impairment in rheumatic and musculoskeletal disease: do psychopathology and sleep quality explain the linkage?Arthr. Care Res.73, 5564. 10.1002/acr.24440

  • 63

    HopkoD. R.RobertsonS. M. C.CarvalhoJ. P. (2009). Sudden gains in depressed cancer patients treated with behavioral activation therapy. Behav. Ther.40, 346356. 10.1016/j.beth.2008.09.001

  • 64

    HughesA. J.DunnK. M.ChaffeeT.BhattaraiJ.BeierM. (2018). Diagnostic and clinical utility of the GAD-2 for screening anxiety symptoms in individuals with multiple sclerosis. Arch. Phys. Med. Rehabil.99, 20452049. 10.1016/j.apmr.2018.05.029

  • 65

    HunterC. L.American Psychological Associationeds. (2009). Integrated Behavioral Health in Primary Care: Step-By-Step Guidance for Assessment and Intervention, 1st Ed. Washington, DC: American Psychological Association.

  • 66

    IrwinM. R.OlmsteadR.CarrollJ. E. (2016). Sleep disturbance, sleep duration, and inflammation: a systematic review and meta-analysis of cohort studies and experimental sleep deprivation. Bio. Psychiatry80, 4052. 10.1016/j.biopsych.2015.05.014

  • 67

    JackowskaM.BrownJ.RonaldsonA.SteptoeA. (2016). The impact of a brief gratitude intervention on subjective well-being, biology and sleep. J. Health Psychol.21, 22072217. 10.1177/1359105315572455

  • 68

    JacksonC. L.PatelS. R.JacksonW. B.LutseyP. L.RedlineS. (2018). Agreement between self-reported and objectively measured sleep duration among white, black, Hispanic, and Chinese adults in the United States: multi-ethnic study of atherosclerosis. Sleep41:zsy057. 10.1093/sleep/zsy057

  • 69

    JamiesonN. J.HansonC. S.JosephsonM. A.GordonE. J.CraigJ. C.HalleckF.et al. (2016). Motivations, challenges, and attitudes to self-management in kidney transplant recipients: a systematic review of qualitative studies. Am. J. Kid. Dis.67, 461478. 10.1053/j.ajkd.2015.07.030

  • 70

    JonkersC. C. M.LamersF.BosmaH.MetsemakersJ. F. M.Van EijkJ. Th. M. (2012). The effectiveness of a minimal psychological intervention on self-management beliefs and behaviors in depressed chronically ill elderly persons: a randomized trial. Int. Psychogeriatr.24, 288297. 10.1017/S1041610211001748

  • 71

    KalmbachD. A.AndersonJ. R.DrakeC. L. (2018). The impact of stress on sleep: pathogenic sleep reactivity as a vulnerability to insomnia and circadian disorders. J. Sleep Res.27:e12710. 10.1111/jsr.12710

  • 72

    KanagasabaiT.ChaputJ.-P. (2017). Sleep duration and the associated cardiometabolic risk scores in adults. Sleep Health3, 195203. 10.1016/j.sleh.2017.03.006

  • 73

    KaringC. (2021). Prevalence and predictors of anxiety, depression and stress among university students during the period of the first lockdown in Germany. J. Affect. Disord. Rep.5:100174. 10.1016/j.jadr.2021.100174

  • 74

    KlibertJ.RochaniH.SamawiH.Leleux-LaBargeK.RyanR. (2019). The impact of an integrated gratitude intervention on positive affect and coping resources. Int. J. Appl. Posit. Psychol.3, 2341. 10.1007/s41042-019-00015-6

  • 75

    KorzhO.TitkovaA.FylenkoY.LavrovaY. (2022). Evaluation of health-promoting self-care behaviors in hypertensive patients with concomitant chronic kidney disease in primary care. Prim. Health Care Res. Dev.23:e48. 10.1017/S1463423622000299

  • 76

    KroenkeK.SpitzerR. L.WilliamsJ. B. W. (2001). The PHQ-9: validity of a brief depression severity measure. J. Gen. Intern. Med.16, 606613. 10.1046/j.1525-1497.2001.016009606.x

  • 77

    KroenkeK.SpitzerR. L.WilliamsJ. B. W. (2003). The patient health questionnaire-2: validity of a two-item depression screener. Med. Care41, 12841292. 10.1097/01.MLR.0000093487.78664.3C

  • 78

    KroenkeK.SpitzerR. L.WilliamsJ. B. W.MonahanP. O.LöweB. (2007). Anxiety disorders in primary care: prevalence, impairment, comorbidity, and detection. Ann. Intern. Med.146:317. 10.7326/0003-4819-146-5-200703060-00004

  • 79

    LégerD.PartinenM.HirshkowitzM.ChokrovertyS.HednerJ. (2010). Characteristics of insomnia in a primary care setting: EQUINOX survey of 5293 insomniacs from 10 countries. Sleep Med.11, 987998. 10.1016/j.sleep.2010.04.019

  • 80

    LeglerS. R.BealeE. E.CelanoC. M.BeachS. R.HealyB. C.HuffmanJ. C. (2019). State gratitude for one's life and health after an acute coronary syndrome: prospective associations with physical activity, medical adherence and re-hospitalizations. J. Posit. Psychol.14, 283291. 10.1080/17439760.2017.1414295

  • 81

    LiuX.CaoH.ZhuH.ZhangH.NiuK.TangN.et al. (2021). Association of chronic diseases with depression, anxiety and stress in Chinese general population: the CHCN-BTH cohort study. J. Affect. Disord.282, 12781287. 10.1016/j.jad.2021.01.040

  • 82

    LoprestiA. L.HoodS. D.DrummondP. D. (2013). A review of lifestyle factors that contribute to important pathways associated with major depression: diet, sleep and exercise. J. Affect. Disord.148, 1227. 10.1016/j.jad.2013.01.014

  • 83

    LorenzT.AlgnerM.BinderB. (2022). A positive psychology resource for students? Evaluation of the effectiveness of the 6 minutes diary in a randomized control trial. Front. Psychol.13:896741. 10.3389/fpsyg.2022.896741

  • 84

    LöweB.DeckerO.MüllerS.BrählerE.SchellbergD.HerzogW.et al. (2008). Validation and standardization of the generalized anxiety disorder screener (GAD-7) in the general population. Med. Care46, 266274. Available online at: http://www.jstor.org/stable/40221654

  • 85

    LöweB.KroenkeK.GräfeK. (2005). Detecting and monitoring depression with a two-item questionnaire (PHQ-2). J. Psychosom. Res.58, 163171. 10.1016/j.jpsychores.2004.09.006

  • 86

    MacGregorK. L.FunderburkJ. S.PigeonW.MaistoS. A. (2012). Evaluation of the PHQ-9 item 3 as a screen for sleep disturbance in primary care. J. Gen. Intern. Med.27, 339344. 10.1007/s11606-011-1884-5

  • 87

    MacKinnonD. P. (2000). Equivalence of the mediation, confounding and suppression effect. Prevent. Sci.1, 173181. 10.1023/A:1026595011371

  • 88

    MacKinnonD. P.FairchildA. J.FritzM. S. (2007). Mediation analysis. Ann. Rev. Psychol.58, 593614. 10.1146/annurev.psych.58.110405.085542

  • 89

    ManeaL.GilbodyS.HewittC.NorthA.PlummerF.RichardsonR.et al. (2016). Identifying depression with the PHQ-2: a diagnostic meta-analysis. J. Affect. Disord.203, 382395. 10.1016/j.jad.2016.06.003

  • 90

    MartinF.WrightH.MoodyL.WhitemanB.McGillionM.ClyneW.et al. (2020). Help to overcome problems effectively for cancer survivors: development and evaluation of a digital self-management program. J. Med. Internet Res.22:e17824. 10.2196/17824

  • 91

    MayD.LitvinB.AllegranteJ. (2024). Behavioral activation, depression, and promotion of health behaviors: a scoping review. Health Educ. Behav.51, 321331. 10.1177/10901981221090157

  • 92

    McCullaghP. (1980). Regression models for ordinal data. J. R. Stat. Soc. Series B Statist. Methodol.42, 109127. 10.1111/j.2517-6161.1980.tb01109.x

  • 93

    McCulloughM. E.EmmonsR. A.TsangJ.-A. (2002). The grateful disposition: a conceptual and empirical topography. J. Person. Soc. Psychol.82, 112127. 10.1037//0022-3514.82.1.112

  • 94

    MenardS. (2004). Six approaches to calculating standardized logistic regression coefficients. Am. Statist.58, 218223. 10.1198/000313004X946

  • 95

    MillsP. J.RedwineL.WilsonK.PungM. A.ChinhK.GreenbergB. H.et al. (2015). The role of gratitude in spiritual well-being in asymptomatic heart failure patients. Spirit. Clin. Pract.2, 517. 10.1037/scp0000050

  • 96

    MorinC. M. (1993). Insomnia: Psychological Assessment and Management. New York, NY: Guilford Press.

  • 97

    NewmanD. B.GordonA. M.MendesW. B. (2021). Comparing daily physiological and psychological benefits of gratitude and optimism using a digital platform. Emotion21, 13571365. 10.1037/emo0001025

  • 98

    NgM.-Y.WongW.-S. (2013). The differential effects of gratitude and sleep on psychological distress in patients with chronic pain. J. Health Psychol.18, 263271. 10.1177/1359105312439733

  • 99

    OhayonM. M.ReynoldsC. F.DauvilliersY. (2013). Excessive sleep duration and quality of life. Ann. Neurol.73, 785794. 10.1002/ana.23818

  • 100

    O'LaughlinK. D.MartinM. J.FerrerE. (2018). Cross-sectional analysis of longitudinal mediation processes. Multivar. Behav. Res.53, 375402. 10.1080/00273171.2018.1454822

  • 101

    PackardA. E. B.EganA. E.Ulrich-LaiY. M. (2016). HPA axis interactions with behavioral systems. Compr. Physiol.6, 18971934. 10.1002/cphy.c150042

  • 102

    Personality and Assessment Group Department of Psychology, University of Zurich. (n.d.). Gratitude Questionnaire–6, German Translation. Available online at: https://www.michael-mccullough.com/research/measures-and-translations (accessed June 18, 2023).

  • 103

    PetersonB.HarrellF. E. (1990). Partial proportional odds models for ordinal response variables. Appl. Stat.39:205. 10.2307/2347760

  • 104

    PlummerF.ManeaL.TrepelD.McMillanD. (2016). Screening for anxiety disorders with the GAD-7 and GAD-2: a systematic review and diagnostic meta-analysis. Gen. Hosp. Psychiatry39, 2431. 10.1016/j.genhosppsych.2015.11.005

  • 105

    PrinceS. A.CardilliL.ReedJ. L.SaundersT. J.KiteC.DouilletteK.et al. (2020). A comparison of self-reported and device measured sedentary behaviour in adults: a systematic review and meta-analysis. Int. J. Behav. Nutr. Phys. Act17:31. 10.1186/s12966-020-00938-3

  • 106

    PrinceS. A.ReidR. D.PipeA. L.McDonnellL. A. (2017). An evaluation of CardioPrevent: a technology-enabled, health-behavior change program for the global reduction of cardiovascular risk. Curr. Opin. Cardiol.32, 580589. 10.1097/HCO.0000000000000438

  • 107

    PuthM.-T.WeckbeckerK.SchmidM.MünsterE. (2017). Prevalence of multimorbidity in Germany: Impact of age and educational level in a cross-sectional study on 19,294 adults. BMC Public Health17:826. 10.1186/s.12889-017-4833-3

  • 108

    R Core Team (2022). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available online at: https://www.R-project.org/

  • 109

    RevelleW. R. (2022). psych: Procedures for Psychological, Psychometric, and Personality Research. Northwestern University, Evanston, Illinois. Available online at: https://CRAN.R-project.org/package=psych

  • 110

    RogersL. Q.CourneyaK. S.OsterR. A.AntonP. M.RobbsR. S.ForeroA.et al. (2017). Physical activity and sleep quality in breast cancer survivors: a randomized trial. Med. Sci. Sports Exer.49, 20092015. 10.1249/MSS.0000000000001327

  • 111

    RoigE. (2017). “Uttering ‘race' in colorblind France and post-racial Germany,” in Rassismuskritik und Widerstandsformen, eds. K. Fereidooni and M. El (Wiesbaden: Springer Fachmedien Wiesbaden), 613627.

  • 112

    RosenstockI. M. (1966). Why people use health services. Milbank Mem. Fund Quart.44, 94127. 10.2307/3348967

  • 113

    RosenstockI. M.StrecherV. J.BeckerM. H. (1988). Social learning theory and the health belief model. Health Educ. Quart.15, 175183. 10.1177/109019818801500203

  • 114

    RossR. A.FosterS. L.IonescuD. F. (2017). The role of chronic stress in anxious depression. Chronic Stress1:247054701668947. 10.1177/2470547016689472

  • 115

    RosseelY. (2012). lavaan: an R Package for structural equation modeling. J. Stat. Soft. 48. 10.18637/jss.v048.i02

  • 116

    RyffC. D.KeyesC. L. M. (1995). The structure of psychological well-being revisited. J. Pers. Soc. Psychol.69, 719727. 10.1037//0022-3514.69.4.719

  • 117

    ScarlettS.NolanH. N.KennyR. A.O'ConnellM. D. L. (2021). Discrepancies in self-reported and actigraphy-based sleep duration are associated with self-reported insomnia symptoms in community-dwelling older adults. Sleep Health7, 8392. 10.1016/j.sleh.2020.06.003

  • 118

    SchoenC.OsbornR.HowS. K. H.DotyM. M.PeughJ. (2008). In chronic condition: experiences of patients with complex health care needs, in eight countries, 2008. Health Affairs27, w1w16. 10.1377/hlthaff.28.1.w1

  • 119

    SchulerM.StrohmayerM.MühligS.SchwaighoferB.WittmannM.FallerH.et al. (2018). Assessment of depression before and after inpatient rehabilitation in COPD patients: psychometric properties of the German version of the Patient Health Questionnaire (PHQ-9/PHQ-2). J. Affect. Disord.232, 268275. 10.1016/j.jad.2018.02.037

  • 120

    SchulteT.HofmeisterD.Mehnert-TheuerkaufA.HartungT.HinzA. (2021). Assessment of sleep problems with the Insomnia Severity Index (ISI) and the sleep item of the Patient Health Questionnaire (PHQ-9) in cancer patients. Supp. Care Cancer29, 73777384. 10.1007/s00520-021-06282-x

  • 121

    SchwarzerR.JerusalemM. (1995). “Generalized self-efficacy scale,” in Measures in Health Psychology: A User's Portfolio, eds. J. Weinman, S. Wright, and M. Johnston (Windsor: NFER-Nelson), 3537.

  • 122

    ShorofskyM.BourbeauJ.KimoffJ.JenR.MalhotraA.AyasN.et al. (2019). Impaired sleep quality in COPD is associated with exacerbations. Chest156, 852863. 10.1016/j.chest.2019.04.132

  • 123

    SiengsukonC. F.AlshehriM.WilliamsC.DrerupM.LynchS. (2020). Feasibility and treatment effect of cognitive behavioral therapy for insomnia in individuals with multiple sclerosis: a pilot randomized controlled trial. Mult. Scler. Relat. Disord.40:101958. 10.1016/j.msard.2020.101958

  • 124

    SimonP. (2012). Collecting ethnic statistics in Europe: a review. Ethnic Racial Stud.35, 13661391. 10.1080/01419870.2011.607507

  • 125

    SiroisF. M. (2001). The Wellness Behaviors Inventory.Ottawa, ON: Carleton University.

  • 126

    SiroisF. M. (2003a). The Control Beliefs Inventory.Ottawa, ON: Carleton University.

  • 127

    SiroisF. M. (2003b). Beyond health locus of control: a multidimensional measure of health-related control beliefs and motivations (PhD diss.). Carleton University, Ottawa, ON, Canada.

  • 128

    SiroisF. M. (2004). Procrastination and intentions to perform health behaviors: the role of self-efficacy and the consideration of future consequences. Pers. Indiv. Diff.37, 115128. 10.1016/j.paid.2003.08.005

  • 129

    SiroisF. M. (2007). “I'll look after my health, later”: a replication and extension of the procrastination–health model with community-dwelling adults. Pers. Indiv. Diff.43, 1526. 10.1016/j.paid.2006.11.003

  • 130

    SiroisF. M. (2015). A self-regulation resource model of self-compassion and health behavior intentions in emerging adults. Prevent. Med. Rep.2, 218222. 10.1016/j.pmedr.2015.03.006

  • 131

    SiroisF. M.KitnerR.HirschJ. K. (2015). Self-compassion, affect, and health-promoting behaviors. Health Psychol.34, 661669. 10.1037/hea0000158

  • 132

    SouthwellS.GouldE. (2017). A randomised wait list-controlled pre–post–follow-up trial of a gratitude diary with a distressed sample. J. Posit. Psychol.12, 579593. 10.1080/17439760.2016.1221127

  • 133

    SpitzerR. L.KroenkeK.WilliamsJ. B. W.LöweB. (2006). A brief measure for assessing generalized anxiety disorder: the GAD-7. Arch. Intern. Med.166:1092. 10.1001/archinte.166.10.1092

  • 134

    StächeleT.VolzH.-P. (2013). Taschenatlas Stress [Pocket Atlas of Stress].Berlin: Aesopus Publishing.

  • 135

    StaplesL. G.DearB. F.GandyM.FogliatiV.FogliatiR.KarinE.et al. (2019). Psychometric properties and clinical utility of brief measures of depression, anxiety, and general distress: the PHQ-2, GAD-2, and K-6. Gen. Hosp. Psychiatry56, 1318. 10.1016/j.genhosppsych.2018.11.003

  • 136

    SuY.WangS.-B.ZhengH.TanW.-Y.LiX.HuangZ.-H.et al. (2021). The role of anxiety and depression in the relationship between physical activity and sleep quality: a serial multiple mediation model. J. Affect. Disord.290, 219226. 10.1016/j.jad.2021.04.047

  • 137

    SwainN.Lennox ThompsonB.GallagherS.PaddisonJ.MercerS. (2020). Gratitude Enhanced Mindfulness (GEM): a pilot study of an internet-delivered programme for self-management of pain and disability in people with arthritis. J. Posit. Psychol.15, 420426. 10.1080/17439760.2019.1627397

  • 138

    ThorsteinssonE. B.BrownR. F.OwensM. T. (2019). Modeling the effects of stress, anxiety, and depression on rumination, sleep, and fatigue in a nonclinical sample. J. Nerv. Ment. Dis.207, 355359. 10.1097/NMD.0000000000000973

  • 139

    UllmannS.StraußB.TieslerF.SchneiderN.GensichenJ.Brenk-FranzK. (2022). Patient characteristics as predictors of sleep disturbances in multimorbid primary care patients—a secondary analysis. J. Psychos. Res.153:110705. 10.1016/j.jpsychores.2021.110705

  • 140

    VallejoM. A.Vallejo-SlockerL.Fernández-AbascalE. G.MañanesG. (2018). Determining factors for stress perception assessed with the Perceived Stress Scale (PSS-4) in Spanish and other European samples. Front. Psychol.9:37. 10.3389/fpsyg.2018.00037

  • 141

    VerniaF.Di RuscioM.CicconeA.ViscidoA.FrieriG.StefanelliG.et al. (2021). Sleep disorders related to nutrition and digestive diseases: a neglected clinical condition. Int. J. Med. Sci.18, 593603. 10.7150/ijms.45512

  • 142

    WalshR. F. L.SmithL. T.KlugmanJ.TitoneM. K.NgT. H.GoelN.et al. (2023). An examination of bidirectional associations between physical activity and mood symptoms among individuals diagnosed and at risk for bipolar spectrum disorders. Behav. Res. Ther.161:104255. 10.1016/j.brat.2023.104255

  • 143

    Weziak-BialowolskaD.BialowolskiP.NiemiecR. M. (2023). Character strengths and health-related quality of life in a large international sample: a cross-sectional analysis. J. Res. Pers.103:104338. 10.1016/j.jrp.2022.104338

  • 144

    WoodA. M.FrohJ. J.GeraghtyA. W. A. (2010). Gratitude and well-being: a review and theoretical integration. Clin. Psychol. Rev.30, 890905. 10.1016/j.cpr.2010.03.005

  • 145

    WoodA. M.JosephS.LinleyP. A. (2007). Coping style as a psychological resource of grateful people. J. Soc. Clin. Psychol.26, 10761093. 10.1521/jscp.2007.26.9.1076

  • 146

    WoodA. M.MaltbyJ.StewartN.JosephS. (2008). Conceptualizing gratitude and appreciation as a unitary personality trait. Pers. Indiv. Diff.44, 621632. 10.1016/j.paid.2007.09.028

  • 147

    WuS. M.AmtmannD. (2013). Psychometric evaluation of the perceived stress scale in multiple sclerosis. ISRN Rehabil. 2013, 19. 10.1155/2013/608356

  • 148

    YeeT. W. (2015). Vector Generalized Linear and Additive Models: With an Implementation in R.New York, NY: Springer.

  • 149

    YuanS. L. K.CoutoL. A.MarquesA. P. (2021). Effects of a six-week mobile app versus paper book intervention on quality of life, symptoms, and self-care in patients with fibromyalgia: a randomized parallel trial. Braz. J. Phys. Ther.25, 428436. 10.1016/j.bjpt.2020.10.003

  • 150

    ZengL.-N.ZongQ.-Q.YangY.ZhangL.XiangY.-F.NgC. H.et al. (2020). Gender difference in the prevalence of insomnia: a meta-analysis of observational studies. Front. Psychiatry11:577429. 10.3389/fpsyt.2020.577429

  • 151

    ZimmermannT.PuschmannE.Van Den BusscheH.WieseB.ErnstA.PorzeltS.et al. (2016). Collaborative nurse-led self-management support for primary care patients with anxiety, depressive or somatic symptoms: cluster-randomised controlled trial (findings of the SMADS study). Int. J. Nurs. Stud.63, 101111. 10.1016/j.ijnurstu.2016.08.007

Summary

Keywords

insomnia, sleep disturbance, health behaviors, health self-efficacy, psychological distress, gratitude, primary care

Citation

Altier HR, Hirsch JK, Weber A, Kohls N, Schelling J, Toussaint LL, Sirois FM and Offenbächer M (2025) Gratitude and sleep disturbance in primary care patients: the mediating roles of health self-efficacy, health behaviors, and psychological distress. Front. Sleep 4:1459854. doi: 10.3389/frsle.2025.1459854

Received

05 July 2024

Accepted

19 March 2025

Published

17 April 2025

Volume

4 - 2025

Edited by

Kristi E. Pruiksma, The University of Texas Health Science Center at San Antonio, United States

Reviewed by

Jessica M. Meers, United States Department of Veterans Affairs, United States

Rebecca Campbell, University of Arizona, United States

Updates

Copyright

*Correspondence: Heather R. Altier Jörg Schelling

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics