Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Identifying pathways to religious service attendance among older adults: A lagged exposure-wide analysis

  • Richard G. Cowden ,

    Roles Conceptualization, Writing – original draft, Writing – review & editing

    rcowden@fas.harvard.edu

    Affiliation Human Flourishing Program, Institute for Quantitative Social Science, Harvard University, Cambridge, Massachusetts, United States of America

  • Julia S. Nakamura,

    Roles Conceptualization, Formal analysis, Writing – review & editing

    Affiliation Department of Psychology, University of British Columbia, Vancouver, BC, Canada

  • Zhuo Job Chen,

    Roles Writing – review & editing

    Affiliation School of Nursing, University of North Carolina at Charlotte, Charlotte, NC, United States of America

  • Brendan Case,

    Roles Writing – review & editing

    Affiliation Human Flourishing Program, Institute for Quantitative Social Science, Harvard University, Cambridge, Massachusetts, United States of America

  • Eric S. Kim,

    Roles Funding acquisition, Project administration, Supervision, Writing – review & editing

    Affiliations Human Flourishing Program, Institute for Quantitative Social Science, Harvard University, Cambridge, Massachusetts, United States of America, Department of Psychology, University of British Columbia, Vancouver, BC, Canada, Lee Kum Sheung Center for Health and Happiness, Harvard T.H. Chan School of Public Health, Boston, MA, United States of America

  • Tyler J. VanderWeele

    Roles Conceptualization, Funding acquisition, Methodology, Supervision, Writing – review & editing

    Affiliations Human Flourishing Program, Institute for Quantitative Social Science, Harvard University, Cambridge, Massachusetts, United States of America, Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, United States of America, Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, United States of America

Abstract

We used prospective data (spanning 8 years) from a national sample of older U.S. adults aged > 50 years (the Health and Retirement Study, N = 13,771) to evaluate potential factors that lead to subsequent religious service attendance. We applied a lagged exposure-wide epidemiologic design and evaluated 60 candidate predictors of regular subsequent religious service attendance. Candidate predictors were drawn from the following domains: health behaviors, physical health, psychological well-being, psychological distress, social factors, and work. After rigorous adjustment for a rich set of potential confounders, we observed modest evidence that changes in some indices of physical health, psychological well-being, psychological distress, and social functioning predicted regular religious service attendance four years later. Our findings suggest that there may be opportunities to support more regular religious service attendance among older adults who positively self-identify with a religious/spiritual tradition (e.g., aid services for those with functional limitations, psychological interventions to increase hope), which could have downstream benefits for various dimensions of well-being in the later years of life.

Introduction

A consistent theme across many cultures and religious traditions is that holistic well-being intersects the mind, body, and spirit [13]. Religion is a multidimensional concept involving beliefs, behaviors, rituals, and ceremonies related to the transcendent (e.g., God, Higher Power), and research has shown that religion may have positive or negative influences on health and well-being [4, 5]. Religious engagement often includes private (e.g., prayer, scripture reading) and public forms of participation (e.g., religious service attendance, study groups). A growing body of robust longitudinal evidence suggests that religious service attendance in particular is associated with indicators of improved physical, behavioral, and psychosocial health and well-being across adulthood, including among older adults [613]. For example, research suggests that more regular religious service attendance is related to lower risk of cause-specific (e.g., suicide) and all-cause mortality, improved health behaviors (e.g., reduced smoking), better psychological well-being (e.g., reduced depression), and higher social well-being (e.g., more social connectedness). Religious service attendance is also associated with virtuous and prosocial actions, such as volunteering and charitable giving [1419], which confer benefits to broader society (e.g., reducing the social and economic burden of substance abuse and depression, providing tangible support to those in need).

Religious service attendance appears to affect well-being more strongly than other aspects of religious affiliation taken in isolation, including private prayer or doctrinal beliefs [2022]. Although most of these findings are drawn from American samples consisting largely of Christians, similar results have been found using data on Christians in highly-secular Denmark [23], Israeli Jews [24], and a diverse mix of Taoists, Buddhists, animists, and Christians in Taiwan [22].

Despite the potential role of religious service attendance in supporting the well-being of individuals and communities, evidence indicates that religious service attendance is declining in rich, highly developed parts of the world (e.g., North America [25]). For example, estimates suggest that the percentage of U.S. adults who do not attend religious services has more than doubled over the last few decades (from 11% in the early 1970s to 26% in 2014), with the most striking declines in religious service attendance occurring after the start of the 21st century [26]. This trend cannot be exclusively accounted for by population replacement, as Twenge and colleagues [26] observed a general pattern of decline in religious orientation among all generations in the U.S. population. Research suggests that this cross-generational trend may to be linked to a decrease in church membership, including a roughly 10% decline in church membership among older generations of U.S. adults (e.g., traditionalists, baby boomers) over the last two decades [27].

Many factors seem to drive religious disaffiliation and declining participation in the U.S. and elsewhere. One standard theory is that religion suffers in developed economies with high levels of education and robust social safety nets, which unsettle traditional certitudes and substitute for the social support historically provided by religious communities [28, 29]. This thesis only imperfectly accounts for trends in the U.S., however, where religious participation has declined fastest among the least educated and least affluent, which is part of a more general trend toward social disaffiliation among the lower classes (in large part due to their declining economic prospects), as evident in declining rates of marriage, unionization, and civic participation more broadly [30]. Over the past 20 years, many religious communities have also been impacted by scandal, notably over sexual abuse within churches and its subsequent cover-up by leaders, which has negatively affected public trust in religious communities and likely contributed to disaffiliation as well [31].

Given the apparent decline in religious service attendance within the U.S. (and other parts of the Western world), research is needed to identify and better understand the factors that might influence religious service attendance. Such evidence could contribute to determining suitable pathways to increase religious service attendance or identifying barriers that impede more regular religious service attendance, particularly among those who positively self-identify with a religious/spiritual tradition and subpopulations (e.g., older adults) for whom the benefits of religious services attendance might be especially pronounced. For example, advancing age is often accompanied by hardships in various domains of life, including physical (e.g., functional limitations), social (e.g., loss of loved ones), and emotional (e.g., grief) afflictions [32]. Religious service attendance can provide older adults with opportunities to access, develop, or acquire a wide range of supportive resources, both spiritual (e.g., closer connection to the divine) and nonspiritual (e.g., social relationships), which they can draw on as they encounter such adversities during their later years of life. A more thorough understanding of the potential pathways and/or barriers to religious service attendance among older adults could be a useful stepping-stone towards future research and applied work aimed at increasing religious service attendance in this subpopulation, which may precipitate additional gains in resources that support adjustment during the later years of life. In the present study, we apply a rigorous analytic template to examine an array of potential predictors of religious service attendance in a large cohort of older U.S. adults.

Prior research that has examined factors leading to an increase in religious service attendance has contributed substantially to the literature, but there are several notable gaps in the existing evidence. First, many prior studies in this area (including those with older adults) have relied on cross-sectional data, and such designs cannot infer causality [7]. For example, cross-sectional data may reveal that depression is associated with a lower likelihood of religious service attendance either because depressed people may stop attending, or because people who cease attending religious services are more likely to become depressed, which generates concerns about potential reverse causation [33]. Additionally, the observed associations in cross-sectional studies may overestimate the actual effect of an exposure on an outcome [34]. Hence, additional longitudinal evidence is needed to enrich our understanding of the strength and directionality of associations between candidate predictors and religious service attendance.

Second, existing longitudinal studies on this topic (including those involving older adults) typically address one or a narrow range of exposure variables [11, 33, 3540]. For example, depression is one of the most frequently studied predictors of religious service attendance, but it has seldom been examined alongside other potential determinants of religious service attendance. In one of the few longitudinal studies that extended beyond psychological predictors of religious service attendance, Ferraro and Kelley-Moore [36] observed that not smoking and more frequent volunteering were associated with increased subsequent religious service attendance. However, none of the physical health indicators (e.g., chronic illness) that they examined predicted change in religious service attendance. Although these prior studies (most of which have principally applied theory-based methods to select exposures) have identified important potential determinants of religious service attendance, a more complete picture of the dynamics that shape religious service attendance among older adults may be acquired by examining a wider range of possible predictors across multiple domains of functioning. Specifically, examining an array of predictors across different domains of functioning (including those that have been considered less frequently in previous work) within the same sample may provide novel insights into the strongest predictors (as compared to many others across domains) of subsequent religious service attendance. Research along these lines could also help to build the kind of foundation needed to develop interventions that address pathways to religious service attendance in a more holistic and integrative way.

Third, longitudinal studies in this area have rarely controlled for prior values of the exposures, which can help reduce concerns about reverse causation [41]. Adjusting for prior values of the exposures can also elucidate how changes in each exposure—much like what one might expect to see from an intervention targeting that exposure—are associated with an outcome of interest. Because many potential determinants of religious service attendance (e.g., depression) may change over time, adjustment for prior values of an exposure can strengthen causal inference [42]. If such analytic decisions can be applied to improve causal inference in observational studies examining predictors of religious service attendance, researchers and practitioners are likely to be equipped with stronger evidence that they could use to make more informed decisions about how to increase religious service attendance in a particular target population and the resources they may need to allocate in order to achieve this objective.

To help enrich the existing empirical literature in this area, the purpose of the current study is to explore the determinants of change in religious service attendance among older U.S. adults. We applied a lagged exposure-wide epidemiologic design, which is a data-driven approach aimed at elucidating promising predictors of a particular outcome that can then be examined further in subsequent research. This analytic approach includes steps that rigorously control for potential confounding and dampen concerns about reverse causation [41], thereby strengthening causal inference. Using a large cohort of older U.S. adults, we test for evidence concerning potential causal effects of 60 candidate determinants of religious service attendance, including indicators of health behaviors, physical health, psychological well-being, psychological distress, social factors, and work.

Method

Study population

We used data from the Health and Retirement Study (HRS), a national panel study of adults aged > 50 years in the U.S. In 2006, about half of HRS respondents completed an enhanced face-to-face interview (when psychosocial data were first collected). The other half of respondents were assessed in 2008. Participants then completed a psychosocial questionnaire (response rates of 88% in 2006 and 84% in 2008), which they mailed to the University of Michigan upon completion [43]. These sub-cohorts alternate reporting on psychosocial factors, with each participant reporting psychosocial data every four years. In this study, data from the 2006 and 2008 sub-cohorts were combined to increase sample size and statistical power. Participants were excluded if they did not report psychosocial data in 2006/2008, because over half of our candidate predictors were psychosocial factors. This resulted in a final sample of 13,771 participants.

This study used data from three time points: (1) covariates were assessed in 2006/2008 (T0), (2) candidate predictors were assessed in 2010/2012 (T1), and (3) the outcome (religious service attendance) was assessed in 2014/2016 (T2). The HRS (see further study documentation on the HRS website: http://hrsonline.isr.umich.edu/) is sponsored by the National Institute on Aging (NIA U01AG009740) and conducted by the University of Michigan [44]. The ethics board at Harvard University exempted the present study from review because it used de-identified and publicly available data.

Measures

Religious service attendance.

Respondents were asked: “About how often have you attended religious services during the past year?” Response options included: > 1x/week, 1x/week, 2-3x/month, ≥ 1x/year, and never. Results for both (1) the binary religious service attendance (0 = < 1x/week, 1 = ≥ 1x/week) and (2) the ordinal religious service attendance variables (1 = rarely or never, 2 = 2-3x/month, 3 = 1x/week, 4 = > 1x/week) are available in the main text (binary variable) and supplementary materials (ordinal variable). We used the binary measure of religious service attendance (< 1x/week or ≥ 1x/week) in our main analyses because religious service attendance at the threshold of ≥ 1x/week has been associated with improved health and well-being [7]. Further, ordinal measures of religious service attendance are more difficult to interpret (as these analyses represent movement between subcategories of religious service attendance), and there was relatively little movement between categories using the ordinal variable in our study (see S1 Table, which shows the changes in religious service attendance from T0 to T2).

Covariates.

We adjusted for numerous covariates assessed at T0, including age (continuous), gender (male, female), race/ethnicity (White, African-American, Hispanic, other), marital status (not married, married), income (< $50,000, $50,000-$74,999, $75,000-$99,999, ≥ $100,000), total wealth (based on quintiles of the score distribution for total wealth in this sample), educational attainment (no degree, GED/high school diploma, ≥ college degree), employment status (no, yes), health insurance (no, yes), geographic region (Northeast, Midwest, South, West), personality traits (openness, conscientiousness, extraversion, agreeableness, neuroticism; continuous), childhood abuse (no, yes), and religious service attendance (as either a binary or ordinal variable to match the coding of the outcome variable). To examine change in each candidate predictor, we additionally adjusted for prior values of all candidate predictors assessed at T0 (except for the depression variable, which was constructed from the depressive symptoms variable).

Predictors.

We evaluated 60 candidate predictors assessed at T1 (2010/2012). These included measures of (1) health behaviors (frequent physical activity, smoking, heavy drinking, sleep problems); (2) physical health (diabetes, hypertension, stroke, cancer, heart disease, lung disease, arthritis, overweight/obese, physical functioning limitations, cognitive impairment, chronic pain, self-rated health, hearing, eyesight); (3) psychological well-being (positive affect, life satisfaction, optimism, purpose in life, personal mastery, health mastery, financial mastery); (4) psychological distress (depression, depressive symptoms, hopelessness, negative affect, perceived constraints, anxiety symptoms, trait anger, state anger, cynical hostility, stressful life events, financial strain, daily discrimination, major discrimination); (5) social factors (living with a spouse/partner, frequency of contact with (i) children, (ii) other family, and (iii) friends, loneliness, closeness with spouse, number of close (i) children, (ii) other family, and (iii) friends, social support from (i) spouse, (ii) children, (iii) other family, and (iv) friends, social strain from (i) spouse, (ii) children, (iii) other family, and (iv) friends, volunteering, helping friends, neighbors, or relatives, social status ladder ranking, change in social status ladder ranking); and (6) work (in the labor force). S1 Appendix and the HRS materials provide further details about each of these variables [43, 45, 46].

Multiple imputation.

There were missing data for the covariates (up to 10.2%), candidate predictors (up to 56.2%), and outcome (33.8%). Depending on the predictor being evaluated, complete-case analyses resulted in a loss of 60.4% to 78.1% of the total sample. Thus, we imputed missing data for the covariates, candidate predictors, and outcome using imputation by chained equations, and five datasets were created. This method tends to be more flexible than other methods for handling missing data [47, 48], and it also helps address problems that arise from attrition [49, 50].

Statistical analysis

Primary analysis.

We used an exposure-wide analytic approach [41], with separate models estimated for each candidate predictor (see S2 Appendix). In our primary analysis, religious service attendance was a binary outcome with a prevalence ≥ 10%. Thus, we used generalized linear models (with a log link and Poisson distribution) to individually regress religious service attendance at T2 (2014/2016) on each candidate predictor at T1 (2010/2012). All continuous predictors were standardized (mean = 0, standard deviation = 1) so that their effect sizes could be interpreted as a standard deviation change in the exposure variable. We marked multiple p-value cut-offs in our tables (including Bonferroni-corrected), as multiple testing practices vary widely and are continuously evolving. We also report exact confidence intervals. Unless otherwise specified, we focus our discussion on the primary analytic results with religious service attendance modeled as a binary outcome in the full sample.

Additional analyses.

We conducted several additional analyses. First, we used E-values to assess the minimum strength of association (on the risk ratio scale) that an unmeasured confounder would need to have with both the exposure and the outcome to explain away the observed association [51]. This helps us evaluate the robustness of the results from the primary analysis to potential unmeasured confounding. Second, we repeated the primary analysis among the subsample of participants who regularly attended religious services (≥ 1x/week) at T1 (S2 Table). This allowed us to identify what factors predict continued religious service attendance, four years later, for people who were already attending religious services regularly at T1. Third, we repeated the primary analysis among the subsample of participants who did not attend religious services regularly (< 1x/week) at T1 (S3 Table). This allowed us to identify what factors predict regular religious service attendance, four years later, for people who did not regularly attend religious services at T1. Fourth, we replicated all analyses using the abovementioned ordinal religious service attendance variable (1 = rarely or never, 2 = 2-3x/month, 3 = 1x/week, 4 = > 1x/week) by applying ordinal logistic regression, although this makes the additional assumption that the odds ratio for the covariates across each of two adjacent ordinal outcomes categories is constant. This assumption allows for greater statistical power and potential additional nuance in the results. Specifically, we repeated these analyses (1) in the entire analytic sample (S4 Table), (2) only among people who attended religious services regularly (≥ 1x/week) at T1 (S5 Table), and (3) only among people who did not attend religious services regularly (< 1x/week) at T1 (S6 Table). Fifth, we reproduced the primary analysis using complete cases to assess the impact of multiple imputation on results (S7 Table).

Results

At T0, participants were on average 69 years old (SD = 10), a majority of whom were women (58%) and married (62%). Table 1 provides the distribution of covariates by religious service attendance at T0.

thumbnail
Table 1. Distribution of participant characteristics (T0) by religious service attendance (T0)a,b,c.

https://doi.org/10.1371/journal.pone.0278178.t001

Primary analysis

Results for the exposure-wide analysis involving candidate predictors of regular subsequent religious service attendance are reported in Table 2. There was little evidence of associations between health behaviors and regular subsequent religious service attendance.

thumbnail
Table 2. Candidate predictors (T1) of subsequent religious service attendance (T2)a,b,c.

https://doi.org/10.1371/journal.pone.0278178.t002

For physical health factors, having physical functioning limitations was associated with a 16% lower likelihood of attending religious services regularly (95% CI: 0.77, 0.92) four years later. For each standard deviation increase in self-rated health, participants were 9% more likely to subsequently attend religious services regularly (95% CI: 1.05, 1.13). However, there was little evidence of associations between other physical health indicators and regular subsequent religious service attendance.

Among psychological well-being factors, participants were more likely to regularly attend religious services four years later for each standard deviation increase in positive affect (RR = 1.06; 95% CI: 1.02, 1.11), life satisfaction (RR = 1.05; 95% CI: 1.02, 1.09), purpose in life (RR = 1.06; 95% CI: 1.00, 1.11), and health mastery (RR = 1.05; 95% CI: 1.01, 1.09). Among psychological distress factors, participants were 6% less likely to regularly attend religious services four years later for each standard deviation increase in hopelessness (95% CI: 0.90, 0.98) and perceived constraints (95% CI: 0.90, 0.98). There was little evidence of associations between other psychological well-being and distress factors and regular subsequent religious service attendance.

For social factors, participants who were in contact with friends 1-2x/month (RR = 1.13; 95% CI: 1.02, 1.26), 1-2x/week (RR = 1.19; 95% CI: 1.08, 1.31), or ≥ 3x/week (RR = 1.19; 95% CI: 1.06, 1.33) were more likely to regularly attend religious services four years later. For each standard deviation change in positive social support from friends (95% CI: 1.00, 1.08) and social strain from friends (95% CI: 1.00, 1.07), participants were 4% more likely to subsequently attend religious services regularly. Participants who volunteered 1–49 hours/year (RR = 1.32; 95% CI: 1.20, 1.45), 50–99 hours/year (RR = 1.35; 95% CI: 1.22, 1.49), 100–199 hours/year (RR = 1.40; 95% CI: 1.26, 1.55), and ≥ 200 hours/year (RR = 1.39; 95% CI: 1.23, 1.57) were more likely to regularly attend religious services four years later. There was little evidence of associations between other social factors and regular subsequent religious service attendance. Work was also not associated with regular subsequent religious service attendance.

Additional analyses

E-values suggested that many of the observed associations were moderately robust to unmeasured confounding (Table 3). To illustrate, for physical functioning limitations, an unmeasured confounder associated with both physical functioning limitations and religious service attendance by risk ratios of 1.66 each (above and beyond the covariates already adjusted for) could explain away the association, but weaker joint confounding associations could not. Further, to shift the confidence interval to include the null, an unmeasured confounder associated with both physical functioning limitations and religious service attendance by risk ratios of 1.39 could suffice, but weaker joint confounding associations could not.

thumbnail
Table 3. Robustness to unmeasured confounding (E-values) for the associations between candidate predictors (T1) and subsequent religious service attendance (T2)a,b,c.

https://doi.org/10.1371/journal.pone.0278178.t003

Results were somewhat different (see S2 and S3 Tables) when comparing candidate predictors for the subsample of participants who attended religious services regularly at T1 vs. those who did not. For example, some indices of physical health (e.g., cancer), psychological well-being (e.g., positive affect), and psychological distress (e.g., hopelessness) were predictive of regular subsequent attendance among participants who did not attend regularly at T1, but not among those who did attend regularly at T1. When religious service attendance was modeled as an ordinal variable, results for the entire sample, among participants who attended religious services regularly at T1, and among those who did not attend religious services regularly at T1 showed some differences from each other as well as from results of analyses in which religious service attendance was modeled as a binary variable (see S4S6 Tables for the entire analytic sample, for those who did regularly attend religious services at T1, and for those who did not regularly attend religious services at T1, respectively). For example, when using the ordinal religious service attendance variable with the full analytic sample (see S4 Table), we observed associations for some predictors (e.g., frequent physical activity, cognitive impairment, optimism, state anger, helping friends, neighbors, or relatives) that were not associated with subsequent religious service in the primary analysis. This is not entirely surprising, given that ordinal logistic regression tends to increase power within analyses (though at the cost of additional assumptions). The odds ratios in this analysis concern the odds ratio for each T1 predictor comparing two adjacent religious service attendance categories at T2 (e.g., rarely or never compared to 2-3x/month, or 2-3x/month compared to 1x/week, or 1x/week compared to > 1x/week). For odds ratios comparing non-adjacent categories, one exponentiates the odds ratio in the table by how many steps the categories are apart. Results of the complete-case analyses were somewhat similar to the results we observed in the primary analysis (see S7 Table).

Discussion

Using longitudinal data from a prospective sample of older U.S. adults, we evaluated potential determinants of regular (≥ 1x/week) religious service attendance across several domains, including health behavior, physical health, psychological well-being, psychological distress, social factors, and work. In the primary analysis, we observed that changes in some physical health (i.e., physical functioning limitations, self-rated health), psychological well-being (i.e., positive affect, life satisfaction, purpose in life, health mastery), psychological distress (i.e., hopelessness, perceived constraints), and social factors (i.e., contact with friends, positive social support from friends, negative social strain from friends, volunteering) were associated with regular subsequent service attendance four years later. Changes in several predictors that have seldom received attention in this literature (e.g., health mastery) were associated with subsequent religious service attendance in our sample of older adults.

Health behaviors and physical health

Our results did not provide substantial evidence of associations between health behaviors (e.g., frequent physical activity, heavy drinking) and regular subsequent religious service attendance. These findings diverge from some prior longitudinal research. For example, status as a smoker has been modestly associated with a lower likelihood of subsequent religious service attendance among younger and older adults [36]. Although the effect estimate we observed for smoking in the full analytic sample was similar to effect estimates reported in previous work, the confidence interval was wide. This is consistent with the possibility of a null or modest effect of smoking on regular service attendance in the older U.S. adult population, which is because levels of smoking were relatively stable over time in this sample. Additional research involving a wide range of health behaviors in younger samples and over longer periods of time could clarify whether changes in health behaviors might lead to more frequent religious service attendance.

Study participants who reported improvements in physical functioning limitations and self-rated health were more likely to regularly attend religious services four years later. Compared to some earlier work involving older US adults [52], our findings suggest that functional disabilities might have longer-term effects on religious service attendance than previously found. However, none of the other indices of physical health status that we evaluated (e.g., hearing, eyesight) were associated with regular subsequent religious service attendance. Taken together, this study’s fundings indicate that difficulties performing more basic physical activities (e.g., getting up from a chair, walking across a room) and a general decline in physical health that typically comes with aging may make it challenging for some older adults to regularly attend religious services in the longer-term. In light of the potential benefits that older adults might glean from “the fellowship and inspirational experiences” (p. 246) that accompany religious service attendance [53], professional care services (e.g., mobility support) and onsite assistance by communities of faith could enable older adults with functional limitations to attend religious services more regularly.

We found little evidence of association between the chronic conditions examined in this study (e.g., diabetes, hypertension) and regular subsequent religious service attendance. These findings resonate with previous longitudinal research [36] that showed no association of religious service attendance with either generalized chronic illness categories (e.g., circulatory conditions) or specific chronic illnesses (e.g., cancer diagnosis). However, one interesting exception emerged from the stratified analyses; cancer diagnosis was associated with a higher likelihood of regular subsequent religious service attendance among participants who did not attend regularly at T1 (see S3 Table). This finding may reflect the stress mobilization effect [54, 55], which refers to the notion that negative life events can trigger religious/spiritual coping responses (e.g., seeking comfort from God by attending religious services). An alternative but complementary explanation is that irregular attenders at T1 who lived with and survived cancer over the four-year study period may have gone through a process of religious/spiritual transformation that is reflected in more frequent faith-based activities, such as regular religious service attendance [56]. Interestingly, no other chronic illnesses that we assessed were associated with subsequent religious service attendance among those who did not attend regularly at T1. This may be because cancer is one of the most feared chronic illnesses among older adults [57], and it is possible that the perceived existential threat of cancer precipitates stronger change in religious/spiritual involvement when compared to other kinds of serious illness. In future research, subgroup analyses stratified by prior attendance levels may be of practical interest in identifying those who might benefit from initiatives to support their spiritual/religious response to cancer and other life-threatening health conditions.

Psychological well-being and psychological distress

Among the indices of psychological well-being that were examined in this study, we found that increases in some aspects of hedonic well-being (i.e., positive affect, life satisfaction) were associated with a higher likelihood of regular religious service attendance at follow-up. These findings are consistent with the idea that positive subjective experiences tend to prompt approach behavior (e.g., exploration of novel situations, participation in activities [58], which might lead some older adults to increase their engagement with the environment by attending religious services on a regular basis. There was also evidence that increases in some aspects of eudaimonic well-being (i.e., purpose in life, health mastery) were associated with regular subsequent religious service attendance. One possible explanation for these findings is that eudaimonic experiences, such as a renewed or enhanced sense of purpose in life, often motivate people to make positive lifestyle choices aimed at nurturing their well-being [59], including an increase in religious participation in view of the psychosocial and spiritual goods it often affords [60, 61]. Furthermore, older adults often shift their focus from the ‘resumé virtues’ that helped them succeed professionally to the ‘eulogy virtues’ by which they hope to be remembered [62]; it would be no surprise if this shift in priorities prompted a heightened interest in one’s spiritual life as well.

Of the indices of psychological distress that were examined, increases in hopelessness and perceived constraints about personal control were associated with a lower likelihood of regular subsequent religious service attendance. These findings coincide with evidence from previous cross-sectional research [63, 64], but it is not possible to discern the direction of causation from such cross-sectional data. Given the afflictions that tend to befall people in the later stages of life (e.g., deterioration of health, loss of spouses and/or friends), our findings resonate with recent calls for public health activities to aim at bolstering the psychological well-being of older adults [65]. For example, several evidenced-based programs oriented towards promoting hope have received empirical support for use in older populations [66], some of which could be adapted for large-scale implementation. Beyond the potential psychological benefits that might accompany such interventions, the findings of this study suggest that improvements in certain aspects of psychological well-being (e.g., hope, internal locus of control) may lead to an increase in regular service attendance. In turn, robust longitudinal evidence has shown that religious service attendance could have benefits for health and well-being in a variety of life domains [79, 11, 12].

None of the other indices of psychological distress were associated with regular subsequent religious service attendance, including depression. This latter finding warrants specific consideration because there has been substantial scholarly interest in examining the effect of depression on religious service attendance, with prior longitudinal research demonstrating mixed results [33, 35]. There could be many methodological reasons underlying the divergent results seen across studies, such as sample composition (e.g., previous research suggests that depression is more strongly associated with religious service attendance at younger ages [37]) and measurement of depression (e.g., diagnostic interview versus self-report assessment). An important contribution of this study is that we estimated the associations of both depression and hopelessness with religious service attendance, and we also simultaneously adjusted for both variables in all models, neither of which has been done in previous longitudinal research. Drawing on theory and some empirical evidence that suggests hopelessness can be a sufficient proximal cause of a subtype of depression that is referred to as hopelessness depression [67, 68], future research should ensure that assessment of depression is not conflated with hopelessness (or lack of hope). For example, the Center for Epidemiologic Studies Depression Scale 10, which has been used in prior research on religious service attendance and depression [11], has an item that captures a lack of hope about the future. If hopelessness (or a lack of hope) is integrated into measures of depression, it is difficult to pinpoint whether depression or hopelessness (or both) contribute to reduced religious service attendance.

Social factors and work

Compared to all other candidate determinants examined in this study, we found that hours spent volunteering was the strongest predictor of regular subsequent religious service attendance in the full analytic sample. Our results correspond with prior longitudinal research that reported a positive association between more frequent voluntary association participation and religious service attendance [18, 36]. Collectively, these findings suggest that interventions and public health campaigns designed to encourage volunteering among older adults could lead to an increase in religious service attendance. Given the potential for volunteering to support healthy aging in multiple areas of life (e.g., physical health, psychological well-being [69], there may be other important benefits that arise from resources that are dedicated toward encouraging volunteering among older adults (e.g., economic value of volunteer activities [70]). The results of this study may also be of interest to communities of faith, as community service or outreach projects that they organize to benefit society at large could also provide them with opportunities to reach new members and further contribute to the spiritual lives of older adults.

Complementing previous longitudinal studies that have focused on aspects of social well-being as outcomes of religious service attendance among older adults (e.g., social support, social integration [7, 71], we found that an increase in both social contact with friends and social support from friends were associated with a higher likelihood of regular subsequent religious service attendance. These findings suggest that friendships during the later years of life could play an important role in whether older adults attend religious services on a regular basis. For example, older adults who are unable to drive often rely on their friends for transportation to and from religious services and other activities that take place within religious communities [72]. Interestingly, there was also evidence indicating that an increase in social strain from friends was associated with regular subsequent religious service attendance. This finding may reflect religious coping, with older adults who experience strained friendships possibly turning to the religious realm as a source of comfort for dealing with relational difficulties [73] or as an alternative resource because the quality of resources in the social domain have diminished [74]. In alignment with some prior longitudinal research involving middle-aged and older adult participants [36], we found little evidence of association between work and subsequent religious service attendance. With some research indicating that entry into the labor force predicts a decline in religious service attendance among younger adult samples [75], religious service attendance could be more strongly affected by employment changes that accompany transitional periods at younger ages (e.g., transition from adolescence to adulthood). Further research is needed to better understand the impact of life course transitions on patterns of religious service attendance within adulthood.

Limitations, strengths, and future research

This study had several limitations. First, the sample included U.S. adults aged 50 years and older, a population for which some determinants of religious service attendance documented in prior research (e.g., living with a school-age child [39]) are largely irrelevant. Hence, it would be worthwhile to replicate this exposure-wide analysis in younger cohorts. It could also be interesting to extend this kind of analysis to other parts of the world (e.g., Africa) where religions such as Christianity and Islam are experiencing rapid growth within comparatively younger populations [76], and to differentiate between those who are religiously affiliated and those who are not. Second, this study addressed a single dimension of religiosity (i.e., religious service attendance) that does not provide a complete account of religious/spiritual life and is more relevant to the practices of some religious traditions than others. Our understanding of what determines religious involvement, including whether specific factors are common predictors of different dimensions of religiosity, would be enhanced by methodologically rigorous research focused on other aspects of religious/spiritual life. Third, four years of follow-up data may not be sufficient to detect effects for antecedents of religious service attendance that vary or accumulate over time. For example, there was little evidence of an association between loneliness and subsequent religious service attendance. If loneliness has a more pronounced effect on religious service attendance earlier in life, the data and analytic approach used in this study would not be able to detect this effect. Fourth, the findings of this study are limited by our use of self-report data. Although many potential determinants were assessed across several domains of human life, future research could build on this study by integrating objective psychological testing and markers of physical health into measurement procedures. Fifth, a concern with using observational data is that estimated effects might be biased by unmeasured confounding (e.g., intrinsic religiosity) or reverse causality. Our analytic design choices attempted to reduce concerns about unmeasured confounding and reverse causality, and E-values for the estimated effects indicated that the observed associations (e.g., volunteering) were at least somewhat robust to unmeasured confounding. However, we cannot rule out the possibility that the observed effects could be explained away by some combination of statistical uncertainty and unmeasured confounding.

Despite selected limitations, this study used a robust analytic approach to adduce evidence of potential causal effects of numerous candidate predictors of subsequent religious service attendance. Other strengths of our study include our ability to evaluate many understudied predictors of religious service attendance and use of a large, national sample of older adults. Although the results that emerged for selected aspects of physical health, psychological well-being, psychological distress, and social factors are certainly of interest, of equal interest is that there was little evidence that other factors, of those examined, were associated with regular subsequent religious service attendance. It may be that religious service attendance is less malleable or changeable later in life than earlier in life, and further research is needed to investigate this possibility.

Conclusion

This prospective cohort study of older U.S. adults provided modest evidence indicating that improvements in some aspects of functioning in different domains of life (i.e., physical psychological, social) were associated with regular service attendance four years later. Our results suggest that there may be opportunities to support more regular religious service attendance among older adults who positively self-identify with a religious/spiritual tradition (e.g., aid services for those with functional limitations, psychological interventions to increase hope), which itself may confer additional benefits for various dimensions of well-being in the later years of life. Notwithstanding the findings of this study, further research is needed to identify and understand the factors that might influence changes in religious service attendance over time and throughout the life course.

Supporting information

S1 Appendix. Assessment of candidate predictors.

https://doi.org/10.1371/journal.pone.0278178.s001

(DOCX)

S2 Appendix. Proof illustrating how adjusting for the prior value of a predictor can help us evaluate how “changes” in the predictor are associated with subsequent religious service attendance.

https://doi.org/10.1371/journal.pone.0278178.s002

(DOCX)

S1 Table. Changes in religious service attendance from T0 to T2.

https://doi.org/10.1371/journal.pone.0278178.s003

(DOCX)

S2 Table. Candidate predictors (T1) of subsequent religious service attendance (T2) among those who regularly attended religious services at T1.

https://doi.org/10.1371/journal.pone.0278178.s004

(DOCX)

S3 Table. Candidate predictors (T1) of subsequent religious service attendance (T2) among those who did not regularly attend religious services at T1.

https://doi.org/10.1371/journal.pone.0278178.s005

(DOCX)

S4 Table. Candidate predictors (T1) of subsequent religious service attendance (T2) with religious service attendance as an ordinal variable.

https://doi.org/10.1371/journal.pone.0278178.s006

(DOCX)

S5 Table. Candidate predictors (T1) of subsequent religious service attendance (T2) among those who regularly attended religious services at T1 (religious service attendance modeled as an ordinal variable).

https://doi.org/10.1371/journal.pone.0278178.s007

(DOCX)

S6 Table. Candidate predictors (T1) of subsequent religious service attendance (T2) among those who did not regularly attend religious services at T1 (religious service attendance modeled as an ordinal variable).

https://doi.org/10.1371/journal.pone.0278178.s008

(DOCX)

S7 Table. Complete-case analyses for candidate predictors (T1) of subsequent religious service attendance (T2).

https://doi.org/10.1371/journal.pone.0278178.s009

(DOCX)

References

  1. 1. Balboni MJ, Puchalski CM, Peteet JR. The relationship between medicine, spirituality and religion: three models for integration. J Relig Health. 2014;53(5):1586–98. pmid:24917445
  2. 2. Mark GT, Lyons AC. Maori healers’ views on wellbeing: the importance of mind, body, spirit, family and land. Soc Sci Med. 2010;70(11):1756–64. pmid:20338680
  3. 3. Juvenal. The satires. Barr W, editor. Oxford University Press; 1992.
  4. 4. Chatters LM. Religion and health: public health research and practice. Annu Rev Public Health. 2000;21:335–67. pmid:10884957
  5. 5. Koenig HG, Al-Zaben F, VanderWeele TJ. Religion and psychiatry: recent developments in research. BJPsych Adv. 2020;26(5):262–72.
  6. 6. Braam AW, Hein E, Deeg DJH, Twisk JWR, Beekman ATF, Van Tilburg W. Religious involvement and 6-year course of depressive symptoms in older Dutch citizens: results from the Longitudinal Aging Study Amsterdam. J Aging Health. 2004;16(4):467–89. pmid:15271266
  7. 7. Chen Y, Kim ES, VanderWeele TJ. Religious-service attendance and subsequent health and well-being throughout adulthood: evidence from three prospective cohorts. Int J Epidemiol. 2020;49(6):2030–40.
  8. 8. Chen Y, Koh HK, Kawachi I, Botticelli M, VanderWeele TJ. Religious service attendance and deaths related to drugs, alcohol, and suicide among US health care professionals. JAMA Psychiatry. 2020;77(7):737–44. pmid:32374360
  9. 9. Chen Y, VanderWeele TJ. Associations of religious upbringing with subsequent health and well-being from adolescence to young adulthood: an outcome-wide analysis. Am J Epidemiol. 2018;187(11):2355–64. pmid:30215663
  10. 10. Hummer RA, Rogers RG, Nam CB, Ellison CG. Religious involvement and U.S. adult mortality. Demography. 1999;36(2):273–85. pmid:10332617
  11. 11. Li S, Stampfer MJ, Williams DR, VanderWeele TJ. Association of religious service attendance with mortality among women. JAMA Intern Med. 2016;176(6):777–85. pmid:27183175
  12. 12. Pawlikowski J, Białowolski P, Węziak-Białowolska D, VanderWeele TJ. Religious service attendance, health behaviors and well-being—an outcome-wide longitudinal analysis. Eur J Public Health. 2019;29(6):1177–83. pmid:31038170
  13. 13. Strawbridge WJ, Shema SJ, Cohen RD, Kaplan GA. Religious attendance increases survival by improving and maintaining good health behaviors, mental health, and social relationships. Ann Behav Med Publ Soc Behav Med. 2001;23(1):68–74. pmid:11302358
  14. 14. Kim YI, Jang SJ. Religious service attendance and volunteering: A growth curve analysis. Nonprofit Volunt Sect Q. 2016;46(2):395–418.
  15. 15. Lunn J, Klay R, Douglass A. Relationships among giving, church attendance, and religious belief: the case of the Presbyterian Church (USA). J Sci Study Relig. 2001;40(4):765–75.
  16. 16. Lyons M, Nivison-Smith I. Religion and giving in Australia. Aust J Soc Issues. 2006;41(4):419–36.
  17. 17. Musick MA, Wilson J. Volunteering and depression: the role of psychological and social resources in different age groups. Soc Sci Med. 2003;56(2):259–69. pmid:12473312
  18. 18. Putnam RD, Campbell DE. American grace: how religion divides and unites us. Simon & Schuster; 2012.
  19. 19. Thoits PA, Hewitt LN. Volunteer work and well-being. J Health Soc Behav. 2001;42(2):115–31. pmid:11467248
  20. 20. Musick MA, House JS, Williams DR. Attendance at religious services and mortality in a national sample. J Health Soc Behav. 2004;45(2):198–213. pmid:15305760
  21. 21. VanderWeele TJ, Yu J, Cozier YC, Wise L, Argentieri MA, Rosenberg L, et al. Attendance at religious services, prayer, religious coping, and religious/spiritual identity as predictors of all-cause mortality in the Black Women’s Health Study. Am J Epidemiol. 2017;185(7):515–22. pmid:28338863
  22. 22. Yeager DM, Glei DA, Au M, Lin HS, Sloan RP, Weinstein M. Religious involvement and health outcomes among older persons in Taiwan. Soc Sci Med. 2006;63(8):2228–41. pmid:16797809
  23. 23. la Cour P, Avlund K, Schultz-Larsen K. Religion and survival in a secular region: a twenty year follow-up of 734 Danish adults born in 1914. Soc Sci Med. 2006;62(1):157–64.
  24. 24. Litwin H. What really matters in the social network–mortality association? a multivariate examination among older Jewish-Israelis. Eur J Ageing. 2007;4(2):71–82. pmid:28794773
  25. 25. Brenner PS. Cross-national trends in religious service attendance. Public Opin Q. 2016;80(2):563–83. pmid:27274579
  26. 26. Twenge JM, Sherman RA, Exline JJ, Grubbs JB. Declines in American adults’ religious participation and beliefs, 1972–2014. SAGE Open. 2016;6(1):1–13.
  27. 27. Jones JM. U.S. church membership falls below majority for first time. Gallup [Internet]. 2021 Mar 29 [cited 2022 Mar 15]; Available from: https://news.gallup.com/poll/341963/church-membership-falls-below-majority-first-time.aspx
  28. 28. Taylor C. A secular age. Harvard University Press; 2007.
  29. 29. Zuckerman P. Society without God: what the least religious nations can tell us about contentment. New York University Press; 2008.
  30. 30. Case A, Deaton A. Deaths of despair and the future of capitalism. Princeton University Press; 2020.
  31. 31. Bottan NL, Perez-Truglia R. Losing my religion: the effects of religious scandals on religious participation and charitable giving. J Public Econ. 2015;129:106–19.
  32. 32. Manning LK, Miles A. Examining the effects of religious attendance on resilience for older adults. J Relig Health. 2018;57(1):191–208. pmid:28744592
  33. 33. Maselko J, Hayward RD, Hanlon A, Buka S, Meador K. Religious service attendance and major depression: a case of reverse causality? Am J Epidemiol. 2012;175(6):576–83. pmid:22350581
  34. 34. VanderWeele TJ. Effects of religious service attendance and religious importance on depression: examining the meta-analytic evidence. Int J Psychol Relig. 2021;31(1):21–6.
  35. 35. Balbuena L, Baetz M, Bowen R. Religious attendance after elevated depressive symptoms: is selection bias at work? PeerJ. 2014;2:e311–e311. pmid:24688885
  36. 36. Ferraro KF, Kelley-Moore JA. Religious seeking among affiliates and non-affiliates: do mental and physical health problems spur religious coping? Rev Relig Res. 2001;42(3):229–51.
  37. 37. Horowitz JL, Garber J. Relation of intelligence and religiosity to depressive disorders in offspring of depressed and nondepressed mothers. J Am Acad Child Adolesc Psychiatry. 2003;42(5):578–86. pmid:12707562
  38. 38. Miller L, Weissman M, Gur M, Greenwald S. Adult religiousness and history of childhood depression: eleven-year follow-up study. J Nerv Ment Dis. 2002;190(2):86–93. pmid:11889361
  39. 39. Schleifer C, Chaves M. Family formation and religious service attendance: untangling marital and parental effects. Sociol Methods Res. 2017;46(1):125–52.
  40. 40. Storm I. Does security increase secularity? evidence from the British Household Panel Survey on the relationship between income and religious service attendance. J Relig Eur. 2017;10(3):328–49.
  41. 41. VanderWeele TJ, Mathur MB, Chen Y. Outcome-wide longitudinal designs for causal inference: a new template for empirical studies. Stat Sci. 2020;35(3):437–66.
  42. 42. VanderWeele TJ, Jackson JW, Li S. Causal inference and longitudinal data: a case study of religion and mental health. Soc Psychiatry Psychiatr Epidemiol. 2016;51(11):1457–66. pmid:27631394
  43. 43. Smith J, Ryan L, Fisher GG, Sonnega A, Weir D. HRS psychosocial and lifestyle questionnaire 2006–2016 [Internet]. Ann Arbor, Michigan: Survey Research Center, University of Michigan; 2017. Available from: https://hrs.isr.umich.edu/sites/default/files/biblio/HRS%202006-2016%20SAQ%20Documentation_07.06.17_0.pdf
  44. 44. Sonnega A, Faul JD, Ofstedal MB, Langa KM, Phillips JW, Weir DR. Cohort profile: the Health and Retirement Study (HRS). Int J Epidemiol. 2014;43(2):576–85. pmid:24671021
  45. 45. Fisher GG, Faul JD, Weir DR, Wallace RB. Documentation of chronic disease measures in the Health and Retirement Study (HRS/AHEAD) [Internet]. Ann Arbor, Michigan: Survey Research Center, University of Michigan; 2005. Available from: https://hrs.isr.umich.edu/sites/default/files/biblio/dr-009.pdf
  46. 46. Jenkins KR, Ofstedal MB, Weir DR. Documentation of health behaviors and risk factors measured in the Health and Retirement Study (HRS/AHEAD) [Internet]. Ann Arbor, Michigan: Survey Research Center, University of Michigan; 2008. Available from: https://hrs.isr.umich.edu/sites/default/files/biblio/dr-010.pdf
  47. 47. Groenwold RHH, Donders ART, Roes KCB, Harrell FEJ, Moons KGM. Dealing with missing outcome data in randomized trials and observational studies. Am J Epidemiol. 2012;175(3):210–7. pmid:22262640
  48. 48. Sterne JAC, White IR, Carlin JB, Spratt M, Royston P, Kenward MG, et al. Multiple imputation for missing data in epidemiological and clinical research: potential and pitfalls. BMJ. 2009;338:b2393. pmid:19564179
  49. 49. van Ginkel JR, Linting M, Rippe RCA, van der Voort A. Rebutting existing misconceptions about multiple imputation as a method for handling missing data. J Pers Assess. 2020;102(3):297–308. pmid:30657714
  50. 50. Harel O, Mitchell EM, Perkins NJ, Cole SR, Tchetgen Tchetgen EJ, Sun B, et al. Multiple imputation for incomplete data in epidemiologic studies. Am J Epidemiol. 2018;187(3):576–84. pmid:29165547
  51. 51. VanderWeele TJ, Ding P. Sensitivity analysis in observational research: introducing the E-value. Ann Intern Med. 2017;167(4):268–74. pmid:28693043
  52. 52. Idler EL, Kasl SV. Religion among disabled and nondisabled persons II: attendance at religious services as a predictor of the course of disability. J Gerontol Ser B. 1997;52B(6):S306–16. pmid:9403524
  53. 53. Roff LL, Klemmack DL, Simon C, Cho GW, Parker MW, Koenig HG, et al. Functional limitations and religious service attendance among African American and White older adults. Health Soc Work. 2006;31(4):246–55. pmid:17176972
  54. 54. Pargament KI. The psychology of religion and coping: Theory, research, practice. Guilford Press; 1997.
  55. 55. Kent BV, Cowden RG, Counted V, Davis EB, Rueger SY, Worthington EL Jr. Do religious/spiritual resources moderate the association between suffering and religious/spiritual struggles? A three-wave longitudinal study of US adults with chronic illness. J Sci Study Relig. 2022;Advance online publication.
  56. 56. Alcorn SR, Balboni MJ, Prigerson HG, Reynolds A, Phelps AC, Wright AA, et al. “If God wanted me yesterday, I wouldn’t be here today”: religious and spiritual themes in patients’ experiences of advanced cancer. J Palliat Med. 2010;13(5):581–8. pmid:20408763
  57. 57. Awang H, Mansor N, Nai Peng T, Nik Osman NA. Understanding ageing: fear of chronic diseases later in life. J Int Med Res. 2017;46(1):175–84. pmid:28760083
  58. 58. Fredrickson BL, Cohn MA. Positive emotions. In: Lewis M, Haviland-Jones JM, Barrett LF, editors. Handbook of emotions. 3rd ed. Guilford Press; 2008. p. 777–96.
  59. 59. Kim ES, Shiba K, Boehm JK, Kubzansky LD. Sense of purpose in life and five health behaviors in older adults. Prev Med. 2020;139:106172. pmid:32593729
  60. 60. Cowden RG, Counted V, Ho MY. Positive psychology and religion/spirituality across cultures in Africa, Asia, and Oceania. In: Davis EB, Worthington EL Jr, Schnitker SA, editors. Handbook of positive psychology, religion, and spirituality. Springer; in press.
  61. 61. Shiba K, Cowden RG, Gonzalez N, Ransome Y, Nakagomi A, Chen Y, et al. Associations of online religious participation during COVID-19 lockdown with subsequent health and well-being among UK adults. Psychol Med. 2022;Advance online publication. pmid:35189993
  62. 62. Brooks D. The road to character. Random House; 2015.
  63. 63. Jackson BR, Bergeman CS. How does religiosity enhance well-being? the role of perceived control. Psychol Relig Spiritual. 2011;3(2):149–61. pmid:24860640
  64. 64. Waters EK, Millard HM, Doyle Z. Losing hope: mental health and religious service non-attendance in Australia. Ment Health Relig Cult. 2015;18(2):114–22.
  65. 65. Long KNG, Kim ES, Chen Y, Wilson MF, Worthington EL Jr, VanderWeele TJ. The role of hope in subsequent health and well-being for older adults: an outcome-wide longitudinal approach. Glob Epidemiol. 2020;2:100018.
  66. 66. Hernandez SC, Overholser JC. A systematic review of interventions for hope/hopelessness in older adults. Clin Gerontol. 2021;44(2):97–111. pmid:31913808
  67. 67. Abramson LY, Metalsky GI, Alloy LB. Hopelessness depression: a theory-based subtype of depression. Psychol Rev. 1989;96(2):358–72.
  68. 68. Haeffel GJ, Hershenberg R, Goodson JT, Hein S, Square A, Grigorenko EL, et al. The hopelessness theory of depression: clinical utility and generalizability. Cogn Ther Res. 2017;41(4):543–55.
  69. 69. Kim ES, Whillans AV, Lee MT, Chen Y, VanderWeele TJ. Volunteering and subsequent health and well-being in older adults: an outcome-wide longitudinal approach. Am J Prev Med. 2020;59(2):176–86. pmid:32536452
  70. 70. Salamon LM, Sokolowski SW, Haddock MA. Measuring the economic value of volunteer work globally: concepts, estimates, and a roadmap to the future. Ann Public Coop Econ. 2011;82(3):217–52.
  71. 71. Hill TD, Bradley CS, Dowd-Arrow B, Burdette AM. Religious attendance and the social support trajectories of older Mexican Americans. J Cross-Cult Gerontol. 2019;34(4):403–16. pmid:31701310
  72. 72. Fields NL, Adorno G, Magruder K, Parekh R, Felderhoff BJ. Age-friendly cities: the role of churches. J Relig Spiritual Aging. 2016;28(3):264–78.
  73. 73. Voytenko V, Pargament KI, Cowden RG, Lemke AW, Kurniati NMT, Bechara AO, et al. Religious coping with interpersonal hurts: psychosocial correlates of the Brief RCOPE in four non-western countries. Psychol Relig Spiritual. 2021;Advance online publication.
  74. 74. Cowden RG, Rueger SY, Davis EB, Counted V, Kent BV, Chen Y, et al. Resource loss, positive religious coping, and suffering during the COVID-19 pandemic: a prospective cohort study of US adults with chronic illness. Ment Health Relig Cult. 2022;25(3):288–304.
  75. 75. Hardie JH, Pearce LD, Denton ML. The dynamics and correlates of religious service attendance in adolescence. Youth Soc. 2016;48(2):151–75. pmid:26900186
  76. 76. Pew Research Center. The changing global religious landscape [Internet]. 2017 [cited 2021 May 21]. Available from: https://www.pewforum.org/2017/04/05/the-changing-global-religious-landscape/