Rui Guo 1, 2 and Yiqin Wang 1 and Hanxia Yan 1 and Jianjun Yan 3 and Fengyin Yuan 4 and Zhaoxia Xu 1 and Guoping Liu 1 and Wenjie Xu 1
Academic Editor:Josiah Poon
1, Laboratory of Information Access and Synthesis of TCM Four Diagnosis, Shanghai Univerisity of Traditional Chinese Medicine, Shanghai 201203, China
2, Center for TCM Information Science and Technology, Shanghai University of Traditional Chinese Medicine, Shanghai 201203, China
3, Center for Mechatronics Engineering, East China University of Science and Technology, Shanghai 200237, China
4, The first Affiliated Hospital, School of clinical Medicine of G.D.P.U, Guangzhou 51000, China
Received 5 September 2014; Accepted 15 January 2015; 9 June 2015
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.
1. Introduction
Traditional Chinese medicine (TCM) is an ancient medical practice system which emphasizes regulating the integrity of the human body and its interrelationship with natural environments [1]. Zheng (meaning syndrome or pattern) is a unique TCM concept. It is the overall physiological and/or pathological pattern of the human body in response to a given internal and external condition, which usually is an abstraction of internal disharmony defined by a comprehensive analysis of the clinical symptoms and signs gathered by a practitioner using inspection, auscultation, olfaction, interrogation, and palpation of the pulses [2]. Chinese practitioners diagnose diseases through "Zheng differentiation." The Zheng differentiation of TCM considers the etiology, location, nature, and condition of a disease during a specific stage of the disease process based on clinical symptoms and signs. Pulse taking is one of the key methods to gather the signs and symptoms of patients by a practitioner. During pulse taking, TCM practitioners place their fingers on the radial artery, from which various physiological and pathological conditions can be detected. Traditional pulse taking has important clinical value on the diagnosis and prognosis of the diseases, especially angiocardiopathy. Accurate pulse taking can only be done by TCM practitioners with years of experience. Therefore, objective and quantitative pulse diagnosis is highly desirable, which is of help to establish the objective and quantitative criteria on disease diagnosis or Zheng differentiation.
Coronary heart disease (CHD) is considered a primary cause of death in developed countries and is predicted to be one of the most common causes of death worldwide by 2020 [3]. Early diagnosis and prevention of CHD are essential and have critical public health implications. Identifying, vascular lesions in early stages to reverse and prevent CHD, stroke, sudden death, and other malignant vascular events are crucial [4]. Previous studies have been conducted to clarify and identify subclinical vascular disease. Thus, a noninvasive, convenient, and efficient method should be developed to detect vascular lesions. In TCM, visceral pathological changes and other information can be obtained by detecting pulses. For instance, pathological changes in CHD are clearly reflected in pulse diagnosis information. Therefore, CHD could be considered a breakthrough point for providing a theoretical basis for investigation the TCM pulse.
Modern medical research has shown that the arterial pulse is caused by heart contraction. The left ventricle ejects blood into the aorta through the aortic valve, causing the velocity, pressure, and diameter in the arterial tree to pulsate [5]. The signal acquired from the radial artery is the comprehensive reflection of the wave form (shape), velocity (fast or slow), period (rhythm), and swing (intensity) of pulse waves in the radial artery [6]. Pulse waves contain human physiological and pathological information. Pulse diagnosis is conducted to investigate this information. Thus, pressure pulse charts, which show pulses in TCM, can be used noninvasively and conveniently to provide insight into visceral diseases, particularly cardiovascular diseases such as CHD. Subjective judgments and description of pulse in TCM rely on the experience of doctors. However, clinical application and development of this technique are restricted. Thus, further study should be conducted using various modern information-processing methods, including time domain analysis [7]; frequency domain analysis, such as the Fourier transform [8]; and combined time-frequency analysis, such as wavelet decomposition [9]. Many quantifiable parameters can be obtained; moreover, fuzzy clustering, the Bayes classifier, support vector machine, artificial neural networks, and other methods can be used to classify and recognize pulse [10-13]. For feature extraction of pulse, although the Fourier transform provides the average distribution of signal energy, this technique fails to characterize time-varying information of signal frequency and cannot describe the time domain local features of signals. For characterizing the nonstationary signals of pulse waves, an extra high-frequency harmonic signal is necessary. However, a correct and reasonable interpretation of the signal cannot be provided because the high-frequency harmonic signal is noninherent [14]. Although wavelet transform is essentially a Fourier transform with an adjustable window, the signal in a wavelet window must be stable and cannot eliminate the limitations of Fourier analysis.
Hilbert-Huang transform (HHT) is a new self-adapting time-frequency analytic method. This method can be used to conduct self-adapting time-frequency decomposition according to local time-varying characteristics and can overcome the defects of insignificant harmonic component showing nonstationary and nonlinear signal in traditional methods. Moreover, this method enables obtaining high time-frequency resolution and sufficient time-frequency aggregation; hence, this technique is suitable for nonstationary and nonlinear signal analysis [15]. HHT is composed of empirical mode decomposition (EMD) and the Hilbert transform. The core of this technique is EMD, which is performed on the basis of the time-scale feature of data. Therefore, EMD is more suitable for nonstationary and nonlinear data processing than the Fourier and wavelet methods depending on the transcendental function based on the decomposition method. We used the random forest algorithm as the classifier. The training and prediction speeds of this algorithm are high, and an internal unbiased estimation of a generalization error can be generated. The interaction between features and their degree of importance can be detected, and the over fitting does not occur. For an unbalanced classified data set, this algorithm can balance the error and can be easily parallelized. Therefore, we analyzed and determined the pulse condition of patients with CHD in this study by using EMD time series analysis method and the random forest recognition algorithm.
2. Clinical Material
TCM pulse refers to the pulse sensed by doctors as they palpating the examinee's radial artery with their fingers. Imitating TCM doctors, measurement equipment (cooperatively developed by our research team and Shanghai Asia-Pacific Computer Co. Ltd) was employed to acquire pulse recordings, which provided the basis for objective pulse analysis.
Pulse recordings used in this study were acquired from 342 volunteers for 60 sec with a sampling rate of 720 Hz. Two groups without respiratory system and nervous system disorders were studied. Each subject was instructed to relax for more than 5 min before pulse was recorded.
Group 1 included 225 inpatients with CHD aged [figure omitted; refer to PDF] years from Longhua Hospital and Shuguang Hospital, which are affiliated to Shanghai University of Traditional Chinese Medicine.
Group 2 included 117 normal subjects, who are selected as control subjects aged [figure omitted; refer to PDF] years. The subjects were players in the "2010 Zhangjiang ball game competition for the elderly" and staff from Shanghai University of Traditional Chinese Medicine. These subjects have no documented history of cardiovascular disorders.
3. Method
3.1. HHT
Huang et al. proposed the HHT method [16]. The proposed technique is based on the intrinsic mode function (IMF) and EMD. EMD involves decomposing a given signal from a small scale to a large scale to obtain the component signal IMF according to local characteristic time scale. The IMF obtained through decomposition must satisfy two conditions: [figure omitted; refer to PDF] for the entire signal length, the numbers of extreme points and zero crossing of an IMF component must be equal to or differ by 1 at most.; [figure omitted; refer to PDF] at any time, the average value of an upper envelope point defined by the maximum and a lower envelope defined by the minimum is 0.
EMD was performed as follows.
[figure omitted; refer to PDF] Three sample interpolation fittings were used for obtaining the upper envelope curves and lower envelope curves of the signal to calculate the average value of the upper and lower envelope curves at each point and, thus, obtain the average curve [figure omitted; refer to PDF] .
[figure omitted; refer to PDF] The average curve minus [figure omitted; refer to PDF] was subtracted from the original signal [figure omitted; refer to PDF] . The first component [figure omitted; refer to PDF] was subsequently obtained.
If [figure omitted; refer to PDF] satisfied two conditions of the IMF, then [figure omitted; refer to PDF] was the IMF component in the first order. Otherwise, the difference value between [figure omitted; refer to PDF] and the other envelope median value was calculated again. The IMF component [figure omitted; refer to PDF] in the first stage can only be obtained if the difference value sequence satisfies the two conditions of the IMF.
[figure omitted; refer to PDF] The component [figure omitted; refer to PDF] was subtracted from the original signal to obtain the residual signal of the original signal [figure omitted; refer to PDF] . The signal [figure omitted; refer to PDF] was redefined as the original signal. Steps [figure omitted; refer to PDF] to [figure omitted; refer to PDF] were repeated and [figure omitted; refer to PDF] IMF components were obtained until [figure omitted; refer to PDF] was converted into monotone function or reached a constant value.
By performing EMD, n IMF components and a residual signal [figure omitted; refer to PDF] can be obtained. Thus, the original signal can be represented using the following equation: [figure omitted; refer to PDF]
Each IMF component [figure omitted; refer to PDF] and its Hilbert transform were used to construct an analytic signal as shown in the following equation: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] is the amplitude function, which shows the instantaneous amplitude energy of the signal at each sampling point and [figure omitted; refer to PDF] is the phase function, which shows the instantaneous phase of a signal at each sampling point; instantaneous frequency [figure omitted; refer to PDF] can be obtained by calculating its derivative. Thus, the amplitude and frequency of the signal are functions of time and are plotted on the time-frequency plane to obtain the Hilbert spectra [figure omitted; refer to PDF] . Hilbert showed the global transformation rule for a signal amplitude with time and frequency conversions; this corresponds to the distribution of signal energy in various characteristic scales (time or frequency) to a certain extent.
3.2. Energy Feature of a Pulse Signal Based on EMD
IMF component was obtained after a signal was decomposed through EMD. The energy of each IMF component was calculated according to the following equation: [figure omitted; refer to PDF]
For subsequent statistical analysis, normalized energy can be obtained according to the following: [figure omitted; refer to PDF]
3.3. Sample Entropy of Pulse Signals Based on EMD
Richman and Moorman [17] developed a new and related complexity measure, sample entropy (SampEn), based on the research of Grassberger, and so forth. The sample entropy calculation of each IMF component is introduced as follows.
[figure omitted; refer to PDF] For a time series (pulse signal) of [figure omitted; refer to PDF] points [figure omitted; refer to PDF] , [figure omitted; refer to PDF] . The time series [figure omitted; refer to PDF] forms a set of m -dimensional vectors [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , where [figure omitted; refer to PDF] is the vector of [figure omitted; refer to PDF] data points from [figure omitted; refer to PDF] to [figure omitted; refer to PDF] .
[figure omitted; refer to PDF] The distance between two such vectors is defined as the maximum difference in their corresponding scalar components: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] .
[figure omitted; refer to PDF] Given [figure omitted; refer to PDF] of [figure omitted; refer to PDF] , for every value from 1 to [figure omitted; refer to PDF] , the number of [figure omitted; refer to PDF] is calculated, and [figure omitted; refer to PDF] to exclude self-matches. The ratio of the number to [figure omitted; refer to PDF] is defined as follows: [figure omitted; refer to PDF] where [figure omitted; refer to PDF] . The average value of [figure omitted; refer to PDF] is defined as follows: [figure omitted; refer to PDF]
[figure omitted; refer to PDF] Given [figure omitted; refer to PDF] of [figure omitted; refer to PDF] , set [figure omitted; refer to PDF] where [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] . The average value of [figure omitted; refer to PDF] is then obtained using the following equation: [figure omitted; refer to PDF]
(5) [figure omitted; refer to PDF] is the probability that two sequences match at [figure omitted; refer to PDF] points, whereas [figure omitted; refer to PDF] is the probability that two sequences match for [figure omitted; refer to PDF] points. The parameter [figure omitted; refer to PDF] is defined as follows: [figure omitted; refer to PDF]
The sample entropy of pulse signal with finite length [figure omitted; refer to PDF] is estimated using the following statistic: [figure omitted; refer to PDF]
The parameters [figure omitted; refer to PDF] and [figure omitted; refer to PDF] are used to estimate SampEn [18]. Pincus suggested that [figure omitted; refer to PDF] , [figure omitted; refer to PDF] to 0.25δ , and [figure omitted; refer to PDF] is the standard deviation of the original signal [figure omitted; refer to PDF] , [figure omitted; refer to PDF] . [figure omitted; refer to PDF] reflects the rate of information production as [figure omitted; refer to PDF] increases from 2 to 3. The higher the SampEn value is, the higher the rate of the information production, indicating that the signal is more complex.
3.4. Random Forest Recognition Method
A random forest [19] is composed of numerous decision trees, which are formed using a stochastic method. Thus, it is also called a random decision tree. Trees in a random forest do not correlate. After test data are used as input in a random forest to classify each decision tree, the category with the highest classification results in all decision trees is selected as the final result. Therefore, a random forest is a classifier that contains multiple-decision trees, and its output category relies on the mode of output categories of individual trees.
A random forest resampling technique uses the bootstrap method, which entails repeatedly and randomly selecting [figure omitted; refer to PDF] samples from the original training sample set [figure omitted; refer to PDF] to generate new training sample sets. Subsequently, [figure omitted; refer to PDF] classification trees are generated according to the bootstrap sample set to construct random forests. The classification results of the new data rely on the score formed by the vote of the classification tree. The algorithm is presented as follows.
[figure omitted; refer to PDF] The original training set is [figure omitted; refer to PDF] , and the bootstrap method is used to randomly select [figure omitted; refer to PDF] new self-help sample sets and to construct [figure omitted; refer to PDF] classification trees. Samples not drawn at each time constitute [figure omitted; refer to PDF] data outside the bag.
[figure omitted; refer to PDF] In total, [figure omitted; refer to PDF] variables are set, and [figure omitted; refer to PDF] variables ( [figure omitted; refer to PDF] ) are randomly selected at each node of each tree. The variable with the greatest classification ability is selected. The variable classification threshold is determined by examining each classification point.
[figure omitted; refer to PDF] Each tree grows to the maximum size, with on pruning.
[figure omitted; refer to PDF] The generated multiple classification trees constitute random forests. New data are obtained and classified according to the random forest classifier. The classification results depend on the number of votes provided by the tree classifiers.
Random forests are an improvement of the decision tree algorithm, in which multiple decision trees are merged. Each tree is established on the basis of an independently extracted sample. All of the trees in the forest are uniformly distributed. Classification error depends on the classification ability of each tree and the correlation among these trees. Feature selection is performed to divide each node by using a stochastic method. The errors generated in various circumstances are then compared. Internal estimation error, or the classification and correlation ability, is detected to determine the number of selected features. The classification capability of a single tree may be low. The most likely classification of a test sample is selected after a high number of decision trees are randomly generated and after statistical analysis is performed according to the classification result of each tree.
The number of decision trees in a random forest in this study was 500, and [figure omitted; refer to PDF] took the mean square root of [figure omitted; refer to PDF] .
4. Results
4.1. Statistical Analysis of Energy and Sample Entropy of Pulse Signal Based on EMD
EMD involves adaptively decomposing a signal frequency into a series of IMFs from a high level to a low level. IMF at each level adaptively showed signal characteristics with various resolutions. The amplitude of the IMF at each level at each time differed. The amplitude showed the strength change in the signal at a modal. Most pulse signals can be decomposed to level 7 or higher through EMD (IMF1 -IM [figure omitted; refer to PDF] , [figure omitted; refer to PDF] ), and only two pulse signals of patients with CHD could be decomposed to level 6 ( [figure omitted; refer to PDF] , [figure omitted; refer to PDF] ). In EMD, the modal energy at high levels was low; their effects on the entire system were weak. Therefore, we analyzed the IMF1 -IMF7 in front of all modals further. The IMF at each level of a patient with CHD and one normal pulse signal after EMD are shown in Figure 1.
Figure 1: IMFs of a pulse graph for a normal subject and a patient with CHD after EMD.
[figure omitted; refer to PDF]
[figure omitted; refer to PDF]
Figure 1 shows that the components, including IMF1 -IMF7 and the residual parameters res., were obtained after EMD of the pulse signal. The frequencies of IMF1 -IMF7 decreased successively, and the amplitudes of IMF1 -IMF7 increased progressively. The differences of IMFs between Figures 1(a) and 1(b) can be observed. For example, Figure 1(b) shows that components of IMF1 -IMF7 had higher morphological variation than those in Figure 1(a), reflecting the irregularity of a normal pulse signal, and had more high-frequency parts especially in IMF3 -IMF7 than those in Figure 1(a). In order to quantitatively describe the differences between the CHD patients and the healthy subjects, we extracted the IMF energy and the IMF sample entropy of the pulse signals to make an analysis.
We observed that the variances were nonhomogeneous in the distribution of the IMF energy and the IMF sample entropy of the pulse signals in the CHD and normal groups by using IBM SPSS20.0 statistical software. Thus, we used a nonparametric test for statistical analysis. For the nonparametric test of independent samples of the two groups, we used the rank sum test method to calculate the statistical difference between the two groups. Table 1 shows the statistical difference in the average rank of IMF energy between the two groups. The average rank with IMF normalized energy in the normal group was significantly greater than that in the CHD group. Table 2 shows the statistical difference in the average rank of the IMF sample entropy between the two groups. Owing to the high-frequency modes [figure omitted; refer to PDF] and [figure omitted; refer to PDF] were caused by interference, the modes [figure omitted; refer to PDF] and [figure omitted; refer to PDF] were discarded without the following analysis. The intermediate-frequency modes, such as [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , [figure omitted; refer to PDF] , and [figure omitted; refer to PDF] in the CHD group were significantly lower than those in the normal group. No statistically significant difference was observed in [figure omitted; refer to PDF] between the two groups.
Table 1: Rank sum test of IMF energy in two groups.
| (Mean rank) | ||||||
Group | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] |
The CHD | 172.10 | 132.50 | 148.68 | 153.93 | 152.81 | 146.98 | 142.89 |
The normal subjects | 250.65* | 240.32* | 209.35* | 199.29* | 201.45* | 212.59* | 220.43* |
[figure omitted; refer to PDF] P < 0.001, versus normal group.
Table 2: Rank sum test of IMF sample entropy in two groups.
| (Mean rank) | ||||||
Group | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] | [figure omitted; refer to PDF] |
The CHD | 180.28 | 188.17 | 156.74 | 143.86 | 150.93 | 158.98 | 164.09 |
The normal subjects | 148.86* | 133.77* | 193.91* | 218.56* | 205.04* | 189.64* | 179.84 |
[figure omitted; refer to PDF] P < 0.001, versus normal group.
4.2. Pulse Recognition Based on the Random Forest Classifier
We classified and recognized the energy and sample entropy characteristics of the IMFs of the two groups of pulses by using random forest classifier, and the recognition results are shown in Table 3.
Table 3: Average recognition rate of different pulse characteristics obtained using random forest (%).
Group | Feature | ||
IMF sample entropy | IMF entropy | Combination of IMF sample entropy and IMF sample entropy | |
The CHD group | 88.32 ± 3.89 | 89.20 ± 6.76 | 95.49 ± 2.78 |
The normal group | 53.48 ± 10.15 | 73.99 ± 15.47 | 80.07 ± 9.13 |
The two groups | 76.35 ± 4.77 | 84.00 ± 5.04 | 90.21 ± 4.09 |
Table 3 shows that the average recognition rate was 76.35%, when we used only the sample entropy of IMFs as the feature vector to recognize pulses in the two groups. Furthermore, the average recognition rate was 84%, when we used only the IMF energy as the feature vector to recognize the pulses in the two groups. If the IMF energy and IMF sample entropy were jointly used as the feature vector to recognize pulses in the two groups, the average recognition rate could reach 90.21%.
5. Discussion
Adaptive HHT exhibits a high time-frequency resolution and local orthogonal self-adaption; moreover, this technique is superior to the wavelet transform and other signal analysis methods. Adaptive HHT can be used to examine instantaneous frequency and energy at multiple scales and is a powerful tool for analyzing biomedical signals [20]. The IMF sequence of pulse signals obtained using EMD is directly separated from the original time series data. The number of decompositions need not be determined in advance without being affected by human factors. Thus, the IMF sequence can effectively reflect the inherent physical characteristics of original data. Its decomposition is objective, inherent, and adaptive. The IMF at each level represents band information with a specific meaning. Table 1 shows that the IMF energy at each level in the CHD group was significantly lower than that in the normal group, which indicates that the energy of the cardiovascular system in patients with CHD was likely lower than that in the normal subjects. Human biological waves such as the pulse wave can directly reflect physical function and health status. Human biological wave energy can be produced through the metabolism of water, air, sunlight, food, and other substances, which are absorbed by molecules and cells of the human body and can be naturally supplied with nutrients and oxygen. High pulse wave energy in the healthy people may result when qi of the viscera and meridians moves harmoniously. As Table 2 shows, the sample entropy of pulse signals in patients with CHD in modes IMF3 -IMF6 , which showed statistically significant differences, was lower than that in the normal subjects. Entropy is the rate at which new information is generated; the entropy value reflects the complexity of a system. A higher entropy value shows that the signal generated by the system is possibly random and irregular and that the system adaptability is stronger. A lower entropy value shows that the signal generated by the system is simpler and regular, and that the system adaptability is lower. The result of Table 2 indicated that the adaptive ability of the cardiovascular system of patients with CHD was likely lower than that of the normal subjects. This finding is consistent with a relatively weaker physiological complexity of the human body in the pathological state and, thus, corresponds to regular pathological states. The energy and sample entropy of IMFs are valuable features for characterizing the pathological state of pulse in patients with CHD. The HHT method can be used to analyze pulse signals of nonstationary dynamic change. This method can sensitively capture the primary features of various pulse signal components with the dynamic changes in time and frequency. The energy and sample entropy of IMF provides a new basis for the feature extraction and pattern recognition of various pulses.
Random forest is a classifier that produces highly accurate with high training and prediction speeds. This classifier can generate an internal unbiased estimation of a generalization error during classification. Furthermore, overfitting never occurs. Random forest can also balance the errors in unbalanced data sets. In this study, we classified and determined the IMF energy and IMF sample entropy of pulse signals by using random forest classifier. The energy of IMFs, sample entropy of IMFs, and their combination were inputted as pulse feature vectors; the corresponding average recognition rates were 84%, 76.35%, and 90.21%, respectively. Compared with the separate use of IMF energy or IMF sample entropy as a feature vector, the combined use of IMF energy and IMF sample entropy as a feature vector improved the classified accuracy for the CHD group and the normal group. Although the sample size in the CHD group was differed from that in the normal group, the unbalanced sample capacity of random forest was superior. Nevertheless, we achieved satisfactory classification results.
6. Conclusion
Pulse diagnosis is a characteristic diagnostic method in TCM. Pulse detection is a noninvasive technique with the simple and easy operation and the stable performance, which does not require expensive equipment. A pressure pulse signal that corresponds to the pulse condition in TCM can be detected noninvasively and conveniently to obtain the pathological and physiological status of the cardiovascular system. We extracted the energy and sample entropy of IMFs of pulse signals as pulse features and used random forests as the classifier of pulse signals. The results illustrate that the proposed methods or pulse-signal processing and classification is effective and efficient. This study offers a new method for developing and promoting a noninvasive pulse diagnostic technique. Furthermore, this research provides objective and quantitative parameters of TCM pulse and lays a foundation for research on objective and quantitative criteria on disease diagnosis or Zheng differentiation.
Acknowledgments
This work was supported by National Natural Science Foundation of China for Youth under Grant no. 81302913 and National Natural Science Foundation of China under Grant nos. 81173199 and 81473594.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
[1] Y. Wang, A. Xu, "Zheng: a systems biology approach to diagnosis and treatments," Science , vol. 346, no. 6216, supplement, pp. S13-S15, 2014.
[2] F. Cheung, "TCM: made in China," Nature , vol. 480, no. 7378, pp. S82-S83, 2011.
[3] E. L. Jones, R. P. Naoumova, "Genetics of amilial combined hyperlipidemia and risk of coronary heart disease," Human Molecular Genetics , vol. 13, pp. 149-160, 2004.
[4] H.-Y. Wang, J. Wang, W.-P. Liu, "Application of early vascular disease detection system in patients with cardiovascular disease," Medical Journal of Chinese People's Health , vol. 18, no. 5, pp. 336-339, 2006.
[5] W. W. Nichols, M. F. O'Rourke McDonald's Blood Flow in Arteries: Theoretic, Exprimental and Clinical Principles , Edward Arnold, London, UK, 2005., 5th.
[6] Z. R. Liu, "Pulse and hemodynamics in traditional Chinese medicine," Chinese Journal of Science , vol. 5, no. 6, pp. 411-413, 1982.
[7] H. X. Yan, Y. Q. Wang, F. F. Li, "Research and application of pulse analysis," Chinese Archives of Traditional Chinese Medicine , vol. 23, no. 1, pp. 129-130, 2005.
[8] B. Jiang, Z. C. Song, "Spectral analysis of pulse signals," Techniques of Automation and Application , vol. 26, no. 8, pp. 38-40, 2007.
[9] L. Chen, L.-J. Yang, "Identifying pulse signals based on fractal theories and wavelet," Aeronautical Computing Technique , vol. 38, no. 4, pp. 41-44, 2008.
[10] J.-J. Shu, Y. Sun, "Developing classification indices for Chinese pulse diagnosis," Complementary Therapies in Medicine , vol. 15, no. 3, pp. 190-198, 2007.
[11] W.-P. Zhang, Y. Zhang, S.-S. Zhang, "Application on wavelet and neural network in the analysis and pattern recognition of the manifestation of the pulse for detection of cerebrovascular disease," Progress in Biomedical Engineering , vol. 29, no. 2, pp. 84-86, 2008.
[12] Q.-L. Guo, K.-Q. Wang, D.-Y. Zhang, N.-M. Li, "A wavelet packet based pulse waveform analysis for cholecystitis and nephrotic syndrome diagnosis," in Proceedings of the International Conference on Wavelet Analysis and Pattern Recognition (ICWAPR '08), pp. 513-517, Hong Kong, August 2008.
[13] H.-Y. Wang, S. Xu, "Automatic pulse recognition method based on bayesian classifier," Chinese Journal of Biomedical Engineering , vol. 28, no. 5, pp. 735-742, 2009.
[14] Z.-D. Zhao, X.-H. Tang, Z.-J. Zhao, M. Pan, Y.-Q. Chen, "Spectrum analysis of heart sound signal based on Hilbert-Huang transform," Chinese Journal of Sensors and Actuators , vol. 18, no. 1, pp. 18-22, 2005.
[15] Y. De-jie, C. Jun-seng, Y. yu Application of Hilbert-Huang Transform in Mechanical Fault Diagnosis , Science Press, Beijing, China, 2006.
[16] N. E. Huang, Z. Shen, S. R. Long, M. C. Wu, H. H. Shih, Q. Zheng, N.-C. Yen, C. C. Tung, H. H. Liu, "The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis," Proceedings of the Royal Society of London , vol. 454, no. 1971, pp. 903-995, 1998.
[17] J. S. Richman, J. R. Moorman, "Physiologica time-series analysis using approximate entropy and sample entropy," The American Journal of Physiology-Heart and Circulatory Physiology , vol. 278, no. 6, pp. H2039-H2049, 2000.
[18] S. M. Pincus, "Assessing serial irregularity and its implications for health," Annals of the New York Academy of Sciences , vol. 954, pp. 245-267, 2001.
[19] L. Breiman, "Random forests," Machine Learning , vol. 45, no. 1, pp. 5-32, 2001.
[20] C.-K. Peng, "Potential HHT Applications in Biomedical Signal Analysis," http://www.fuentek.com/technologies/hht.htm
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2015 Rui Guo 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.
Abstract
Objective. This research provides objective and quantitative parameters of the traditional Chinese medicine (TCM) pulse conditions for distinguishing between patients with the coronary heart disease (CHD) and normal people by using the proposed classification approach based on Hilbert-Huang transform (HHT) and random forest. Methods. The energy and the sample entropy features were extracted by applying the HHT to TCM pulse by treating these pulse signals as time series. By using the random forest classifier, the extracted two types of features and their combination were, respectively, used as input data to establish classification model. Results. Statistical results showed that there were significant differences in the pulse energy and sample entropy between the CHD group and the normal group. Moreover, the energy features, sample entropy features, and their combination were inputted as pulse feature vectors; the corresponding average recognition rates were 84%, 76.35%, and 90.21%, respectively. Conclusion. The proposed approach could be appropriately used to analyze pulses of patients with CHD, which can lay a foundation for research on objective and quantitative criteria on disease diagnosis or Zheng differentiation.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer