Hostname: page-component-848d4c4894-4hhp2 Total loading time: 0 Render date: 2024-05-31T16:33:02.844Z Has data issue: false hasContentIssue false

Associations of online religious participation during COVID-19 lockdown with subsequent health and well-being among UK adults

Published online by Cambridge University Press:  22 February 2022

Koichiro Shiba*
Affiliation:
Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA Human Flourishing Program, Institute for Quantitative Social Science, Harvard University, Cambridge, MA, USA
Richard G. Cowden
Affiliation:
Human Flourishing Program, Institute for Quantitative Social Science, Harvard University, Cambridge, MA, USA
Natasha Gonzalez
Affiliation:
Independent Researcher, Madrid, Spain
Yusuf Ransome
Affiliation:
Department of Social and Behavioral Sciences, Yale University School of Public Health, New Haven, CT, USA
Atsushi Nakagomi
Affiliation:
Department of Social Preventive Medical Sciences, Center for Preventive Medical Sciences, Chiba University, Chiba, Japan
Ying Chen
Affiliation:
Human Flourishing Program, Institute for Quantitative Social Science, Harvard University, Cambridge, MA, USA
Matthew T. Lee
Affiliation:
Human Flourishing Program, Institute for Quantitative Social Science, Harvard University, Cambridge, MA, USA
Tyler J. VanderWeele
Affiliation:
Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA, USA Human Flourishing Program, Institute for Quantitative Social Science, Harvard University, Cambridge, MA, USA Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
Daisy Fancourt
Affiliation:
Department of Behavioural Science and Health, University College London, London, UK
*
Author for correspondence: Koichiro Shiba, E-mail: shiba_k@g.harvard.edu
Rights & Permissions [Opens in a new window]

Abstract

Background

In-person religious service attendance has been linked to favorable health and well-being outcomes. However, little research has examined whether online religious participation improves these outcomes, especially when in-person attendance is suspended.

Methods

Using longitudinal data of 8951 UK adults, this study prospectively examined the association between frequency of online religious participation during the stringent lockdown in the UK (23 March –13 May 2020) and 21 indicators of psychological well-being, social well-being, pro-social/altruistic behaviors, psychological distress, and health behaviors. All analyses adjusted for baseline socio-demographic characteristics, pre-pandemic in-person religious service attendance, and prior values of the outcome variables whenever data were available. Bonferroni correction was used to correct for multiple testing.

Results

Individuals with online religious participation of ≥1/week (v. those with no participation at all) during the lockdown had a lower prevalence of thoughts of self-harm in week 20 (odds ratio 0.24; 95% CI 0.09–0.62). Online religious participation of <1/week (v. no participation) was associated with higher life satisfaction (standardized β = 0.25; 0.11–0.39) and happiness (standardized β = 0.25; 0.08–0.42). However, there was little evidence for the associations between online religious participation and all other outcomes (e.g. depressive symptoms and anxiety).

Conclusions

There was evidence that online religious participation during the lockdown was associated with some subsequent health and well-being outcomes. Future studies should examine mechanisms underlying the inconsistent results for online v. in-person religious service attendance and also use data from non-pandemic situations.

Type
Original Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2022. Published by Cambridge University Press

Introduction

The goal of public health is to promote comprehensive health, including psychological, social, and spiritual well-being (Grad, Reference Grad2002; UN General Assembly, 2015; VanderWeele, Reference VanderWeele2017a). Religion is principally oriented toward the promotion of spiritual well-being but can enhance other dimensions of health and well-being (Idler, Reference Idler2014; Ransome, Reference Ransome2020; VanderWeele, Reference VanderWeele and Balboni2017c). An integral component of many religious traditions is religious service attendance. Approximately 40% of the world's population attends religious services at least weekly (Pew, Reference Pew2017). Growing evidence suggests that religious service attendance is associated with indices of favorable health and well-being, including reduced mortality and unhealthy behaviors (e.g. smoking), lower levels of psychological distress (e.g. depressive symptoms), and higher levels of psychosocial well-being (e.g. life satisfaction, social connectedness) across adulthood (Chen, Kim, & VanderWeele, Reference Chen, Kim and VanderWeele2021; Chen, Koh, Kawachi, Botticelli, & VanderWeele, Reference Chen, Koh, Kawachi, Botticelli and VanderWeele2020; Chen & VanderWeele, Reference Chen and VanderWeele2018; Li, Stampfer, Williams, & VanderWeele, Reference Li, Stampfer, Williams and VanderWeele2016; Pawlikowski, Białowolski, Węziak-Białowolska, & VanderWeele, Reference Pawlikowski, Białowolski, Węziak-Białowolska and VanderWeele2019).

