PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Med Care. Author manuscript; available in PMC 2010 February 22.
Published in final edited form as:
PMCID: PMC2825686
NIHMSID: NIHMS175740

Organizational Determinants of Hospital End-of-Life Treatment Intensity

Caroline Y. Lin, MD,* Max H. Farrell,* Judith R. Lave, PhD,* Derek C. Angus, MD, MPH,* and Amber E. Barnato, MD, MPH, MS*

Abstract

Background

There is substantial hospital-level variation in end-of-life (EOL) treatment intensity.

Objective

To explore the association between organizational factors and EOL treatment intensity in Pennsylvania (PA) hospitals.

Research Design

Cross-sectional mixed-mode survey of Chief Nursing Officers of PA hospitals linked to hospital-level measures of EOL treatment intensity calculated from PA Health Care Cost Containment Council (PHC4) hospital discharge data.

Hospitals

One hundred sixty-four hospitals, of which 124 (76%) responded to the survey.

Measures

The dependent variable was an index of hospital EOL treatment intensity; the independent variables included administrative data-derived structural and market characteristics and 29 survey-derived hospital or ICU programs, policies, or practices.

Results

In models restricted to independent variables drawn from administrative sources (available for all 164 hospitals), bed size (P < 0.001), proportion of admissions among black patients (P < 0.001), and county-wide hospital market competitiveness (Herfindahl-Hirschman index) (P = 0.001) were independently associated with greater EOL treatment intensity (adjusted R2 = 0.5136). In models that additionally included hospital programs, policies, and practices (available for 124 hospitals), only an ICU long length of stay review committee (P = 0.03) was independently associated with greater EOL treatment intensity (adjusted R2 = 0.5357).

Conclusions

Information about hospital and ICU programs, policies, and practices believed relevant to the treatment of patients near the end of life offers little additional explanatory power in understanding hospital-level variation in EOL treatment intensity than administratively-derived variables alone. Future studies should explore the contribution of more difficult to measure social norms in shaping hospital practice patterns.

Keywords: intensive care units, terminal care, life support care, hospitals, organization and administration

Nearly 40% of all deaths nationwide occur in the acute care setting and 20% involve the use of intensive care services.1 The use of acute care hospital and intensive care services for patients in their last 12 months of life consumes one-quarter of total Medicare spending for inpatient care.2 Wide variation in expenditures and in the use of intensive care and specific life-sustaining procedures at the end of life between hospital referral regions raises concerns about efficiency.3,4 Regional variation in end-of-life inpatient service use and expenditures have been attributed in part to physician, hospital, and ICU bed supply,3 which may translate into different local provider practice norms5 but are not necessarily reflective of local patient preferences.6 Less is known about the factors that influence variation in treatment intensity between hospitals, although bed supply and local practice norms7 undoubtedly play a role. Hospitals are perhaps more attractive than hospital referral regions for study because they are potentially accountable organizations with a more proximate locus of control over decisions to use intensive treatments.

The purpose of the current study was to explore organizational determinants of hospital-level variation in end-of-life treatment intensity. We used a newly developed index of hospital end-of-life treatment intensity that quantifies the use of ICU and life-sustaining treatments among patients at a high probability of death, rather than among those retrospectively identified as decedents,4 to address biases introduced by the decedent follow-back method.8 Further, we augmented traditional hospital characteristics—such as size, teaching status, and case mix or illness severity—with information about the hospital programs, policies, and practices that may have direct and indirect impacts on rates of ICU and life-sustaining treatment use. A better understanding of the organizational factors associated with end-of-life treatment intensity could be used to better align patient preferences with end-of-life treatment.9,10

Methods

Study Design

We conducted a retrospective analysis of hospital treatment patterns among patients admitted to Pennsylvania acute care hospitals, using Pennsylvania Health Care Cost Containment Council (PHC4) hospital discharge data between April 1, 2001 and March 31, 2005. We chose PHC4 data because it contains a predicted probability of in-hospital mortality calculated from key clinical findings abstracted from the medical chart during the first 48 hours of admission.11 We then linked hospital-level data from PHC4 to a cross-sectional survey of the respective hospital Chief Nursing Officers (CNOs), fielded between June 1, 2005 and May 31, 2006. The survey elicited information about 15 hospital-wide (eg, a palliative care consult service, a hospitalist program) and 14 ICU-specific (eg, daily multidisciplinary rounds, symptom management protocols) programs, policies, and practices that were hypothesized to directly or indirectly impact end-of-life treatment intensity.

