Skip to main content

Multi-scale symbolic entropy analysis provides prognostic prediction in patients receiving extracorporeal life support



Extracorporeal life support (ECLS) can temporarily support cardiopulmonary function, and is occasionally used in resuscitation. Multi-scale entropy (MSE) derived from heart rate variability (HRV) is a powerful tool in outcome prediction of patients with cardiovascular diseases. Multi-scale symbolic entropy analysis (MSsE), a new method derived from MSE, mitigates the effect of arrhythmia on analysis. The objective is to evaluate the prognostic value of MSsE in patients receiving ECLS. The primary outcome is death or urgent transplantation during the index admission.


Fifty-seven patients receiving ECLS less than 24 hours and 23 control subjects were enrolled. Digital 24-hour Holter electrocardiograms were recorded and three MSsE parameters (slope 5, Area 6–20, Area 6–40) associated with the multiscale correlation and complexity of heart beat fluctuation were calculated.


Patients receiving ECLS had significantly lower value of slope 5, area 6 to 20, and area 6 to 40 than control subjects. During the follow-up period, 29 patients met primary outcome. Age, slope 5, Area 6 to 20, Area 6 to 40, acute physiology and chronic health evaluation II score, multiple organ dysfunction score (MODS), logistic organ dysfunction score (LODS), and myocardial infarction history were significantly associated with primary outcome. Slope 5 showed the greatest discriminatory power. In a net reclassification improvement model, slope 5 significantly improved the predictive power of LODS; Area 6 to 20 and Area 6 to 40 significantly improved the predictive power in MODS. In an integrated discrimination improvement model, slope 5 added significantly to the prediction power of each clinical parameter. Area 6 to 20 and Area 6 to 40 significantly improved the predictive power in sequential organ failure assessment.


MSsE provides additional prognostic information in patients receiving ECLS.


In recent years, extracorporeal life support (ECLS) has been increasingly used as a life-saving intervention for a variety of critically ill patients [1]. Thus, ECLS has been used as cardiac and/or pulmonary support in various clinical settings, such as fulminant myocarditis, bridge-to-heart transplantation, severe respiratory failure, cardiogenic shock after cardiac surgery, assistance for cardiopulmonary resuscitation (CPR), and septic shock [2]-[6].

However, the mortality associated with ECLS remains high, and not all critically ill patients will benefit from it [2],[4]. The outcome for the ECLS recipient is influenced not only by patient characteristics (disease severity, type of illness, other organ support) [4],[7],[8], but also by procedural complications related to ECLS [9]. Furthermore, ECLS is a resource-intensive procedure with high cost [10]. Issues about suitability of patients and, more importantly, when to cease ECLS, are particularly important. Therefore, it is important to identify patients who are likely or unlikely to benefit from use of this high-risk, high-cost treatment. Several commonly used scoring systems in the ICU, such as the acute physiology and chronic health evaluation score (APACHE), multiple organ dysfunction score (MODS), sequential organ failure assessment (SOFA), and their successors, have been evaluated in ECLS recipients [8],[11]. However, the predictive power among scoring systems varies among studies [8],[11]-[13].

Analysis of the variation of heart rate dynamics, also known as heart rate variability (HRV), is commonly used to assess autonomic function in human studies [14],[15] due to its simplicity, noninvasive character, and low cost. Its application to risk stratification of patients with cardiovascular disease has been documented to be independent of conventional clinical parameters [16],[17]. In recent years newer methods of calculating and expressing HRV based on nonlinear and non stationary signal modeling have been developed and successfully applied [18]-[21]. Compared to traditional linear HRV parameters, the nonlinear metrics showed better prediction power for cardiovascular events in several studies [22],[23]. One nonlinear method, multiscale entropy (MSE) analysis, was developed to quantify heterogeneous complexity. MSE extends the traditional entropy algorithm to quantify the information richness over multiple time scales that operate in physiological systems [18]-[20]. In a previous study, MSE provided the best prognostic prediction in patients with heart failure [23]. The utility of MSE extends beyond patients with cardiovascular disease. Thus, MSE also predicts the outcome of patients with severe trauma requiring ICU admission across the diverse spectrum of traumatic injury [24].

Unpredictable yet frequent ectopic beats are common in critical illness. These ectopic beats introduce large artifacts in the calculation of MSE [19]. The conventional approach to using MSE and other measures of HRV is to visually scan, identify and reject those beats, using interpolation procedures to smooth the time series of interbeat intervals. As critically ill patients so frequently have such arrhythmias, visual or computational pre-processing to reject and smooth not only raises the calculation complexity but also can introduce spurious trends or fluctuations into the signals that could limit clinical usefulness. In this paper we introduce a new method, termed multiscale symbolic entropy analysis (MSsE), which is derived from MSE and mitigates the effect of arrhythmia on HRV analysis.

In the current study, we hypothesized that MSsE could yield a prognostic marker in patients receiving ECLS. The aims of this study were 1) to assess the prognostic significance of parameters derived from MSsE; 2) to compare MSsE parameters to conventional clinical parameters; and 3) to evaluate the effect of combining MSsE parameters with conventional clinical parameters.

Material and methods

Setting and population