The coronavirus disease 2019 (COVID-19) pandemic disrupted traditional forms of religious engagement; in-person religious services were often included in the non-essential activities suspended by public health measures (e.g. stay-at-home orders) to limit transmission of infection (Oxholm, Rivera, Schirrman, & Hoverd, Reference Oxholm, Rivera, Schirrman and Hoverd2021; Pew, Reference Pew2020; Schuchat & CDC COVID-19 Response Team, Reference Schuchat2020). Following public health experts' recommendations, many religious communities shifted congregational activities from in-person to online (World Health Organization, 2020). Although online religious participation could have health and well-being benefits comparable to those from in-person religious service attendance as there are shared features (e.g. content of services and opportunity for reflection), other important health-promoting features of in-person religious service attendance (e.g. communal religious practices and casual communication with fellow adherents that involves bodily co-presence) may not be adequately replaced by online participation (VanderWeele, Reference VanderWeele2020). Given the mass transition of faith-based communities to online forms and the potential continuation of online or hybrid worship during and after the COVID-19 pandemic, research is needed to identify the health-related implications of online religious participation. However, no study to date has rigorously examined the potential effects of online religious participation on health and well-being outcomes using longitudinal data.

In this large prospective study of UK adults, we examined the associations of online religious participation during a period of stringent lockdown in the UK (23 March to 13 May 2020) – when in-person religious services were suspended – with an array of subsequent health and well-being outcomes including psychological well-being, social well-being, pro-social/altruistic behaviors, psychological distress, and health behaviors.

Methods

Data

We used data from the UCL COVID-19 Social Study, a large prospective panel study on psychological and social experiences of UK adults (≥18 years old) during the COVID-19 pandemic. Multiple non-probability sampling strategies were used to obtain a non-random but well-stratified sample, including those from lower socioeconomic backgrounds and vulnerable groups (e.g. older adults, adults with pre-existing mental health conditions). All data were weighted to proportions of gender, age, ethnicity, education, and country of living obtained from the Office for National Statistics. Further details on sampling and weighting are available elsewhere (CSSUserGuide, 2021; Fancourt, Steptoe, & Bu, Reference Fancourt, Steptoe and Bu2021).

The study began on 21 March 2020 and collected data weekly from online participants. For this analysis, we used data from the participants recruited from 21 to 27 March 2020 (week 1: n = 28 847). Of those, 11 424 individuals completed the religion module administered in week 14 (20–26 June 2020) and week 15 (27 June–3 July 2020), which assessed the frequency of online religious participation. We excluded participants with missing survey weights (n = 49) and those who did not participate in the week 20 survey (1–7 August 2020), from which the outcome data were drawn (n = 2424). The final analytic sample consisted of 8951 individuals. Figure 1 illustrates the timeline of data collection and key dates for the COVID-19 pandemic in the UK.

Fig. 1. Timelines of data collection and key dates for the COVID-19 pandemic in the UK.

Measures

Online religious participation

The religion module was administered in weeks 14 and 15. Participants retrospectively reported the frequency of online religious service attendance during lockdown by answering the question, ‘During lockdown have you engaged in any online or digital religious activities such as watching streamed services, watching video prayers, joining online faith discussion groups, or receiving other digital spiritual support?’ from the options: 1 = more than once a week, 2 = once a week, 3 = two or three times a month, 4 = one or more times a year, and 5 = not at all. Consistent with previous research involving in-person service attendance (Chen et al., Reference Chen, Koh, Kawachi, Botticelli and VanderWeele2020), responses were grouped into three categories: ‘not at all’, ‘<1/week’, or ‘≥1/week’.

Outcomes

To reduce the possibility of recall bias for the online religious participation variable, we used the outcome data from week 20, so there was an approximately 1-month interval between the religion module and outcome assessment. We examined 21 indicators of health and well-being as outcomes, including psychological well-being (life satisfaction, happiness, and meaning), social well-being (social support and loneliness), pro-social/altruistic behaviors (volunteering, caring, compliance with social isolation, and compliance with social distancing), psychological distress (depressive symptoms, anxiety, number of minor stressors, number of major stressors, and thoughts of self-harm), and health behaviors (unhealthy change in smoking, unhealthy change in alcohol drinking, unhealthy change in diet, gentle physical activity, high-intensity physical activity, exercising at home, and good sleep). Online Supplementary Table S1 provides further details about each outcome measurement.

Covariates

All covariates were taken from week 1 (21–27 March 2020); the strict lockdown in the UK began on 23 March 2020. These covariates included sociodemographic factors (age, gender, race/ethnicity, living alone, education, employment, key worker role, low income), number of self-reported health conditions, frequency of service attendance before the pandemic, compliance with not leaving home, social relationships (frequency of meeting up with people in usual life, number of close friends), personality (neuroticism, extraversion, openness, agreeableness, conscientiousness), and baseline health behaviors (current smoking status, number of alcoholic drinks in the past week). To reduce the possibility of reverse causation (i.e. health and well-being affecting participation in online religious activities), we also adjusted for prior values of the outcomes. Specifically, adjustments were made for all week 20 outcome measures using the data provided in week 1, except for happiness, compliance with social distancing, thoughts of self-harm, and unhealthy change in diet, which were not measured in week 1 (CSSUserGuide, 2021).

Statistical analysis

We used an outcome-wide analytic approach, which enables holistic assessment of the impact of a single exposure on a range of outcomes and has several methodological advantages (e.g. being less susceptible to p-hacking and publication bias) (VanderWeele, Reference VanderWeele2017b; VanderWeele, Mathur, & Chen, Reference VanderWeele, Mathur and Chen2020). We used separate regression models to regress outcomes on online religious participation during lockdown, adjusting for both covariates measured at week 1 and prior outcome values wherever data were available. We used different models depending on the type of outcome: (1) linear regression for continuous outcomes (life satisfaction, happiness, meaning, social support, loneliness, compliance with social isolation, depressive symptoms, anxiety, number of minor and major stressors), (2) Poisson regression for non-rare binary outcomes with a prevalence ≥10% (compliance with social distancing, no unhealthy change in drinking, no unhealthy change in diet, gentle physical activity, high-intensity physical activity, exercising at home, good sleep), and (3) logistic regression for rare binary outcomes with a prevalence <10% (volunteering, caring, no unhealthy change in smoking behaviors, thoughts of self-harm). All continuous outcomes were standardized (mean = 0, s.d. = 1), so the effect estimates can be interpreted as standard deviation changes in the outcome variable. Poisson regression models for non-rare binary outcomes estimate risk ratios, and logistic regression models for rare binary outcomes estimate odds ratios approximating risk ratios. We used Bonferroni correction to account for multiple testing.

We performed three sensitivity analyses. First, to evaluate the robustness of our effect estimates to unmeasured confounding, we calculated E values for each exposure–outcome association (VanderWeele & Ding, Reference VanderWeele and Ding2017). E values quantify the minimum strength of association on the risk ratio scale that an unmeasured confounder would need to have with both the exposure and outcome, above and beyond the adjusted covariates, to explain away the observed association. Second, we performed subgroup analysis to examine whether associations between online religious participation and outcomes differ by frequency of in-person service attendance before the COVID-19 pandemic (‘Not at all’ and ‘<1/week or more often’). Third, we examined outcomes from week 16 (one week after the religion module) to examine the robustness of the results to the timing of outcome assessment.

We used multiple imputation by chained equations to impute missing data on all variables, using the mice R package (van Buuren & Groothuis-Oudshoorn, Reference van Buuren and Groothuis-Oudshoorn2011). After generating five imputed datasets, we performed the analyses described above using each imputed dataset and combined the results across imputations. All analyses were conducted in R, version 3.6.0.

Results

Table 1 shows baseline characteristics of the weighted study sample according to the frequency of online religious participation during the lockdown. Compared to those with no participation in online religious activities (n = 7795), individuals with participation of ⩾1/week (n = 542) reported higher psychological well-being (i.e. life satisfaction and meaning) at the beginning of the lockdown (week 1), but no apparent difference in psychological distress, such as depressive symptoms. On the other hand, individuals with participation of <1/week (n = 613) reported the lowest psychological well-being (e.g. mean life satisfaction = 5.32 for participation of <1/week v. 5.49 for no participation and 5.95 for participation of ≥1/week) and the highest baseline loneliness and psychological distress (e.g. mean depressive symptoms = 7.1 for participation of <1/week v. 6.1 for both no participation and participation of ≥1/week).

Table 1. Weighted sample characteristics at baseline by categories of online religious participation during the lockdown (n = 8951)

s.d., standard deviation.

Lockdown was imposed in the UK on 23 March 2020, and we used data from week 1 of the survey (21–27 March 2020). Data were weighted to the proportions of gender, age, ethnicity, education, and country of living obtained from the Office for National Statistics. Columns do not necessarily sum up exactly to 8951 because of rounding after weighting.

Table 2 shows estimated β coefficients (continuous outcomes), risk ratios (non-rare binary outcomes), and odds ratios (rare binary outcomes) for the online religious participation categories, adjusting for covariates. In week 20, those who participated in online religious activities <1/week had higher life satisfaction (standardized β = 0.25; 95% CI 0.11–0.39) and happiness (standardized β = 0.25; 0.08–0.42) than those with no participation. These associations remained below the p = 0.05 threshold after accounting for multiple testing via Bonferroni correction. Participation of ≥1/week (v. no participation at all) was associated with decreased thoughts of self-harm (odds ratio 0.24; 0.09–0.62). There was modest evidence of associations between participation of ≥1/week and increased compliance with social distancing (risk ratio = 1.08; 1.01–1.16), no unhealthy change in alcohol drinking in the past week (risk ratio = 1.08; 1.01–1.15), and increased loneliness (standardized β = 0.14; 0.01–0.26). However, none of these associations were below the p = 0.05 after Bonferroni correction. We found little evidence of associations between online religious participation of either frequency and other health and well-being outcomes.

Table 2. Online religious participation during lockdown and subsequent health and well-being in the UK (n = 8951)a

CI, confidence interval; RR, risk ratio; OR, odds ratio.

a Online religious participation during lockdown (23 March–13 May 2020) was assessed in the religion module conducted in week 14 (20–26 June 2020) and week 15 (27 June–3 July 2020). Outcomes were assessed in week 20 (1–7 August 2020). Covariates were measured at the beginning of the lockdown (week 1, 21–27 March 2020). The analytic sample was restricted to those who had participated in the religion module (weeks 14–15) and completed the survey in both weeks 1 and 20. Multiple imputation was performed to impute missing data on the covariates and the outcomes.

b All continuous outcomes (life satisfaction, happiness, meaning, social support, loneliness, compliance with social isolation, depressive symptoms, anxiety, and number of minor and major stressors) were standardized (mean = 0, standard deviation, 1), and β was the standardized effect size. The estimates for the outcomes of volunteering, caring, no unhealthy change in smoking behaviors, and thoughts of self-harm were odds ratios estimated via weighted logistic regression; these outcomes were rare (prevalence <10%), so the odds ratios would approximate the risk ratios. The estimates for other dichotomized outcomes (compliance with social distancing, no unhealthy change in drinking, no unhealthy change in diet, gentle physical activity, high-intensity physical activity, and good sleep) were risk ratios estimated via weighted Poisson regression.

c All models were controlled for participants' sociodemographic characteristics (age, gender, race, living alone, education, employment, any key worker role, and low income), health (number of health conditions, current smoking, and number of alcohol drinks in the past week), pre-pandemic religious service attendance, not leaving home, social relationships (frequency of meeting up with people in usual life and number of close friends), personality (neuroticism, extraversion, openness, agreeableness, and conscientiousness), and the prior outcome values wherever data were available. Specifically, the adjustment of the prior values were made for the following outcomes simultaneously: psychological well-being (life satisfaction and meaning), social well-being (social support and loneliness), pro-social/altruistic behaviors (volunteering, caring, and compliance with social isolation), psychological distress (depressive symptoms, anxiety, and number of major/minor stressors), health behaviors (change in smoking, change in alcohol drinking, gentle/high-intensity physical activity, exercising at home, and sleep). Data were weighted to the proportions of gender, age, ethnicity, education, and country of living obtained from the Office for National Statistics.

*p < 0.05 before Bonferroni correction; **p < 0.01 before Bonferroni correction; ***p < 0.05 after Bonferroni correction (the p value cutoff for Bonferroni correction is p = 0.05/21 outcomes = p < 0.0024).

The calculated E values (Table 3) suggested that some observed associations between online religious participation and subsequent well-being were moderately robust to an unmeasured confounder. For example, for the association between online religious participation of <1/week and life satisfaction (standardized β = 0.25), an unmeasured confounder would need to be associated with both the exposure and outcome – above and beyond the adjusted covariates – by a risk ratio of 1.82 to fully explain away the observed association and by 1.45-fold to shift the CI to include the null value. As shown in online Supplementary Table S2, the conditional associations of the observed covariates with outcomes were generally smaller than the E values for online religious participation, even for covariates with particularly strong associations with an outcome. For example, the risk ratio for the conditional association between life satisfaction and its prior value from week 1 was 1.17, whereas the E value for online participation of <1/week was 1.82.

Table 3. Robustness to unmeasured confounding (E values) of associations between online religious participation during lockdown and subsequent health and well-being (n = 8951)

CI, confidence interval.

a See VanderWeele and Ding (Reference VanderWeele and Ding2017) for the formula for calculating E values.

b E values for effect estimates are 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, above and beyond the measured covariates, to fully explain away the observed association of each level of online religious participation during lockdown (reference: ‘Not at all’) with various outcomes as shown in Table 2.

c E values for the 95% CI limit closest to the null denote 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, above and beyond the measured covariates, to shift the 95% CI to include the null value.

Subgroup analysis stratified by pre-pandemic in-person service attendance (online Supplementary Table S3) showed similar point estimates. However, some associations became null after stratifying (e.g. online participation of <1/week and happiness among those without pre-pandemic in-person service attendance). When we assessed outcomes from an earlier time point (week 16), we found similar trends for most outcomes (relative to week 20) but no associations with life satisfaction and happiness (online Supplementary Table S4).

Discussion

In this longitudinal study of UK adults, we examined the associations of online religious participation during the initial strict lockdown due to the COVID-19 pandemic with subsequent health and well-being. Our main findings are threefold. First, the most frequent online religious participation (≥1/week) was associated with a lower prevalence of thoughts of self-harm compared to no online religious participation. Second, we found an association with higher psychological well-being (i.e. life satisfaction and happiness) for intermittent participation (<1/week), but not for most frequent participation (≥1/week), compared to no participation. Third, except for these two findings, there was little evidence that online religious participation was associated with other indices of subsequent health and well-being.

The observed association between the most frequent online participation and lower prevalence of thoughts of self-harm is consistent with prior evidence on in-person service attendance. For example, previous US-based studies found associations with lower rates of suicide and deaths of despair (deaths from suicides and drugs/alcohol use) for in-person participation of ≥1/week but not participation of <1/week (Chen et al., Reference Chen, Koh, Kawachi, Botticelli and VanderWeele2020; VanderWeele, Li, Tsai, & Kawachi, Reference VanderWeele, Li, Tsai and Kawachi2016). These associations could be explained by strong beliefs against self-harm and suicide that are part of religious values (Dervic et al., Reference Dervic, Oquendo, Grunebaum, Ellis, Burke and Mann2004; Koenig, Koenig, King, & Carson, Reference Koenig, Koenig, King and Carson2012). Many religious traditions – including Christianity, which constitutes 38% of the entire UK population and 80% of religious individuals in the UK and 83% of this study sample (Curtice, Clery, Perry, Phillips, & Rahim, Reference Curtice, Clery, Perry, Phillips and Rahim2019) – explicitly prohibit self-injury and suicide because life is considered a gift from God and worth protection and care. Scholars have also hypothesized that religious service attendance could lower suicide risk by promoting hope, purpose, and meaning (Idler, Reference Idler2014); however, our outcome-wide analysis did not provide evidence for the association between online religious participation and a single-item measure of meaning.

The inverse U-shaped relationship between frequencies of online religious participation and life satisfaction and happiness was inconsistent with evidence that the more frequent in-person service attendance is increasingly associated with greater psychological well-being (Chen et al., Reference Chen, Kim and VanderWeele2021; Chen & VanderWeele, Reference Chen and VanderWeele2018). There are at least three possible explanations for this inconsistency. First, of the individuals with the most frequent online religious participation (≥1/week), 69.3% had attended in-person religious services ≥1/week before the pandemic. These individuals with high levels of pre-pandemic religious engagement (v. individuals with less frequent or no pre-pandemic in-person service attendance) may have lost more resources due to the suspended in-person services during the lockdown, such as social support from a congregation. Online religious participation may not have sufficiently replaced resources to fulfill the increased need for those who participated most frequently during the lockdown (Holmgreen, Tirone, Gerhart, & Hobfoll, Reference Holmgreen, Tirone, Gerhart, Hobfoll, Cooper and Quick2017). This theorizing is supported by the increased loneliness among people with the most frequent online participation (Table 2). Second, as shown in Table 1, individuals with online participation of <1/week had poorer baseline well-being at the beginning of the lockdown; hence, they may have had more room for improving life satisfaction and happiness. Third, evidence suggests that moderate social contact via telephone or video may be beneficial, but high levels of virtual communication can result in greater distress (Bu, Steptoe, Mak, & Fancourt, Reference Bu, Steptoe, Mak and Fancourt2021; Twenge, Spitzberg, & Campbell, Reference Twenge, Spitzberg and Campbell2019).

Despite observed associations with decreased thoughts of self-harm and increased life satisfaction and happiness, online religious participation was not associated with other outcomes, which contrasts with the evidence linking in-person service attendance to a range of health and well-being outcomes (Chen et al., Reference Chen, Kim and VanderWeele2021; Chen & VanderWeele, Reference Chen and VanderWeele2018). We provide four potential explanations for our results. First, online communication cannot replicate some physical features of in-person social interactions (e.g. touching, hugging), which could reduce stress (Cohen, Janicki-Deverts, Turner, & Doyle, Reference Cohen, Janicki-Deverts, Turner and Doyle2015; Jakubiak & Feeney, Reference Jakubiak and Feeney2017). Second, online religious participation is unlikely to provide the same benefits as in-person services because sacred places often play a central role in inner religious/spiritual experiences (Counted, Neff, Captari, & Cowden, Reference Counted, Neff, Captari and Cowden2020; Fancourt & Steptoe, Reference Fancourt and Steptoe2019; Mazumdar & Mazumdar, Reference Mazumdar and Mazumdar2004). Third, online services may not have provided the same quality of faith teachings and communal religious/spiritual experiences as in-person services. Fourth, evidence suggests that social support tends to be less beneficial when provided by people experiencing the same stressors (Rueger, Malecki, Pyun, Aycock, & Coyle, Reference Rueger, Malecki, Pyun, Aycock and Coyle2016). The helpfulness of online religious participation in our study may have been affected by the widespread impact of the COVID-19 pandemic. Indeed, the sensitivity analyses examining outcomes from earlier time points – times when stressors were likely more severe for many because some lockdown regulations were still imposed – indicated that the associations between online religious participation and psychological well-being were null.

The frequency of pre-pandemic religious service attendance in this study population was somewhat low (22.8% reported participation), given that 38% of the UK population identified themselves as Christians. However, the proportion of people who attend religious services in the UK has been reported to be smaller than the proportion of religious individuals in the same population, according to the British Social Attitudes survey – a survey of a representative sample of adults ≥18 years old in the UK. In 2018, approximately 11% of the UK population (v. 5.7% in our sample) attended a religious service at least once a week, 20% (v. 17.1% in our sample) reported more infrequent participation, and 66% reported no attendance (v. 77.2% in our sample) (Curtice et al., Reference Curtice, Clery, Perry, Phillips and Rahim2019). Because the UK's religious population has been decreasing over time, the even more infrequent self-reported religious service attendance in our study sample may be due to the time difference between the British Social Attitudes survey (2018) and our study (2020).

There were notable changes in the frequency when comparing pre-pandemic in-person service attendance and online attendance during the lockdown. Specifically, there was an increase in never attendees (77.2% before the pandemic to 87.1% during the lockdown) and frequent attenders (5.7–6.9%), while infrequent attendees decreased (17.1–6.1%). Although speculative, there are three potential explanations. First, the shift from in-person religious services to online services can cause both an increase and decrease in the frequency of service attendance. For example, the online format might have facilitated attendance among some people by removing barriers (e.g. not having to spend time to go to churches). At the same time, some people may have become less interested in service attendance because the online format does not come with the joy of social interactions with peers that they used to gain from in-person attendance or has some technical complexities (e.g. people not knowing how to attend services via the Internet). This might explain the increase in most frequent participation and no participation and decrease in infrequent participation. Secondly, some individuals who infrequently or never attended religious services before the pandemic may have started to seek opportunities to cope with the hardship/anxiety they experienced during the pandemic, which may explain the increase in most frequent service attendance during the lockdown. This is a form of reverse causation (i.e. health/well-being leading to more religious participation rather than service attendance causally affecting people's health/well-being), which we partly addressed by adjusting for outcome values at the beginning of the lockdown in the outcome-wide regression analysis. Third, while the infrequent attendance category (<1/week) can include a wide range of frequencies (ranging from, e.g. once every other week to once per year), online attendance was evaluated for a narrower timeframe of only 2 months (i.e. the hard lockdown between March and May 2020). Therefore, the category of infrequent online attendance during the lockdown could only catch a much smaller range of frequency (the least possible frequency was once over the 2 months). In other words, many individuals who would have been categorized as infrequent attenders if a wider timeframe were given (e.g. once per year) ended up being categorized as never attenders.

We note six limitations of the current study. First, frequency of online religious participation during the strict lockdown was retrospectively reported approximately 5 weeks after the lockdown ended on 13 May 2020, making the measurement susceptible to recall bias. We addressed this issue, at least partly, by obtaining the primary outcome information from the study wave conducted one month after measuring online religious participation. Second, we cannot preclude the possibility of unmeasured confounders; however, we leveraged the panel structure of the data and adjusted for a rich set of covariates, including prior levels of outcomes and pre-pandemic in-person religious service attendance. Moreover, the sensitivity analysis suggested robustness of the observed associations to residual confounding, at least for thoughts of self-harm as an outcome. Third, the earliest experience of online religious participation (i.e. exposure) may have preceded the pandemic or the covariate measurement in the first week of the hard lockdown. If one participated in online religious service before the pandemic, there will be residual confounding because we did not measure pre-pandemic online service attendance; however, we adjusted for a series of covariates, including pre-pandemic in-person attendance, which is likely correlated with online attendance strongly. If the covariate measurement was after the exposure, the models effectively adjust for mediators rather than confounders, and our results would be conservative estimates of the total effects of online religious participation. Fourth, some associations between online religious participation and the outcomes may be due to post-lockdown in-person service attendance rather than online participation because the outcomes in the main analysis were taken after lockdown regulations were revoked almost entirely. Fifth, there was attrition from the original participants in week 1 (n = 28 847) to the analytic sample (n = 8951). We conducted an ad hoc analysis comparing weighted characteristics in week 1 of these two samples and confirmed that the two samples had comparable distributions of the observed covariates (online Supplementary Table S5). However, we cannot eliminate the possibility that unmeasured factors might be associated with both the attrition and the outcomes. If the attrition is also associated with the exposure (i.e. online religious participation), the attrition may result in selection bias and induce non-causal associations (Hernán, Hernández-Díaz, & Robins, Reference Hernán, Hernández-Díaz and Robins2004). Lastly, the generalizability and transportability of the findings are limited. The sample was from an Internet-based study and not representative of the UK general population, although we used survey weights to make the estimates applicable to the wider population. Moreover, because the sample was from the UK population, 83% of the religious participants were Christian; thus, the results may not be transportable to other populations or religious groups. The results should be interpreted in the context of the pandemic and resulting lockdown because online religious participation may have differential impacts during non-crisis times when societies as a whole experience less distress.

In conclusion, there was only limited evidence that online religious participation during the lockdown promoted health and well-being among adults in the UK. To explore whether online religious participation has the potential to promote the same health and well-being benefits as in-person participation, future studies should investigate mechanisms underlying the weaker associations with the outcomes for online religious participation compared to prior findings regarding in-person service attendance. Research with data from non-pandemic situations is also needed to examine the benefits of online religious participation under less stressful conditions.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291722000551

Acknowledgements

The researchers are grateful for the support of a number of organizations with their recruitment efforts including the UKRI Mental Health Networks, Find Out Now, UCL BioResource, HealthWise Wales, SEO Works, FieldworkHub, and Optimal Workshop.

Financial support

This work was supported by the Nuffield Foundation [grant number WEL/FR-000022583]; the MARCH Mental Health Network funded by the Cross-Disciplinary Mental Health Network Plus initiative supported by UK Research and Innovation [grant number ES/S002588/1]; Wellcome Trust [grant number 221400/Z/20/Z, 205407/Z/16/Z]; the Wellbeing for Planet Earth Foundation; and the John Templeton Foundation [grant number 61075].

Conflict of interest

None.

Ethical standards

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008.

References

Bu, F., Steptoe, A., Mak, H., & Fancourt, D. (2021). Time use and mental health in UK adults during an 11-week COVID-19 lockdown: A panel analysis. The British Journal of Psychiatry, 219(4), 551556. doi: 10.1192/bjp.2021.44.CrossRefGoogle ScholarPubMed
Chen, Y., Kim, E. S., & VanderWeele, T. J. (2021). Religious-service attendance and subsequent health and well-being throughout adulthood: Evidence from three prospective cohorts. International Journal of Epidemiology, 49(6), 20302040.CrossRefGoogle ScholarPubMed
Chen, Y., Koh, H. K., Kawachi, I., Botticelli, M., & VanderWeele, T. J. (2020). Religious service attendance and deaths related to drugs, alcohol, and suicide among US health care professionals. JAMA Psychiatry, 77(7), 737744.CrossRefGoogle ScholarPubMed
Chen, Y., & VanderWeele, T. J. (2018). Associations of religious upbringing with subsequent health and well-being from adolescence to young adulthood: An outcome-wide analysis. American Journal of Epidemiology, 187(11), 23552364.CrossRefGoogle ScholarPubMed
Cohen, S., Janicki-Deverts, D., Turner, R. B., & Doyle, W. J. (2015). Does hugging provide stress-buffering social support? A study of susceptibility to upper respiratory infection and illness. Psychological Science, 26(2), 135147.CrossRefGoogle ScholarPubMed
Counted, V., Neff, M. A., Captari, L. E., & Cowden, R. G. (2020). Transcending place attachment disruptions during a public health crisis: Spiritual struggles. Resilience, and Transformation. Journal of Psychology and Christianity, 39(4), 276287.Google Scholar
CSSUserGuide. (2021). Retrieved July 1, 2021, from https://github.com/UCL-BSH/CSSUserGuide.Google Scholar
Curtice, J., Clery, E., Perry, J., Phillips, M., & Rahim, N. (2019). British social attitudes: The 36th report, London. Retrieved from The National Centre for Social Research website https://www.bsa.natcen.ac.uk/media/39363/bsa_36.pdf.Google Scholar
Dervic, K., Oquendo, M. A., Grunebaum, M. F., Ellis, S., Burke, A. K., & Mann, J. J. (2004). Religious affiliation and suicide attempt. The American Journal of Psychiatry, 161(12), 23032308.CrossRefGoogle ScholarPubMed
Fancourt, D., & Steptoe, A. (2019). Present in body or just in mind: Differences in social presence and emotion regulation in live vs. virtual singing experiences. Frontiers in Psychology, 10, 778.CrossRefGoogle ScholarPubMed
Fancourt, D., Steptoe, A., & Bu, F. (2021). Trajectories of anxiety and depressive symptoms during enforced isolation due to COVID-19 in England: A longitudinal observational study. The Lancet. Psychiatry, 8(2), 141149.CrossRefGoogle Scholar
Grad, F. P. (2002). Review of the preamble of the constitution of the World Health Organization. Bulletin of the World Health Organization, 80(12), 981984.Google Scholar
Hernán, M. A., Hernández-Díaz, S., & Robins, J. M. (2004). A structural approach to selection bias. Epidemiology, 15(5), 615.CrossRefGoogle ScholarPubMed
Holmgreen, L., Tirone, V., Gerhart, J., & Hobfoll, S. E. (2017). Conservation of resources theory. In Cooper, C., & Quick, J. C. (Eds.), The handbook of stress and health: A guide to research and practice (pp. 443457). Cambridge, MA, USA: Elsevier Academic Press.CrossRefGoogle Scholar
Idler, E. L. (2014). Religion as a social determinant of public health. Oxford, United Kingdom: Oxford University Press.CrossRefGoogle Scholar
Jakubiak, B. K., & Feeney, B. C. (2017). Affectionate touch to promote relational, psychological, and physical well-being in adulthood: A theoretical model and review of the research. Personality and Social Psychology Review, 21(3), 228252.CrossRefGoogle ScholarPubMed
Koenig, H., Koenig, H. G., King, D., & Carson, V. B. (2012). Handbook of religion and health. New York, NY, USA: Oxford University Press.Google Scholar
Li, S., Stampfer, M. J., Williams, D. R., & VanderWeele, T. J. (2016). Association of religious service attendance with mortality among women. JAMA Internal Medicine, 176(6), 777785.CrossRefGoogle ScholarPubMed
Mazumdar, S., & Mazumdar, S. (2004). Religion and place attachment: A study of sacred places. Journal of Environmental Psychology, 24(3), 385397.CrossRefGoogle Scholar
Oxholm, T., Rivera, C., Schirrman, K., & Hoverd, W. J. (2021). New Zealand religious community responses to COVID-19 while under level 4 lockdown. Journal of Religion and Health, 60(1), 1633.CrossRefGoogle ScholarPubMed
Pawlikowski, J., Białowolski, P., Węziak-Białowolska, D., & VanderWeele, T. J. (2019). Religious service attendance, health behaviors and well-being-an outcome-wide longitudinal analysis. European Journal of Public Health, 29(6), 11771183.CrossRefGoogle ScholarPubMed
Pew, R. C. (2017). The changing global religious landscape. Retrieved April 16, 2021, from https://www.pewforum.org/2017/04/05/the-changing-global-religious-landscape/.Google Scholar
Pew, R. C. (2020). Americans oppose religious exemptions from coronavirus-related restrictions. Retrieved from https://www.pewforum.org/2020/08/07/americans-oppose-religious-exemptions-from-coronavirus-related-restrictions/.Google Scholar
Ransome, Y. (2020). Religion, spirituality, and health: New considerations for epidemiology. American Journal of Epidemiology, 189(8), 755758.CrossRefGoogle ScholarPubMed
Rueger, S. Y., Malecki, C. K., Pyun, Y., Aycock, C., & Coyle, S. (2016). A meta-analytic review of the association between perceived social support and depression in childhood and adolescence. Psychological Bulletin, 142(10), 10171067.CrossRefGoogle ScholarPubMed
Schuchat, A., & CDC COVID-19 Response Team. (2020). Public health response to the initiation and spread of pandemic COVID-19 in the United States, February 24–April 21, 2020. MMWR. Morbidity and Mortality Weekly Report, 69(18), 551556.CrossRefGoogle Scholar
Twenge, J. M., Spitzberg, B. H., & Campbell, W. K. (2019). Less in-person social interaction with peers among U.S. adolescents in the 21st century and links to loneliness. Journal of Social and Personal Relationships, 36(6), 18921913.CrossRefGoogle Scholar
UN General Assembly. (2015). Transforming our world: The 2030 agenda for sustainable development, A/RES/70/1. Retrieved from https://sustainabledevelopment.un.org/content/documents/21252030%20Agenda%20for%20Sustainable%20Development%20web.pdf.Google Scholar
van Buuren, S., & Groothuis-Oudshoorn, K. (2011). mice: Multivariate imputation by chained equations in R. Journal of Statistical Software, 45(3), 167.Google Scholar
VanderWeele, T. J. (2017a). On the promotion of human flourishing. Proceedings of the National Academy of Sciences of the USA, 114(31), 81488156.CrossRefGoogle ScholarPubMed
VanderWeele, T. J. (2017b). Outcome-wide epidemiology. Epidemiology, 28(3), 399402.CrossRefGoogle ScholarPubMed
VanderWeele, T. J. (2017c). Religion and health: A synthesis. In Balboni, M. J. (Ed.), Spirituality and religion within the culture of medicine: From evidence to practice (Vol. 419, pp. 357401). New York, NY, USA: Oxford University Press, xii.Google Scholar
VanderWeele, T. J. (2020). Love of neighbor during a pandemic: Navigating the competing goods of religious gatherings and physical health. Journal of Religion and Health, 59(5), 21962202.CrossRefGoogle ScholarPubMed
VanderWeele, T. J., & Ding, P. (2017). Sensitivity analysis in observational research: Introducing the E-value. Annals of Internal Medicine, 167(4), 268274.CrossRefGoogle ScholarPubMed
VanderWeele, T. J., Li, S., Tsai, A. C., & Kawachi, I. (2016). Association between religious service attendance and lower suicide rates among US women. JAMA Psychiatry, 73(8), 845851.CrossRefGoogle ScholarPubMed
VanderWeele, T. J., Mathur, M. B., & Chen, Y. (2020). Outcome-wide longitudinal designs for causal inference: A new template for empirical studies. Statistical Science, 35(3), 437466.CrossRefGoogle Scholar
World Health Organization. (2020). Practical considerations and recommendations for religious leaders and faith-based communities in the context of COVID-19: interim guidance, 7 April 2020. Retrieved from World Health Organization website https://apps.who.int/iris/handle/10665/331707.Google Scholar
Figure 0

Fig. 1. Timelines of data collection and key dates for the COVID-19 pandemic in the UK.

Figure 1

Table 1. Weighted sample characteristics at baseline by categories of online religious participation during the lockdown (n = 8951)

Figure 2

Table 2. Online religious participation during lockdown and subsequent health and well-being in the UK (n = 8951)a

Figure 3

Table 3. Robustness to unmeasured confounding (E values) of associations between online religious participation during lockdown and subsequent health and well-being (n = 8951)

Supplementary material: File

Shiba et al. supplementary material

Shiba et al. supplementary material

Download Shiba et al. supplementary material(File)
File 50.8 KB