| HOME | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|---|
| ||||||||||||||||||||||||||||||||
RESEARCH ARTICLE |
a Department of Epidemiology and Public Health, Yale University School of Medicine, New Haven, Connecticut
Peter H. Van Ness, Department of Epidemiology and Public Health, Yale University School of Medicine, 60 College Street, New Haven, CT 06510 E-mail: peter.vanness{at}yale.edu.
Decision Editor: Fredric D. Wolinsky, PhD
| Abstract |
|---|
|
|
|---|
Methods. The study population is a racially and religiously diverse sample of 2,812 community-dwelling men and women from New Haven, Connecticut; it was representative of New Haven in 1982, but not of the United States then or now. The primary study hypothesis proposes that attendance at religious services as measured in 1982 contributes to the prediction of lower levels of cognitive dysfunction in 1985 and 1988. A multivariate logistic regression analysis controlling for seven sociodemographic, four behavioral, and seven biomedical covariables is conducted.
Results. Study results show an inverse association between religious attendance in 1982 and cognitive dysfunction in 1985 (odds ratio = 0.64; 95% confidence interval = 0.49, 0.85). Religious attendance as measured in 1982 no longer contributes to the prediction of cognitive dysfunction in 1988.
Discussion. The short-term nature of the effect is partially explained by differential mortality. Persons with infrequent religious service attendance and with high levels of cognitive dysfunction in 1982 were more likely to die in the period from 1985 to 1988.
Religious beliefs and practices have been the subject of biostatistical and epidemiological research since the origins of these disciplines. In 1872, Francis Galton, an originator of regression analysis, published "Statistical Inquiries into the Efficacy of Prayer." He claimed to find little evidence of prayer's efficacy, for instance, in averting disease and death; his conclusions, however, are not without bias. From only anecdotal evidence, he opined that "there is reason to expect devout and superstitious men should be unreasonable" and that "praying people are not practical" (Galton 1872
). In 1897, Émile Durkheim included religious affiliation as an explanatory factor in his statistical examination of the occurrence of suicide in Western Europe. While acknowledging a protective association of Roman Catholicism and Judaism with suicide, Durkheim explained this influence by the lack of "a spirit of free inquiry" in these traditional and communal religions, and, in turn, he associated this lack with cultural "inferiority" (Durkheim 1897
, Durkheim 1966
).
In recent decades, there has been an increase in the number of epidemiological studies considering religion as an exposure factor (Koenig, McCullough, and Larson 2001
; Larson, Swyers, and McCullough 1997
; Levin and Schiller 1987
). Popular texts have claimed generally salutary health effects for religion (Benson 1997
; Koenig 1999
). High rates of religious belief and activity among elderly populations have attracted researchers' attention (Idler 1994
). Few studies, however, have explicitly looked at religion in relation to the outcomes of cognitive dysfunction and dementia. These are important health problems, especially for elderly populations. U.S. prevalence studies of dementia report rates from 4.5 to 8.2% among Americans 65 and older; rates are thought to double for approximately every 5 years of aging beyond 65 (Hendrie 1998
). Given the increasing proportion of older people in the United States, concern is growing about the prevalence of Alzheimer's disease and other forms of dementia. In fact, some commentators say that fear of Alzheimer's disease is a health risk for some older Americans (Centofani 1998
). Thus, attention to this topic is timely. By addressing the outcome of cognitive dysfunction, this prospective cohort study will fill a gap in epidemiological research regarding religion.
Two studies using the same data set as is analyzed here have motivated the hypotheses that will be tested. From both a cross-sectional and longitudinal perspective Idler and Kasl 1997a
have shown that attendance at religious services is associated with lower levels of functional disability in this elderly population (Idler and Kasl 1997b
). The current study uses very similar measures of religiousness, as did Idler and Kasl. One purpose of the study is to determine whether the association of religion and good health in this study population extends from the domain of physical functioning to cognitive functioning, and, if so, whether religious attendance, as opposed to some other measure of religiousness, is responsible for the association. A more recent longitudinal study has shown that social disengagement is a risk factor for cognitive decline in this elderly cohort; lack of attendance at religious services is one of the six components contributing to the authors' scale of social disengagement (Bassuk, Glass, and Berkman 1999
). Thus, a second purpose of the current study is to control for social engagement to establish whether or not the relation of religiousness to cognitive dysfunction is independent of a measure of social engagement.
The primary hypothesis of the study is that attendance at religious services in 1982 contributes to the prediction of levels of cognitive dysfunction in 1985 and 1988. On the basis of previous studies of this population that have claimed salutary associations for religious attendancepreeminently the two studies cited previouslythe direction of the association is proposed to be inverse: higher levels of religious attendance will be associated with lower levels of cognitive dysfunction. Religion is assumed to be a sociobehavioral factor similar to social support and education, factors that themselves are positively associated with several health benefits, including better cognitive functioning (Stephens, Kinney, Norris, and Ritchie 1987
; Weintraub, Powell, and Whitla 1994
). No specific mechanism for a possible protective impact of religious attendance on cognitive functioning is hypothesized; however, the proposal by Bassuk and colleagues 1999
that social engagement provides a "mental stimulation" that inhibits cognitive decline seems plausible and relevant. Attendance at religious services is potentially very stimulating in that it characteristically involves several of the domains that psychologists say comprise human intelligence (e.g., the verbal, musical, and emotional domains), and it engages them in both perceptual and productive modalities (e.g., worshippers listen to sermons and recite prayers, they hear anthems and sing hymns, and so forth; Gardner 1983
; Goleman 1995
; Sternberg 1985
).
A second hypothesis contends that most of the inverse association of religious attendance and cognitive dysfunction remains after controlling for social engagement of other types. If confirmed, this hypothesis provides evidence that religious attendance is a religious form of social engagement that is statistically independent of other forms of social engagement and thereby deserves the specific attention it has been given in recent social epidemiological literature. Finally, it is hypothesized that religious attendance is a more influential predictor of cognitive dysfunction than an alternative subjective measurement of religious identity; this hypothesis is motivated by the work of Idler and Kasl 1997a
who found this to be the case for the outcome of functional disability. Also, feelings of closeness and comfort with God are probably not accompanied by mental stimulation to the same extent as are attendance and participation in religious services. These feelings are rarely accompanied with vivid auditions and visions, and they are characteristically more perceptual than productive experiences.
| Methods |
|---|
|
|
|---|
Exposure
Two separate exposure measures are used. The more objective exposure variable consists of responses to a question about the frequency of the respondent's attendance at religious services. Six answers are offered ranging from "more than once a week" attendance to "never or almost never" attend religious services. Responses are dichotomized as persons attending religious services once a week or more frequently versus persons attending religious services less frequently than once a week. This point of bifurcation reflects a long-standing convention in the epidemiological literature on religion (Comstock and Partridge 1972
; Strawbridge, Cohen, Shema, and Kaplan 1997
).
The second more subjective exposure measure gauges the extent of the respondent's religious self-identity. It is composed of responses to two questions. The first offers five options ranging from being "deeply religious" to being "against religion." The second inquires about the extent to which the respondent experiences religion to be a source of strength and comfort, offering three options ranging from "a great deal" to "none." Responses to these two questions are combined into a single scale and dichotomized as those persons who understand themselves to be deeply religious and greatly comforted and strengthened by religion versus those who have less strong religious identifications and experiences. The epidemiological literature contains no established convention about where to dichotomize such responses. The dichotomization used here divides the population more evenly than alternative options.
Outcome
The primary outcome measure consists of the number of erroneous responses to all 10 questions in Pfeiffer's Short Portable Mental Status Questionnaire (SPMSQ). This questionnaire is designed to test whether respondents can orient themselves effectively in regard to their temporal, physical, and social environments (Pfeiffer 1975
). Questions ask, for instance, about the date and day of the week. Others ask about the current and most recent past President of the United States or about simple arithmetic. Refusals were counted as errors in accordance with Pfeiffer's instructions. Respondents who have missing data (i.e., who gave neither a correct answer, an incorrect answer, nor an explicit refusal to answer) for two or more of the 10 questions on the questionnaire were counted as having a missing total test score. Respondents with one or two missing answers had their number of errors mathematically adjusted to be proportionally represented on a 010 scale of errors. One motivation for this way of handling the data was that some questions did not apply to some respondents at the time of data collection; for instance, requests to recall telephone numbers did not apply to persons lacking home telephones.
Pfeiffer grouped the SPMSQ scores into four categories, ranging from "Intact Intellectual Functioning" to "Severe Intellectual Impairment." The questionnaire was validated in two elderly cohorts: one drawn from institutional living situations and the other from clinic clientele. Test-retest correlations for two groups tested were .82 and .83. There was 92% agreement between SPMSQ results and clinical diagnoses when the SPMSQ indicated high impairment and 82% agreement when mild or no impairment was indicated (Pfeiffer 1975
, p. 439).
When used in the healthier Yale Health and Aging Project cohort, researchers have changed Pfeiffer's grouping to a trichotomy in which there is an intact category (01 errors), a mild dysfunction category (23 errors), and a severe dysfunction category (410 errors; Bassuk et al. 1999
; Moritz, Kasl, and Berkman 1989
, Moritz, Kasl, and Berkman 1995
). Bassuk and colleagues 1999
report that, in one validation study, "sensitivity and specificity in identifying moderate or severe impairment were 85% and 96% respectively" (Field 1993
). In subsequent waves of data collection results from Folstein's Mini-Mental State Examination were also available as a way of confirming the results of the SPMSQ; in data analysis, answers to this examination were converted into a scale of errors to conform to the scoring of the Pfeiffer questionnaire (Folstein, Folstein, and McHugh 1975
). Proxy responses were not accepted to questions comprising the Pfeiffer and Folstein examinations.
Covariables
Covariables are of three types: sociodemographic covariables of religious denominational preference, age, gender, race, income, education, and marital status; behavioral covariables of alcohol drinking, physical activity, smoking, and social engagement; and biomedical covariables of baseline cognitive dysfunction, depression, diabetes, functional disability, heart disease, hypertension, and stroke. Linearity test showed that age was related to the outcome variable in a linear fashion and so was entered into multivariate models as a continuous variable (mean = 74.6 and standard deviation = 6.9). Education was defined by a scale ranging from 0 (no formal schooling) to 17+ (some graduate education); it also had a linear relation to the outcome and was entered into multivariate models as continuous (mean = 9.0 and standard deviation = 3.8). Gender was defined a dichotomous categorical variable coded sequentially for males and females. Race was defined in terms of White and non-White racial categories. Income was categorized into three levels. Marital status recorded whether or not a respondent was married at the time of the interview.
A respondent was considered to be an alcohol drinker if he or she drank more than 6.8 ounces of ethyl alcohol (as contained in beer, wine, or liquor) within the last month. Physical activity was measured on a dichotomized scale reflecting whether people never (1), sometimes (2), or often (3) engaged in sports, walking, gardening, or physical exercises. The resulting scale with range from 4 to 12 was dichotomized to yield two groups of most approximately equal size (i.e., at score greater or equal to 7). Smokers were categorized as never, past, and current smokers. The social engagement variable was adapted from the social disengagement scale created by Bassuk and colleagues deleting the items that are either accounted for by the religious attendance variable or are used as covariables (i.e., marital status). The scoring in this scale represents engagement rather than disengagement and more directly reflects the original interview questions.
Respondents were defined as depressed if they had a score equal to or greater than 16 on the Center for Epidemiological Studies Depression scale. The medical covariables of diabetes, heart disease, hypertension, and stroke were defined by responses to questions as to whether the respondent had ever been definitely or possibly diagnosed with these conditions. People were defined as functionally disabled if they had a disability score of 1 or more on a seven-item activities of daily living scale; respondents were counted as disabled if they were unable to perform a given activity or if they needed assistance in doing it.
Previous studies have shown that cognitive dysfunction tends to increase with age in an elderly population and decrease with level of education even though cognitive function measures like the ones mentioned previously are not primarily intended to measure factual knowledge or general intelligence (Hendrie 1998
; Pfeiffer 1975
). Both age and education have been shown to be associated with religiousness and so are potential confounders (Ehmann 1999
). Functional disability is an especially important potential confounder: in unadjusted analyses, religious service attendees as a group may have relatively little cognitive dysfunction because good physical health enables them physically to attend services more readily than others.
Statistical Analysis
Statistical analysis sought to avoid categorization of continuous variables if possible. When categorization was necessary, the least arbitrary approach was chosen. The distribution of the outcome variable of cognitive dysfunction does not follow a normal distribution. More respondents have no or few errors than many errors, and there are a considerable number of zero values (see Fig. 1). Following a strategy advocated by Chang and Pocock 2000
for "analyzing data with clumping at zero," we fit a proportional odds model (i.e., an ordinal polytomous logistic regression model) using the trichotomy cited previously, except that an additional cut point between 0 and 1 was introduced to see if 0 responses differed from others. The proportional odds assumption was not met, suggesting that increments in effect were not uniform across the scale and that linear regression was not an appropriate model for this data.
|
The selection process used in the fitting pruned multiple logistic regression models was manual backward elimination. Models with all relevant variables were fit first and then variables with the highest Wald test p values were removed one by one. Religious predictor variables were forced into the model, and covariables were retained whose observed levels of significance fell below 0.25. A modified version of the social disengagement scale of Bassuk and colleagues was included in the main models to clearly show the result of controlling for this factor.
The first longitudinal model was first fit having 1982 data for religious predictor variables and relevant covariables, and 1985 data for the cognitive dysfunction outcome. Included in this model was a variable indicating the 1982 baseline cognitive dysfunction status of all respondents. The continuous version of the cognitive dysfunction variable was used for this purpose. A second longitudinal regression model was fit to see if 1982 religiousness predicts cognitive dysfunction in 1988. This model was the same as the first longitudinal model, except that 1988 cognitive function data was used.
All bivariate and multivariate analyses were weighted to account for oversampling and nonrandom patterns of refusal to participate in the study. Also, weights were included to make the sample representative of the New Haven population in 1982. Complex sample survey methods were used for variance estimation of bivariate and multivariate models to account for the stratified and clustered character of the sample. Univariate analyses were conducted with the SAS software package, version 6.12 (SAS Institute, Cary, NC; SAS/STAT User's Guide, 1999
). To account for the stratified and clustered nature of the sample design SUDAAN software, version 7.11 (Research Triangle Institute, Research Triangle Park, NC), was used in SAS-callable form for the bivariate and multivariate models (Shah, Barnwell, and Bieler 1997
). S-Plus 2000 was used to create graphs (MathSoft, Inc., Seattle, WA; S-Plus 2000 User's Guide, 1999
).
| Results |
|---|
|
|
|---|
|
Bivariate Analyses
Bivariate analyses provide an unadjusted estimate of the magnitude and direction of the relationships between independent and dependent variables (Table 2 ). The dichotomized versions of the two main religion predictors are only modestly correlated. Idler and Kasl 1997b
(p. S309) have shown, however, that the religious attendance and religious identity variables are each well correlated with subsequent measures of the same concept. Most striking are the differences between the relationships of religious attendance and religious identity variables to the cognitive dysfunction outcomes. Religious attendance has an odds ratio (OR) less than 1 in the first two waves of data. Quite differently, the religious identity variable consistently has an OR greater than 1. The OR for religious attendance is only modestly reduced below the null value. Education, income, and the social variables of marital status and social engagement have ORs that are farther below 1, indicating stronger associations. Age, of course, has an elevated OR. Of the biomedical variables, functional disability, stroke, and depression are the greatest possible risk factors.
|
Multivariate Analyses
The two longitudinal models present rather different results. The 3-year model shows an inverse association (OR = .64; 95% confidence interval [CI] = 0.490.85) that is statistically significant. This association, however, disappears at 6 years (OR = 1.00; 95% CI = 0.711.41). Similar results occur when the Folstein measure (similarly dichotomized to yield to maximally equal groups) is used as the outcome variable instead of the Pfeiffer scale (respectively, OR = 0.72; 95% CI = 0.55 0.95 and OR = 1.11; 95% CI = 0.801.55). The data suggest that it is not a long-lasting relationship. A similar but less dramatic attenuation of relative risk is evident for the income, social engagement, and education variables. The deletion of the statistically nonsignificant variables did not much change the OR and CI of the religious attendance variable. Also, when the social engagement variable was deleted from the 19821985 model, the OR for the religious attendance variable was only slightly changed (OR = 0.62; 95% CI = 0.48 0.81). This indicates that controlling for social engagement eliminated only a very small portion of the inverse association of religious attendance with cognitive dysfunction.
The logistic regression models for the religious identity variable show a different story. When the religious attendance variable is replaced by the religious identity variable in the multivariate models reported in Table 3 , religious identity is shown to have a nonsignificantly elevated OR (respectively, OR = 1.06; 95% CI = 0.801.41 and OR = 1.13; 95% CI = 0.801.61). Only the 1982 cross-sectional model has an OR less than 1 (OR = 0.76; 95% CI = 0.561.02). Here, the addition of covariables moves the unadjusted OR for the religious identity variable from above the null value to below it. The deletion of the social engagement scale from these religious identity multivariate models introduces even less change than for the religious attendance models. Finally, when the religious identity variable is added to the 1985 religious model reported in Table 3 , the OR for the religious attendance variable is slightly reduced (OR = 0.62; 95% CI = 0.470.82) and that of the religious identity variable is slightly elevated (OR = 1.17; 95% CI = 0.881.54). The p value of the Wald test for religious identity is 0.267, and so the religious identity variable was not included in the pruned models.
|
|
|
| Discussion |
|---|
|
|
|---|
2 = 44.0, p < .0001) and among those with baseline cognitive dysfunction (by 1988, 44.1% were dead, whereas only 29.0% of persons without dysfunction died,
2 = 67.1, p < .0001). Thus, those infrequent religious attendees in 1982 who became cognitively impaired by 1985 were more likely to die by 1988. If they had survived then, the 1982 baseline religious attendance could have shown an association through 1988. This interpretation is particularly appropriate when the cut point is set at higher values such as 3+ or 4+ errors, because among those thus defined as cognitively impaired the proportion that dies by 1988 is greater (respectively, 49.4% and 53.5%). The pattern of increasing ORs for higher value cut points in the 1988 model supports this point.
The religious identity variable contributes little to the prediction of levels of cognitive dysfunction. This confirms the secondary hypothesis that it is a less influential predictor of cognitive dysfunction than religious attendance. These results agree well with the two studies of the same data set that motivated this analysis. For instance, Idler and Kasl 1997a
, Idler and Kasl 1997b
found that religious attendance, but not religious identity, made a significant contribution to the prediction of future levels of functional disability.
Implications
Idler and Kasl 1997b
argue that higher religious attendance predicts lower levels of functional disability. Their basic argument can be extended to the relationship between religious attendance and future levels of cognitive dysfunction. The longitudinal and multivariate regression model results reported here are not supportive of the explanation that religious attendance "selects" lower levels of cognitive dysfunction in the sense that people with little cognitive dysfunction attend religious services more frequently than others because of their higher level of functioning. In this longitudinal study, the measurement of religious attendance in 1982 occurred before the measurement of the cognitive dysfunction outcome in 1985 and so this later level of cognitive dysfunction could not be influencing the previously measured level of religious attendance; also, in the multivariate analysis, baseline cognitive dysfunction was included as a control variable. Results are, however, compatible with the claim that frequent participation in religious services promotes lower levels of future cognitive dysfunction.
Although Idler and Kasl 1997b
(p. S315) argue for an "independent role" for religious attendance in predicting level of physical disability, Bassuk and colleagues 1999
make a somewhat different claim. Speaking of the six types of social engagement contributing to their index (and of religious attendance as one of the six), they emphasize that "no particular type of social contact is essential." Any social contact that elicits a palpable degree of "mental stimulation" may inhibit cognitive decline in old age, "possibly by maintaining a critical density of neocortical synapses" (Bassuk et al. 1999
). After controlling for social engagement, the inverse association of religious attendance and cognitive dysfunction remains. Although there is no direct evidence for it from this study beyond the meaning of statistically controlling for a variable, we believe that religious attendance is a religious form of social engagement that makes a contribution that is unique in the sense that its inverse association with cognitive functioning would remain even as one increases social engagement in other, nonreligious domains of social life. We cannot reject, however, either with additional empirical results or additional statistical controls, the argument that religious attendance is merely a form of social engagement. That is an interpretation of the very meaning of what the religious attendance variable measures. Now, and in these data, it is not a testable proposition.
The implications of this study for the understanding of religiousness are several. First, this study demonstrates that religiousness is not an undifferentiated single phenomenon. The dichotomized religious attendance and religious identity variables are not highly correlated, and they yield quite divergent ORs in logistic regression models of their relationship with cognitive dysfunction. In two further respects, they are quite different. The OR of the religious attendance variable is farther below 1 in the 19821985 model than in the 1982 cross-sectional model; for the religious identity model, just the opposite is true. Also, when the religious attendance and religious identity variables are put in the same logistic regression model, the religious attendance OR is slightly lower than when alone, whereas the religious identity OR is higher than when in a similar regression model alone. These results support the need to develop multidimensional instruments for measuring spirituality and religiousness (Fetzer Institute/National Institute of Aging Working 1999
). A second implication of this study is that the magnitude of the associations of religiousness with health outcomes may be sensitive to demographic characteristics of sex, race, age, and income.
Limitations
Concepts can be informative only to the extent that they are measured precisely and comprehensively. The measurement of religiousness in this study addressed only two dimensions of religiousness, and whereas attendance at religious services was measured, the nature of the respondents' participation in the services attended was not. Also, the measurement of cognitive dysfunction was made with screening instruments that are not capable of diagnostic sensitivity. Thus, the relationship between religious attendance and cognitive dysfunction was not fully specifiable. Even though the study was longitudinal and religious attendance was measured before cognitive dysfunction, one cannot on this basis alone infer a causal relationship from their inverse association. Finally, the results of the study are limited to the population studied and cannot be readily generalizable to other settings and populations.
| Acknowledgments |
|---|
We acknowledge with appreciation the assistance of Jerome Abramson, Peter Charpentier, Evelyne Gahbauer, Thomas Gill, Christian Miner, and Viola Vaccarino.
Received for publication January 23, 2001. Accepted for publication August 30, 2001.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. D. Hill Religious Involvement and Healthy Cognitive Aging: Patterns, Explanations, and Future Directions J. Gerontol. A Biol. Sci. Med. Sci., May 1, 2008; 63(5): 478 - 479. [Full Text] [PDF] |
||||
![]() |
T. D. Hill, A. M. Burdette, J. L. Angel, and R. J. Angel Religious Attendance and Cognitive Functioning Among Older Mexican Americans J. Gerontol. B. Psychol. Sci. Soc. Sci., January 1, 2006; 61(1): P3 - P9. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||
| HOME | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|---|