Lung protective mechanical ventilation and two year survival in patients with acute lung injury: prospective cohort studyBMJ 2012; 344 doi: http://dx.doi.org/10.1136/bmj.e2124 (Published 05 April 2012) Cite this as: BMJ 2012;344:e2124
- Dale M Needham, associate professor12,
- Elizabeth Colantuoni, assistant scientist3,
- Pedro A Mendez-Tellez, assistant professor4,
- Victor D Dinglas, research programme supervisor1,
- Jonathan E Sevransky, assistant professor1,
- Cheryl R Dennison Himmelfarb, associate professor6,
- Sanjay V Desai, assistant professor1,
- Carl Shanholtz, associate professor7,
- Roy G Brower, professor1,
- Peter J Pronovost, professor456
- 1Outcomes After Critical Illness and Surgery Group, Division of Pulmonary and Critical Care Medicine, Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
- 2Department of Physical Medicine and Rehabilitation, Johns Hopkins University School of Medicine
- 3Department of Biostatistics, Johns Hopkins University Bloomberg School of Public Health, Baltimore, MD, USA
- 4Department of Anesthesiology and Critical Care Medicine, Johns Hopkins University School of Medicine
- 5Department of Health Policy and Management, Johns Hopkins University Bloomberg School of Public Health, Baltimore, MD, USA
- 6Johns Hopkins University School of Nursing, Baltimore, MD, USA
- 7Division of Pulmonary and Critical Care Medicine, University of Maryland, Baltimore, MD, USA
- Correspondence to: D M Needham
- Accepted 27 February 2012
Objective To evaluate the association of volume limited and pressure limited (lung protective) mechanical ventilation with two year survival in patients with acute lung injury.
Design Prospective cohort study.
Setting 13 intensive care units at four hospitals in Baltimore, Maryland, USA.
Participants 485 consecutive mechanically ventilated patients with acute lung injury.
Main outcome measure Two year survival after onset of acute lung injury.
Results 485 patients contributed data for 6240 eligible ventilator settings, as measured twice daily (median of eight eligible ventilator settings per patient; 41% of which adhered to lung protective ventilation). Of these patients, 311 (64%) died within two years. After adjusting for the total duration of ventilation and other relevant covariates, each additional ventilator setting adherent to lung protective ventilation was associated with a 3% decrease in the risk of mortality over two years (hazard ratio 0.97, 95% confidence interval 0.95 to 0.99, P=0.002). Compared with no adherence, the estimated absolute risk reduction in two year mortality for a prototypical patient with 50% adherence to lung protective ventilation was 4.0% (0.8% to 7.2%, P=0.012) and with 100% adherence was 7.8% (1.6% to 14.0%, P=0.011).
Conclusions Lung protective mechanical ventilation was associated with a substantial long term survival benefit for patients with acute lung injury. Greater use of lung protective ventilation in routine clinical practice could reduce long term mortality in patients with acute lung injury.
Trial registration Clinicaltrials.gov NCT00300248.
Survivors of severe critical illness, such as acute lung injury, and its more severe subset, acute respiratory distress syndrome (ARDS), commonly experience increased mortality and morbidity in the months and years after hospital discharge.1 2 3 Compared with age matched and sex matched controls, patients discharged from intensive care are two to five times more likely to die during three to 15 years’ follow-up.3 Few interventions have been evaluated for improving this increased long term mortality.4
Randomised trials and meta-analyses have shown that use of volume limited and pressure limited mechanical ventilation (lung protective ventilation) in patients with acute lung injury substantially decreases short term mortality.5 6 7 8 9 A randomised trial of lung protective ventilation carried out by the ARDS Network 8 found an 8.8% absolute reduction in short term mortality. This trial evaluated a ventilator tidal volume of 6 mL/kg predicted body weight (calculated on the basis of a patient’s sex and height8) and a plateau pressure (airway pressure measured after a 0.5 second end inspiratory pause) of ≤30 cm of water compared with a tidal volume of 12 mL/kg predicted body weight and a plateau pressure of ≤50 cm of water. Understanding the effect of lung protective ventilation on long term survival is important,4 5 especially since critical care interventions with a mortality benefit at hospital discharge may not have a sustained benefit thereafter.10 Lung protective ventilation may have long term benefits owing to its reductions in inflammation, ventilator induced lung injury, and duration of non-pulmonary organ failure in the setting of the intensive care unit.8 11
We evaluated the association of lung protective ventilation and average ventilator tidal volume with two year survival in patients with acute lung injury. We carried out this evaluation using a prospective cohort study design to understand the epidemiology and effects of lung protective ventilation as part of routine medical care and to avoid the ethical issues associated with randomly allocating patients to higher tidal volumes and pressures as part of a randomised trial.
From 2004-7 we enrolled into a prospective cohort study 520 consecutive mechanically ventilated patients with acute lung injury, diagnosed according to the American-European Consensus criteria (box).12 Patients were recruited from 13 medical, surgical, and trauma intensive care units at four academic teaching hospitals in Baltimore, Maryland, United States. Intensive care units that specialised in patients with neurological conditions, and patients with acute lung injury with primary neurological disease or head trauma were not eligible. We excluded patients who had pre-existing comorbid illness with a life expectancy of six or less months (for example, metastatic cancer), pre-existing cognitive impairment or communication/language barriers, no fixed address for follow-up, been transferred from another hospital and had pre-existing acute lung injury of more than 24 hours’ duration, been mechanically ventilated for more than five days before onset of acute lung injury, previous lung resection, and a physician order for no escalation of care in the intensive care unit (for example, no vasopressors) at onset of acute lung injury. We obtained informed consent for prospective follow-up in writing after patients regained decision making capacity (or through proxies if patients remained incapable of decision making).13
American-European Consensus Conference criteria for acute lung injury and acute respiratory distress syndrome12
Timing: acute onset
Oxygenation: PaO2:FiO2 ≤300 mg Hg (acute lung injury) or ≤200 mg Hg (acute respiratory distress syndrome)
Frontal chest radiography findings: bilateral infiltrates
No clinical evidence of left atrial hypertension or pulmonary artery wedge pressure ≤18 mm Hg
Assessment of primary outcome: mortality
We prospectively evaluated patients at 3, 6, 12, and 24 months after onset of acute lung injury. Mortality status and date of death were obtained from family members and subsequently verified using a commercial version of the social security death master file14 or other publicly accessible records of death. Any patient for whom we could not determine mortality status was censored at the last day known to be alive.
Assessment of primary exposure: lung protective mechanical ventilation
After onset of acute lung injury, we recorded mechanical ventilator settings twice daily at 6 am and 6 pm for the duration of the patient’s ventilation at the study site hospital. A ventilator setting was defined as being “eligible” for lung protective ventilation if mechanical ventilation was delivered through an endotracheal or tracheostomy tube with a fraction of inspired oxygen (FiO2) ≥0.50 or positive end expiratory pressure >5 cm of water. We chose this definition of eligible because it approximated the threshold at which the ARDS Network ventilation protocol permitted trials of spontaneous breathing in which tidal volume and plateau pressure were not specified and controlled, and discontinuation of mechanical ventilation could be evaluated.
A ventilator setting was defined as “adherent” to lung protective ventilation if it was an eligible setting that met two criteria: tidal volume ≤6.5 mL/kg predicted body weight (the threshold used in the ARDS Network tidal volume trial to designate study site adherence to the goal tidal volume of 6.0 mL/kg predicted body weight) and a plateau pressure of ≤30 cm of water (based on documentation of respiratory therapy; in pressure regulated modes of ventilation, in which plateau pressure is not measured, we used peak pressure or the sum of positive end expiratory pressure and the prescribed increment in inspiratory pressure). We considered ventilator settings that used a mode of mechanical ventilation for which tidal volume is not routinely measured (that is, high frequency oscillation or airway pressure release modes) to be not eligible for assessing adherence to lung protective ventilation but we still included them in the statistical analysis using the cumulative number of each of these ventilator modes as time varying covariates. We excluded from analyses those patients without at least one eligible ventilator setting before death or discontinuation of mechanical ventilation because adherence to lung protective ventilation is undefined in this patient subset.
Assessment of baseline and time varying covariates
In this study we measured covariates potentially associated with ventilator settings or mortality, based on previous literature (for example,15) and a systematic process.16 We evaluated 10 baseline variables: age, sex, body mass index category, comorbidity (Charlson index17), severity of illness within 24 hours of admission to an intensive care unit (acute physiology and chronic health evaluation II score),18 primary risk factor for acute lung injury, type of intensive care unit (medical versus surgical), patient location before admission to intensive care (for example, emergency department), year of enrolment to study (first, second, or third), and study site identifier (1, 2, 3, or 4). Throughout the duration of ventilation we obtained time varying covariates from the medical record twice daily at 6 am and 6 pm, when recorded: positive end expiratory pressure, pressure of arterial oxygen (PaO2), FiO2, arterial pH, and actual respiratory rate. From the data we calculated and included in the model a twice daily measure of static compliance of the respiratory system.19 Additional time varying covariates were measured once daily in the intensive care unit, including organ dysfunction (sequential organ failure assessment score20), sedation and delirium status (Richmond agitation and sedation scale,21 and confusion assessment method for the intensive care unit22), dose of systemic corticosteroids and neuromuscular blocking agents, and net fluid balance (total fluid input minus total fluid output).
We used descriptive statistics to summarise the data and compared these using Wilcoxon rank-sum and Fisher’s exact tests, as appropriate. A Kaplan-Meier curve was used to describe the survival experience over the two years of follow-up.
A multivariable Cox proportional hazards regression model with time varying covariates was used to estimate the hazard ratios of mortality as a function of lung protective ventilation (the primary exposure variable) and all 27 covariates chosen a priori based on clinical relevance. The primary exposure was modelled as the time varying number of ventilator settings adherent to lung protective ventilation, with adjustment for the total duration of mechanical ventilation. When a similar magnitude of hazard ratios existed between categories we collapsed the categorical variables to binary variables. We examined Martingale residual plots to confirm appropriate modelling of continuous covariates, with no departures from linearity detected. Analysis of variance inflation factors confirmed the absence of any important multicollinearity in the full multivariable model. The Cox model allowed for separate baseline hazard functions for mortality during stay in the intensive care unit compared with after stay, as defined by a time varying indicator of discharge from the intensive care unit. We assessed the proportional hazard assumption for the full model using a trend test based on the Schoenfeld residuals23 and graphical displays of the cumulative Schoenfeld residuals versus event times for each explanatory variable, indicating no important violation of this assumption.24
To illustrate the effect of lung protective ventilation on mortality, we used the multivariable Cox model to predict the survival function and cumulative risk of mortality at 1, 6, 12, 18, and 24 months after onset of acute lung injury for a prototypical patient who survived for at least four days of eligible mechanical ventilation and had median values for all continuous covariates and mode values for all binary covariates.
For days in the intensive care unit during which a sedation or delirium assessment was not completed, we made a conservative assumption that coma (Richmond agitation and sedation scale score −4 or −5) and delirium were not present. For missing data on plateau pressure we used multiple imputation25 (with five datasets), as in previous research using ARDS Network data.26 As a secondary analysis, we repeated all analyses using a subset of the entire dataset excluding patients with complete missing ventilator data required for measuring the primary exposure to compare with results from the primary analysis using multiple imputation for missing data.
We carried out additional a priori analyses. To determine if the findings were robust to the tidal volume and plateau pressure thresholds that defined an adherent ventilator setting, we reran the multivariable Cox model using different thresholds for adherence both for tidal volume (≤6.0 to ≤9.0 mL/kg predicted body weight) and for plateau pressure (≤25 to ≤45 cm of water). In addition, we explored the association of mean tidal volume (analysed as a time varying variable over the duration of mechanical ventilation) with two year survival, modelling mean tidal volume as both a three level categorical exposure variable (<6.5, 6.5-8.5, and >8.5 mL/kg predicted body weight) and a continuous exposure variable. Using the same multivariable regression model, we also evaluated whether mean tidal volume, considered as a continuous exposure variable without accounting for plateau pressure, had a linear association with two year survival. In this evaluation, we calculated the estimated change in hazard of mortality across mean tidal volume between the multivariable Cox model assuming a linear relation with survival compared with assuming a flexible, non-linear relation achieved by adding a cubic spline (with three degrees of freedom) into the model.
We also carried out a post hoc propensity score analysis approach for non-binary treatments27 to compare with the primary analysis. Separate propensity scores were created using Poisson regression models for each of the number of eligible ventilator settings and the number of adherent ventilator settings using the previously described covariates for the multivariable regression models. Data were stratified into nine strata based on three levels of each propensity score, with a Cox model fit within each stratum with further adjustment for a subset of the relevant covariates. We calculated the overall propensity score adjusted hazard ratio and associated variance as a weighted average of the strata specific estimates.
We calculated that we would require a sample size of 520 patients to show a relative hazard of 0.70 for two year survival comparing two equally sized groups of patients (those receiving greater than compared with less than the median proportion of eligible ventilation days using lung protective ventilation) with an α of 0.05 and statistical power of more than 85%. Statistical significance was defined as a two sided P<0.05. All analyses were completed using R statistical software (version 2.11.1).
Overall, 35 of 520 patients were excluded from analysis: 32 (6%) had no eligible ventilator settings and three (<1%) had missing data on height required for calculating predicted body weight and adherence to lung protective ventilation (fig 1⇓). Thus data from 485 patients, with 6240 eligible ventilator settings and 12 202 total ventilator settings, were available for analysis. For the secondary analysis of complete data, 30 (6%) of the 485 eligible patients were excluded owing to missing ventilator variables required for measuring the primary exposure, leaving 455 patients with 4938 eligible ventilator settings and 10 321 total ventilator settings.
Tables 1⇓ and 2⇓ present the patients’ characteristics, ventilator variables, and other covariates related to the intensive care unit, by adherence to lung protective ventilation and by two year mortality status. Of 485 patients, 274 (57%) were male and the median age (years) was 53 (interquartile range 42-63). Patients were predominantly admitted to a medical (compared with surgical) intensive care unit, with a median severity of illness (acute physiology and chronic health evaluation II) score of 27 (interquartile range 20-33, table 1). The median number of eligible ventilator settings per patient was 8 (interquartile range 3-15). The median duration (days) of mechanical ventilation was 9 (interquartile range 5-17), stay in the intensive care unit was 13 (8-22), and hospital stay was 21 (13-36).
Of 6240 eligible ventilator settings, 2548 (41%) were adherent to lung protective ventilation, with no significant temporal trends in adherence observed over the three year study period. Of the 485 patients, 417 (86%) received 50% or less and 68 (14%) received more than 50% of their twice daily ventilator settings adherent to lung protective ventilation. A total of 180 (37%) patients never received lung protective ventilation during any of their twice daily ventilator settings. Patients who were ever (versus never) exposed to lung protective ventilation were younger, and a significantly greater proportion were male and received neuromuscular blockers (table 1).
Overall, 311 (64%) of the 485 patients died during the two years after onset of acute lung injury. Mortality status at two years was not available for 29 (6.0%) of the 485 patients, and these patients were censored in the analysis. The mortality rate increased over time, particularly over the first year of follow-up, from 44% at 30 days to 52% at 90 days and 62% at one year (fig 2⇓). Those who were alive after two years were significantly younger with lower comorbidity, severity of illness, and organ failure scores (table 1).
During the two year follow-up, survival was independently associated with several covariates: younger age; lower comorbidity, organ failure scores, cumulative fluid balance, and days with neuromuscular blockers; and a greater number of high frequency oscillation or airway pressure release ventilator settings (table 3⇓). After adjusting for the total duration of mechanical ventilation and all other covariates, for each additional adherent ventilator setting measured twice daily, the risk of mortality over two years decreased by 3% (hazard ratio 0.97, 95% confidence interval 0.95 to 0.99, P=0.002). For a prototypical patient, this finding equates to an estimated absolute risk reduction in two year mortality of 4.0% (95% confidence interval 0.8% to 7.2%, P=0.012) for 50% ventilator adherence and 7.8% (1.6% to 14.0%, P=0.011) for 100% adherence compared with a 49.7% baseline mortality under the assumption of no adherence to lung protective ventilation (fig 3⇓).
The results of the primary analysis of 485 patients were robust to varying both the tidal volume and the plateau pressure thresholds used to define adherence to lung protective ventilation. In addition, the primary results using multiple imputation for missing data were similar when replicated in the secondary analysis of 455 patients, which excluded those with complete missing data for the primary exposure, with a hazard ratio for each additional adherent ventilator setting of 0.96 (95% confidence interval 0.94 to 0.99, P=0.001). Finally, results of the post hoc propensity score analyses were consistent with results from the primary analysis.
Compared with a mean tidal volume <6.5 mL/kg predicted body weight, the adjusted hazard ratios for two year mortality for a mean tidal volume of 6.5 to 8.5 mL/kg predicted body weight was 1.59 (1.19 to 2.14, P=0.001) and for >8.5 mL/kg predicted body weight was 1.97 (1.23 to 3.16, P=0.004). Moreover, when tidal volume was modelled as a continuous variable, no evidence supported a non-linear relation between two year survival and mean tidal volume (fig 4⇓; P=0.182 for the non-linear terms of a cubic spline model), with an adjusted hazard ratio of 1.18 (1.07 to 1.31, P=0.001), indicating an 18% relative increase in mortality for each 1 mL/kg predicted body weight increase in average tidal volume.
In this multisite, prospective cohort study of patients with acute lung injury, adherence to lung protective ventilation was associated with improved two year survival, with an absolute reduction in mortality of about 4% and 8% when 50% and 100% adherence, respectively, were compared with no adherence. These findings were robust to differing thresholds of tidal volume and plateau pressure used to define adherence to lung protective ventilation. Moreover, average tidal volume showed an independent linear relation with two year survival, with an 18% relative increase in the risk of mortality over two years for each 1 mL/kg predicted body weight increase in average tidal volume over the duration of mechanical ventilation. Thus, within the setting of routine clinical practice, use of lung protective ventilation was strongly associated with a substantial, long term survival benefit for patients with acute lung injury.
Comparison with other studies
No existing studies of lung protective ventilation have been designed to evaluate long term survival in patients.6 The high two year mortality observed in this study is consistent with previous research showing increased long term mortality for patients after discharge from an intensive care unit versus age and sex matched population controls, and with the known risks for long term mortality, including high comorbidity, organ failure, and severity of illness,3 28 29 which were observed in this patient sample.
This observational study allowed evaluation of lung protective ventilation within a setting of routine clinical practice in four teaching hospitals. Only 41% of all eligible ventilator settings were adherent to lung protective ventilation, and when twice daily collected ventilator settings were evaluated 37% of patients never received lung protective ventilation. These results represent an improvement from previous studies30 31 32 and from usual care practice before randomisation in the original ARDS Network trial in which mean baseline tidal volume was 10 mL/kg predicted body weight.8 There are many known reasons for non-adherence to lung protective ventilation31 33 34 35 and still substantial room for improvement with use of this life saving intervention.6 36 37 Rigorous knowledge translation research, aimed at improving the implementation of clinical research into practice, is needed to maximise the public’s return on investment from clinical and preclinical research that established the short term efficacy of lung protective ventilation.37 These knowledge translation efforts are especially important given the ageing population and the higher incidence and mortality associated with acute lung injury and mechanical ventilation with increasing age.38 39 40
In this study, the duration of neuromuscular blockade (but not cumulative dose) was associated with an increased risk of mortality over two years. These findings contrast with the 90 day survival benefit of early neuromuscular blockade in severe acute respiratory distress syndrome from a recent randomised trial.41 This difference in findings is possibly due, in part, to noticeably different use of these drugs in this study compared with the randomised trial protocol.
This study also found a two year survival benefit with use of high frequency oscillation and airway pressure release modes of ventilation in selected patients, with adjusted hazard ratios for each additional high frequency oscillation or airway pressure release ventilator setting (measured twice daily) of 0.93 (95% confidence interval 0.88 to 0.97, P=0.002) and 0.98 (0.96 to 0.99, P=0.005), respectively. However, the study was not specifically designed and analysed to evaluate these alternative modes of ventilation, which have not been compared with conventional lung protective ventilation in any large scale randomised trial.6 Hence the results for these alternative modes of ventilation must be interpreted with caution.
Limitations of the study
This study has potential limitations. Firstly, as this study was observational we cannot prove causation between lung protective ventilation and increased two year survival since both measured and unmeasured characteristics differed between patient groups receiving and not receiving lung protective ventilation. Given the decrease in short term mortality established in previous studies,5 6 7 8 9 however, a randomised trial to evaluate the long term effects of lung protective ventilation would not be ethical. Causality is plausible given the known short term benefits of lung protective ventilation6 7 8 9 11 and the dose-response effect of lung protective ventilation observed in this study. Secondly, the plateau pressure required for determining adherence to lung protective ventilation was not always recorded in the medical record, occurring in 21% of the eligible ventilator settings, with a median number of missing values across all 485 patients of 1 (interquartile range 0-3). This is not an uncommon finding in routine clinical practice, with plateau pressure, peak end expiratory pressure, or tidal volume missing in about 40% of patients in both a population based observational study42 and in the ARDS Network trial before randomisation of the groups.8 These missing data have the potential to bias results; however, both the primary analysis with imputation of missing data and a secondary complete case analysis showed similar results, which provide some reassurance. Greater efforts to routinely measure and record plateau pressure and tidal volume in mL/kg predicted body weight, and to explicitly consider these ventilator variables during bedside rounds, may improve adherence to lung protective ventilation.33 36 Finally, only teaching hospitals were evaluated and exclusion criteria were applied in selecting eligible patients for this study; hence generalisability of the findings may be limited. However, four hospitals with 13 intensive care units participated, providing variation in the medical care provided to participants, and the number of exclusion criteria were limited (compared with randomised trials of lung protective ventilation) which aids in generalisability.
Lung protective mechanical ventilation is associated with a substantial survival benefit for patients with acute lung injury over two year’s follow-up. Average tidal volume showed a linear relation with two year survival, such that even a relatively small decrease in average tidal volume over stay in the intensive care unit was independently associated with an important decrease in risk of mortality. Given the study’s findings that patients with acute lung injury often did not receive lung protective ventilation, greater efforts to implement lung protective ventilation in routine clinical practice should be undertaken to reduce patients’ long term mortality.
What is already known on this topic
Survivors of critical illness, such as those with acute lung injury, are at increased risk of mortality for years after hospital discharge
Volume limited and pressure limited (lung protective) mechanical ventilation decreases short term mortality in patients with acute lung injury, but its association with longer term survival has not been evaluated
What this study adds
Lung protective mechanical ventilation was associated with a substantial long term survival benefit for patients with acute lung injury
Lung protective ventilation is under-utilised in routine clinical practice, necessitating further knowledge translation research to improve uptake of this evidence based therapy
Cite this as: BMJ 2012;344:e2124
We thank the patients who participated in the study and the research assistants and coordinators who assisted with data collection and management for the study, including Nardos Belayneh, Rachel Bell, Kim Boucher, Abdulla Damluji, Carinda Feild, Thelma Harrington, Praveen Kondreddi, Robert Lord, Frances Magliacane, Stacey Murray, Kim Nguyen, Susanne Prassl, Arabela Sampaio, Kristin Sepulveda, Shabana Shahid, Faisal Siddiqi, Michelle Silas, and Jennifer Titus. JES is currently assistant professor, Division of Pulmonary, Allergy and Critical Care Medicine, Emory University School of Medicine, Atlanta, Georgia.
Contributors: DMN, PAM-T, JES, CRDH, CS, and PJP contributed to the conception and design of the study. DMN, EC, PAM-T, VDD, JES, CS, RGB, and PJP analysed and interpreted the data. DMN drafted the article and all other authors critically revised it for important intellectual content. All authors gave final approval of the manuscript version to be published. DMN and EC are responsible for the overall content as guarantors.
Funding: This research was supported by the National Institutes of Health (Acute Lung Injury Specialized Centers of Clinically Oriented Research grant No P050 HL 73994). DNM received support from a clinician-scientist award from the Canadian Institutes of Health Research. JES received support from the National Institutes of Health (NIGMS K-23 GMO7-1399). All researchers are independent of the funding bodies. The funding bodies had no role in the study design; in the collection, analysis, and interpretation of data; in the writing of the report; and in the decision to submit the manuscript for publication.
Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf (available on request from the corresponding author) and declare: no support from any organisation for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; and no other relationships or activities that could appear to have influenced the submitted work.
Ethical approval: This study was approved by the institutional review boards of Johns Hopkins University (No NA_00041630) and all other participating study sites.
Data sharing: No additional data available.
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.