Hospital Intensity

The primary dependent variable was an index of each hospital's “end-of-life” treatment intensity. As we describe in detail elsewhere,4 the index is a factor score of 6 case-mix standardized (observed-to-expected) treatment ratios for admissions over age 65 with a high probability of dying. We defined high probability of dying as those in the 95th percentile of predicted risk of death upon admission, corresponding to a probability of 0.21 or higher (mean = 0.41, 25th percentile = 0.26, 75th percentile = 0.51). We calculated an end-of-life intensity index for hospitals with at least 50 decedents and 50 admissions with a high probability of dying during the study period. The 6 standardized ratios included in the index were empirical Bayes' shrinkage estimates12 of: (1) the ICU admission rate, (2) ICU length of stay, (3) intubation and mechanical ventilation rate, (4) hemodialysis rate, (5) tracheostomy rate, and (6) gastrostomy rate during the 4 study years. These 6 measures were chosen using clinical judgment and empirical data reduction techniques among several intensive, life-sustaining treatments. Among acute care PA hospitals in our study, the treatment intensity index scores range from −2.08 to 3.12 (increasing values indicating higher treatment intensity), with a mean of −0.04, median of −0.03, and standard deviation of 0.99.

Hospital Characteristics

Structural and Market Characteristics

We obtained or calculated hospital characteristics during the study period from PHC4 administrative and discharge data, from Centers for Medicare and Medicaid (CMS) administrative data, and US Department of Agriculture rural-urban continuum codes. These characteristics included bed size, resident-to-bed ratio, mean illness severity, composite quality score, proportion of admissions among black patients, urbanicity, and the Herfindahl-Hirschman index (a 0–1 measure of area-level competitiveness based on countywide hospital market share; higher numbers reflect less competition [1 is an area monopoly]).13

Hospital and ICU Programs, Policies, and Practices

We identified the candidate hospital-wide and ICU programs, policies, and practices through a structured review of the literature and semi-structured interviews with over 100 informants from 14 Pennsylvania hospitals. On the basis of these semi-structured interviews, we identified the CNO/Vice President of Patient Care Services as the single best informant about these hospital programs, policies, and practices.

We asked 2 national experts in end-of-life care to review the survey instrument for face validity before field testing and revised the instrument accordingly. We then field tested the revised web-based survey with 10 Pennsylvania hospital CNOs and debriefed the 6 respondents by phone after survey completion, revising the user interface and response options in response to their feedback.

We fielded the final survey to the CNOs of all Pennsylvania acute care hospitals between June 2005 and May 2006 using a mixed-mode (self-administered web-based followed by computer-assisted telephone interview for nonresponders). Each respondent received a letter of invitation signed by one of the authors (AEB) describing the study and providing a link to the web-based survey, a username, and password. This invitation was accompanied by a letter endorsing the study signed by the Pennsylvania Secretaries of Health and Aging.

CNOs who did not complete the self-administered web-based survey after 1 round of mail invitations and 2 follow-up telephone or e-mail reminders were contacted and scheduled to complete a telephone interview. We offered no financial incentive for participation.

Statistical Analyses

We calculated descriptive statistics of the administratively-derived hospital characteristics for all hospitals in the sample, and compared these observable characteristics between survey nonresponder and responder hospitals using 2-sample Student t tests and χ2 tests as appropriate. To explore organizational determinants of end-of-life treatment intensity, we performed univariable tests of association with 2-sample Student t tests and χ2 tests on the entire hospital sample for administratively-derived independent variables (structural and market characteristics) and for the subsample of survey respondents for the survey-derived independent variables (programs, policies, practices). We then implemented 2 multivariable linear regression models including variables significant at the P < 0.1 level in univariable tests: a model restricted to the administratively-derived dependent variables available for all hospitals in the sample and a model that additionally included survey-derived variables in the subsample of survey respondents. We performed all analyses with STATA 10.0 (StataCorp, College Station, TX).

Institutional Review Board Approval

This study was approved by the University of Pittsburgh Institutional Review Board and deemed exempt from the requirement of written informed consent.