This prospective study was conducted between March 2008 and March 2010 at the National Taiwan University Hospital, which is an ECLS referral center and performs approximately 85 extracorporeal life-support procedures a year [2],[6]. Patients were eligible for the present study if they were 18 years or older and had received ECLS for circulatory or respiratory failure that required mechanical support [3]. The decision to use ECLS was made by experienced intensive care specialists or cardiac surgeons. The subjects were enrolled within the first 24 hours after ECLS implantation. Demographic, clinical features and outcomes of patients were recorded. We also recorded the usage of an intra-aortic balloon pump (IABP) or CPR during hospitalization. Catecholamine dose was evaluated by the inotrope equivalent method, that is, calculated as:

mg / kg / min = dopamine + dobutamine + 100 * epinephrine + 100 * norepinephrine + 100 * isoproterenol + 15 * milrinone )


The APACHE II, SOFA, MODS, logistic organ dysfunction sore (LODS) at ICU admission were calculated according to previous publications [11],[26]-[29]. The primary endpoint was death or urgent cardiac transplantation during the index admission. We followed the patients until discharge or death from the index admission. Healthy volunteers with normal cardiac function evaluated by echocardiography and without cardiovascular disease were enrolled as the control group. The institutional review board of National Taiwan University Hospital approved the present study and informed consent was given by each patient’s family member in the ECLS group due to the unconsciousness of patients and each subject in the control group. Informed consent in the ECLS group was obtained from patients’ family members in the order of spouse, sons or daughters, parents, grandchildren, grandparents, siblings, aunts, uncles, nieces, and nephews.

Holter recording

All enrolled subjects were placed on standard ambulatory (Holter) electrocardiogram (ECG) recorders for 24-hour recording. The ECG signals were sampled at 250 Hz and stored on an SD card for subsequent offline analysis.

MSsE calculation

The MSE has proven to be an effective tool in exploring the characteristics of heart rate dynamics and can predict important clinical outcomes [18],[23],[30],[31]. The original MSE comprises two steps: 1) coarse-graining the signals using different time scales; 2) quantifying the degree of irregularity in each coarse-grained time series using sample entropy. However, the major challenge of applying MSE to severe illness patients is the frequent ectopic beats that can reduce the reliability of the MSE results [19]. To address the issue of ectopic beats, we have adopted the intuitive idea that the coarse grained time sequence will be better reconstructed with the median value rather than mean value over non-overlapping windows. In addition, a further step proposed to attenuate the influence from the ectopic beats is to use the sign of the coarse-grained signal to measure its entropy. Dynamics of heart rate after coarse-graining in different scales were therefore denoted as the increase (+) or decrease (-) sign in forming coarse-grained binary sequences. Intuitively, by focusing only on the direction while ignoring the amplitude of the change, the effect of these outliers could be both localized and diminished.

In order to compute the entropy of a binary sequence, the sequence should be divided into subsets consisting of L consecutive binary bits. The detailed procedure for calculating the entropy of the coarse-grained binary sequences can be found in the caption of Figure 1. We plotted the entropy as a function of scale, which provided the quantitation of structure richness embedded in heartbeat fluctuations over different timescales. This symbolic dynamics approach, now termed multiscale symbolic entropy analysis (MSsE), can provide similar information as the original MSE in quantifying the complexity of the signal while being far less sensitive to the unpredictable yet frequent appearance of ectopic beats. According to previous results that applied MSE to the heartbeat recordings of healthy young subjects [18],[19], the sample entropies of heartbeat fluctuations in individuals resembles those of Brownian motions over short timescales (<5 heartbeats) while they behave more like those of pink (1/f) noise over longer timescales. Such a crossover phenomenon is believed to be related to different physiological processes influencing cardiac dynamics over different timescale regions, that is, over short timescales respirations dominantly entrain heartbeat fluctuations (respiratory sinus arrhythmia) via its influence on the parasympathetic nervous system; and over longer timescales (>5 heartbeats), heartbeat fluctuations are influenced by many different physiological functions such as baroreflex and circadian rhythm [32] that lead to complex fluctuation patterns similar to 1/f noise. These findings suggest that we should select Slope5 (slope 5), Area6–20 (area 6 to 20), and Area6–40 (area 6 to 40) as the indices with physiological meaning to profile the MSsE curve (see Figure 2).

Figure 1

The first operation of MSsE is coarse graining. We divide the time series {x1, x2 … xM} into non-overlapping boxes of size N (scale "N"). The median of each local box is taken, thus producing a new time series y 1 N , y 2 N y k N . The second stage in forming the sign time series is to transform the coarse-grained series into yet another new series by taking the directions in its change. We measure the change against a threshold value and acquire the sign series b 1 N , b 2 N b k N , where b i N is either +1 or -1, depending on whether the corresponding y i N is increasing or decreasing. To quantify the complexity of the sign sequence, we sort all sequences into categories of sub-sets consist of L consecutive binary bits (L-bit; L = 8 in this study). The probability distribution of all patterns of sub-sets is recorded. To avoid over-counting similar patterns, the data sequence of total length L should be divided into multiple m-dimensional vectors; each consists of m consecutive bits {(b1, b2, … bm); (b2, b3, … bm +1); …}. The conditional probability is determined numerically by the ratio of number of each paired vectors which are of exactly same binary codes for dimension "m+1" to the number for the identical vectors of dimension "m". By identifying the patterns of the same conditional probability, it allows us to rank the m-bit patterns according to the information they imply (large rank number means lower conditional probability). The expectation value of the rank conceptually indicates the degree of uncertainty.

