|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Epidemiology and Outcomes |

,



,
* Renal, Electrolyte and Hypertension Division of the Department of Medicine,
Center for Clinical Epidemiology and Biostatistics, and
Department of Biostatistics and Epidemiology, University of Pennsylvania School of Medicine, Philadelphia, Pennsylvania, and
Winthrop University Hospital, Mineola, New York
Correspondence: Dr. Harold I. Feldman, Center for Clinical Epidemiology and Biostatistics, 923 Blockley Hall, 423 Guardian Drive, Philadelphia, Pennsylvania 19104. Phone: 215-898-0901; Fax: 215-898-0643; E-mail: hfeldman{at}mail.med.upenn.edu
| Abstract |
|---|
|
|
|---|
Design, settings, participants, & measurements: A retrospective cohort of 34,963 Fresenius Medical care dialysis patients from 1996 was assembled. Hemoglobin variability, absolute hemoglobin level, and temporal hemoglobin trend were measured over rolling 6-mo exposure windows. Their association with mortality was estimated using history-adjusted marginal structural analysis that adjusts for time-dependent confounding by applying weights to observations inversely related to the predictability of observed levels of hemoglobin.
Results: In the primary analysis, each g/dl increase in hemoglobin variability was associated with an adjusted hazard ratio (HR) [95% confidence interval (CI)] for all-cause mortality of 1.93 (1.20 to 3.10). Neither higher absolute hemoglobin level nor increasing hemoglobin trend were significantly associated with mortality; adjusted HR (95% CI) 0.85 (0.64 to 1.11) and 0.60 (0.25 to 1.45), respectively.
Conclusions: Marginal structural analysis demonstrates that hemoglobin variability is associated with increased mortality among chronic hemodialysis patients, and that this effect is more pronounced than appreciated using standard statistical techniques that do not take time-dependent confounding into account.
| Introduction |
|---|
|
|
|---|
As is the case in all nonrandomized research, the observed association between Hgb-Var and death was subject to potential confounding. Confounding occurs when a variable (or variables) affects the likelihood of exposure (e.g. Hgb-Var) and of outcome (e.g. death), but does not serve as an intermediate between exposure and outcome. Standard methods of statistical adjustment are often sufficient to adjust for these confounders when they do not vary over time compared with their baseline values.
Unfortunately, there are instances when application of standard statistical techniques may yield biased associations, such as in the presence of time-dependent confounding. Time-dependent confounders (TDC) are variables that have two different roles that are difficult to distinguish: (1) they influence subsequent exposure and outcome like confounders that are fixed over time, and (2) they also serve as intermediates (or links) between exposure and outcome. Although it is appropriate to adjust for their role as confounders, it is inappropriate to adjust for their role as intermediates on the causal pathway from exposure to outcome. If standard statistical techniques are used to adjust for TDC, they adjust for both roles simultaneously and, therefore, may underestimate true associations owing to the adjustment for the role as a causal pathway factor. For example, if Hgb-Var is both influenced by comorbid conditions and promotes these same conditions as part of the mechanism through which it causes poor clinical outcomes, simple adjustment for these comorbidities would lead to biased estimates of the Hgb-Var outcome relationship. Figure 1 serves to illustrate one potential TDC of the association between Hgb-Var and death, namely intravenous iron administration, which may independently be associated with Hgb-Var and mortality (Figure 1A), and may also serve on the causal pathway (Figure 1B).
|
Considering that the prior association may have been biased by the presence of time-dependent confounders, we have reanalyzed our data using HA-MSM to better estimate the causal association between Hgb-Var and mortality.
| Materials and Methods |
|---|
|
|
|---|
Demographic data on age, gender, race, duration of ESRD and diabetic status (adult, juvenile, none) were collected at baseline. Hospitalization data and data on infections were not available. Laboratory values and medication exposure data, including hemoglobin, urea reduction ratio, Kt/V, albumin, bicarbonate, aspartate aminotransferase, calcium, phosphate, intact parathyroid hormone, ferritin, transferrin saturation, iron level, intravenous iron dose, and erythropoietin dose were updated monthly. Data on activated vitamin D treatment were not available.
The primary exposures of interest were Hgb-Var, absolute level of hemoglobin, and temporal trend in hemoglobin. To derive these values, a best-fit (ordinary least squares) line of hemoglobin over time was plotted for each subject: absolute hemoglobin was the intercept of this line; temporal hemoglobin trend was the slope of this line; Hgb-Var was the residual SD (the SD of the individual distances between observed hemoglobin values and the line) (5). Hgb-Var cannot be defined at a single point in time. To balance between maximizing accuracy of Hgb-Var measurements (which is enhanced by longer-exposure windows) and minimizing delay between Hgb-Var exposure and outcome (which is enhanced by narrower windows), we defined hemoglobin parameters over 6-mo windows.
The outcome of interest was death from any cause. Subjects were followed from the first date of their qualifying enrollment period (January 1, 1996 or July 1, 1996); follow-up time began on the first day of the tenth month (allowing for a 3-mo covariate assessment period followed by a 6-mo hemoglobin exposure window), and continued until death or censoring occurred. There were a total of 527,967 patient-months of potential follow-up. Censoring criteria included transfer of care, renal transplantation, or study end: September 30, 1998.
All analyses were performed in SAS 9.1 (Cary, NC).
Outcome Model
The primary outcome was based on a logistic regression analog of the history-adjusted marginal structural Cox proportional hazards model (8,11). To fit this model, each subject-month was considered independently. (Robust variance estimates were used to account for nonindependence of observations within subject.) For each subject-month (t), survival was predicted based on the hemoglobin parameters (Hgb-Var, absolute level of hemoglobin, and hemoglobin trend) calculated over the preceding 6-mo window (t – 6 through t – 1). Adjustment was made for covariates measured over the 3-mo period preceding this window (months t – 9 through t – 7) (Figure 2). Covariates were considered before the hemoglobin exposure window because prior values of these covariates could be confounders, but could not serve as causal pathway intermediates. Concurrent values of covariates (those from months t – 6 to t – 1), which may have either "confounding" or "pathway intermediate" effects were adjusted for by application of weights (see Appendix). Analysis was restricted to subject-months for which all associated hemoglobin and covariate data were available (260,028 and 390,621 follow-up months for weighted and unweighted models, respectively). To facilitate clinical interpretation, HR are expressed per g/dl (Hgb-Var and absolute hemoglobin) or g/dl/mo (hemoglobin trend). (Given that 1 g/dl has different meaning for each parameter, caution is advised in contrasting the relative effects of each parameter on mortality.)
|
To examine the effects of marginal structural analysis, we fit unweighted models (otherwise analogous to the HA-MSM) as comparators.
| Results |
|---|
|
|
|---|
During follow-up, 9160 subjects (26.2%) died, 9965 (28.5%) transferred care or were lost to follow-up, and 15,838 (45.3%) survived to study end.
Hgb-Var was distributed with a mean (SD) of 0.60 (0.33) g/dl and median [(interquartile range (IQR)] of 0.53 (0.36 to 0.76) g/dl. Absolute level of hemoglobin was distributed with mean (SD) 10.25 (1.25) g/dl. Temporal trend in hemoglobin was distributed with a mean (SD) of 0.02 (0.29) g/dl/mo (5). [To facilitate clinical interpretation, HR (below) are reported in clinical units (g/dl or g/dl/mo). A difference of 1 g/dl in Hgb-Var represents a greater relative difference than a difference of 1 g/dl in absolute hemoglobin level, owing to the tighter distribution (smaller SD) of the former.]
Association between Hemoglobin Parameters and Death
In the marginal structural analysis, each g/dl rise in Hgb-Var was associated with an adjusted HR (95% CI) for death of 1.93 (1.20 to 3.10). Using an unweighted model, which does not address time-dependent confounding, the HR (95% CI) was 1.34 (1.24 to 1.44; Table 1).
|
In the marginal structural analysis, each g/dl/mo increase in temporal hemoglobin trend was associated with a HR (95% CI) for death of 0.60 (0.25 to 1.45). Using an unweighted mode, the HR (95% CI) was 0.28 (0.24 to 0.32).
| Discussion |
|---|
|
|
|---|
Several reports have demonstrated substantial Hgb-Var among dialysis patients (1–4). This observed variability may derive from patterns of care, and from comorbidities that influence sensitivity to these therapies. Only one previous report has examined the association between Hgb-Var and mortality, demonstrating a 33% increase in mortality per g/dl increase in Hgb-Var (5). (These findings were very similar to those yielded by the unweighted model in this study, differing only because of use of rolling versus static exposure windows.) However, findings from that study were likely biased by the effects of time-dependent confounding, explaining the discrepancy between its findings and those in study presented here.
Other than in chronic diseases such as kidney failure, hemoglobin levels are typically maintained within a narrow range to ensure consistent oxygen delivery to peripheral tissues. Repeated oscillations in hemoglobin may result in interruption of tissue oxygen delivery, with resulting ischemic tissue damage and inappropriate activation of cardiovascular compensatory mechanisms such as cardiac myocyte hypertrophy and left ventricular hypertrophy and dilation (24–26). In addition, hemoglobin variation has been shown to promote autonomic dysfunction, which in turn predisposes patients to sudden death (27,28).
History-adjusted marginal structural models use inverse probability of exposure weighting (or a stabilized analog) to adjust for the confounding effects of TDC, but not their pathway intermediate effects (8,10,11). In the absence of unmeasured confounders, weighting in this manner eliminates bias due to time-dependent confounding. Unfortunately, their application simultaneously creates imprecision (i.e. widens CI), as is evident comparing the CI width between marginal structural and unweighted analyses.
Comparing findings from HA-MSM and unweighted models, the association between Hgb-Var and death appears potentiated: HR (95% CI) 1.93 (1.20 to 3.10) versus 1.34 (1.24 to 1.44), respectively. Intuitively, this makes sense considering that patients with greater comorbid disease burdens are expected to have more variable response to exogenous erythropoietin (i.e. greater Hgb-Var), and also to be at increased risk of death. Conversely, the protective effects of absolute hemoglobin level and temporal hemoglobin trend become attenuated as weighting increases. Again, this makes intuitive sense considering that healthier patients are expected to have greater hemoglobin levels and more positive hemoglobin trends, and be at decreased risk of death.
Reassurance as to the effectiveness of HA-MSM to appropriately control for TDC can be found by analogy. Observational studies have demonstrated a significant and potent association between absolute hemoglobin level and mortality, (1,14–17) a finding not borne out by randomized trials (18,19). In examining our data, the association between higher hemoglobin levels and mortality becomes attenuated upon use of marginal structural analysis, and no longer reaches statistical significance, consistent with the results of randomized trials.
Other examples in the literature where the related marginal structural models have been applied include analysis of the Multicenter AIDS Cohort study, in which marginal structural analysis confirmed that zidovudine treatment was associated with improved survival among HIV-positive men, whereas standard statistical methods suggested it was harmful (11). Similarly, in a post hoc analysis of the Physicians Health Study, marginal structural analysis demonstrated that aspirin use was more beneficial in preventing cardiovascular mortality than appreciated by traditional statistical techniques (20). In the renal literature, marginal structural models have been used to confirm a survival benefit associated with activated vitamin D treatment of chronic dialysis patients (21), to demonstrate the safety of intravenous iron administration, (22) and to examine the effects of carnitine supplementation on hospitalization in dialysis patients (23).
Several limitations of the study presented here should be noted. First, the success of HA-MSM to estimate causal associations is contingent upon the untestable assumption that no residual confounding is present; this is especially true of "endogenous" variables such as hemoglobin that are not under direct human control (8,11). Through inverse probability of exposure (IPE) weighting, we attempted to account for many candidate TDC, but the possibility exists that there were others for which we could not fully account. This limitation is particularly relevant here given the limited data available on comorbid disease, infection, or hospitalization.
Second, as with all observational studies, there was the opportunity for information and selection biases. The latter is particularly concerning given the survival requirement imposed by our metric of Hgb-Var.
Third, the cohort we studied is nearly a decade old, and there have been temporal changes in the management of anemia in dialysis patients since that time. Although we know of no reason why the fundamental relationship between Hgb-Var and absolute level of hemoglobin and death would be different now as compared with a decade ago, we cannot exclude this possibility.
Fourth, although the population studied was geographically and demographically diverse, all subjects were treated in one network of dialysis providers. Generalization of results to other populations or settings should be undertaken cautiously. In particular, although the results presented here suggest that reducing Hgb-Var might be beneficial, it is premature to conclude that specific approaches to doing this would lead to reductions in mortality.
Finally, it should be noted that no statistical technique is a valid substitute for well done randomized clinical trials when such trials are feasible. However, history-adjusted marginal structural modeling (as well as the related marginal structural modeling) is one potential analytical technique that can be applied to observational data when randomization is infeasible.
In conclusion, this study demonstrates that hemoglobin variability is associated with increased mortality among chronic hemodialysis patients using marginal structural analysis, which controls for time-dependent confounding in a fashion superior to traditional multivariable statistical analyses. Furthermore, this association was more pronounced than those demonstrated using methods that do not take time-dependent confounding into account. If confirmed in other observational studies that also permit generalization beyond the population studied here, hemoglobin variability may ultimately become an important therapeutic target. Before this can occur, clinical trials will need to demonstrate the superiority of anemia management strategies that reduce hemoglobin variability. Studies such as this provide further impetus to develop and implement such trials.
| Appendix |
|---|
|
|
|---|
9.7 g/dl, 9.7 to 10.6 g/dl, 10.6 to 11.3 g/dl, and >11.3 g/dl. We selected polychotomous logistic regression rather than dichotomous logistic regression or ordinal logistic regression because the former does not provide for sufficient resolution of hemoglobin level, whereas the latter makes stronger assumptions about the form of the association between predictors and hemoglobin level that may not apply. The predictor variables in this model were laboratory measures and medication exposure from the previous 3 mo, and hemoglobin concentration over each of the prior 6 mo. We used stabilized weights in place of raw weights, because they generally produce more stable estimates (or "estimates with higher precision") (8,11). Stabilized weights are a ratio: the denominator is precisely the probability of exposure outlined above; the numerator is the probability of exposure as predicted by the prior 6-mo hemoglobin concentration and laboratory and medication data from the 3 mo before the hemoglobin exposure window.
Censoring Weights
Discontinued or censored follow-up occurred in our study cohort at the time of death, at the end of the study, or when a subject was lost to follow-up. The latter censoring event may have been associated with the exposure or outcome (mortality) (i.e. informative censoring), thereby, creating opportunities for bias. This was minimized by estimating a stabilized censoring weight for each subject using models analogous to those for stabilized IPE weights, except where the outcome for each month was censoring (1 = yes, 0 = no) (8,11). The weights that we applied in the outcome analysis were the product of the stabilized IPE times censoring weight (8,11).
| Disclosures |
|---|
|
|
|---|
| Acknowledgments |
|---|
| Footnotes |
|---|
Received October 8, 2007. Accepted February 1, 2008.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. J. Arneson, D. Zaun, Y. Peng, C. A. Solid, S. Dunning, and D. T. Gilbertson Comparison of methodologies to characterize haemoglobin variability in the US Medicare haemodialysis population Nephrol. Dial. Transplant., May 1, 2009; 24(5): 1378 - 1383. [Full Text] [PDF] |
||||
![]() |
A. L. M. de Francisco, P. Stenvinkel, and S. Vaulont Inflammation and its impact on anaemia in chronic kidney disease: from haemoglobin variability to hyporesponsiveness NDT Plus, January 1, 2009; 2(suppl_1): i18 - i26. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. M. Brunelli, K. E. Lynch, E. D. Ankers, M. M. Joffe, W. Yang, R. I. Thadhani, and H. I. Feldman Association of Hemoglobin Variability and Mortality among Contemporary Incident Hemodialysis Patients Clin. J. Am. Soc. Nephrol., November 1, 2008; 3(6): 1733 - 1740. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |