Journal Home
Search for

Volume 89, Issue 8, Pages 1482-1491 (August 2008)


View previous. 13 of 40 View next.

A Prospective Study of Health and Risk of Mortality After Spinal Cord Injury

James S. Krause, PhDaCorresponding Author Informationemail address, Rickey E. Carter, PhDb, E. Elisabeth Pickelsimer, DAb, Dulaney Wilson, MSPHb

Abstract 

Krause JS, Carter RE, Pickelsimer EE, Wilson D. A prospective study of health and risk of mortality after spinal cord injury.

Objective

To test hypothesized relationships between multiple health parameters and mortality among persons with spinal cord injury (SCI) while controlling for variations in biographical and injury characteristics.

Design

Prospective cohort study with health data collected in late 1997 and early 1998 and mortality status ascertained in December 2005.

Setting

A large rehabilitation hospital in the Southeastern United States.

Participants

A total of 1389 adults with traumatic SCI, at least 1 year postinjury.

Interventions

Not applicable.

Main Outcome Measures

The primary outcome was time from survey to mortality (or time of censoring). Mortality status was determined using the National Death Index and the Social Security Death Index. There were 225 deaths (16.2%) by December 31, 2005.

Results

Cox proportional hazards modeling identified several significant health predictors of mortality status, while controlling for biographic and injury factors. Two sets of analyses were conducted—the first identifying the significance of a single variable of interest and the second analysis building a comprehensive model based on an optimal group of variables. Multiple types of health conditions were associated with mortality. The best set of health predictors included probable major depression, surgeries to repair pressure ulcers, fractures and/or amputations, symptoms of infections, and days hospitalized. Inclusion of these variables, along with a general health rating, improved prediction of survival compared with biographic and injury variables alone, because the pseudo R2 increased from .12 to .18 and the concordance from .730 to .776.

Conclusions

In addition to secondary conditions that have been the traditional focus of prevention efforts (eg, pressure ulcers, urinary tract infections), amputations, fractures, and depressive symptoms were associated with higher risk for mortality; however, further research is needed to identify the association of specific conditions with causes of death and to determine whether interventions can modify these conditions and ultimately improve survival.

Article Outline

Abstract

Causes of Death

Health Factors and Mortality

Purpose

Methods

Participants

Prospective Data Collection Procedures

Determination of Mortality Status

Measures

Data Analyses

Results

Discussion

Study Implications

Study Limitations

Future Research

Conclusions

Acknowledgment

References

Copyright

ALTHOUGH THERE HAS BEEN a trend toward substantial improvement in year 1 survival rates among persons with SCI, there has been no improvement beyond that in the past 25 years, and persons with SCI remain at high risk for premature death.1, 2, 3 Enhancing longevity is not solely dependent on new medical breakthroughs. It may also result from new knowledge from epidemiologic research that helps professionals better apply existing techniques by identifying and targeting people at highest risk and by a better understanding of the risks and protective factors leading to early mortality.

Our understanding of the epidemiology of mortality after SCI remains rather limited. We know that the chance of dying is greatest in the first year after injury.4, 5, 6 The risk is greatest for those with a more severe injury and multiple trauma at time of injury4 and for elderly persons who are likely to have compromised health and a greater number of comorbidities.6 Evidence from several studies, including a study spanning 50 years of 3179 participants, clearly indicates that a higher neurologic level and completeness of SCI and older age at injury are related to an elevated risk of mortality.6, 7, 8, 9, 10, 11, 12 Additional strong predictors of premature mortality include ventilator dependency,12, 13 preinjury cardiovascular disease, alcohol or substance abuse, psychiatric diagnosis, and being female.9

Contrary to findings of Lidal et al,9 when compared by age and injury severity, the MSCIS found that women have slightly higher survival rates, yet the difference was not as great as between men and women in the general population.5 The MSCIS further reported that being white5 and being married14 were determinants of survival.

Causes of Death 

return to Article Outline

Investigators attempting to identify health factors that may precipitate mortality have generally relied on retrospective studies or studies of causes of death, rather than prospective cohort studies collecting uniform data on health outcomes and assessing their later relationship with mortality. Whiteneck et al15 used life-table techniques to determine survival rate, causes of death, and psychosocial outcomes of persons with SCI in 2 British hospitals. The most prevalent causes of death were genitourinary problems (primarily renal failure), followed by cardiovascular and respiratory diseases. In a study of 195 persons experiencing new SCI between 1955 and 1994, Soden et al16 reported incidences of death from septicemia, pneumonia and influenza, diseases of the urinary system, and suicide were significantly greater than in the general population.

In an analysis of underlying and contributing causes of death, DeVivo and Stover17 classified a total of 1403 deaths recorded on death certificates into categories suggested by the National Center for Health Statistics. They noted that pneumonia and influenza were the primary cause of death (17.7%), followed by nonischemic heart disease (16.5%) and septicemia (12%). They calculated SMRs to compare causes of deaths of persons with SCI to the ratio of deaths expected in the general population for all causes. The highest SMRs were for septicemia (64.2 times the general population), disease of the pulmonary circulation (47.1), pneumonia and influenza (35.6), symptoms and ill-defined conditions (13.8), and diseases of the urinary system (10.9). It is noteworthy that almost all causes of death were higher among the SCI sample, with the exception of homicide, legal intervention, and cancer.

A study of 361 men with SCI age 23 to 87 years who had been injured 1 to 56.5 years assessed the relationship between comorbid medical conditions and other health-related factors to mortality.18 The NDI provided cause-specific mortality over a median of 55.6 months (range, 0.33−74.4) through December 31, 2000. Mortality was elevated (SMR=1.47) compared with US rates. Diabetes (RR=2.62), heart disease (RR=3.66), reduced pulmonary function, and smoking were major risk factors for death. The most common underlying and contributing causes of death were diseases of the circulatory system (40%) and of the respiratory system (24%).18 In a less detailed follow-up, DeVivo et al2 noted that respiratory-related cause of death occurring after the first anniversary of injury was the only cause for which the relative odds increased meaningfully during the latest period (76% increase over 1988–1992 compared with all other causes).

CVD is the most common underlying cause of death noted on death certificates of persons with or without SCI. Long-term SCI or being nonambulatory is a risk factor for a greater prevalence of CVD, obesity, lipid disorders, metabolic syndrome, and diabetes—all factors related to increased mortality.19

Health Factors and Mortality 

In the aforementioned studies, the health conditions leading to mortality were directly identified from causes of death mostly listed on death certificates. In contrast, there has been limited research that has independently and prospectively measured secondary conditions and other indicators of morbidity and linked them with eventual mortality. We therefore know very little about how the development of a specific health condition may increase the likelihood of mortality or how subtle indicators of change in health may also signal risk of early death. Nevertheless, the success of interventions to enhance longevity is dependent on an understanding of these factors.

Early studies of the association of multiple variables with mortality suggested the psychosocial and vocational outcomes were more highly correlated with mortality than were general health indicators.20, 21 Because the time between collection of the prospective data and determination of mortality status was 11 years and 15 years, respectively, health status could have changed substantially in the years after data collection, as was suggested by a study that used a shorter follow-up interval (4y) but had identified only 22 deceased participants.22 A relatively recent report from the MSCIS in the United States14 found general health ratings were associated with an elevated risk of mortality, but the data on more specific conditions were limited.

Purpose 

The purpose of this study was to identify which health outcomes are associated with elevated risk of mortality after SCI. By focusing on health factors and mortality, this study evaluated the primary component of a hierarchical empirical risk model for mortality23 while controlling for the independent effects of biographical and injury factors. There are 2 related objectives: (1) to evaluate the extent to which health outcomes, as measured by general indices, indicators of health impact, and secondary conditions, are associated with hazard of mortality after SCI; and (2) to evaluate the extent to which accounting for different sets of predictors enhances our prediction of mortality, moving progressively from injury severity to a comprehensive model that includes all the health factors available in the study.

Our hypotheses were (1) when statistically controlling for biographical and injury characteristics, several types of health factors will be associated with hazard of mortality, including indices of general health, health impact, and secondary conditions; and (2) when building an optimal risk model for mortality, inclusion of health factors will enhance our prediction of hazard for mortality above and beyond that of biographical and injury factors alone.

Methods 

return to Article Outline

Participants 

After obtaining approval from the institutional review board, participants were identified from 3 types of records of a specialty hospital in the Southeastern United States: (1) MSCIS patient database, (2) MSCIS registry, and (3) outpatient directory. To be included in the study, participants were adults with traumatic SCI occurring at least 1 year prior to the study that resulted in some residual impairment (those with no neurologic deficits were excluded). Of an original cohort of 1929 patients, 1389 participated (response rate, 72%).

Seventy-four percent were men, and 75% were white. The overall breakdown by sex and race and ethnicity was 54.7% white men, 19.3% white women, 20.2% minority men, and 5.8% minority women. Of the minority participants, 87.7% were African American. Average age at time of injury was 31.8 years (IQR, 20.6–39.8y). At the time of prospective data collection, participants' mean age was 40.7 years (IQR, 30.3–48.9y). They had been injured a mean of 9.8 years (IQR, 3.5–12.2y), and the average education was 13.1 years. The primary etiology was vehicular crashes (50.9%), followed by falls or flying objects (17.3%), acts of violence (12.7%), sports (11.9%), and other (7.2%). Fifty-five percent reported cervical injuries. Just over a fifth (21.4%) reported functional movement below level of injury sufficient to assist ambulation.

Prospective Data Collection Procedures 

Letters announcing the study were sent to all prospective participants, followed by a mailed survey 4 to 5 weeks later. All introductory letters followed U.S. Postal Service standard guidelines to identify changes of address whenever possible. Two additional mail-out attempts and a follow-up phone call were made to nonresponders. Prospective data collection was initiated in July 1997 and ended in April 1998, although nearly all data had been collected by the end of December 1997. To elicit the highest possible response rates, participants were offered a $20 stipend and were made eligible for drawings totaling $1500.

Determination of Mortality Status 

Mortality status was determined an average of 8 years after the primary data collection (December 2005). Persons not found to be deceased by December 31, 2005, were assumed to be alive. To ensure a comprehensive determination of mortality, we elected to use both the NDI of the National Center for Health Statistics24 and the SSDI of the Social Security Administration.25 Although neither data source is reported to be 100% inclusive,26, 27, 28 researchers often use the NDI and the SSDI to determine mortality. NDI provides access to information on decedents through a centralized computer index by searching death records provided by state offices. NDI provides information on the name of the state where death occurred, the death certificate number, and the date of death for each subject searched and identified. NDI death records are available approximately 16 months after the conclusion of a given year (ie, NDI records were not available for 2004–2005). Among persons in the MSCIS whose full name and Social Security number were known, sensitivity and specificity of using the SSDI were 92.4% and 99.5%, respectively.5 In a later study of 371 known decedents, the SSDI correctly identified almost 87% of decedents age 50 years and older, 69.8% of decedents age 30 to 49 years, and only 34.6% of decedents younger than age 30 years. Factors that influenced sensitivity were older age, male sex, white race, married status, and survival greater than 1 year after injury.29

Measures 

Several instruments were used to measure the health factors and secondary conditions, including subsets of items from the BRFSS30 and LSQ-R.31 The OAHMQ32 was used to measure depressive symptoms. Selected items were also developed for the study.

Several types of health outcomes were measured that fall into the following broad categories: (1) general health indicators, (2) treatment variables, and (3) secondary conditions. The secondary condition variables fell into 3 general classes: symptoms of SCI complications (particularly symptoms of infections, eg, urinary tract complications, sweats, chills), irreversible conditions (eg, amputations), and reoccurring conditions (eg, pressure ulcers, injuries, depression). To minimize the potential for an inflated type I error rate, the biologic plausibility of each health variable was considered prior to analysis. Whereas the first 2 classes of variables, general health indicators and treatment variables, are indicators of overall health, secondary condition variables are more viable targets for intervention.

Table 1 summarizes the variables to be used in the study. The first variable (self-health rating) may be best conceptualized as the summation of health, not linked to any specific secondary condition, event, or health variable. It was measured by a single BRFSS item: self-rating current health. Two other health impact variables were measured by 2 additional BRFSS items: days in poor physical health and days in poor mental health, each of which required participants to indicate the number of days out of the past 30 for which their physical or mental health was not good. The CDC has used these items extensively and has existing data available for normative purposes. Three treatment variables were used from the LSQ-R. Each of the 3 LSQ-R items requires participants to indicate over the previous 12 months the number of nonroutine physician office visits, (2) the number of hospitalizations, and (3) the number of days hospitalized.31 Several types of secondary conditions were assessed. Symptoms of secondary conditions were measured by items reflecting complications of SCI. These included fevers, sweats and chills, and UTIs (self-report based on diagnoses by physicians). Participants were asked to indicate the number of occurrences of each of these in the past year using ordinal rankings (0, 1–2, 3–6, 7–12, ≥13). A summary measure (summated score) was created. Irreversible secondary conditions reflected historical events that could not be reversed (eg, upper-limb or lower-limb amputation), events that led to a major nonrecurrent complication (upper-extremity or lower-extremity fracture, DVT), or the accumulation of conditions that, although in theory reversible, in practice rarely occur (eg, upper-limb or lower-limb contractures, curvature of the spine). There were 8 items, each of which asked participants whether the condition had ever occurred. The items included (1) blood clots (DVT), (2) amputation of a lower extremity (leg, foot), (3) amputation of an upper extremity (hand, arm), (4) curvature of the spine (more difficulty sitting up straight), (5) broken bone in a lower extremity (leg, foot), (6) broken bone in an upper extremity (hand, arm), (7) contracture in hip, knee, or ankle, and (8) contracture in elbow or shoulder.

Table 1.

Variables Used in the Data Analysis, Type of Variable, and the Instrument From Which They Were Selected

Biologic DomainType of VariableInstrument
Variable
A. General health indicators
General health ratingOrdinal: scale ABRFSS
No. days poor health over past 30 daysContinuousBRFSS
No. days poor mental health over past 30 daysContinuousBRFSS
B. Treatments
Nonroutine physician visits over past yearOrdinal: scale BLSQ-R
Hospitalizations (any reason) over past yearOrdinal: scale BLSQ-R
Days spent hospitalized over past yearOrdinal: scale CLSQ-R
C. Secondary conditions
C1. Symptoms of SCI complications (in past year)
Stomach pain and distention (bloating)Ordinal: scale DSCI health survey
Bowel accidentsOrdinal: scale DSCI health survey
Rectal bleedingOrdinal: scale DSCI health survey
Urine leaking or accidentsOrdinal: scale DSCI health survey
FeversOrdinal: scale DSCI health survey
Sweats and chillsOrdinal: scale DSCI health survey
UTIs (must be diagnosed by a doctor)Ordinal: scale DSCI health survey
No. of UTI-like symptomsContinuousDerived as sum of fevers, sweats/chills, and diagnosed UTIs
C2. Irreversible conditions
Report fractures to lower or upper extremities since injuryBinarySCI health survey
Report amputations to lower or upper extremities since injuryBinarySCI health survey
Fracture or amputation since injuryBinaryDerived indicator of any fracture or amputation since injury
Reported blood clots, curvature of the spine, or contracture since injuryBinaryDerived from items on SCI health survey
C3. Recurrent conditions
Pressure ulcers
Since injury, how many pressure ulcers were severe enough to require medical attentionOrdinal: scale BSCI health survey
Pressure ulcer that required attentionBinaryDerived as having at least 1 pressure ulcer
Since injury, how many surgeries to repair pressure ulcersOrdinal: scale BSCI health survey
Since injury, what best describes pressure ulcer historyCategorical: scale ESCI health survey
History of frequent pressure ulcersBinaryDerived indicator representing people who always seem to have ulcers, often requiring surgery or hospitalization
Subsequent injuries
In past year, no. of injuries that required medical careOrdinal: scale BSCI health survey
Since injury, no. of injuries that required medical careOrdinal: scale BSCI health survey
Since injury, no. of times hospitalized because of injuriesOrdinal: scale BSCI health survey
Depressive symptoms
Probable major depressionBinary: scale FOAHMQ

Scale A: poor (1), fair (2), good (3), very good (4), excellent (5).

Scale B: 0 (0), 1 (1), … , 9 (9), 10+ (10).

Scale C: 0 (0), 1 (1), … , 29 (29), 30+ (30).

Scale D: How many times did you have the following problems? 0, 1−2 (1.5), 3−6 (4.5), 7−12 (9.5), 13+ (14).

Scale E: A, Have never had any pressure ulcers; B, had pressure ulcers immediately after injury, but rarely, if at all, since initial rehabilitation hospitalization; C, get pressure ulcers about every couple years; D, get at least 1 pressure ulcer a year; and E, always seem to have ulcers, often requiring surgery or hospitalization.

Scale F: OAHMI scores: normal depressive symptomatology (OAHMI scores <6); clinically significant symptomatology (OAHMI score range, 6 to <11); and probable major depression (OAHMI score range, 11−22). Binary depression variable of OAHMI scores ≥11.

For each of the following scales, the value provided in parentheses is the ordinal score used for the analysis. If no score is given (and the type of variable is categorical), the corresponding indicator (ie, dummy variable) was constructed.

In contrast, recurrent secondary conditions may have repeated cycles. We assessed 3 types of conditions: pressure ulcers, subsequent injuries, and depressive symptoms. Pressure ulcers were defined as “open sores in pressure areas, such as your tailbone, ischium, heel, elbows. They are usually caused by pressure but may also be caused by friction or shearing (rubbing), moisture, burns, or falls.”33(p1259) Participants were to indicate (1) the number of pressure ulcers they had in the past year, (2) the number of days they were forced to reduce their sitting time as the result of a pressure ulcer over the past year (they were presented 8 multiple-choice categories), (3) whether they had a pressure ulcer at time of the study, (4) the number of surgeries since the onset of their SCI to heal pressure ulcers, and (5) which of 5 categories best described their pressure ulcer history (categories were described along a continuum from nonrecurrent to recurrent).

Developed in collaboration with the National Center for Injury Prevention of the CDC, subsequent injuries were defined in the following manner: “The following questions relate to INJURIES, including broken bones, burns, or lacerations (cuts). Injuries happen as the result of some type of mishap or event, such as a fall, collision, motor vehicle wreck, or act of violence.”34(p1504) Participants were asked the number of times in the past year they had been injured seriously enough to receive medical care in a clinic, emergency department, or hospital and whether any of these injuries had been serious enough to require treatment in a hospital for at least 1 night. Subsequent injuries were reported in 19% of participants over a 12-month period, with 27% of participants who reported at least 1 injury also reporting 1 or more injury-related hospitalizations.34

The 22-item OAHMQ32 was used to measure depressive symptoms. The OAHMQ was developed to measure depression in older adults and among people with physical disabilities by including few items that reflect physical or vegetative symptomatology, because these types of items often invalidate commonly used measures of depression and other clinical syndromes. Scores 11 or higher were considered to indicate probable major depression. Kemp and Adams32 reported test-retest correlations were found to be .87 (P<.001) and α values were .93 (P<.001). In terms of validity, sensitivity and specificity were reported to be .93 and .87, respectively. The OAHMQ was also compared to 2 other depression scales with established validity and reliability, the Geriatric Depression Scale and the Symptom Checklist−90−Revised. Correlation for the OAHMQ and both measures was .70 (P<.001). The OAHMQ has been used in previous research with SCI.35

