|
|
||||||||
1From the Odum Institute, University of North Carolina, Chapel Hill, North Carolina; and the 2Department of Epidemiology and Public Health, and the 3Bascom Palmer Eye Institute, University of Miami School of Medicine, Miami, Florida.
| Abstract |
|---|
|
|
|---|
METHODS. The National Health Interview Survey (NHIS) is a population-based annual survey designed to be representative of the U.S. civilian noninstitutionalized population. The National Death Index of 135,581 NHIS adult participants, 18 years of age and older, from 1986 to 1996 provided the mortality linkage through 2002. A generalized linear structural equation model (GSEM) with latent variable was used to estimate the results of a system of equations with various outcomes. Standard errors and test statistics were corrected for weighting, clustering, and stratification.
RESULTS. VI affects mortality, when direct adjustment was made for the covariates. Severe VI increases the hazard rate by a factor of 1.28 (95% CI: 1.07–1.53) compared with no VI, and some VI increases the hazard by a factor of 1.13 (95% CI: 1.07–1.20). VI also affects mortality indirectly through self-rated health and disability. The total effects (direct effects plus mediated effects) on the hazard of mortality of severe VI and some VI relative to no VI are hazard ratio (HR) 1.54 (95% CI: 1.28–1.86) and HR 1.23 (95% CI: 1.16–1.31), respectively.
CONCLUSIONS. In addition to the direct link between VI and mortality, the effects of VI on general health and disability contribute to an increased risk of death. Ignoring the latter may lead to an underestimation of the substantive impact of VI on mortality.
One challenge faced by investigators studying the relationship between VI and mortality is the selection of covariates, given that VI broadly affects psychosocial functioning. For example, a relationship has been found between VI and self-rated health,5 21 32 which in turn has been a predictor of mortality in some studies.33 34 The inclusion of self-rated health as a covariate in mortality models may therefore lead to an underestimate of the mortality risk attributable to VI.
The purpose of this study was to examine the direct and indirect (mediated) effects of VI on mortality risk while controlling for the presence of other chronic conditions and other health indicators. Using structural equation modeling and controlling for relevant covariates such as age, sex, and educational attainment, we estimated direct effects on mortality risk and tested the hypotheses that VI has an indirect adverse effect on mortality via self-rated health and disability (see Fig. 1 ).
|
| Materials and Methods |
|---|
|
|
|---|
Variables
Participants were asked to indicate whether they were blind in one or both eyes or had any other trouble seeing with one or both eyes, even when wearing glasses. Visual impairment (VI) was subsequently coded as (1) severe VI: blind in both eyes; (2) some VI: VI in both eyes, blind in one eye, visually impaired in the other eye, or blind or visually impaired in one eye only, with the other eye having good vision or not mentioned; and (3) no VI: no impairment in either eye. For some respondents, proxy reporters answered the VI questions.
Disability was measured as a latent variable with two indicators: (1) the number of days in bed in the past 2 weeks and (2) the number of reduced-activity days in the past 2 weeks. Self-reported health status was measured on a five-point ordinal scale: (1) excellent, (2) very good, (3) good, (4) fair, and (5) poor.
Covariates include age in years, sex (female reference), three-category racial identity (white [reference], black, and other races), marital status (not married [reference] and married), education level (less than high school [reference], high school, and more than high school), and the number of nonocular health conditions. Respondents listed their health conditions, which were subsequently classified according to an International Classification of Diseases 9th revision (ICD-9) code based on follow-up questions detailing the conditions. A three-category nonocular health condition variable was then created as follows: no conditions (reference), one condition, and two or more conditions.
Analysis
Structural equation modeling (SEM) is a general modeling framework that incorporates many common statistical methods, including regression, analysis of variance (ANOVA), confirmatory factor analysis, and simultaneous equations.35 SEM has advanced to include the simultaneous estimation of generalized linear equations (GSEM), such as logistic, Poisson, and Cox proportional hazards regression models.36 GSEM offers several advantages over traditional analytic methodology. First, as demonstrated in the present analysis, it allows for the estimation of multiple equations simultaneously, so that associations between multiple predictor and outcome variables can be assessed in the same model—even when the distribution of outcome measures varies from dichotomous (e.g., VI), to ordinal (e.g., self-rated health), to Poisson (e.g., disability indicators), to continuous, to time-dependent (e.g., mortality) events. Second, constructs such as disability can be estimated net of random measurement error. Theoretically meaningful constructs can be developed by using latent variables (or factors), which are unobserved variables that are indirectly measured by multiple observed variables through a system of equations. Heuristically, the multiple observed variables are optimally combined into a composite representing the latent variable. This method has great potential for improving the measurement quality of health data collected in surveys, given the ability to adjust for random measurement error. Third, GSEM provides a powerful tool for the assessment of mediation effects (indirect pathways). Mediation is estimated and tested in a single step, with potentially more statistical power than traditional multistep methods.37 Finally, GSEM software (e.g., M-Plus38 ) can incorporate sample weights and the complex sample survey design (clustering and stratification) into the analysis. This advance permits the appropriate application of GSEM to complex sample survey data including the NHIS.
The GSEM depicted in Figure 1 was fit to the data with the aforementioned additional covariate controls.36 The equations for the VI and self-rated health outcomes are logistic regressions, that for the disability latent variable outcome is ordinary linear regression, and that for the mortality outcome is a Cox proportional hazards regression. The disability latent factor is continuous, but the paths linking disability to the variables days in bed and days of restricted-activity are Poisson. The disability latent variable combines these variables while removing the random measurement error from each. All equations are estimated simultaneously by using a weighted maximum-likelihood estimator,39 with standard errors that are corrected for all features of the complex sample features of the NHIS.40
First, the model was estimated without the mortality outcome, to obtain traditional SEM estimated fit statistics for the model without mortality (model 1).35 41 42 Second, the mortality outcome was added, and estimates were provided from this model (model 2). Third, model 2 was re-estimated, treating the ordinal self-rated health variable as continuous to calculate indirect effects on mortality through self-rated health. Results from this model were identical with the results in model 2, in which self-rated health was treated as ordinal with the exception of a slight change in the effects of nonocular conditions on disability and the change from logit to linear parameters for the self-rated health outcome. Finally, the mortality outcome is rescaled from time intervals in days to time intervals in years, to obtain baseline hazards estimates (model 3). Results from model 3 were the same as those from model 2 for all outcomes except mortality. Results from the mortality equations were substantively the same in size and significance of effects.
Indirect or mediation effects are calculated by multiplying the two parameters involved in the mediation relationship.35 43 For example, to obtain the effect of VI on mortality through the disability mechanism, the raw coefficients for the effect of VI on disability are multiplied by the effect of disability on mortality. The new parameter is exponentiated to obtain the hazard ratio (HR) for the indirect effect. Total effects are the effects of the independent variable (VI) on mortality through all pathways, including the direct and all mediation pathways after adjustment for model covariates. Total effects are calculated by summing the raw coefficients for the direct effects together with the indirect effects. Standard errors for both indirect and total effects were obtained by using the delta method.44
Descriptive and model-based analyses were completed with adjustments for sample weights and design effects (SUDAAN 9.0, 200445 ; Research Triangle Institute, Research Triangle Park, NC, and M-Plus 4.2138 statistical packages, respectively).
| Results |
|---|
|
|
|---|
The SEM model without the mortality linkage (model 1) had an excellent model fit. The comparative fit index (CFI) and Tucker-Lewis index (TLI) are 1.00 and 0.998, respectively, where a value of 1.00 indicates excellent fit.35 The root-mean-square-error of approximation (RMSEA) is 0.006, where 0.00 indicates excellent fit.46 47 The
2 test of model fit is significant, indicating a poor fit (52.960, df = 8; P < 0.001), where a nonsignificant
2 is desired. However, with more than 135,000 cases, the power to detect minor departures from the data is high.48
Reliability estimates for the disability indicators are obtained from model 1. Both indicators are very reliable, with estimates of 0.94 and 0.99 for the number of bed days and number of restricted-activity days, respectively.
Tables 1 2 3 present the estimates when the mortality outcome is added to the model (model 2). Table 1 gives the direct effects of all variables on mortality. Both some and severe VI significantly affects mortality with HRs 1.13 (95% CI: 1.07–1.20) and 1.28 (95% CI: 1.07–1.53), respectively. Note that the direct effects from the GSEM are very similar to the effects one would obtain from a standard, single-equation Cox regression model that includes self-rated health and disability as independent control variables, which are HR 1.14 (95% CI: 1.07–1.21) and HR 1.30 (95% CI: 1.09–1.56) for some VI and severe VI, respectively.
|
|
|
Table 2 presents the indirect effects of VI on mortality through the self-rated health variable and the disability variable, after adjustment for model covariates. The indirect effects are calculated as the product of the two direct effects involved in the mediation relationship. For example, the log-hazard effect of VI on disability multiplied by the log-hazard effect of disability on mortality represents the indirect log-hazard effect of VI on mortality through disability. Table 2 also gives the total effect of VI on mortality, by combining the direct effect and both of the mediation (indirect) effects, where total effects are calculated as the sum of the log-hazard direct and indirect effects.
The indirect effect of VI on mortality through self-rated health increases the hazard by a factor of 1.07 (95% CI: 1.06–1.08) and severe VI increases the hazard by a factor of 1.17 (95% CI: 1.13–1.21). The total effects of VI on mortality are the direct effects plus the mediation effects through both self-rated health and disability. The total effects are: HR 1.23 (95% CI: 1.16–1.31) for some VI and HR 1.54 (95% CI: 1.28–1.86) for severe VI. The total effects can be compared to the direct effects only, which are relatively smaller and are comparable to what is obtained in a traditional Cox proportional hazards model.
Another way to view the degree to which VI impacts mortality through multiple pathways is to consider the difference in the baseline hazard functions for the three levels of VI. Figure 2A portrays the line graph of the baseline hazard function for a standard Cox proportional hazards model (albeit with a latent variable predictor) controlling for covariates. Figure 2B portrays the line graph of the baseline hazard function for the full GSEM model accounting for both the direct and indirect impacts of VI and adjusting for the covariates. The baseline represents the hazard function for a 50-year-old, white, married male in good health, with no nonocular conditions and 1 standard deviation below average disability. Comparison of Figures 2A and 2B reveals that the impact of VI on all-cause mortality is much greater when the estimated effects of VI operating through the health and disability mediators are considered.
|
| Discussion |
|---|
|
|
|---|
GSEM also allows for the measurement of constructs by using latent variables (or factors). Latent variables lead to testable hypotheses that help in the understanding of associations between VI and mortality. For example, the concept of frailty, loosely defined as aging-associated declines in organ systems, is associated with an increased risk of mortality.49 50 Presently, no standardized measure of frailty exists that is typically operationalized as an index including components such as low grip strength, low levels of physical activity, and low peak expiratory flow.49 50 Some investigators also include VI levels in the frailty index. These indexes have also been shown to be predictive of mortality.51 52 SEM could be used to identify a common set of measures that load optimally on a frailty construct (which could be estimated net of random measurement error) and then determine the extent to which frailty is predictive of mortality. SEM therefore, can be used by investigators to test specifically whether VI should be considered a key frailty indicator or whether VI influences mortality risk through a pathway independent of frailty.
The inclusion of frailty in mortality models that include assessment of other ocular conditions such as cataract, glaucoma, and age-related macular degeneration and mortality risk may provide a theoretically plausible explanation for these previously documented associations. Furthermore, SEM provides a very flexible analytic platform to examine simultaneously other proposed hypotheses that may account for associations between ocular diseases and mortality risk, such as the cardiovascular risk factor hypothesis, which posits that risk factors common to both cardiovascular mortality risk and the development of ocular conditions (e.g., smoking) may partially explain associations between ocular disease and increased mortality risk.53 54 55 As demonstrated in the present analysis SEM provides a flexible analytic approach to the testing of potential mediation effects such as between smoking, the development of select ocular diseases, and mortality risk.
The NHIS lacks routine assessment of important risk factors, such as depression, as well as a comprehensive assessment of smoking status, which would be an important covariate to include in any mortality analysis. Current smoking status was assessed in a subset of NHIS participants in years 1987, 1988, and 1990 to 1994 (n = 40,211).56 We reran our final SEM models in this subset and found only modest reductions in our HR estimates for VI when compared with the results based on our full sample. The total HR estimate for all pathways for some VI is reduced from 1.23 (Table 2) to 1.16 (95% CI: 1.05–1.28) when smoking is accounted for in the model. Corresponding HR estimates for severe VI are 1.54 and 1.48 (95% CI: 1.15–1.89).
Another potential limitation of the NHIS is that proxy information was gathered for household members who were not present for the interview. In this sample, 37% of the respondents were proxy reporters. We reran the full GSEM model with proxy as a covariate, and results for all the VI effects were nearly identical with those reported for this study. We also evaluated results for the proxy and self-report populations separately and found that most results were very similar. There are no statistically significant differences in the effects of VI on self-rated health and disability for the two types of reporters. However, there is a significant difference in the direct effect of severe VI on mortality for self-reporters versus proxy-reporters (z = 2.38, P = 0.017) where the effect is stronger in the case of proxy reporting, perhaps because self-reporters are somewhat more likely to report some VI (OR: 1.14; 95% CI: 1.06–1.23) and less likely to report severe VI (OR: 0.62; 95% CI: 0.47–0.83). The direct effect of some VI on mortality is similar for both types of reporters. Finally, the NHIS is a cross-sectional survey with mortality linkage. No repeat household interviews were completed that would have permitted assessment of change in VI and mortality risk.
Inherent in any self-reported health survey is the presence of misclassification of VI status. Data on the validity of the two VI questions used in the NHIS are not available; however, validation studies of other global self-reported VI measures correlate modestly with clinical indicators of VI. For example, overall VI assessed by either one or two items within the National Eye Institute Visual Function Questionnaire (NEI-VFQ) correlates significantly with clinically assessed visual acuity (range of correlations, 0.65–0.68).57 58
Despite these limitations, the NHIS has several important strengths, including a large sample size representative of the U.S. noninstitutionalized household population, excellent household response rates (95%–98%), and a mortality linkage that is 96% complete, resulting in the identification of more than 17,842 deaths in our analysis. The large sample size permitted the identification of more than 385 adults with severe VI—a sample size that is impossible to obtain in smaller population-based clinical ocular epidemiology studies.
To summarize, we found that reported VI was predictive of an increased risk of mortality even after control for sociodemographic measures, chronic conditions, self-reported health, and disability. Visual impairment had both direct and indirect effects on mortality, operating in part through influence on self-rated health and disability. Investigators interested in associations between VI and mortality are encouraged to test their models by using GSEM, which will permit a more complete examination of pathways through which VI increases mortality risk.
| Footnotes |
|---|
Submitted for publication January 2, 2008; revised March 7, 2008; accepted June 4, 2008.
Disclosure: S.L. Christ, None; D.J. Lee, None; B.L. Lam, None; D.D. Zheng, None; K.L. Arheart, None
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Sharon L. Christ, Odum Institute, 06 Manning Hall, CB #3355, UNC at Chapel Hill, Chapel Hill, NC 27599-3355; slchrist{at}email.unc.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. J. Karpa, P. Mitchell, K. Beath, E. Rochtchina, R. G. Cumming, and J. J. Wang Direct and Indirect Effects of Visual Impairment on Mortality Risk in Older Persons: The Blue Mountains Eye Study Arch Ophthalmol, October 1, 2009; 127(10): 1347 - 1353. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Factor-Litvak and A. Sher Invited Commentary: Coming Out of the Box Am. J. Epidemiol., May 15, 2009; 169(10): 1179 - 1181. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |