Research Article  Open Access
RongGuan Yeh, ChungWu Lin, Maysam F. Abbod, JiannShing Shieh, "TwoDimensional Matrix Algorithm Using Detrended Fluctuation Analysis to Distinguish Burkitt and Diffuse Large BCell Lymphoma", Computational and Mathematical Methods in Medicine, vol. 2012, Article ID 947191, 8 pages, 2012. https://doi.org/10.1155/2012/947191
TwoDimensional Matrix Algorithm Using Detrended Fluctuation Analysis to Distinguish Burkitt and Diffuse Large BCell Lymphoma
Abstract
A detrended fluctuation analysis (DFA) method is applied to image analysis. The 2dimensional (2D) DFA algorithms is proposed for recharacterizing images of lymph sections. Due to Burkitt lymphoma (BL) and diffuse large Bcell lymphoma (DLBCL), there is a significant different 5year survival rates after multiagent chemotherapy. Therefore, distinguishing the difference between BL and DLBCL is very important. In this study, eighteen BL images were classified as group A, which have one to five cytogenetic changes. Ten BL images were classified as group B, which have more than five cytogenetic changes. Both groups A and B BLs are aggressive lymphomas, which grow very fast and require more intensive chemotherapy. Finally, ten DLBCL images were classified as group C. The shortterm correlation exponent α1 values of DFA of groups A, B, and C were , and , respectively. It was found that α1 value of BL image was significantly lower () than DLBCL. However, there is no difference between the groups A and B BLs. Hence, it can be concluded that α1 value based on DFA statistics concept can clearly distinguish BL and DLBCL image.
1. Introduction
Natural phenomena almost are random, nonlinear, nonstationary, disordered, and uncertain systems. It is difficult to follow the traditional Newtonian rules to be completely forecast or control. Therefore, opening up fractal methods to investigate complex, rough, fragment shape, and noninteger dimension naturally objects (trees, coastlines, clouds, and mountains, etc.) is necessity. In the 1960s, the mathematician Mandelbrot had been indicated the fractal objects whose complex geometry cannot be characterized by an integral dimension [1]. This phenomenon is often expressed by spatial or timedomain statistical scaling laws and is mainly characterized by the powerlaw behavior of realworld physical systems. It represents fractal applies to objects in space or fluctuations in time that possess a form of selfsimilarity [2]. The object has selfsimilarity means; the variant of object expressed same qualitatively, irrespective reduction or magnification of the object. Selfsimilarity is one way to calculate fractal dimension. For example, one can subdivide a line segment into selfsimilar intervals, each with the same length, and each of which can be magnified by a factor of to yield the original segment [3]. Due to the fractal geometry having an approximately copy of the whole, the fractal dimension is consistent over a wide range of scales, which is known as scale invariance [4]. This property provided a useful measurement of complexity object.
DFA method was developed from a modified root mean square analysis of a random walk to exclude the local trend induced by characteristic time scales from the fluctuations of the multicomponent systems and get a longrange correlation [5–7]. It was originally a method to measure scale invariant behavior developed by Peng et al. [6] that evaluated trends of all sizes in the presence or absence of fractal correlation properties of time series data [8, 9]. This method has been applied to heart rate dynamics such as autonomic nervous system [10], congestive heart failure [8], dilated cardiomyopathy [11], ventricular fibrillation [12], and other physiological nonstationary time series systems (DNA sequences [13], neuron spiking [14, 15], human gait analysis [16], electroencephalogram (EEG) in sleep [17–20], stock returns [21], periodic trends [22], estimating dependence [23], etc.). Experience has shown that monodimensional detrended fluctuation analysis (DFA) used in the scaling analysis of fractal time series is accurate and easy to implement regardless in longterm and shortterm time scale series [24, 25]. In recent years, there are some modified DFA method researches that are proposed such as generalized the monodimensional DFA and multifractal detrended fluctuation analysis (MFDFA) to higherdimensional versions and derived multifractal detrended crosscorrelation analysis method to investigate the multifractal behaviors in the powerlaw cross correlations between two time series or higherdimensional quantities recorded [26–28]. The multifractal detrended crosscorrelation analysis based on DFA (MFXDFA) [27] is actually a multifractal generalization of the detrended crosscorrelation analysis (DCCA) [29], which has other variants such as the multifractal detrended crosscorrelation analysis based on DMA (MFXDMA) [30]. Those study results validated well for distinguishing fractal/multifractal properties of synthetic surfaces (including fractional Brownian and multifractal surfaces), one/twodimensional cross correlation of two financial time series, and linear/nonlinear correlation analysis of traffic time series (to find the cross correlation of traffic flow and volume data).
Although there are many varieties of malignant lymphomas, one of them is aggressive Bcell lymphoma. Diffuse large Bcell lymphoma (DLBCL) is the largest category of aggressive Bcell lymphomas. Less than 50% of patients can be cured by combination chemotherapy [31]. DLBCL has two important subgroups, which are germinal center Bcelllike (GCB) and activated Bcelllike (ABC) lymphoma. In medicine, cDNA microarrays method can successfully use to distinguish GCB and ABC DLBCL. The advantage of distinguish GCB and ABC DLBCL subgroups has significantly different 5year survival rates after multiagent chemotherapy (GCB over 60%) [32, 33]. A similar situation exists between Burkitt lymphoma (BL) and DLBCL. Both lymphomas were all classified as aggressive Bcell nonHodgkin’s lymphoma in the World Health Organization [34]. Therefore, how to distinguish the difference between BL and DLBCL is a challenge, as the two diseases require different treatment and have different cure rate. Existing diagnosis and classification between BL and DLBCL evaluated their morphologic, immunophenotypic, and cytogenetic features and clinical outcomes [35, 36]. Recently, the Cui et al. study [37] had successful applied nonmedical methods (i.e., statistical and engineering methods, linguistic analysis, and ensembled artificial neural networks) to classify two types of GCB and ABC DLBCL. Because fractal temporal process may generate fluctuations on different area scales that are statistically selfsimilarity [38], therefore, the same concept of fractal temporal process and the statistically selfsimilarity of cell image are used as shown in Figure 1 because the lymphoma cells exist big and small cells at the same time which can easily display statistical selfsimilarity characteristics. In this paper, a nonmedical method/twodimensional (2D) algorithms of DFA has been proposed based on the original design method concepts. The proposed method was used to recharacterize the images of lymph sections. It is anticipated that 2D DFA could be helpful to distinguish BL and DLBCL section images.
2. Material and Methods
2.1. Material
A total of 38 lymph section images cataloged into 3 lymphoma groups were used in the classification as shown in Table 1. Eighteen BL images were classified as group A, which have one to five cytogenetic changes. Ten BL images were classified as group B, which have more than five cytogenetic changes. Both group A and B BLs are high grade aggressive lymphomas, which grow very fast and require more intensive chemotherapy. Finally, ten DLBCL images were grouped as C. Some images of healthy cell, BL and DLBCL (4080 × 3072 pixels), are shown in Figure 2.

(a)
(b)
(c)
(d)
2.2. TwoDimensional Analysis Algorithm of DFA
The algorithm of monodimensional DFA method described in [6] quantifies fractallike correlation properties by calculating the scaling property of the rootmeansquare fluctuation of the integrated and detrended time series data. To illustrate the DFA algorithm, a time series signal (with samples) is used and analyzed as in the following equations: where is the th sample of a signal; is the average of overall signal; and is the value of the th sample of the integrated time series. Then, the fluctuation of integrated and detrended time series for a given window with scale of is calculated by where is the fluctuation of an integrated time series for a given window with scale of , and is the th point on the trend derived using a predetermined window with scale of . A straight line of versus plot indicates the presence of power law (fractal) correlation between scales and the fluctuations of the detrended time series. The slope of loglog plot is defined the DFA scaling exponent α. For the α exponent, Peng et al. indicate it as an indicator that describes the “roughness” of the original time series: the larger the value of α, the smoother the time series [6].
As the application of DFA in image analysis, twodimensional algorithm of DFA should be refined in both integration and detrending processes. Considering the integration in both dimensions, the formula of integration should be as the following steps.
Step 1. The color image should be grayscale processing at first. Consider a selfsimilar surface, which is denoted by a twodimensional array , where and . Content of is the pixel value of surface image. Like monodimensional DFA method, the first integrated equation is as follows: for , , where is the value of the pixel on the integrated image, is the th pixel on the th column, is the average of the th column, is the th pixel on the th row, and is the average of the th row.
Step 2. The surface is partitioned into disjoint square segments of the same size , where , , and min . can be denoted by , segments such that for , , and , where and .
Step 3. Use least square method to calculate the trend of matrix expressed as . The trend matrix is the mathematical model of linear regression. The trend decomposes into row and column, two directions.
Step 4. The mean squared error (MSE) represents the fluctuation of an integrated and detrended segment image with scale of as for value expected to be minimized.
Step 5. The overall detrended fluctuation is calculated by averaging overall the segments, that is Substituting (4) into (5) to get (6) Then Moreover, the least square plans with size of were used to fit the th segment of local trends in images. Then, the th trend can be removed from the integrated image to derive the th fluctuation. DFA scaling exponent is defined as the powerlaw correlation between the using scales and the derived fluctuations. Finally, the DFA scaling exponent α by the straight line of versus plot is obtained.
2.3. Statistical Analysis
Values were expressed as means ± SD. Data were analyzed by oneway analysis of variance (ANOVA) (SigmaStat statistical software, Jandel Scientific, San Rafael, CA). The Tukey test was conducted for multiple comparisons when the null hypothesis was not applicable for the same group. It was also used for all pairwise comparisons of the mean responses to the different treatment groups. Differences were considered significant at a value of .
3. Results and Discussion
3.1. Simulation
24 images ( pixels, same as lymph section images scale) were used to check the proposed twodimensional analysis DFA method for image analysis. Different shape, color, and size images are used. Some images of simulation are shown in Figure 3. The image shape has circle and square, two types, which color has white, black, red, green, blue, and mix color, six modes, and size has large and small, two categories. The scaling exponent α can be estimated by a linear fit on the loglog plot of versus . The monodimensional DFA method has three types α value. It represents the correlation properties of the statistically selfsimilarity. Reflecting on twodimensional DFA method, the global scaling exponent α value was calculated within the range of between and , the shortterm correlation exponent α1 was calculated within the range between and [7, 39, 40], and the longterm correlation exponent α2 was calculated within the range of between and . Plotted versus picture by the DFA analysis method of red small circle image was shown in Figure 4. End of the plotted curve is displayed nearly the horizon line. This trend is similar as other 23 simulation images. Because this curve trend, not displaying the line slope style, is unable to calculate the α and α2 values, we selected α1 to calculate total 24 simulation images and compared between the different shape, color, and size as shown in Table 2. The results show that different shapes, colors, and size images have different α1 values. The 24 simulation images show three cases. First case is similar shape with same color which is all the large size α1 value bigger than the small size. The small size shape means the numbers are more, and the shape area and length have larger values and complexity which are the same α exponent roughness concept of monodimensional DFA method from Peng et al. [6]. It is found that complex images have lower value of α1. Second case is similar size (i.e., the square diagonal length equal circle diameter) and color situation where the square shape α1 value is greater than or equal (only white large image) the circle shape α1 value. This observation corresponds the larger area image (i.e., circle area) and has lower α1 value as initially observed. Third case is similar shape and size situation where α1 values of different colors almost have no significant differences. This means the image color influence is very low when using 2D DFA method to analysis image characteristic. Therefore, when clinical images are influenced by light, the resulting images which are different in colors and cause a change in the characteristic value can be reduced. In summary, the index of the α1 value has distinctive capability and consistency in image analysis.

(a)
(b)
(c)
(d)
3.2. Healthy Cell, BL and DLBCL Lymphoma Image Calculated Results
Figure 5 shows a healthy cell image plotted using versus picture based on the DFA analysis method. Figure 6 shows a BL image of group B. All the simulation images show that the lymphoma images have similar curve trend. Therefore the DFA method was used to calculate α1 value for a healthy cell image for a total of 38 lymphoma images classified into three groups. The α1 value for healthy cell of Figure 2(a) is 0.50 which is due to the clarity of healthy cell images. The mean values of the shortterm correlation exponent α1 of groups A, B, and C are , , and , respectively, as shown in Table 3. It is observed that the healthy cell α1 value was different from lymphoma cells. Furthermore, the α1 value of BL image (both BLs) was significantly lower () than DLBCL image (including groups A and C, groups B and C) and has no difference between both BLs (groups A and B). Hence, it could be concluded that α1 values based on DFA statistics concept could clearly distinguish pathologic states between BL and DLBCL images.
 
Values are expressed as mean ± standard deviation. was considered statistically significant difference using the ANOVA method. ^{ a} for group A versus group C comparison using the Tukey test. ^{ b} for group B versus group C comparison using the Tukey test. 
4. Conclusions
DFA method was utilized to measure the scale invariant behavior that evaluates trends in the presence or absence of fractal correlation properties of time series data [8, 9]. In this study, 2D DFA method was derived and used to explore the trend of fractal images where the α1 value can be easily identified.
In this paper, the DFA method has been applied to image analysis. The twodimensional matrix algorithm of DFA, both integration and detrending processes, was used for time series data field. The DFA method has been used to investigate the characteristic of different type of simulated and lymphoma image. The lymphoma images test results show that the shortterm correlation exponent α1 value of DFA obtained from BL and DLBCL have statistical significant difference. This result is very encouraging, which α1 value could be an index, to help the doctor for distinguishing between BL and DLBCL.
However, the authors had been testing the matrix performance of a twodimensional image as in equation of [26]. The 2D DFA results calculation of the matrix performance are very time consuming. For example, one image (4080 × 3072 pixels) takes over three days. This is disadvantageous for the realtime requirements. Therefore, the 1D DFA concept is being used to solve the aforementioned image’s problem. We assume the lymphoma cell shape, color, and size were influenced by the row and column cells before it. By this way, the verification results of calculation time were great reduced for the 2D DFA, about 3 hours for one image. Fortunately, statistical analysis of calculation results can distinguish between BL and DLBCL lymphoma. However, groups A and B in BL still cannot be classified in this study. Therefore, further investigations are needed to improve the sensitivity and specificity of classification. Results of this study can be compared to 2D DFA algorithm in [26] or 2DDMA algorithm in [41, 42] to investigate how widely this method can be applied in clinical analysis.
Acknowledgments
The authors wish to thank the National Science Council (NSC) of Taiwan (Grant no. NSC 992622E155007CC2) for supporting this research. This research was also supported by the Center for Dynamical Biomarkers and Translational Medicine, National Central University, Taiwan, which is also sponsored by National Science Council (NSC1002911I008001).
References
 B. Mandelbrot, “The variation of certain speculative prices,” The Journal of Business, vol. 36, no. 4, pp. 394–419, 1963. View at: Publisher Site  Google Scholar
 R. Lopes and N. Betrouni, “Fractal and multifractal analysis: a review,” Medical Image Analysis, vol. 13, no. 4, pp. 634–649, 2009. View at: Publisher Site  Google Scholar
 W. Klonowski, “Signal and image analysis using chaos theory and fractal geometry,” Machine Graphics and Vision Journal, vol. 9, pp. 403–432, 2000. View at: Google Scholar
 O. Zmeškal, M. Veselý, M. Nežádal, and M. Buchníček, “Fractal Analysis of Image Structures,” Harmonic and Fractal Image Analysis, pp. 3–5, 2001. View at: Google Scholar
 C. K. Peng, S. V. Buldyrev, S. Havlin, M. Simons, H. E. Stanley, and A. L. Goldberger, “Mosaic organization of DNA nucleotides,” Physical Review E, vol. 49, no. 2, pp. 1685–1689, 1994. View at: Publisher Site  Google Scholar
 C. K. Peng, S. Havlin, H. E. Stanley, and A. L. Goldberger, “Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series,” Chaos, vol. 5, no. 1, pp. 82–87, 1995. View at: Google Scholar
 R. G. Yeh, J. S. Shieh, G. Y. Chen, and C. D. Kuo, “Detrended fluctuation analysis of shortterm heart rate variability in late pregnant women,” Autonomic Neuroscience, vol. 150, no. 12, pp. 122–126, 2009. View at: Publisher Site  Google Scholar
 H. E. Stanley, L. A. N. Amaral, A. L. Goldberger, S. Havlin, P. C. Ivanov, and C. K. Peng, “Statistical physics and physiology: monofractal and multifractal approaches,” Physica A, vol. 270, no. 1, pp. 309–324, 1999. View at: Publisher Site  Google Scholar
 R. G. Yeh, G. Y. Chen, J. S. Shieh, and C. D. Kuo, “Parameter investigation of detrended fluctuation analysis for shortterm human heart rate variability,” Journal of Medical and Biological Engineering, vol. 30, no. 5, pp. 277–282, 2010. View at: Publisher Site  Google Scholar
 T. Nakamura, H. Horio, S. Miyashita, Y. Chiba, and S. Sato, “Identification of development and autonomic nerve activity from heart rate variability in preterm infants,” BioSystems, vol. 79, no. 1–3, pp. 117–124, 2005. View at: Publisher Site  Google Scholar
 N. G. Mahon, A. E. Hedman, M. Padula et al., “Fractal correlation properties of RR interval dynamics in asymptomatic relatives of patients with dilated cardiomyopathy,” European Journal of Heart Failure, vol. 4, no. 2, pp. 151–158, 2002. View at: Publisher Site  Google Scholar
 T. H. Mäkikallio, J. Koistinen, L. Jordaens et al., “Heart rate dynamics before spontaneous onset of ventricular fibrillation in patients with healed myocardial infarcts,” American Journal of Cardiology, vol. 83, no. 6, pp. 880–884, 1999. View at: Publisher Site  Google Scholar
 C. K. Peng, S. V. Buldyrev, A. L. Goldberger et al., “Longrange correlations in nucleotide sequences,” Nature, vol. 356, no. 6365, pp. 168–170, 1992. View at: Publisher Site  Google Scholar
 S. Bahar, J. W. Kantelhardt, A. Neiman et al., “Longrange temporal anticorrelations in paddlefish electroreceptors,” Europhysics Letters, vol. 56, no. 3, pp. 454–460, 2001. View at: Publisher Site  Google Scholar
 S. Blesić, S. Milošević, D. Stratimirović, and M. Ljubisavljević, “Detrended fluctuation analysis of time series of a firing fusimotor neuron,” Physica A, vol. 268, no. 34, pp. 275–282, 1999. View at: Publisher Site  Google Scholar
 J. M. Hausdorff, S. L. Mitchell, R. Firtion et al., “Altered fractal dynamics of gait: reduced strideinterval correlations with aging and Huntington's disease,” Journal of Applied Physiology, vol. 82, no. 1, pp. 262–269, 1997. View at: Google Scholar
 J. M. Lee, D. J. Kim, I. Y. Kim, K. S. Park, and S. I. Kim, “Detrended fluctuation analysis of EEG in sleep apnea using MIT/BIH polysomnography data,” Computers in Biology and Medicine, vol. 32, no. 1, pp. 37–47, 2002. View at: Publisher Site  Google Scholar
 J. M. Lee, D. J. Kim, I. Y. Kim, K. Suk Park, and S. I. Kim, “Nonlinearanalysis of human sleep EEG using detrended fluctuation analysis,” Medical Engineering and Physics, vol. 26, no. 9, pp. 773–776, 2004. View at: Publisher Site  Google Scholar
 M. Staudacher, S. Telser, A. Amann, H. Hinterhuber, and M. RitschMarte, “A new method for changepoint detection developed for online analysis of the heart beat variability during sleep,” Physica A, vol. 349, no. 34, pp. 582–596, 2005. View at: Publisher Site  Google Scholar
 C. J. Stam, T. Montez, B. F. Jones et al., “Disturbed fluctuations of resting state EEG synchronization in Alzheimer's disease,” Clinical Neurophysiology, vol. 116, no. 3, pp. 708–715, 2005. View at: Publisher Site  Google Scholar
 P. GrauCarles, “Longrange powerlaw correlations in stock returns,” Physica A, vol. 299, no. 34, pp. 521–527, 2001. View at: Publisher Site  Google Scholar
 R. Nagarajan and R. G. Kavasseri, “Minimizing the effect of periodic and quasiperiodic trends in detrended fluctuation analysis,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 777–784, 2005. View at: Publisher Site  Google Scholar
 R. Weron, “Estimating longrange dependence: finite sample properties and confidence intervals,” Physica A, vol. 312, no. 12, pp. 285–299, 2002. View at: Publisher Site  Google Scholar
 R. G. Yeh, J. S. Shieh, Y. Y. Han, Y. J. Wang, and S. C. Tseng, “Detrended fluctuation analyses of shortterm heart rate variability in surgical intensive care units,” Biomedical Engineering, vol. 18, no. 2, pp. 67–72, 2006. View at: Publisher Site  Google Scholar
 R. G. Yeh, Y. Y. Han, J. S. Shieh, Y. J. Wang, S. C. Tseng, and Y. C. Fu, “Nonrandomness index applied for heart rate variability in surgical intensive care units using frequency and rank order statistics,” Biomedical Engineering, vol. 19, pp. 303–311, 2007. View at: Publisher Site  Google Scholar
 G. F. Gu and W. X. Zhou, “Detrended fluctuation analysis for fractals and multifractals in higher dimensions,” Physical Review E, vol. 74, no. 6, Article ID 061104, 2006. View at: Publisher Site  Google Scholar
 W. X. Zhou, “Multifractal detrended crosscorrelation analysis for two nonstationary signals,” Physical Review E, vol. 77, no. 6, Article ID 066211, 2008. View at: Publisher Site  Google Scholar
 K. Dong and P. Shang, “Statistical properties of detrended crosscorrelation analysis,” Journal of Beijing Jiaotong University, vol. 34, no. 6, pp. 64–67, 2010. View at: Google Scholar
 B. Podobnik and H. E. Stanley, “Detrended crosscorrelation analysis: a new method for analyzing two nonstationary time series,” Physical Review Letters, vol. 100, no. 8, Article ID 084102, 2008. View at: Publisher Site  Google Scholar
 Z. Q. Jiang and W. X. Zhou, “Multifractal detrending movingaverage crosscorrelation analysis,” Physical Review E, vol. 84, no. 1, Article ID 016106, 2011. View at: Publisher Site  Google Scholar
 T. Ando, M. Suguro, T. Hanai, T. Kobayashi, H. Honda, and M. Seto, “Fuzzy neural network applied to gene expression profiling for predicting the prognosis of diffuse large Bcell Lymphoma,” Japanese Journal of Cancer Research, vol. 93, no. 11, pp. 1207–1212, 2002. View at: Google Scholar
 G. Wright, B. Tan, A. Rosenwald, E. H. Hurt, A. Wiestner, and L. M. Staudt, “A gene expressionbased method to diagnose clinically distinct subgroups of diffuse large B cell lymphoma,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 17, pp. 9991–9996, 2003. View at: Publisher Site  Google Scholar
 C. P. Hans, D. D. Weisenburger, T. C. Greiner et al., “Confirmation of the molecular classification of diffuse large Bcell lymphoma by immunohistochemistry using a tissue microarray,” Blood, vol. 103, no. 1, pp. 275–282, 2004. View at: Publisher Site  Google Scholar
 C. Bellan, L. Stefano, D. F. Giulia, E. A. Rogena, and L. Lorenzo, “Burkitt lymphoma versus diffuse large Bcell lymphoma: a practical approach,” Hematological Oncology, vol. 28, no. 2, pp. 53–56, 2010. View at: Publisher Site  Google Scholar
 N. Nakamura, H. Nakamine, J. I. Tamaru et al., “The distinction between Burkitt lymphoma and diffuse large Bcell lymphoma with cmyc rearrangement,” Modern Pathology, vol. 15, no. 7, pp. 771–776, 2002. View at: Publisher Site  Google Scholar
 X. F. Zhao, A. Hassan, A. Perry, Y. Ning, S. A. Stass, and L. P. Dehner, “CMYC rearrangements are frequent in aggressive mature Bcell lymphoma with atypical morphology,” International Journal of Clinical and Experimental Pathology, vol. 1, no. 1, pp. 65–74, 2008. View at: Google Scholar
 X. Cui, C. W. Lin, M. F. Abbod, Q. Liu, and J. S. Shieh, “Diffuse large Bcell lymphoma classification using linguistic analysis and ensembled artificial neural networks,” Journal of the Taiwan Institute of Chemical Engineers, vol. 43, no. 1, pp. 15–23, 2012. View at: Publisher Site  Google Scholar
 A. L. Goldberger, L. A. N. Amaral, J. M. Hausdorff, P. C. Ivanov, C. K. Peng, and H. E. Stanley, “Fractal dynamics in physiology: alterations with disease and aging,” Proceedings of the National Academy of Sciences of the United States of America, vol. 99, no. 1, pp. 2466–2472, 2002. View at: Publisher Site  Google Scholar
 T. H. Mäkikallio, T. Seppänen, K. E. J. Airaksinen et al., “Dynamic analysis of heart rate may predict subsequent ventricular tachycardia after myocardial infarction,” American Journal of Cardiology, vol. 80, no. 6, pp. 779–783, 1997. View at: Publisher Site  Google Scholar
 T. H. Mäkikallio, T. Ristimäe, K. E. J. Airaksinen, C. K. Peng, A. L. Goldberger, and H. V. Huikuri, “Heart rate dynamics in patients with stable angina pectoris and utility of fractal and complexity measures,” American Journal of Cardiology, vol. 81, no. 1, pp. 27–31, 1998. View at: Publisher Site  Google Scholar
 A. Carbone, “Algorithm to estimate the Hurst exponent of highdimensional fractals,” Physical Review E, vol. 76, no. 5, Article ID 056703, 2007. View at: Publisher Site  Google Scholar
 G. F. Gu and W. X. Zhou, “Detrending moving average algorithm for multifractals,” Physical Review E, vol. 82, no. 1, Article ID 011136, 2010. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2012 RongGuan Yeh et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.