Data Analyses 

A 3-stage hierarchical strategy to model building was employed to identify the association of each health variable with mortality and to define an optimal set of health predictors of mortality. Cox proportional hazards modeling was used with the number of days between the survey and event (ie, mortality) as the dependent variable. The censoring date was December 31, 2005, the date up to which mortality status could be verified.

During the first stage of analysis, a base model consisting of biographic and injury characteristics, including sex, race (white-minority), age, years lived since injury, and injury severity, was specified (table 1 summarizes the variables and their coding). Injury severity was defined by a combination of 2 items that were used to generate 5 categories that were similar, but not equivalent, to those frequently reported in the SCI mortality literature. The groups were defined as a combination of injury level and neurologic completeness of injury, with injury level broken down into 3 categories (C1-4, C5-8, T1-S5) among those who did not report motor functional recovery sufficient for ambulation. Among those who had recovery sufficient for ambulation, we further broke this group into those with cervical injuries and those with noncervical injuries. Table 2 summarizes the breakdown of participants according to this categorization.

Table 2.

Breakdown as a Function of Injury Level and Neurologic Completeness and Ability to Ambulate

CategoryNo.Percent
C1-4, nonfunctional17613.2
C5-8, nonfunctional40730.5
Noncervical, nonfunctional46935.2
Cervical, ambulatory14210.7
Noncervical, ambulatory13910.4
Unclassified (not used in analysis)56

Percentages exclude the 56 observations not classified.

The second stage of the analysis focused on adding single, biologically plausible health outcome variables to the base model, thereby screening each of these potential predictors for inclusion in the final stage model. All variables significant at the α equal to .10 level of significance were considered for subsequent modeling.36 Once variables had been screened within a biologic domain (eg, recurrent secondary conditions), those identified were then assessed for multicolinearity within a biologic domain, and backward selection was used to identify the key predictor variables for each biologic domain. The variables identified by using this process are summarized in table 3.

Table 3.

Cox Proportional Hazards Analysis for the Preliminary and Full Models

Biologic DomainStage 2 ResultsSelected for Stage 3Stage 3 Results
VariableHR95% CIPHR95% CIP
Core adjustment variables
Functional status
C1-4, nonfunctionalNANANAYes4.251.88−9.64<.01
C5-8, nonfunctionalNANANAYes2.871.30−6.34<.01
Noncervical, nonfunctionalNANANAYes2.581.17−5.68.02
Cervical, ambulatoryNANANAYes1.290.50−3.33.60
Noncervical, ambulatory (referent)NANANAYes1.00NDND
Caucasian vs non-CaucasianNANANAYes1.140.82−1.60.43
SexNANANAYes1.220.87−1.73.25
Age at injuryNANANAYes1.061.05−1.07<.01
Years since injuryNANANAYes1.041.02−1.06<.01
A. General health indicators
General health rating0.620.53−0.72<.01No§NDNDND
Poor health days1.021.01−1.04<.01No§NDNDND
Poor mental health days1.031.01−1.04<.01No§NDNDND
B. Treatments
Nonroutine doctor visits1.081.04−1.12<.01YesNSNDND
Hospitalizations1.161.09−1.23<.01NoNDNDND
Days spent in hospital1.031.02−1.05<.01Yes1.021.00−1.03.03
C1. Symptoms for SCI complications
Stomach pain/distention1.021.00−1.05.05NoNDNDND
Bowel accidents0.990.96−1.03.71NoNDNDND
Rectal bleeding0.990.96−1.03.59NoNDNDND
Urine leaking or accidents0.990.97−1.02.61NoNDNDND
Fevers1.051.02−1.08<.01NoNDNDND
Sweats and chills1.051.02−1.08<.01NoNDNDND
UTIs1.071.03−1.10<.01NoNDNDND
No. of infection symptoms1.031.02−1.04<.01Yes1.021.00−1.03.03
C2. Irreversible conditions
Fractures to extremities1.501.05−2.15.03NoNDNDND
Amputations to extremities3.692.14−6.37<.01NoNDNDND
Fracture or amputation since injury1.811.31−2.52<.01Yes1.781.27−2.51<.01
Other irreversible secondary conditions1.240.94−1.64.13NoNDNDND
Pressure ulcers
No. of ulcers needing medical care since injury1.151.09−1.21<.01NoNDNDND
At least 1 ulcer needing care1.921.43−2.62<.01YesNSNDND
Surgeries to repair ulcers1.201.13−1.27<.01Yes1.161.08−1.23<.01
History of frequent ulcers3.682.18−6.19<.01NoNDNDND
Subsequent injuries
Injuries in the past year needing medical care1.131.00−1.26.04NoNDNDND
Injuries since injury needing medical care1.181.10−1.26<.01YesNSNDND
Hospitalized for injury since injury1.161.09−1.25<.01NoNDNDND
Depression
Probable major depression2.211.66−2.94<.01Yes1.861.36−2.53<.01

Abbreviations: CI, confidence interval; NA, not applicable; ND, no data to enter into cell; NS, variable was eliminated from the stage 3 model using backward selection.

HRs from a Cox model predicting the time, measured in days, from survey to either mortality or censoring (12/31/05) and adjusted for core (biographic and injury) variables listed in the first biologic domain.

HRs for the stage 3 final model are adjusted for all other variables in which an estimate is provided. Dashes indicate that the variable was not considered for stage 3 analysis.

HRs for core adjustment variables not shown for stage 2 results because there is variability of the estimates depending on which of the candidate variables is used in the analysis.

§

General health indicators, while significant predictors in stage 2, were not considered for stage 3 of the analysis. They were used, however, for sensitivity analyses and alternative model specifications.

The final stage of the analysis formulated a Cox proportional hazards model that consisted of the base model in addition to the variables identified in stage 2 of the analysis. After assessment for multicolinearity, backward elimination was used to identify the final fitted model. The general health indicators domain was not considered for the backward elimination; however, the general health rating and physical health rating from the BRFSS were included in alternative (exploratory) models. We chose to focus on secondary conditions because they have more direct implications for prevention.

The proportional hazards assumption was assessed using the Schoenfeld residuals37 and found to be tenable. The fit of the model was assessed using the likelihood ratio test and the c statistic.38 The likelihood ratio test was used to calculate the Nagelkerke pseudo R2.39 The value of the c statistic is closely related to the area under a receiver operating characteristic curve and is interpretable as the probability that the cases (ie, deaths) have higher risks as measured by the linear component of the regression model. Accordingly, a value of 0.5 is for chance prediction, and the discrimination of the model is improved as the c value approaches 1.0.40, 41 Once the final model was determined, all pairwise interaction terms of the secondary conditions were included in a new model to assess goodness of fit further. A Wald linear contrast indicated these interaction terms were not needed in the model (P>.50), and accordingly, these interaction terms were removed. All model building was conducted using the SAS.a The validation of the proportional hazards assumption and the estimation of the c statistic were performed using Stata.b

Results 

return to Article Outline

Of the 1389 participants originally studied, missing data related to the health indicators, to the uncertainty of the injury severity, and to an unknown time of SCI onset reduced the sample size for statistical modeling. The final statistical model consisted of 1265 observations and 188 observed deaths. Of the variables listed in table 1, many were found to be independent predictors of mortality after adjustment for general biographic and injury variables (see table 3). Complications related to SCI were among the predictors identified for further modeling. In particular, an increased number of UTIs, sweats and chills, and fevers were all associated with an increased risk of mortality. The summated rating of these symptoms (mean ± SD, 10±10) remained highly statistically significant, and for a 1 SD increase, the adjusted hazard of mortality increased by 34% (HR=1.34, P<.01). This finding persisted in the final model with negligible attenuation of the magnitude of the HR for the 1 SD change (HR=1.33, P=.03) when all other variables were held constant in the final model.

Two classifications (or biologic domains) consisting of treatment-related variables were studied. The first classification consisted of nonspecific measures of physician and hospitalization utilization. The second classification focused on subsequent injury-related physician and hospitalization utilization. While all 3 variables contained in the subsequent injury classification were significant predictors at stage 2 of the analysis, none persisted to the final model. The more general classification did yield a predictor that remained in the final model. In particular, the number of hospital bed days a participant had over the past year (mean, 4.4±8.7) was a significant predictor in both stages 2 and 3 of the analysis. For each SD increase in the number of hospital bed days (and when all other variables were held constant), the adjusted hazard for mortality increased by 15% (HR=1.13, P=.03).

A participant's pressure ulcer history, recent and since time of injury, was an important prognostic class of variables. Variables representing the frequency of pressure ulcers as well as the severity of the pressure ulcers were identified for further modeling. In addition, people who classified themselves as “Always seem to have sores, often requiring surgery or hospitalization” were at a much greater risk (HR=2.5, P<.01) of mortality than the other 4 classifications of pressure ulcer history; however, this variable was eliminated from consideration by the backward selection within the biological domain that was included during stage 2 of the analysis. The number of surgeries to repair pressure ulcers was a statistically significant predictor during both stages 2 and 3 of the analysis. The final model estimated that for every additional ulcer that required surgery, the risk for mortality increased by 10% (HR=1.16, P<.01) when adjusted for all other variables in the model.

There was mixed evidence for the importance of irreversible symptoms. A history of an amputation or fracture since injury predicted mortality, yet a history of other irreversible symptoms (eg, DVT, curvatures of the spine, contractures) was not predictive of mortality. While a history of an amputation alone had the largest HR in stage 2, only 2% of the sample had an amputation at the time of the study. It was for this reason that a derived indicator representing injury to the extremities (either amputation or fracture to an extremity) was considered for the final stage of modeling. This approach attenuated the magnitude of the association of amputations alone, but because 18% of the sample had either an amputation or a fracture, the variable was conducive to regression modeling strategies. In particular, the hazard for mortality was 78% higher (HR=1.78, P<.01) in those participants who reported having either a fracture or amputation in the final model.

The subsequent injuries domain did not yield a predictor for the final model; however, all 3 variables were associated with mortality in stage 2 of the modeling. After backward elimination with domain, only the number of injuries since the time of SCI was considered for stage 3. This predictor was the last variable eliminated during stage 3.

The importance of probable major depression persisted in the final model. The hazard for mortality was 86% higher (HR=1.86, P<.01) for those meeting the definition of probable major depression (scores ≥11) than for those not meeting the definition. While not considered for the final model, the stage 2 results for the number of poor mental health days over the past 30 days (mean, 6.2±8.6d) supported this finding. In particular, for each SD increase in the number of poor mental health days, the hazard increased by a factor of 23% when controlled only for the core adjustment variables. In addition, as a participant's perception of his/her general health decreased (lower score) or as the number of poor health days increased, the risk of mortality increased (see table 3, stage 2 results). While only the diagnosis of probable major depression was considered for stage 3 of the analysis, the other general physical health variables were evaluated further in sensitivity analyses.

Table 4 presents 5 models developed through model building exercises. The first model includes only the injury severity variable broken down into 5 categories related to neurologic level of injury and neurologic completeness of injury. This is the best single indicator of the SCI per se with mortality. The pseudo R2 was only .02 and the c statistic was only .578. The addition of biographic variables and age at injury onset substantially increased both indicators (R2=.121, c=.730). Model 3 is the final model determined during stage 3 of the analysis. The model provided a high degree of discrimination (c=.771), which was substantially higher than the first model with injury severity alone (increase of .193 in the c statistic), but only modestly higher than the full biographic injury model alone (increase of .041 in the c statistic). The pseudo R2 value (R2=.167) represented a more substantial increase than the biographic injury model (increase of .046). The fourth model, adding the 2 general health variables (overall rating of health and number of poor health days), further improved the pseudo R2 value slightly (from .167 to .178), but with a negligible increase in discrimination (c increased from .771 to .776).

Table 4.

Summary of the Pseudo R2, R Change, c Statistic, and Change in the c Statistic for Each of 4 Stages of Model Building and an Alternative Model

ModelDescriptionR2Changec StatisticChange
1Injury severity only.016ND.578ND
2Core biographic and injury variables.121.105.730.152
3Full final stage 3 model.167.046.771.041
4Alternative model that includes the general health ratings.178.011.776.005
ATreatment and secondary conditions.075.059.676.098

Abbreviation: ND, no data to enter into cell.

1. Five category breakdown (C1-4, C5-8, noncervical, and so forth).

2. Add core biographic variables includes Caucasian race, male sex, and age at injury.

3. Final model. Exclusive of general health indicators.

4. Alternative model that includes general health ratings.

A. Alternative model that included treatments and secondary conditions, but not biographic or injured variables. The change figures for this model are compared with the first model.

Fig 1, Fig 2 graphically illustrate the estimated survival curves for the final model. As illustrated in figure 1, the rate of mortality in the C1-4, nonfunctional injury classification was substantially higher than in the other 4 injury classifications.


View full-size image.

Fig 1. Estimated survival curves for the 5 functional injury classifications at the mean level of all other variables included in the final model.



View full-size image.

Fig 2. Estimated survival curves for persons free of health complications (Good Health) relative to persons with numerous complications (Poor Health). Good health was defined as nondepressed, 0 days hospitalized over the past year, free of symptoms of infections, no broken bones or amputations, and no surgical repairs of ulcers. Poor health was defined as probable major depression, with 13.1+1 SD days hospitalized, 20.0+1 SD occurrences of infectionlike symptoms, a broken bone or amputation, and 2.1+1 SD surgeries to repair pressure ulcers. All other variables included in the model were assumed to be at the mean level.


Figure 2 illustrates the impact of the health indicators on survival for 2 general classes of people. The first class presented is a good health class in which people were free of health complications (ie, not depressed, no ulcers, and so forth). The second class was defined to represent people who had a complication value of 1 SD above the mean for each respective health indicator.

Discussion 

return to Article Outline

The unique contribution of this study is the identification of the association of multiple health risk factors for mortality after SCI using prospective data rather than causes of death. This allowed us to test 2 study hypotheses relating health factors for mortality and the additional contributions of health factors to the overall efficacy of prediction of mortality.

The first hypothesis was confirmed, because multiple health factors were found to be associated with mortality. The final model included 4 primary types of secondary conditions including symptoms of secondary conditions (particularly infections), irreversible conditions (amputations and fractures), surgeries for pressure ulcers, and the presence of probable major depression. A treatment variable, number of days hospitalized, was also included in the final model. The diversity of these factors speaks to the importance of overall health promotion in multiple areas of life. It must be noted that, although these were the factors that were included in the final model, many other conditions were significantly related to mortality. The final model simply incorporates those that predict the greatest amount of unique variance. Other factors, including subsequent injuries, were also related to increased risk for mortality but shared variance with other predictors of the equation. Similarly, although we selected 1 predictor from each secondary condition, other parameters were also important. For instance, although we ultimately retained the number of surgeries for pressure ulcers in the final equation, recurrent pressure ulcers and pressure ulcers within the past year were also related to early mortality during preliminary analyses. In sum, diverse secondary conditions should serve as markers for increased risk of early mortality.

The second hypothesis, that inclusion of health factors would enhance prediction of hazard for mortality above and beyond that of biographic and injury factors alone, was also supported, although the increase in prediction varied between indicators. The pseudo R2 estimate was enhanced by the addition of the health predictors, increasing from .121 to .167, a 38% increase. However, the discrimination only increased from .730 to .771. Inclusion of general health ratings further increased the pseudo R2 estimate to .178 and the discrimination to .776.

The overall magnitude of the pseudo R2 suggests there is a great deal that one needs to learn about prediction of mortality; however, until more mortality events are observed, this prediction is not likely to improve dramatically. The modest improvement in discrimination also suggests substantial limitations of our ability to enhance classification or prediction above and beyond that of the biographic and injury predictors with current methodologies that use snapshot measures of predictors (ie, a single point in time) when those predictors change over time. This is consistent with other research from the MSCIS in which 4 different models were tested and the highest generalized R2 was only .131.22

One particularly interesting finding is that the simplest model based in functional status alone accounted for minimal variation in mortality (pseudo R2=.016, c=.578), and the true enhancement and prediction came from the biographic variables, which included age at injury onset. When we used an alternative model that was restricted to significant health variables (hospitalizations, secondary conditions) and compared this with the base model from injury severity alone, the alternative model accounted for a greater amount of variance (the pseudo R2 for functional status alone was .016, compared with a pseudo R2 of .075 for the alternative model, which included hospitalizations and secondary conditions). Similarly, the alternative model had a higher discrimination (c=.676) than the model with functional status alone (c=.578). Therefore, even though functional status (ie, injury severity) is the most commonly used parameter in tables related to life expectancy after SCI other than chronologic age, by itself it is not as strong a predictor of mortality as overall health.

Study Implications 

Diminished life expectancy after SCI is one of the most well established findings in SCI research. Yet how often does increased risk of mortality enter into routine outpatient assessments, patient conferences, or physician or health provider discussions with people with SCI? The current study does not address this question, yet it certainly indicates that this is a question that needs to be addressed. How can we intervene to reduce the risk of mortality if it is not even a topic discussed with our patients with SCI?

We have sufficient research to direct interventions toward early identification of risk on the basis of the presence of multiple health indicators. Pressure ulcers, UTIs, amputations or fractures, and the development of depressive disorders are all predictive of early mortality, and the presence of these conditions should serve as a red flag for intervention. Some interventions can and should be incorporated into routine assessments, the minimum of which is to inform patients with SCI regularly of the association between their health status and early mortality so they have heightened awareness of risk.

Although we have not demonstrated a causative relationship (only an association) between these factors and mortality, secondary conditions with the strongest associations with early mortality are logical targets for intervention. One type of screening that appears to be particularly important is that of depressive symptoms and the possibility of a major depressive disorder. Because most people who find their way to outpatient clinics present with problems, health care providers and outpatient clinics may assume that depressive symptoms are a natural comorbidity of health complications. A minimum of brief screening is required, sometimes followed up by a full assessment by a licensed psychologist, to rule out more dangerous depressive diagnoses. Because the final model suggests that a broad constellation of risk factors, reflective of diverse underlying constructs, are predictive of mortality, a more holistic approach to promote better health may be needed to reduce mortality.

Study Limitations 

This study had several limitations that limit the generalizability of the findings. First, we have established an association of health outcomes only with mortality, not causality. Second, the data were heavily left censored. Participants entered into the analysis an average of just under 10 years after SCI onset, so a substantial number of potential participants may have died prior to initiation of the data collection. Similar to the most recent MSCIS analysis,22 there were no data in the first postinjury year when many factors such as ventilator dependency, injury severity, and age exert their greatest influence on the likelihood of mortality. Third, NDI data were not available from January 1 to December 31, 2005, so there may be a small number of deaths that were not detected because of the 92.4% sensitivity of the SSDI among deceased persons.14 Fourth, we did not collect data on ventilator status, which is highly correlated with probability of mortality, particularly during the first few years after SCI onset.22 However, we do not anticipate a significant number of these participants given the left censoring (ie, because they have the shortest longevity, and many would have died before enrollment).

There were additional factors that limit the extent to which our design had the power to evaluate fully the relationships between health factors in mortality. First, the sample size was smaller in scope than that used with studies from the model SCI systems in the United States (these studies have focused primarily on biographic and injury characteristics). Although this limitation is balanced with the more detailed data on risk and protective health factors, larger sample sizes would have greater power to identify significant predictors of mortality. Second, recall bias may have influenced the accuracy of certain reports of health conditions during the prospective data collection. Inaccuracy in self-reported health would be an error in measurement and would weaken our ability to identify relationships between health and mortality. Last, risk and protective factors were only measured once with as much as an 8-year lag between collection and determination to mortality status. If data on these variables were collected more recently, they would serve as more precise predictors of mortality. For instance, probable major depression was a prominent predictor of mortality, yet these symptoms may represent chronic dysphoria consistent with a major depressive disorder (highly stable) or transient symptoms more related to situational circumstances. Similarly, the absence of these at one point in time does not preclude their later development (ie, a portion of participants reporting no symptoms would have developed symptoms at a later point in time). Therefore, it may be best to consider the current estimates as conservative.

Future Research 

Three types of research are suggested by the current findings. First, given that this study was based on a multistage model of prediction of mortality, and health and secondary conditions are only 1 step of factors from the model, future research should investigate additional risk and protective factors including behaviors, psychologic characteristics, and environmental factors. Investigation of these factors is particularly important because identification of risk prior to the development of complications will enhance prevention strategies. Second, we need to identify risk and protective predictors for the health factors that related to early mortality. This essentially speaks to the same issue: identification of risk and protective factors early in the chain of events leading to early mortality. The third recommendation for future research is both the most important and most challenging: development and evaluation of interventions to increase longevity after SCI. Our findings clearly may serve as at least a partial foundation for the development of interventions in that they identify people at high risk for early mortality and the health conditions that lead to the elevated risk of mortality. It is only through continued epidemiologic and clinical intervention research that long-term survival will likely be enhanced.

Conclusions 

return to Article Outline

Using a prospective design, we identified several types of secondary health conditions that were associated with hazard for mortality after SCI including surgeries for pressure ulcers, infection symptoms, UTIs, irreversible conditions (amputations and fractures), and the presence of probable major depression. Consideration of these factors led to substantial enhancements in prediction and classification beyond injury severity alone but only modest enhancements in prediction above a base model that included additional biographic and injury factors including age at injury. Although not demonstrating causation, the results suggest the need for prevention strategies that target people who present with multiple secondary conditions, including such basic intervention strategies as providing patients with concrete information about the association of their health with the likelihood of early mortality.

Suppliers

Acknowledgment 

return to Article Outline

The opinions here are those of the grantee and do not necessarily reflect those of the funding agencies.

References 

return to Article Outline

1. 1National Spinal Cord Injury Statistical Center. Annual statistical report. Birmingham: Univ Alabama; 2006;.

2. 2DeVivo MJ, Krause JS, Lammertse DP. Recent trends in mortality and causes of death among persons with spinal cord injury. Arch Phys Med Rehabil. 1999;80:1411–1419. Abstract | Full-Text PDF (1119 KB) | CrossRef

3. 3Strauss DJ, DeVivo MJ, Paculdo DR, Shavelle RM. Trends in life expectancy after spinal cord injury. Arch Phys Med Rehabil. 2006;87:1079–1085. Abstract | Full Text | Full-Text PDF (139 KB) | CrossRef

4. 4Burney RE, Maio RF, Maynard F, Karunas R. Incidence, characteristics, and outcome of spinal cord injury at trauma centers in North America. Arch Surg. 1993;128:596–599. MEDLINE

5. 5DeVivo MJ, Stover SL, Black KJ. Prognostic factors for 12 year survival after spinal cord injury. Arch Phys Med Rehabil. 1992;73:156–162. MEDLINE

6. 6Alander DH, Parker JD, Stauffer ES. Intermediate-term outcome of cervical spinal cord-injured patients older than 50 years of age. Spine. 1997;22:1189–1192. MEDLINE | CrossRef

7. 7Catz A, Thaleisnik M, Fishel B, et al. Survival following spinal cord injury in Israel. Spinal Cord. 2002;40:595–598. MEDLINE | CrossRef

8. 8Frankel HL, Coll JR, Charlifue SW, et al. Long-term survival in spinal cord injury: a fifty year investigation. Spinal Cord. 1998;36:266–274. MEDLINE

9. 9Lidal I, Snekkevik H, Aamodt G, Hjeltnes N, Biering-Sorensen F, Stanghelle J. Mortality after spinal cord injury in Norway. J Rehabil Med. 2007;39:145–151. MEDLINE | CrossRef

10. 10O'Connor PJ. Survival after spinal cord injury in Australia. Arch Phys Med Rehabil. 2005;86:37–47. Abstract | Full Text | Full-Text PDF (172 KB) | CrossRef

11. 11Prakash M, Raxwal V, Froelicher VF, et al. Electrocardiographic findings in patients with chronic spinal cord injury. Am J Phys Med Rehabil. 2002;81:601–608. MEDLINE | CrossRef

12. 12Shavelle RM, DeVivo MJ, Strauss DJ, Paculdo DR, Lammertse DP, Day SM. Long-term survival of persons ventilator dependent after spinal cord injury. J Spinal Cord Med. 2006;29:511–519. MEDLINE

13. 13DeVivo MJ, Ivie CS. Life expectancy of ventilator-dependent persons with spinal cord injuries. Chest. 1995;108:226–232. MEDLINE | CrossRef

14. 14Krause J, DeVivo M, Jackson A. Health status, community integration, and economic risk factors for mortality after spinal cord injury. Arch Phys Med Rehabil. 2004;85:1764–1773. Abstract | Full Text | Full-Text PDF (129 KB) | CrossRef

15. 15Whiteneck GG, Charlifue SW, Frankel HL, et al. Mortality, morbidity, and psychosocial outcomes of persons spinal cord injured more than 20 years ago. Paraplegia. 1992;30:617–630. MEDLINE

16. 16Soden RJ, Walsh J, Middleton JW, Craven ML, Rutkowski SB, Yeo JD. Causes of death after spinal cord injury. Spinal Cord. 2000;38:604–610. MEDLINE

17. 17DeVivo MJ, Stover SL. Long term survival and causes of death. In:  Stover SL,  DeLisa JA,  Whiteneck GG editor. Spinal cord injury: clinical outcomes from the model systems. Gaithersburg: Aspen; 1995;p. 289–316.

18. 18Garshick E, Kelley A, Cohen SA, et al. A prospective assessment of mortality in chronic spinal cord injury. Spinal Cord. 2005;43:408–416. MEDLINE | CrossRef

19. 19Myers J, Lee M, Kiralti J. Cardiovascular disease in spinal cord injury: an overview of prevalence, risk, evaluation, and management. Am J Phys Med Rehabil. 2007;86:142–152. MEDLINE | CrossRef

20. 20Krause JS. Survival following spinal cord injury: a fifteen-year prospective study. Rehabil Psychol. 1991;36:89–98.

21. 21Krause JS, Crewe NM. Prediction of long-term survival of persons with spinal cord injury: an 11 year prospective study. Rehabil Psychol. 1987;32:205–213.

22. 22Krause J, Kjorsvig J. Mortality after spinal cord injury: a four-year prospective study. Arch Phys Med Rehabil. 1992;73:558–563. MEDLINE

23. 23Krause JS. Secondary conditions and spinal cord injury: a model for prediction and prevention. Top Spinal Cord Inj Rehabil. 1996;2:217–227.

24. 24National Center for Health Statistics. National Death Index. http://www.cdc.gov/nchs/ndi.htmAccessed September 28, 2007.

25. 25Social Security Death Index (SSDI). http://www.ancestry.comAccessed September 28, 2007.

26. 26Foster S, Schubauer-Berigan M, Waters K. The specificity of the National Death Index and Social Security Administration Death Master File when information on Social Security number is lacking [abstract]. Am J Epidemiol. 2000;151(Suppl):S43.

27. 27Schnorr T, Steenland K. Underascertainment of deaths using the Social Security Administration Death Master File [abstract]. Am J Epidemiol. 1995;141(Suppl):S58.

28. 28Sesso HD, Paffenbarger RS, Lee IM. Comparison of National Death Index and World Wide Web death searches. Am J Epidemiol. 2000;152:107–111. MEDLINE | CrossRef

29. 29DeVivo M, Underhill A, Fine P. Accuracy of World-Wide-Web searches for persons with traumatic brain injury. Brain Inj. 2004;18:1155–1162. MEDLINE | CrossRef

30. 30Powell-Griner E, Anderson JE, Murphy W. State- and sex-specific prevalence of selected characteristics—Behavioral Risk Factor Surveillance System, 1994 and 1995. MMWR CDC Surveill Summ. 1997;46(3):1–31. MEDLINE

31. 31Krause JS. Dimensions of subjective well-being after spinal cord injury: an empirical analysis by gender and race/ethnicity. Arch Phys Med Rehabil. 1998;79:900–909. Abstract | Full-Text PDF (1385 KB) | CrossRef

32. 32Kemp BJ, Adams BM. The Older Adult Health and Mood Questionnaire: a measure of geriatric depressive disorder. J Geriatr Psychiatry Neurol. 1995;8:162–167. MEDLINE

33. 33Krause JS, Broderick L. Patterns of recurrent pressure ulcers after spinal cord injury: identification of risk and protective factors 5 or more years after onset. Arch Phys Med Rehabil. 2004;85:1257–1264. Abstract | Full Text | Full-Text PDF (97 KB) | CrossRef

34. 34Krause JS. Factors associated with risk for subsequent injuries after the onset of traumatic spinal cord injury. Arch Phys Med Rehabil. 2004;85:1503–1508. Abstract | Full Text | Full-Text PDF (68 KB) | CrossRef

35. 35Krause JS, Kemp BJ, Coker JL. Depression after spinal cord injury: relation to gender, ethnicity, aging, and socioeconomic indicators. Arch Phys Med Rehabil. 2000;81:1099–1109. Abstract | Full Text | Full-Text PDF (95 KB) | CrossRef

36. 36Hosmer DW, Lemeshow S. Applied logistic regression. 2nd ed. New York: John Wiley & Sons; 2000;.

37. 37Schoenfeld D. Chi-squared goodness-of-fit tests for the proportional hazards regression model. Biometrika. 1980;67:145–153. CrossRef

38. 38Harrell FE, Lee KL, Mark DB. Multivariate prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med. 1996;15:361–387. MEDLINE | CrossRef

39. 39Nagelkerke NJ. A note on a general definition of the coefficient of determination. Biometrika. 1991;78:691–692. CrossRef

40. 40Harrell FE. Evaluating the yield of medical tests. JAMA. 1982;247:2543–2546. MEDLINE

41. 41Harrell FE. Regression modeling strategies with applications to linear models, logistic regression, and survival analysis. New York: Springer-Verlag; 2001;.

a College of Health Professions, Medical University of South Carolina, Charleston, SC

b Department of Biostatistics, Bioinformatics, and Epidemiology, Medical University of South Carolina, Charleston, SC.

Corresponding Author InformationCorrespondence to James S. Krause, PhD, Dept of Rehabilitation Science, College of Health Professions, Medical University of South Carolina, 77 President St, Ste 117, PO Box 250700, Charleston, SC 29425

 Supported by the National Institute for Disability and Rehabilitation Research (grant no. H133G030117), Model Spinal Cord Injury Systems (grant no. H133N000005), and the National Institutes of Health (grant no. 1R01 NS 48117-01 A1).

No commercial party having a direct financial interest in the results of the research supporting this article has or will confer a benefit upon the authors or upon any organization with which the authors are associated.

Reprints are not available from the author.

 It was beyond our scope to address the other 3 components of the model: psychologic proxy variables, environmental factors, and risk and protective behaviors.

 Convention has been to use 4 groups based on the breakdown according to the 3 levels for ASIA grades A to C, with a single group denoting ASIA grade D regardless of injury level. We have used ambulatory status in lieu of ASIA grades that are not available and have broken this down according to cervical and noncervical injuries.

 DVT is technically reversible (the condition is corrected over time) but is not likely to be recurrent.

a Version 9.1.3; SAS Institute Inc, 100 SAS Campus Dr, Cary, NC 27513.

b Version 9.2; StataCorp, 4905 Lakeway Dr, College Station, TX 77845.

PII: S0003-9993(08)00322-5

doi:10.1016/j.apmr.2007.11.062


View previous. 13 of 40 View next.