Figure 2

Quantification of multiscale symbolic entropy (MSsE): summation of the entropy over different scales can quantify the complexity over certain timescales. However, typical profile of MSsE in extracorporeal life support patients showed a crossover phenomenon around scale 5. Three parameters of the MSsE were assessed: (1) the linear-fitted slope between scales 1 to 5; (2) complexity between intermediate scales (Area 6-20); and (3) the overall complexity (Area 6-40).

We also accessed slope 5 from a shorter time interval of the ECG recording: slope 5 for the first hour and the first 2 hours of ECLS implantation.

Statistical analysis

Categorical variables were analyzed using the chi-square test or Fisher’s exact test. Continuous variables were represented as mean value ± SD and the normality of those variables was evaluated by using Kolmogorov-Smirnov test. Student’s t-test was applied to the between-group comparison. The maximal hazards ratio and independent correlation of variables with event status (mortality or urgent heart transplantation) was determined by Cox regression analysis. Harrell’s C-statistics (the probability of concordance for any two randomly chosen subjects) were calculated as a measure of a model’s ability to discriminate between patients meeting and not meeting a primary outcome [33]-[36]. More specifically:

C=P Z i > Z j | D i = 1 , D j = 0 ,

where Zi, Zj are model-based risks (that is, linear predictors) and Di, Dj are event indicators for two subjects (1 = patients meeting the primary outcome; 0 = patients not meeting the primary outcome). The receiver-operator characteristic (ROC) curve was determined by the logistic regression model. We performed C-statistics to describe discrimination of the baseline model by clinical severity scores and the model that included selected MSsE parameters [34],[36],[37].

Net reclassification improvement (NRI) and integrated discrimination improvement (IDI) modeling were performed to assess the improvement of the prediction using two different logistic regression models [34], with 0.2 and 0.4 used as cutoff points. All statistical analyses were performed using R software [38], version 2.15.2. Statistical significance was set at P <0.05.


A total of 57 patients (44 male) who were receiving ECLS, and 23 control subjects were enrolled in this study. Of the patients 51 received veno-arterial ECLS and 6 received veno-venous ECLS. The demographic data and data on the MSsE parameters are shown in Table 1. The values of MSsE are also shown in Figure 3. Patients receiving ECLS had significantly lower value of slope 5, area 6 to 20, and area 6 to 40.

Table 1 Clinical and MSsE parameters of all subjects
Figure 3

Multiscale symbolic entropy (MSsE) analysis of heart rate dynamics from control subjects, patients surviving after extracorporeal life support (ECLS), and non-surviving ECLS patients. The values are represented as mean ± standard error. The MSsE profile of ECLS patients who survived differed from that of non-survivors.

Among patients receiving ECLS, 35 patients received ECLS for cardiac indications; including 2 patients with acute myocarditis, 9 patients with acute myocardial infarction, 7 patients with intractable heart failure, 6 patients with poor hemodynamic status after cardiac surgery, 8 patients with poor hemodynamic status due to frequent or intractable ventricular arrhythmia, 2 patients with out-of-hospital cardiac arrest and one patient due to massive pulmonary embolism. A total of 22 patients received ECLS for non-cardiac indications; including 2 patients with thyroid storm, 9 patients with acute respiratory distress syndrome, 7 patients with septic shock, one patient with carbon monoxide intoxication, one patient with electrical injury, and two patients with shock due to massive bleeding.

During the index hospitalization, 26 patients died and 3 patients received urgent cardiac transplantation. The other 28 patients were successfully weaned from ECLS and discharged from the hospital alive. In an analysis of baseline characteristics (Table 2), survivors without urgent cardiac transplantation were younger, had longer period of ECLS implantation, lower APACHE score, lower LODS, and lower MODS. With respect to the MSsE parameters (Table 3), survivors without urgent cardiac transplantation had a significantly higher value of slope 5, area 6 to 20, and area 6 to 40.

Table 2 Clinical and laboratory parameters of ECLS recipients
Table 3 MSsE parameters of ECLS recipients

In Cox regression survival analysis (Table 4), the significant factors were age, prior myocardial infarct (MI) history, APACHE score, LODS, MODS, Slope 5, area 6 to 20, and area 6 to 40. Slope 5 had the largest concordance value among MSsE parameters; the LODS had the largest concordance value among clinical severity scores (Table 5). After combining MSsE parameters with measures of clinical severity, concordance values increased. Among the MSsE parameters, slope 5 added the largest improvement.

Table 4 Cox regression analysis for prediction of the primary endpoint using a single variable in ECLS recipients
Table 5 Concordance, AUC, NRI, and IDI model among clinical and MSsE parameters

In ROC analyses (Table 5), slope 5 had the largest area under the curve (AUC) among MSsE parameters; the LODS had the largest AUC among clinical severity scores. When combining MSsE parameters with measures of clinical severity, AUC improved outcome prediction. Among the MSsE parameters, slope 5 added the largest improvement. For example (see Figure 4), the AUC for SOFA was 0.637. After adding slope 5, the AUC for the new model improved to 0.749.

Figure 4

Receiver-operator characteristic curves for models with (solid line) and without slope 5 (broken line). The area under the curve (AUC) for sequential organ failure assessment SOFA was 0.637. After adding slope 5, the AUC of the new model improved to 0.749.

Although the increase in AUC after adding MSsE parameters on conventional risk scores was not statistically significant, that does not answer the question of whether the addition of MSsE metrics improves the predictive power of conventional risk scores. To answer this question, we considered whether reclassification with the addition of MSsE improved predictions for the individual patients. We used two strategies: NRI and IDI. In general, NRI offers a simple intuitive way of quantifying improvement offered by new markers. As described by Pencina and colleagues, `The NRI focuses on reclassification tables constructed separately for participants with and without events, and quantifies the correct movement in categories—“upwards” for events and “downwards” for non-events’ [37]. In the NRI model, slope 5 significantly improved the predictive power of LODS; area 6 to 20 and area 6 to 40 significantly improved the predictive power in MODS. A further improvement on NRI examines all potential cutoff values in the ROC analysis and examines the integral of sensitivity verses the integral of one minus specificity. This is a continuous version of NRI. In IDI model, slope 5 added significantly to the prediction power of each clinical parameter. Area 6 to 20 and area 6 to 40 significantly improved the predictive power in SOFA. The comparison among slope 5 from three time intervals (1 hour, 2 hours, 24 hours) of the ECG recording is shown in Additional file 1.


The major findings of this study were: 1) Among the various clinical and MSsE parameters, slope 5 showed the greatest power in prognostic prediction. 2) MSE parameters (especially slope 5) added the prognostic prediction power in clinical parameters.

The increase in number and severity of critically ill patients worldwide demands that critical care professionals make prudent and objective decisions about the allocation of costly and risky treatments. Such treatments include ECLS, which currently affords survival to only about half the adult patients in which it is initiated. However, due to its high cost [10], issues about the timing to cease ECLS are particularly important. Therefore, it is important to identify patients who are likely to benefit from ECLS. Physiology-based risk-classification tools are therefore needed to support decisions for or against continuous ECLS usage. Although HRV metrics have proved useful in predicting outcomes in some populations, frequent arrhythmias have limited their use in patients with acute and severe critical illness. In this paper, we report the development and application of a new family of HRV metrics based on symbolic entropies that extend the utility of existing nonlinear measures and demonstrate that three members of this new family (slope 5, area 6–20, area 6–40) have the potential to improve conventional clinical risk-stratification tools in determining likelihood of ECLS success.

In our study, the primary outcome rate is quite high (45.6% for mortality, 50.9% for mortality + urgent transplantation). The reasons are probably due to the high disease severity in our participants. Sixteen patients received CPR before ECLS implantation. In our previous experience, the 30-day survival rate of ECLS assisted CPR is 33.9% [2]. Although the mortality rate is still high, however, the outcome of ECLS assisted CPR is much better than conventional CPR [2]. Due to the high mortality and cost of ECLS therapy, effective risk-classification tools are needed to support decisions for or against ECLS initiation and continuation. In a study in Norway, the median estimated cost for the ECMO procedure was 62,545 USD, and the mean estimated total hospital costs is 191,436 USD [10]. In that study, mean duration of an ECLS procedure is 9.5 days (range: 4–23 days). In our study, the average duration of ECLS procedures is 11.69 days in patients met primary outcome. Early risk-classification tools may offer clinician and patient additional insight into their shared decision to use complicated and costly therapies.

In our study, new HRV metrics based on symbolic entropies appeared to enhance risk-classification. These MSsE parameters had good performance in outcome prediction when using alone or combined with clinical parameters. Among the MSsE parameters, slope 5 added the largest improvement. In each clinical parameter, after adding slope 5, the increase of AUC in new model is around 0.06-0.07. The increase exists in all three clinical parameters. There are advantageous features of this tool. First, the cost of measuring MSsE parameters is low and the method can be integrated into ICU data systems. Second, it is measured at bedside, and therefore does not require potentially dangerous transport to an imaging or procedural location. Third, neither radiation nor pharmaceuticals nor additional devices are required.

There are other studies using MSE analysis for acute illness. In trauma patients admitted to ICU, reduced MSE is significantly associated with increasing mortality, and is independent to clinical parameters [39]. In another study, decreased approximate entropy is associated with mortality in trauma patients independent of Glasgow coma score or injury severity score [40]. In a study of combat casualties in an emergency department in Iraq, the MSE index was significantly decreased in patients who required lifesaving intervention compared to those who did not [41]. This evidence shows the potentials of physiology-based risk-classification tools in patients with acute illness. Furthermore, our study demonstrates the additive effects of combination of physiology-based risk-classification tools and clinical parameters.

Interconnectedness of physiological mechanisms is a crucial feature to the output of a network with feedback interactions [42] and demonstrates as multiple time scale correlation for the output fluctuations [43],[44]. Multiple time scale correlation traditionally can be assessed by fractal analysis such as detrended fluctuation analysis (DFA) [45]-[47] i.e., the random time series (no correlation) will exhibit the ~0.5 scaling exponent in the fluctuations vs. time scales for DFA analysis. On the other hand, MSE offer an alternative to verify the existence of multiscale correlation by the slope of the scale function of entropy. Mathematically, a negative slope of MSE analysis indicates the trivial correlation (i.e., simple oscillation, a repetitive pattern of an increase follow by a decrease) or random fluctuations [19]. Therefore, it can be a sign of the loss of feedback regulation. The present study reveals that, in addition to the overall reduction of entropy at different scales (i.e., area 6–20, area 6–40), the patients with death or urgent transplantation show very negative slope (i.e. slope 5) in MSsE analysis. This finding provides preliminary evidence that loss of feedback control could be important features of severe outcome.

In comparison of the outcome predictive power among slope 5 collected from three time interval (from first one hour, from first two hour, and from 24 hours) of ECG recording, generally, the slope 5 from 24 hours have better performance than slope 5 for shorter time ECG recording. However, the slope 5 form one or two hour recording (especially slope 5 form first two hours) still have the ability to predict outcome and also add the performance of clinical parameters. Therefore, parameters form short time ECG recording also have potential for outcome prediction in ECLS patients. Further larger study is need to compare to outcome prediction power of MSsE parameters form short and long time ECG recording in the future.

Our study has several important limitations. It is a single-center study that focuses on the relatively closed population of Taiwan. The study is based on a 24 hour data collection period, an interval that is both long and arbitrary. Our study lumps patients having primarily cardiac disease with patients who have a spectrum of other diseases. This “lumping” convolves and to some extent confuses expression of primary heart pathology with effects of extra-cardiac pathologies upon the heart. The particular entropy metrics that proved useful, while pre-selected for this study, are still somewhat arbitrary. For all of these reasons, our findings must be considered provisional and this study requires replication in another venue to test the reproducibility and consistency of our findings.

These limitations notwithstanding, our study suggests that the addition of metrics derived from non-linear symbolic analysis of heart rate variability improves the ability of conventional risk-stratification tools to create accurate predictions about the outcomes of individual patients within a complicated and very sick population. This is important for two reasons. First, it adds confidence that the “arrhythmia problem in HRV analysis”--the appearance of confounding ectopic beats in critically ill patients--can be addressed objectively and automatically, thus making HRV metrics generally more useful in critically ill patients. Second, it suggests that a more general problem of situation awareness--specifically the ability to make reliable predictions about the outcome of different treatment options--may be improved with the addition of multiscale metrics of physiologic responsiveness to more traditional static indices of physiologic derangement.


Patients receiving ECLS remain high mortality rate. MSsE offers additional insight into the prognosis of patients considered for ECLS.

Key messages

Extracorporeal life support (ECLS) is a resource-intensive procedure with high cost. Issues about suitability of patients and, more importantly, when to cease ECLS are particularly important.

Multi-scale entropy (MSE) derived from heart rate variability (HRV) is a powerful tool in outcome prediction of patients with cardiovascular diseases. Multi-scale symbolic entropy analysis (MSsE), a new method derived from MSE, mitigates the effect of arrhythmia on analysis.

In this study, we found several MSsE parameters provide additional prognostic information in patients receiving ECLS. Slope 5 is the best parameter.

Authors’ contributions

YHL participated in the analysis of clinical data and drafted the manuscript. HCH participated in the collection of clinical data. YCC participated in the analysis of Holter data. LC participated in the analysis of Holter data. MTL participated in the design of the study and drafted the manuscript. LYDL participated in the analysis of clinical data and performed the statistical analysis. YSC participated in the analysis of clinical data. WJK participated in the study design and clinical data collection. YLH participated in the study design and drafted the manuscript. MFC participated in the study design. CKP critically revised the manuscript. TGB critically revised the manuscript. All authors read and approved the final manuscript.

Additional file



acute physiology and chronic health evaluation


area under curve


cardiopulmonary resuscitation


detrended fluctuation analysis




extracorporeal life support


heart rate variability


integrated discrimination improvement


logistic organ dysfunction score


multiple organ dysfunction score


multiscale entropy


multi-scale symbolic entropy


net reclassification improvement


receiver-operator characteristic


sequential organ failure assessment


  1. 1.

    Ainslie PN, Lewis NC, Eves ND: Letter by Ainslie et al. regarding article, “Endogenous circadian rhythm in vasovagal response to head-up tilt”. Circulation. 2011, 124: e358-10.1161/CIRCULATIONAHA.111.032227. author reply e359

    Article  Google Scholar 

  2. 2.

    Chen YS, Lin JW, Yu HY, Ko WJ, Jerng JS, Chang WT, Chen WJ, Huang SC, Chi NH, Wang CH, Chen LC, Tsai PR, Wang SS, Hwang JJ, Lin FY: Cardiopulmonary resuscitation with assisted extracorporeal life-support versus conventional cardiopulmonary resuscitation in adults with in-hospital cardiac arrest: an observational study and propensity analysis. Lancet. 2008, 372: 554-561. 10.1016/S0140-6736(08)60958-7.

    Article  Google Scholar 

  3. 3.

    Sun HY, Ko WJ, Tsai PR, Sun CC, Chang YY, Lee CW, Chen YC: Infections occurring during extracorporeal membrane oxygenation use in adult patients. J Thorac Cardiovasc Surg. 2010, 140: 1125-1132. 10.1016/j.jtcvs.2010.07.017.

    Article  Google Scholar 

  4. 4.

    Pappalardo F, Pieri M, Greco T, Patroniti N, Pesenti A, Arcadipane A, Ranieri VM, Gattinoni L, Landoni G, Holzgraefe B, Beutel G, Zangrillo A: Predicting mortality risk in patients undergoing venovenous ECMO for ARDS due to influenza A (H1N1) pneumonia: the ECMOnet score. Intensive Care Med. 2013, 39: 275-281. 10.1007/s00134-012-2747-1.

    Article  Google Scholar 

  5. 5.

    Hsu KH, Chi NH, Yu HY, Wang CH, Huang SC, Wang SS, Ko WJ, Chen YS: Extracorporeal membranous oxygenation support for acute fulminant myocarditis: analysis of a single center’s experience. Eur J Cardiothorac Surg. 2011, 40: 682-688.

    Google Scholar 

  6. 6.

    Huang CT, Tsai YJ, Tsai PR, Ko WJ: Extracorporeal membrane oxygenation resuscitation in adult patients with refractory septic shock. J Thorac Cardiovasc Surg. 2013, 146: 1041-1046. 10.1016/j.jtcvs.2012.08.022.

    Article  Google Scholar 

  7. 7.

    Wu MY, Lin PJ, Lee MY, Tsai FC, Chu JJ, Chang YS, Haung YK, Liu KS: Using extracorporeal life support to resuscitate adult postcardiotomy cardiogenic shock: treatment strategies and predictors of short-term and midterm survival. Resuscitation. 2010, 81: 1111-1116. 10.1016/j.resuscitation.2010.04.031.

    Article  Google Scholar 

  8. 8.

    Wu VC, Tsai HB, Yeh YC, Huang TM, Lin YF, Chou NK, Chen YS, Han YY, Chou A, Lin YH, Wu MS, Lin SL, Chen YM, Tsai PR, Ko WJ, Wu KD: Patients supported by extracorporeal membrane oxygenation and acute dialysis: acute physiology and chronic health evaluation score in predicting hospital mortality. Artif Organs. 2010, 34: 828-835. 10.1111/j.1525-1594.2009.00920.x.

    Article  Google Scholar 

  9. 9.

    Huang PM, Ko WJ, Tsai PR, Kuo SW, Hsu HH, Chen JS, Lee JM, Lee YC: Aggressive management of massive hemothorax in patients on extracorporeal membrane oxygenation. Asian J Surg. 2012, 35: 16-22. 10.1016/j.asjsur.2012.04.003.

    Article  Google Scholar 

  10. 10.

    Mishra V, Svennevig JL, Bugge JF, Andresen S, Mathisen A, Karlsen H, Khushi I, Hagen TP: Cost of extracorporeal membrane oxygenation: evidence from the Rikshospitalet University Hospital, Oslo, Norway. Eur J Cardiothorac Surg. 2010, 37: 339-342. 10.1016/j.ejcts.2009.10.008.

    Article  Google Scholar 

  11. 11.

    Tsai CW, Lin YF, Wu VC, Chu TS, Chen YM, Hu FC, Wu KD, Ko WJ, Group NS: SAPS 3 at dialysis commencement is predictive of hospital mortality in patients supported by extracorporeal membrane oxygenation and acute dialysis. Eur J Cardiothorac Surg. 2008, 34: 1158-1164. 10.1016/j.ejcts.2008.07.025.

    Article  Google Scholar 

  12. 12.

    Sun CK, Lin YC, Yuen CM, Chua S, Chang LT, Sheu JJ, Lee FY, Fu M, Leu S, Yip HK: Enhanced protection against pulmonary hypertension with sildenafil and endothelial progenitor cell in rats. Int J Cardiol. 2012, 162: 45-58. 10.1016/j.ijcard.2011.05.002.

    Article  Google Scholar 

  13. 13.

    Chang WW, Tsai FC, Tsai TY, Chang CH, Jenq CC, Chang MY, Tian YC, Hung CC, Fang JT, Yang CW, Chen YC: Predictors of mortality in patients successfully weaned from extracorporeal membrane oxygenation. PLoS One. 2012, 7: e42687-10.1371/journal.pone.0042687.

    CAS  Article  Google Scholar 

  14. 14.

    Heart rate variability: standards of measurement, physiological interpretation and clinical use. Circulation. 1996, 93: 1043-1065. 10.1161/01.CIR.93.5.1043.

  15. 15.

    Umetani K, Singer DH, McCraty R, Atkinson M: Twenty-four hour time domain heart rate variability and heart rate: relations to age and gender over nine decades. J Am Coll Cardiol. 1998, 31: 593-601. 10.1016/S0735-1097(97)00554-8.

    CAS  Article  Google Scholar 

  16. 16.

    La Rovere MT, Pinna GD, Maestri R, Mortara A, Capomolla S, Febo O, Ferrari R, Franchini M, Gnemmi M, Opasich C, Riccardi PG, Traversi E, Cobelli F: Short-term heart rate variability strongly predicts sudden cardiac death in chronic heart failure patients. Circulation. 2003, 107: 565-570. 10.1161/01.CIR.0000047275.25795.17.

    Article  Google Scholar 

  17. 17.

    Tsuji H, Larson MG, Venditti FJ, Manders ES, Evans JC, Feldman CL, Levy D: Impact of reduced heart rate variability on risk for cardiac events. The Framingham Heart Study. Circulation. 1996, 94: 2850-2855. 10.1161/01.CIR.94.11.2850.

    CAS  Article  Google Scholar 

  18. 18.

    Costa M, Goldberger AL, Peng CK: Multiscale entropy analysis of complex physiologic time series. Phys Rev Lett. 2002, 89: 068102-10.1103/PhysRevLett.89.068102.

    Article  Google Scholar 

  19. 19.

    Costa M, Goldberger AL, Peng CK: Multiscale entropy analysis of biological signals. Phys Rev E Stat Nonlin Soft Matter Phys. 2005, 71: 021906-10.1103/PhysRevE.71.021906.

    Article  Google Scholar 

  20. 20.

    Peng CK, Costa M, Goldberger AL: Adaptive Data Analysis of Complex Fluctuations in Physiologic Time Series. Adv Adapt Data Anal. 2009, 1: 61-70. 10.1142/S1793536909000035.

    Article  Google Scholar 

  21. 21.

    Seely AJ, Macklem PT: Complex systems and the technology of variability analysis. Crit Care. 2004, 8: R367-R384. 10.1186/cc2948.

    Article  Google Scholar 

  22. 22.

    Makikallio TH, Huikuri HV, Hintze U, Videbaek J, Mitrani RD, Castellanos A, Myerburg RJ, Moller M, Group DS: Fractal analysis and time- and frequency-domain measures of heart rate variability as predictors of mortality in patients with heart failure. Am J Cardiol. 2001, 87: 178-182. 10.1016/S0002-9149(00)01312-6.

    CAS  Article  Google Scholar 

  23. 23.

    Ho YL, Lin C, Lin YH, Lo MT: The prognostic value of non-linear analysis of heart rate variability in patients with congestive heart failure-a pilot study of multiscale entropy. PLoS One. 2011, 6: e18699-10.1371/journal.pone.0018699.

    CAS  Article  Google Scholar 

  24. 24.

    Riordan WP, Norris PR, Jenkins JM, Morris JA: Early loss of heart rate complexity predicts mortality regardless of mechanism, anatomic location, or severity of injury in 2178 trauma patients. J Surg Res. 2009, 156: 283-289. 10.1016/j.jss.2009.03.086.

    Article  Google Scholar 

  25. 25.

    Chen YS, Ko WJ, Lin FY, Huang SC, Chou TF, Chou NK, Hsu RB, Wang SS, Chu SH: Preliminary result of an algorithm to select proper ventricular assist devices for high-risk patients with extracorporeal membrane oxygenation support. J Heart Lung Transplant. 2001, 20: 850-857. 10.1016/S1053-2498(01)00267-4.

    CAS  Article  Google Scholar 

  26. 26.

    Le Gall JR, Klar J, Lemeshow S, Saulnier F, Alberti C, Artigas A, Teres D: The Logistic Organ Dysfunction system. A new way to assess organ dysfunction in the intensive care unit. ICU Scoring Group. JAMA. 1996, 276: 802-810. 10.1001/jama.1996.03540100046027.

    CAS  Article  Google Scholar 

  27. 27.

    Knaus WA, Draper EA, Wagner DP, Zimmerman JE: APACHE II: a severity of disease classification system. Crit Care Med. 1985, 13: 818-829. 10.1097/00003246-198510000-00009.

    CAS  Article  Google Scholar 

  28. 28.

    Marshall JC, Cook DJ, Christou NV, Bernard GR, Sprung CL, Sibbald WJ: Multiple organ dysfunction score: a reliable descriptor of a complex clinical outcome. Crit Care Med. 1995, 23: 1638-1652. 10.1097/00003246-199510000-00007.

    CAS  Article  Google Scholar 

  29. 29.

    Vincent JL, de Mendonca A, Cantraine F, Moreno R, Takala J, Suter PM, Sprung CL, Colardyn F, Blecher S: Use of the SOFA score to assess the incidence of organ dysfunction/failure in intensive care units: results of a multicenter, prospective study. Working group on “sepsis-related problems” of the European Society of Intensive Care Medicine. Crit Care Med. 1998, 26: 1793-1800. 10.1097/00003246-199811000-00016.

    CAS  Article  Google Scholar 

  30. 30.

    Chialvo DR: Physiology: unhealthy surprises. Nature. 2002, 419: 263-10.1038/419263a.

    CAS  Article  Google Scholar 

  31. 31.

    Yuan HK, Lin C, Tsai PH, Chang FC, Lin KP, Hu HH, Su MC, Lo MT: Acute increase of complexity in the neurocardiovascular dynamics following carotid stenting. Acta Neurol Scand. 2011, 123: 187-192. 10.1111/j.1600-0404.2010.01384.x.

    Article  Google Scholar 

  32. 32.

    Hu K, Scheer FA, Laker M, Smales C, Shea SA: Endogenous circadian rhythm in vasovagal response to head-up tilt. Circulation. 2011, 123: 961-970. 10.1161/CIRCULATIONAHA.110.943019.

    Article  Google Scholar 

  33. 33.

    Obuchowski NA: ROC analysis. AJR Am J Roentgenol. 2005, 184: 364-372. 10.2214/ajr.184.2.01840364.

    Article  Google Scholar 

  34. 34.

    Steyerberg EW, Vickers AJ, Cook NR, Gerds T, Gonen M, Obuchowski N, Pencina MJ, Kattan MW: Assessing the performance of prediction models: a framework for traditional and novel measures. Epidemiology. 2010, 21: 128-138. 10.1097/EDE.0b013e3181c30fb2.

    Article  Google Scholar 

  35. 35.

    Harrell FE, Lee KL, Mark DB: Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med. 1996, 15: 361-387. 10.1002/(SICI)1097-0258(19960229)15:4<361::AID-SIM168>3.0.CO;2-4.

    Article  Google Scholar 

  36. 36.

    Pencina MJ, D’Agostino RB: Overall C as a measure of discrimination in survival analysis: model specific population value and confidence interval estimation. Stat Med. 2004, 23: 2109-2123. 10.1002/sim.1802.

    Article  Google Scholar 

  37. 37.

    Pencina MJ, D’Agostino RB, D’Agostino RB, Vasan RS: Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med. 2008, 27: 157-172. 10.1002/sim.2929. discussion 207–112

    Article  Google Scholar 

  38. 38.

    The R Project for Statistical Computing. In [], []

  39. 39.

    Norris PR, Canter JA, Jenkins JM, Moore JH, Williams AE, Morris JA: Personalized medicine: genetic variation and loss of physiologic complexity are associated with mortality in 644 trauma patients. Ann Surg. 2009, 250: 524-530.

    Google Scholar 

  40. 40.

    Batchinsky AI, Cancio LC, Salinas J, Kuusela T, Cooke WH, Wang JJ, Boehme M, Convertino VA, Holcomb JB: Prehospital loss of R-to-R interval complexity is associated with mortality in trauma patients. J Trauma. 2007, 63: 512-518. 10.1097/TA.0b013e318142d2f0.

    Article  Google Scholar 

  41. 41.

    Cancio LC, Batchinsky AI, Baker WL, Necsoiu C, Salinas J, Goldberger AL, Costa MD: Combat casualties undergoing lifesaving interventions have decreased heart rate complexity at multiple time scales. J Crit Care. 2013, 28: 1093-1098. 10.1016/j.jcrc.2013.08.022.

    Article  Google Scholar 

  42. 42.

    Buchman TG: The community of the self. Nature. 2002, 420: 246-251. 10.1038/nature01260.

    CAS  Article  Google Scholar 

  43. 43.

    Goldberger AL, Amaral LA, Hausdorff JM, Ivanov P, Peng CK, Stanley HE: Fractal dynamics in physiology: alterations with disease and aging. Proc Natl Acad Sci USA. 2002, 99: 2466-2472. 10.1073/pnas.012579499.

    Article  Google Scholar 

  44. 44.

    Goldberger AL: Non-linear dynamics for clinicians: chaos theory, fractals, and complexity at the bedside. Lancet. 1996, 347: 1312-1314. 10.1016/S0140-6736(96)90948-4.

    CAS  Article  Google Scholar 

  45. 45.

    Peng CK, Buldyrev SV, Hausdorff JM, Havlin S, Mietus JE, Simons M, Stanley HE, Goldberger AL: Non-equilibrium dynamics as an indispensable characteristic of a healthy biological system. Integr Physiol Behav Sci. 1994, 29: 283-293. 10.1007/BF02691332.

    CAS  Article  Google Scholar 

  46. 46.

    Peng CK, Havlin S, Hausdorff JM, Mietus JE, Stanley HE, Goldberger AL: Fractal mechanisms and heart rate dynamics. Long-range correlations and their breakdown with disease. J Electrocardiol. 1995, 28: 59-65. 10.1016/S0022-0736(95)80017-4.

    Article  Google Scholar 

  47. 47.

    Hu K, Van Someren EJ, Shea SA, Scheer FA: Reduction of scale invariance of activity fluctuations with aging and Alzheimer’s disease: Involvement of the circadian pacemaker. Proc Natl Acad Sci USA. 2009, 106: 2490-2494. 10.1073/pnas.0806087106.

    CAS  Article  Google Scholar 

Download references


The authors would like to thank the staff of General Clinical Trial and Research in National Taiwan University Hospital for their great support. This study was supported by National Taiwan University Hospital (grants NTU 102-S2096), National Taiwan University (National Taiwan University Cutting-Edge Steering Research 10R71608-1), Ministry of Science and Technology (MOST 103-2220-E-002 -011, MOST 103-2221-E-008-006-MY3), Ministry of Science and Technology support for the Center for Dynamical Biomarkers and Translational Medicine, National Central University, Taiwan (NSC 102-2911-I-008-001), and the joint foundations (Grant No 103CGH-NCU-A1, VGHUST103-G1-3-3). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information



Corresponding authors

Correspondence to Men-Tzung Lo or Yi-Lwun Ho.

Additional information

Competing interests

All authors declare that they have no competing interests.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Lin, Y., Huang, H., Chang, Y. et al. Multi-scale symbolic entropy analysis provides prognostic prediction in patients receiving extracorporeal life support. Crit Care 18, 548 (2014).

Download citation


  • Heart Rate Variability
  • Sequential Organ Failure Assessment
  • Detrended Fluctuation Analysis
  • Extracorporeal Life Support
  • Integrate Discrimination Improvement