Results

Sample Characteristics

There were 188 adult acute care Pennsylvania hospitals reporting data to PHC4 between April 1, 2001 and March 31, 2005. Fourteen of these hospitals closed, became nonacute care facilities, or fully merged with other hospitals (thereby losing their unique PHC4 identifier) before fielding of the survey (June 2005–May 2006). Of the remaining 174 hospitals, 10 did not meet eligibility criteria for calculating the intensity index (having at least 50 decedents or 50 admissions with a high probability of dying). Of the remaining 164 eligible acute care Pennsylvania hospitals, 124 (76%) completed the survey. We present descriptive statistics of the overall sample, survey respondents, and nonrespondents in Table 1.

TABLE 1
Characteristics of Study Hospitals

Nonrespondent hospitals were in more competitive markets than respondent hospitals (Herfindahl-Hirschman index 0.330 vs. 0.432, respectively; P = 0.040), but were otherwise similar with respect to end-of-life treatment intensity, bed size, teaching status, measures of illness severity (mean probability of death upon admission and percentage of admissions comprised of transfers from outside facilities), mean compliance with 10 core measures of quality promulgated by the “Hospital Compare” program, proportion of admissions among black patients, and urban location (Table 1).

Hospital Characteristics and Intensity

In univariate analyses among all study hospitals (N = 164), the administratively-derived structural and market characteristics associated with greater end-of-life treatment intensity at the P < 0.1 level included bed size, resident-to-bed ratio, percent of admissions comprised of transfers from other facilities, percent of admissions among black patients, urban location, and higher market competition (Table 2).

TABLE 2
Univariable and Multivariable Associations Between Hospital Structural/Market Characteristics and “End-of-Life” Treatment Intensity*

In multivariable linear regression including only those characteristics that were significant at the P < 0.1 level in univariate testing, only bed size, percent of admissions among black patients, and higher market competition were independently associated with greater end-of-life treatment intensity (adjusted R2 for model = 0.5136). The remaining 3 variables were no longer significant: resident-to-bed ratio, percent of admissions comprised of transfers from other hospitals, and urban location (Table 2).

We present the prevalence of hospital and ICU programs, policies, and practices hypothesized to impact end-of-life treatment intensity reported among respondent hospitals in Table 3. In univariate analyses among the 124 survey respondent hospitals, the hospital and ICU programs, policies, and practices associated with greater end-of-life treatment intensity at the P < 0.1 level included having a training program in either palliative care, pastoral care, or critical care; an ethics consult service; daily multidisciplinary ICU rounds; a private ICU conference room for family meetings; an ICU long length of stay review committee; regularly scheduled ICU family meetings with an attending physician; and an ICU clinical protocol for withdrawing or withholding life-sustaining treatment (Table 4).

TABLE 3
Prevalence of Hospital and ICU Programs, Policies, and Practices (N = 124 Hospitals)
TABLE 4
Univariable Associations Between Hospital and ICU Programs, Policies, and Practices and “End-of-Life” Treatment Intensity* (N = 124 Hospitals)

In multivariable modeling among the 124 hospitals with both administratively and survey-derived data, the variables that remained independently associated with greater end-of-life treatment intensity at the P < 0.05 level included hospital bed size, percent of admissions comprised of black patients, area-level competition, and the presence of an ICU long length of stay review committee (Table 5). These variables explained just over half of the total variation in end-of-life treatment intensity (adjusted R2 0.5357).

TABLE 5
Parsimonious Multivariable Model of “End-of-Life” Treatment Intensity, Including Hospital Structural/Market Characteristics and Hospital/ICU Organizational Characteristics (N = 124 hospitals)*

Discussion

In this study of the organizational predictors of hospital end-of-life treatment intensity in 164 acute care Pennsylvania hospitals, we found that administratively-derived structural and market characteristics, including larger bed size, greater proportion of admissions among black patients, and higher area-level hospital market competition, were more powerful predictors of end-of-life treatment intensity than specific hospital programs, policies, or practices, with the exception of an ICU long length of stay review committee. Despite our effort to capture a comprehensive set of explanatory variables, almost half of hospital-level variation in end-of-life treatment intensity remained unexplained.

We found that hospitals with most of the studied programs, policies, and practices had higher mean hospital end-of-life treatment intensity indices than hospitals without such organizational initiatives, although most of these did not reach statistical significance. This paradoxical relationship between organizational initiatives that we hypothesized would be inversely related to intensity may reflect “reverse causation” in our cross-sectional analysis. Specifically, these initiatives may have been implemented in response to perceptions of higher than average use of ICU and life-sustaining treatments among patients at the end of life. Indeed, a longitudinal—not cross-sectional—design that could assess utilization changes associated with policy implementation would have been optimal. The finding that most programs, policies, and practices significantly associated with intensity in univariable tests were no longer significant in the multivariable model that included structural and market characteristics suggests that our study could have been underpowered or the associations could simply have been spurious due to confounding by structural characteristics. However, a recent national study demonstrated an inverse cross-sectional relationship between state-level penetration of one kind of program—hospital-based palliative care—and total end-of-life Medicare spending14; national hospital-level analyses are lacking.

The association between the concentration of black patients in a region,15 hospital,16,17 or nursing home18 and greater end-of-life treatment intensity is not new. It is not entirely clear whether this observation reflects black patients' preferences for these life-sustaining treatments1923 or some other phenomenon. For example, the higher observed end-of-life ICU use among blacks is attributable to black patients' disproportionate use of hospitals with greater end-of-life ICU use for both black and white patients.16,17

Our study also revealed a significant positive relationship between larger hospitals and end-of-life treatment intensity. This relationship was not attributable independently to tertiary care status (measured by transfers) or to teaching status (measured by the resident-to-bed ratio). It is possible that local practice norms play a significant role in this relationship, in addition to the ease of access to beds and other technologies in larger facilities. Area-level bed supply, as distinct from institution-specific supply, has a well-documented association with utilization and spending.2428 Variation in hospice use, on the other hand, is not explained by area-level market characteristics.29 This may be because hospice use is more dependent on patient preferences, demographics, and local social norms, and may have very low levels of penetration in certain counties.

In our study, higher market competition was associated with higher end-of-life intensity. This relationship existed even after controlling for bed size, a proxy for technological resources. It is possible that more competitive environments lead to higher intensity treatments because of the incentive to acquire and offer more technology and resources to compete for highly remunerated cases—such as transplant and coronary-artery bypass grafting. This may, in turn, result in greater use of those resources for patients at a high risk for dying.

There are several limitations to our study. The validity and reliability of the survey responses provided by a single informant, the CNO, was not assessed and is likely imperfect. Although the CNO is likely the single best informant, she/he may not have been knowledgeable about all aspects of hospital and ICU policy. Reassuringly, 84% reported that they felt “confident” or “very confident” in their answers to the survey questions.

Our study is unable to establish cause-effect relationships between organizational initiatives and end-of-life treatment intensity in this cross-sectional study. We did attempt to ascertain the direction of the relationship by asking CNOs in our survey when each initiative was implemented. However, many of the responses to timing of implementation were “don't know” (on average 50%–75% of responses), leading to inadequate sample sizes for analysis.

Likewise, the timing of the survey did not coincide exactly with the timing of the selected PHC4 discharge data. The dates of the hospital discharge data used to develop estimates of each hospital's treatment intensity (2001–2005) abutted, but did not overlap, the dates of the survey (2005–2006). Each hospital's practice pattern was stable over the 4 years of available data, but we are nonetheless unable to attribute temporal associations or direction of effect.

Although it was a strength to compare hospitals in 1 state with a uniform regulatory, legal, and financing environment, our study may have been underpowered and is not generalizable to other states. The response rate to our survey was impressive at 76%, but nonetheless reflects only 124 hospitals. And, based on the comparison of observable characteristics, survey nonrespondent hospitals were more likely to be in higher competition markets, confirming a likely response bias.

Unlike outcomes such as end-of-life treatment costs, our end-of-life treatment intensity index itself does not have a natural metric, given that it is a factor score based on 6 underlying standardized utilization ratios. It is not strictly possible to explain that a difference in the index of, for example, 1 point is directly related to X% fewer intubations or ICU admissions. Thus, it cannot be used to compare utilization of individual life-sustaining treatments across hospitals.

Another potential weakness is that our end-of-life treatment intensity index was calculated for all admissions, rather than just among ICU admissions, which may have diluted the “true” contribution of ICU-specific practices on intensity. However, because the treatments included in the index were ICU use (admission and length of stay) and those commonly delivered only in the ICU (mechanical ventilation), we believe that our hospital-specific measure is closely tied to ICU-specific treatments. However, among respondent hospitals with more than 1 ICU (N = 105), we asked the CNO to give us information about programs, policies, and practices implemented in the ICU that takes care of the “majority of chronically ill elders” in their hospitals (typically the MICU or a mixed medical-surgical ICU). We did not ask about the programs, policies, and practices in the other ICUs. For these 105 hospitals, there may be some misspecification—the hospital's overall intensity index is affected by practice in an ICU that we have no information about and which may have different practices than the “sampled” ICU.

Finally, it is always possible that unmeasured case-mix variation contributes to the higher end-of-life treatment intensity seen at larger hospitals. However, we did account for as much case-mix variation as possible by including the mean “predicted risk of death on admission” among the elderly (mean illness severity)—derived from key clinical factors on admission—as an independent variable. We did not find a relationship between the hospital's mean predicted risk of death and intensity in univariable analysis. In addition to this variable, each hospital's transfer-in rates and resident-to-bed ratio were eligible for model inclusion—both of which may measure staffing/equipment/capacity for taking care of sick patients. Although each variable was significantly associated with intensity in univariable analyses, neither was retained in the parsimonious multivariable model. However, hospital bed size and share of black patients may be capturing some of this “resource” construct.

In summary, the presence of the 29 hospital programs, policies, and practices that we assessed provide only minimal additional explanation of the large variation in hospital end-of-life treatment intensity in Pennsylvania acute care hospitals. Although several structural and market characteristics were associated with end-of-life treatment intensity, the mechanism of these effects remains opaque, and almost half of the variation remains unexplained. Future studies should explore the contribution of more difficult-to-measure social norms in shaping hospital practice patterns.

Acknowledgments

The authors thank the contributions of Joanne Lynn and Judith Nelson in reviewing the survey instrument, Pennsylvania Secretary of Aging Nora Dowd Eisenhower and Secretary of Health Calvin B. Johnson in providing a letter of support for the project, and research assistant Cleve Heller in collecting the data.

Supported by NIH grant K08 AG021921 (to A.E.B. PI), with additional support from P01 AG019783 (to Jonathan S. Skinner PI) and 1UL1 RR024153 (Steven E. Reis PI).

References

1. Angus DC, Barnato AE, Linde-Zwirble WT, et al. Use of intensive care at the end of life in the United States: an epidemiologic study. Crit Care Med. 2004;32:638–643. [PubMed]
2. Barnato AE, McClellan MB, Kagay CR, et al. Trends in inpatient treatment intensity among Medicare beneficiaries at the end of life. Health Serv Res. 2004;39:363–375. [PMC free article] [PubMed]
3. Wennberg JE, Fisher ES, Goodman DC, et al. Tracking the Care of Patients with Severe Chronic Illness: The Dartmouth Atlas of Health Care 2008. 978-0-9815862-0-5The Dartmouth Institute for Health Policy and Clinical Practice. 2008 Report No.
4. Barnato AE, Farrell MH, Chang CC, et al. Development and validation of measures of hospital “end-of-life” treatment intensity. Med Care. In press. [PMC free article] [PubMed]
5. Sirovich BE, Gottlieb DJ, Welch HG, et al. Variation in the tendency of primary care physicians to intervene. Arch Intern Med. 2005;165:2252–2256. [PubMed]
6. Barnato AE, Herndon MB, Anthony DL, et al. Patient preferences do not explain regional variation in end-of-life treatment. J Gen Intern Med. 2006;21(S4):98.
7. Barnato AE, Bost JE, Farrell MH, et al. Relationship between staff perceptions of hospital norms and hospital-level end-of-life treatment intensity. J Palliat Med. 2007;10:1093–1100. [PMC free article] [PubMed]
8. Bach PB, Schrag D, Begg CB. Resurrecting treatment histories of dead patients: a study design that should be laid to rest. JAMA. 2004;292:2765–2770. [PubMed]
9. SUPPORT Investigators. A controlled trial to improve care for seriously ill hospitalized patients. The study to understand prognoses and preferences for outcomes and risks of treatments (SUPPORT) JAMA. 1995;274:1591–1598. [PubMed]
10. Pritchard RS, Fisher ES, Teno JM, et al. Influence of patient preferences and local health system characteristics on the place of death. SUPPORT Investigators. Study to Understand Prognoses and Preferences for Risks and Outcomes Of Treatment. J Am Geriatr Soc. 1998;46:1242–1250. [PubMed]
11. Pine M, Jordan HS, Elixhauser A, et al. Enhancement of claims data to improve risk adjustment of hospital mortality. JAMA. 2007;297:71–76. [PubMed]
12. Ash AS, Shwartz M, Pekoz EA. Chapter 12: Comparing outcomes across providers. In: Iezzoni LI, editor. Risk Adjustment for Measuring Health Outcomes. Vol. 3. Chicago, IL: Academy Health; 2005. pp. 297–333.
13. Phibbs CS, Robinson JC. A variable-radius measure of local hospital market structure. Health Serv Res. 1993;28:313–324. [PMC free article] [PubMed]
14. Goldsmith B, Dietrich J, Du Q, et al. Variability in Access to Hospital Palliative Care in the United States. J Palliat Med. 2008;11:1094–1102. [PMC free article] [PubMed]
15. Baicker K, Chandra A, Skinner JS, et al. Who you are and where you live: how race and geography affect the treatment of medicare beneficiaries. Health Aff (Millwood) 2004;(Suppl Web Exclusives):VAR33–44. [PubMed]
16. Barnato AE, Berhane Z, Weissfeld LA, et al. Racial variation in end-of-life intensive care use: a race or hospital effect? Health Serv Res. 2006;41:2219–2237. [PMC free article] [PubMed]
17. Barnato AE, Chang CC, Saynina O, et al. Influence of race on inpatient treatment intensity at the end of life. J Gen Intern Med. 2007;22:338–345. [PMC free article] [PubMed]
18. Mitchell SL, Teno JM, Roy J, et al. Clinical and organizational factors associated with feeding tube use among nursing home residents with advanced cognitive impairment. JAMA. 2003;290:73–80. [PubMed]
19. Mebane EW, Oman RF, Kroonen LT, et al. The influence of physician race, age, and gender on physician attitudes toward advance care directives and preferences for end-of-life decision-making. J Am Geriatr Soc. 1999;47:579–591. [PubMed]
20. Diringer MN, Edwards DF, Aiyagari V, et al. Factors associated with withdrawal of mechanical ventilation in a neurology/neurosurgery intensive care unit. Crit Care Med. 2001;29:1792–1797. [PubMed]
21. Garrett JM, Harris RP, Norburn JK, et al. Life-sustaining treatments during terminal illness: who wants what? J Gen Intern Med. 1993;8:361–368. [PubMed]
22. McKinley ED, Garrett JM, Evans AT, et al. Differences in end-of-life decision making among black and white ambulatory cancer patients. J Gen Intern Med. 1996;11:651–656. [PubMed]
23. O'Brien LA, Grisso JA, Maislin G, et al. Nursing home residents' preferences for life-sustaining treatments. JAMA. 1995;274:1775–1779. [PubMed]
24. Kane RL, Friedman B. State variations in Medicare expenditures. Am J Public Health. 1997;87:1611–1619. [PubMed]
25. Fisher ES, Wennberg JE, Stukel TA, et al. Associations among hospital capacity, utilization, and mortality of US Medicare beneficiaries, controlling for sociodemographic factors. Health Serv Res. 2000;34:1351–1362. [PMC free article] [PubMed]
26. Halpern NA, Pastores SM, Thaler HT, et al. Changes in critical care beds and occupancy in the United States 1985 to 2000: differences attributable to hospital size. Crit Care Med. 2006;34:2105–2112. [PubMed]
27. Decker FH. Dying in a nursing home: the role of local bed supply in nursing home discharges. J Aging Health. 2008;20:66–88. [PubMed]
28. Wennberg JE, Cooper MM. The Dartmouth Atlas of Health Care. Chicago, IL: American Hospital Publishing; 1998.
29. Iwashyna TJ, Chang VW, Zhang JX, et al. The lack of effect of market structure on hospice use. Health Serv Res. 2002;37:1531–1551. [PMC free article] [PubMed]