About this Journal Submit a Manuscript Table of Contents
Evidence-Based Complementary and Alternative Medicine
Volume 2012 (2012), Article ID 563958, 11 pages
http://dx.doi.org/10.1155/2012/563958
Research Article

Novel Diagnostic Model for the Deficient and Excess Pulse Qualities

Division of Constitutional Medicine Research, Korea Institute of Oriental Medicine, Daejeon 305-811, Republic of Korea

Received 5 May 2011; Revised 8 June 2011; Accepted 16 June 2011

Academic Editor: Carlo Ventura

Copyright © 2012 Jaeuk U. Kim 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

The deficient and excess pulse qualities (DEPs) are the two representatives of the deficiency and excess syndromes, respectively. Despite its importance in the objectification of pulse diagnosis, a reliable classification model for the DEPs has not been reported to date. In this work, we propose a classification method for the DEPs based on a clinical study. First, through factor analysis and Fisher's discriminant analysis, we show that all the pulse amplitudes obtained at various applied pressures at Chon, Gwan, and Cheok contribute on equal orders of magnitude in the determination of the DEPs. Then, we discuss that the pulse pressure or the average pulse amplitude is appropriate for describing the collective behaviors of the pulse amplitudes and a simple and reliable classification can be constructed from either quantity. Finally, we propose an enhanced classification model that combines the two complementary variables sequentially.

1. Introduction

Pulse wave is a pressure wave propagating through the arterial system, generated by the periodic contraction and relaxation of the heart, and its characteristics are influenced by the compliance of the vascular system, blood viscosity, and the functions of major organs. By diagnosing the pulse, trained practitioners can gather elaborate physiological and pathological information on the cardiovascular system, organ functions, patients’ constitution, emotional conditions, behavioral patterns, and previous illness, as well as body’s homeostatic balance [15].

Pulse diagnosis has been considered a core component of diagnostics in Oriental medicine for thousands of years. In contemporary Oriental medicine, pulse diagnosis is made dominantly at three adjacent positions along the radial artery in both wrists. A palpation position called Gwan is located on the radial artery closest to the styloid process. Chon is about 10 mm distal from Gwan and Cheok about 10 to 15 mm proximal from Gwan [6]. To diagnose the pulse, an Oriental medical doctor (OMD) places the index, middle, and ring fingers, respectively, at Chon, Gwan, and Cheok and applies varying pressure simultaneously or sequentially to determine the pulse qualities.

A recent survey indicates that about 22% of OMDs rely on pulse diagnosis as the primary diagnostic method, which is next to the inquiry (38%) and observation (27%). About 71% of the survey participants asserted that pulse diagnosis was in their diagnoses [7]. Despite its importance and frequent use in clinics, pulse diagnosis has been criticized for the lack of scientific evidence and for the manual palpation and subjective interpretation of pulse qualities. To provide sound scientific evidence and to overcome the experiential boundary of pulse diagnosis, it is essential to develop objective techniques with standardized protocols to obtain pulse signals, and interpret them into pulse qualities defined in terms of a few measureable physical parameters [1, 8, 9].

With advances in fabrication technology for pulse-taking devices [1012], progress has been made on the quantification and objectification of pulse diagnosis. The physiological characteristics of the pulse at the three aforementioned palpation positions have been shown to differ, which implies that the pulse at each position conveys different clinical information [10]. According to the theory of correspondence between palpation positions and organs, the pulse at the left Gwan conveys the heart functioning [1]. Based on this theory, Huang et al. studied the characteristics of the pulse at the left Gwan and reported that the spectral energy of the pulse in the fourth to sixth harmonics was overly damped in palpitation patients compared to normal subjects [11]. On the other hand, Liu et al. found that, in the pulse measured at the left Chon, Zen meditation induces more elastic pulse waveforms which might indicate improved performance of the cardiovascular system [12]. A report asserted that heat stress reduces the radial augmentation index (AIr) and cold stress increases AIr [13]. Recently, several publications report on the technical improvement of signal processing for the pulse waveform analysis [1418].

A pulse analyzer may replace the OMDs’ pulse diagnosis by fingers if it is capable of analyzing the characteristics of the pulse in terms of fundamental physical parameters such as depth, width, length, force, rhythm, contour, speed, and rate [1, 8]. For this purpose, it requires acquisition of the pulse waveform at different hold-down pressures (equivalently, applied pressures), containing two-dimensional spatial distribution of the pulse amplitude along and across the axis of the radial artery. So far, most works on the radial pulse are limited to pattern classification and feature extractions of the pulse waveform obtained at the optimal applied pressure, aiming to distinguish abnormal pulses from normal pulses. As such, to develop desirable pulse analyzers, more extensive studies are needed.

There are proposals on how to interpret classical pulse qualities in terms of machine appropriate physical parameters [1, 9]. For instance, some researchers attempted to classify a few pulse qualities that can be identified by pattern recognition [19, 20]. Particularly, Zhang et al. developed two effective pattern classification algorithms to distinguish five different pulse patterns of moderate, smooth, taut, hollow, and unsmooth pulses.

More recently, a novel diagnostic algorithm to distinguish a deep-lying pulse (sunken pulse) from a superficial pulse (floating pulse) was proposed and validated clinically by the authors [21, 22]. For this purpose, we introduced a normalized coefficient that changes monotonically from 0 to 1; as the pulse amplitude becomes larger at heavy-applied pressures compared to light-applied pressures, the coefficient lies closer to 1 [21]. The floating pulse and sunken pulse are the two pulse qualities representing the pulse depth and they belong to the four principal pulses together with the rapid pulse and slow pulse. Therefore, development of an effective algorithm which classifies a pulse according to its depth is a major achievement.

Another principal pulse parameter is the pulse force or equivalently the pulse power. In this work, we study the pulse classification method according to its force or power. Strictly speaking, no pulse quality is defined only in terms of the force of the pulse [1]. However, the pulse force is the most crucial parameter that determines excess/deficient syndromes and is therefore of great clinical importance. Most pulse qualities that are too weak or excessively strong in its pulsation strength are indicative of the deficiency syndrome or the excess syndrome, respectively. In this study, we consider the deficient/excess pulse to be the representative of forceless/forceful pulse qualities [23]. To develop an objective and reliable classification model, firstly, we delineate samples with deficient and excess pulse qualities (DEPs) based on OMDs pulse diagnoses. By using statistical methods, such as factor analysis and Fisher’s discriminant analysis, we examine some candidate variables that contribute to the OMDs’ conclusions of the DEPs. Finally, we propose a simple but efficient classification model which best explains OMDs’ diagnostic results.

2. Quantification of the Deficient and Excess Pulse Qualities

Pulse force is a complex parameter determined by the interplay between several variables such as the amplitude of cardiac contraction, volume of blood flow, and the tensile compliance of the arterial wall. A forceful pulse is defined as having large pulse amplitude over a range of the hold-down pressures, while a forceless pulse is defined as one with small pulse amplitude (Figure 1). Forceful pulse qualities include the excess (Shi), long (Chang), flooding (Hong), tight (Jin), and wiry (Xian) pulses, and forceless qualities include the deficient (Xu), weak (Ruo), faint (Wei), scattered (San), and soft (Ru) pulses [1]. Excess pulse (Shi Mai) and deficient pulse (Xu Mai) are the representatives of the forceful pulse qualities and forceless pulse qualities, respectively.

563958.fig.001
Figure 1: Illustration of forceful (excess) qualities versus forceless (deficient) qualities in the P-H plane.

The excess pulse is felt strong at all depths from superficial to deep level, felt wide, and it is felt forceful at more than one palpation positions and the pulse is stretched beyond Cheok and/or beyond Chon positions (Figure 1). By combining with other pulse parameters, it can be distinguished from similar pulse qualities. For instance, the flooding pulse has additional properties such as floating and wide, and the wiry pulse to be less forceful and narrower than the excess and flooding pulses. The excess pulse occurs when excess perverse heat is accumulated in the three heaters in the body. Clinically, probable symptoms include insanity, mania, qi pain, yang toxins, vomiting, and other similar symptoms or it may indicate simple accumulation of perverse yang. Depending on the palpation positions, it is likely to indicate food accumulation (Gwan), constipation (Cheok) due to bound heat in the stomach (Gwan) and the intestines (Cheok), and headache, fever, sore throat, stiffness at the root of tongue, or stiffness in the chest and diaphragm (Chon) [2].

On the other hand, the deficient pulse lacks pulsation intensity. It is felt either weak through the entire range of pressure, or it is easily perceived with light pressure and ceases to be felt with heavy pressure due to arterial occlusion under heavy hold-down pressure (Figure 1). It can be distinguished from other forceless pulse qualities if other pulse parameters are additionally considered. For instance, the weak pulse is felt at a deeper level. The deficient pulse usually indicates deficiency in both qi and blood. Likely clinical symptoms include lethargy, shortness of breath, spontaneous sweating, pale complexion, low voice, dizziness, and pale tongue [2].

3. Subjects and Methods

3.1. Study Subjects

The study was approved by the ethics committees of the Korea Institute of Oriental Medicine, and informed written consent for the study was obtained from all subjects prior to study entry (I0903-01-02). Out of hundreds of healthy volunteers in their 20s with no vascular deformity on the radial artery, one hundred subjects were chosen by pairs of OMDs as appropriate candidates for diagnosing the pulse either with deficient or excess pulse qualities. The basic physiological data of the subjects are summarized as number or as mean ± SD in Table 1.

tab1
Table 1: Basic physiological data of the subjects.
3.2. Study Design

The study design is shown in Figure 2. On each day of study, according to a given schedule, different pairs of OMDs from an OMD pool with eleven practitioners with more than five years of clinical experience were assigned to make pulse diagnosis. The paired OMDs were asked to independently examine each subject. To avoid confusion when characteristic pulse feelings were different between the left pulse and the right pulse, only the left wrist of each subject was used for pulse diagnosis. Of the 100 subjects chosen by the paired OMDs as appropriate candidates for pulse diagnosis for the DEPs, diagnoses on 70 subjects were concordant between the OMD pairs. Using a pulse-taking device, the pulse waveforms were obtained at the three palpation positions of Chon, Gwan, and Cheok in the left arms. We analyzed the pulse waveforms in the 70 subjects diagnosed either with the deficient pulse or with the excess pulse and attempted to develop a classification model for the DEPs which best explains OMDs’ diagnoses.

563958.fig.002
Figure 2: Study design and the outcomes of OMDs' pulse diagnoses.
3.3. Pulse Waveform Acquisition

Pulse waveform was obtained by 3D MAC (Daeyomedi Co., Korea) which was commercially available and approved by the Korea Food and Drug Administration (KFDA). The 3D MAC operates with the applanation tonometry method to apply pressure and acquire pulse waveforms at the traditional palpation positions of Chon, Gwan, and Cheok (Figure 3). The device uses a motor-actuated pressure sensor, which contains 5 sensing elements arrayed crosswise within 10×10 mm2. Each sensing element is a piezo-resistive sensor of size about 2×3 mm2. Figure 3(a) illustrates a pulse-taking operation, with its operation procedure outlined in Figure 3(b). As shown in Figure 3(b), after an operator places the sensor at the proximity of Gwan, by an automated algorithm, the device fine-tunes the sensor location and measures pulse waveforms with varying applied pressures at five discrete pressure steps, for example, at 𝑃1, 𝑃2, 𝑃3, 𝑃4, and 𝑃5. Each pressure step is maintained constant for five seconds. After measuring the pulse waveform at Gwan, the sensor moves towards Cheok and Chon to repeat the measurement. In this study, the hold-down pressure at each pressure step was maintained at 𝑃1=37±4, 𝑃2=73±5, 𝑃3=109±5, 𝑃4=143±7, and 𝑃5=184±7 mmHg on average (±standard deviation) and the device estimated that Cheok and Chon were located about 10 mm away from Gwan along the radial artery which were comparable to the average palpation positions by OMDs for pulse diagnosis (Figure 3) [6]. Pulse pressures measured by the 3D MAC were in the acceptable range of repeatability (within about 11% of coefficient of variation).

fig3
Figure 3: (a) Illustration of a pulse-taking operation by 3D MAC and (b) outline of the automated pulse-taking procedure by 3D MAC.
3.4. Signal Processing

Figure 4 illustrates data processing towards the pulse classification algorithm. Raw data (top panel) contained noise due to breathing, uncontrolled movement of subject’s arm, and so forth. Therefore, it required preprocessing to remove noise and to align baseline followed by period segmentation and averaging (second panel). We used a nearest neighbor interpolation technique to remove abrupt signal variation, and the 5th order polynomial approximation and subsequent spline interpolation to remove baseline wander. Finally, as outlined in the bottom panel in Figure 4, for each applied pressure step 𝑃𝑗 (𝑗=1, 2, 3, 4, and 5), we calculated the maximum amplitude of the pulse waveform 𝐻𝑗 which we call the pulse amplitude at given pressure step 𝑃𝑗, which would be used to find the classification method.

fig4
Figure 4: Outline of data manipulation from raw data (top panel) to feature extraction (bottom panel).
3.5. Pulse Amplitudes

In previous subsections, we discussed discrete pressure steps and pulse amplitudes (𝑃𝑗, 𝐻𝑗) without distinguishing different palpation positions. As will be discussed below, however, to determine the DEPs, OMDs use all the pulse amplitudes at Chon, Gwan, and Cheok. To differentiate between different palpation positions, as in Figure 5, we introduce location label so that 𝐻𝑖𝑗 indicates the pulse amplitude at 𝑖th palpation position (𝑖=1, 2, and 3 for Chon, Gwan, and Cheok, resp.) and at 𝑗th pressure step (𝑗=1, 2, 3, 4, and 5 from the lightest to the heaviest applied pressure).

fig5
Figure 5: An example of pulse amplitude (𝐻) versus applied pressure (𝑃) at Chon, Gwan, and Cheok. To distinguish pulse amplitudes at different palpation positions, locational label was newly introduced to indicate that 𝐻𝑖𝑗 is the pulse amplitude at 𝑖th (𝑖=1, 2, and 3) palpation position and at 𝑗th (𝑗=1, 2, 3, 4, and 5) pressure step.

There are some pulse quantities derived from the pulse amplitudes that are potentially relevant in determining the DEPs. The first such quantity is the pulse pressure (PP). The PP, which is known to be an important indicator in predicting coronary heart disease particularly in the middle-aged and the elderly [24], is defined as the difference between systolic blood pressure and diastolic blood pressure in a cardiac cycle. The PP is equivalent to the maximum amplitude among the pulse amplitudes at various applied pressure steps, that is, PP𝑖=𝐻𝑖maxmax(𝐻𝑖1,𝐻𝑖2,𝐻𝑖3,𝐻𝑖4,𝐻𝑖5), where max() returns the largest value among () [25]. The average pulse pressure over Chon, Gwan, and Cheok is then given by PP𝑖=𝐻𝑖max(𝐻1max+𝐻2max+𝐻3max)/3. There may be other pulse quantities to be used relevantly in determining the DEPs. Such a quantity may be the maximum pulse pressure among the three palpation positions, which is defined by PPmax=max(𝐻1max,𝐻2max,𝐻3max). The mean pulse amplitude (MPA) can also contribute in a major way in the determination of the DEPs, MPA𝑖=𝐻avg𝑖(𝐻𝑖1+𝐻𝑖2+𝐻𝑖3+𝐻𝑖4+𝐻𝑖5)/5. Over the three palpation positions, we define the average and maximum of the MPAi by MPA𝑖=(𝐻avg1+𝐻avg2+𝐻avg3)/3 and MPAmax=max(𝐻avg1,𝐻avg2,𝐻avg3), respectively.

3.6. Statistical Method

Statistical analyses were performed by using SPSS version 14.0 (SPSS Inc., USA) and MATLAB version 7. x (Mathworks Inc., USA). Student’s 𝑡-test was performed to compare means of selected continuous variables in the deficient pulse group and excess pulse group. Factor analysis was used to identify groups of highly correlated variables and their relationship to the target variable. We used Fisher’s discriminant analysis to determine the best concordance with the diagnostic results of OMDs. To examine the quality of concordance between algorithmic predictions and OMDs’ diagnoses, we additionally calculated the Matthews correlation coefficient.

4. Results and Discussion

4.1. Diagnoses by Paired OMDs

One hundred subjects were selected by paired OMDs to be included in this study. Among them, the diagnoses on 70 subjects (70%) were concordant between paired OMDs, while the diagnoses of the remaining 30 subjects (30%) were divergent. Among the concordant cases, 26 subjects (37%) were diagnosed with deficient pulses and 44 subjects (63%) with excess pulses.

Table 2 shows that the accuracy (alternatively, concordant diagnoses) between OMD1 and OMD2 was 70% (70 agreements/100 simultaneous diagnoses), and the Matthews correlation coefficient (MCC) was 0.38. The MCC is regarded as one of the best measures of the quality of binary classifications, particularly when the two classes are of very different sizes [26, 27]. An accuracy of about 70% and MCC of about 0.4 are indicative of moderate concordance between OMDs’ diagnoses, noting that the diagnoses of Table 2 were not made by any fixed pairs of OMDs, but by cyclically paired OMDs among a pool of 11 OMDs on each day of study.

tab2
Table 2: Concordance of the pulse diagnoses between paired OMDs.
4.2. Characteristics of the Deficient Pulse Group and Excess Pulse Group

In total, 70 subjects were concordantly diagnosed with deficient or excess pulses by paired OMDs. Figure 6 summarizes the means and standard deviations (SD) of some relevant physiological quantities for each pulse group, stratified by gender. In addition, a Student’s 𝑡-test was applied between the two pulse groups. For the entire cohort, BMI, systolic (BPsystole) average blood pressure BP, the difference between systolic and diastolic blood pressure (ΔBP), and 𝐻𝑖max were significantly different between the two pulse groups at a significance level of 0.05, while diastolic blood pressure and heart rate remained nonsignificantly different (not shown in the figure). Here, as introduced in the previous section, 𝐻𝑖max is the approximate pulse pressure at ith palpation position measured by the pulse-taking device (3-D MAC) and 𝐻𝑖max is the average over Chon, Gwan, and Cheok in the left wrist [28], while ΔBP is the pulse pressure measured at the right brachial artery by a commercial sphygmomanometer (FT-750(R), Jawon Medical, Korea).

fig6
Figure 6: Characteristics of subjects diagnosed with the DEPs stratified by gender; (a) for the entire sample size, (b) for male subjects, and (c) for female subjects. Data presented are the mean ± SD, and the 𝑃 value from Student’s two sample 𝑡-test. *𝑃<0.05, **𝑃<0.005; N.S: nonsignificant. Abbreviated: “BPsystolic” = systolic blood pressure, “BP” = average blood pressure, “ΔBP” = BPsystolicBPdiastolic, and “𝐻𝑖max” stands for the maximum pulse amplitude at the 𝑖th palpation position whose unit was determined by the manufacturer (arb), and “𝐻𝑖max” is the average of 𝐻𝑖max over 𝑖=1, 2, and 3.

Among the five statistically significant quantities, 𝐻𝑖max were the most significantly different between the two pulse groups, which implies that, among the compared quantities, the average pulse pressure over the three palpation positions is the most appropriate choice as the decision parameter for the DEPs. The reduced significance in all variables in each gender is mostly due to reduced sample size; 𝐻𝑖max remained significantly different by gender, while such statistical differences were diminished in other variables.

The pulse pressure is known to be different throughout the large artery tree, while the mean arterial pressure remains constant and the diastolic pressure does not change substantially throughout the large artery tree [25]. In addition, the pulse pressures and other pulse parameters at the right and left arms usually show distinctive characters [29]. Therefore, it is reasonable to have discrepancy in the significance level between the pulse pressures at the left radial artery (𝐻𝑖max) and at the right brachial artery (ΔBP). On the other hand, it has been shown that the BMI is marginally valid in distinguishing the deficient pulse group from the excess pulse group. It implies that an obese individual is more likely to have an excess pulse than a thin individual, which is in agreement with clinical experience.

Since OMDs palpated the radial artery to diagnose the DEPs, among physiological quantities found at various locations along the large artery tree, a properly defined quantity on the radial artery is expected to show the best correlation with OMDs’ diagnostic result. Therefore, it is intuitively correct that the pulse pressure measured at the radial artery (𝐻𝑖max) was the most appropriate parameter in distinguishing the deficient pulse group from the excess pulse group. More importantly, it implies that the pulse pressure 𝐻𝑖max or similarly defined pulse quantities on the radial artery contribute significantly to the determination of the DEPs. In the following, to determine which pulse quantities contribute significantly to the OMDs’ diagnoses of the DEPs, we perform factor analysis accompanied by Fisher’s discriminant analysis.

4.3. Factor Analysis

As described previously, we observed that the average pulse pressure 𝐻𝑖max is one of the most promising candidates as the decision variable for the DEPs. Since 𝐻𝑖max is derived from the pulse amplitudes 𝐻𝑖𝑗, it is desirable to scrutinize 𝐻𝑖𝑗 to determine which of them are important and in which combination it best explains the decision rule of OMDs’ diagnoses of the DEPs. For this purpose, we performed factor analysis of 15 pulse amplitudes from 𝐻11 to 𝐻35. Factor analysis attempts to describe the covariance relationship among many variables in terms of a few underlying latent quantities, called factors [30]. It reduces attribute space from a larger number of variables to a smaller number of factors. To find factors, we followed the principal component method. For the rotation of variables, we followed the varimax procedure. Kaiser-Meyer-Olkin (KMO) measure of sampling adequacy was 0.672, and Bartlett’s test of sphericity was valid with 𝑃<0.01 of significance level, both of which indicate the appropriateness of factor analysis. By factor analysis, we obtained 5 significant factors with eigenvalues larger than one, which accounted for about 80% of the variance as in Table 3. Note that an eigenvalue represents the amount of variance associated with the factor, and therefore only factors with variance greater than one are expected to show better performance than an initial individual variable, that is, 𝐻𝑖𝑗, whose variance is normalized to one [30].

tab3
Table 3: Factor analysis of the 15 pulse amplitudes.

As in Table 3, we determined the 5 most relevant factors which accounts for about 80% of the variance, in which all 15 pulse amplitudes contribute once and only once with likely weight; factor loadings are in equal orders of magnitude, ranging between 0.58 and 0.93 [31]. The two most contributing factors account for about 54% of the variance, in which the pulse amplitudes at all applied pressures at Cheok (𝐻31 to 𝐻35) and the pulse amplitudes at light-applied pressures at Gwan (𝐻21 and 𝐻22) are almost equally involved with similar factor loadings. In summary, all 15 pulse amplitudes were found to be important in explaining the total variance, and the pulse amplitudes with light-applied pressures and heavy-applied pressures were grouped into different factors, implying possibly different roles in pulse classifications between the two pulse amplitude groups with light-applied pressures and heavy-applied pressures.

4.4. Fisher’s Discriminant Analysis with the Five Factors Obtained from Factor Analysis

With the five factors listed in Table 3, we continued to perform Fisher’s discriminant analysis to determine a discriminant function for the DEPs with reference to the OMDs’ diagnoses [30]. The value of Box’s 𝑀 test was 12.29 (𝑃 value = 0.737), which satisfies the assumption of equal variance. The discriminant function was found to be significant with a Wilks’ Lambda of 0.793 (𝑃 value = 0.010). The standardized canonical discriminant function coefficients for the 5 factors are listed in Table 4, and the classification results are shown in Table 5. The magnitude of a coefficient indicates the contribution weight of the given factor. Table 4 shows that factors 𝑓1, 𝑓4, and 𝑓5 contribute almost equally, and the contribution of 𝑓2 and 𝑓3 are, respectively, about a half and a fifth of the major factors which are in the same order of magnitude. It is worth mentioning that no factor alone governs the behavior of the discriminant function but all the participating factors contribute evenly.

tab4
Table 4: The standardized canonical discriminant function coefficients for the five factors in Table 3.
tab5
Table 5: Pulse classification with the 5 factors obtained from factor analysis into the DEPs.

By combining the results in Tables 3 and 4, we notice that all the pulse amplitudes at light- and heavy-applied pressures at Chon, Gwan, and Cheok contribute on equal orders of magnitude to the classification of the DEPs. More rigorously speaking, the most contributing factors of 𝑓1, 𝑓4, and 𝑓5 are comprised of the pulse amplitudes at light-applied pressures at Chon (𝐻11 and 𝐻12) and the pulse amplitudes at heavy-applied pressures at Gwan (𝐻23, 𝐻24, and 𝐻25) and Cheok (𝐻33, 𝐻34, and 𝐻35). It is clinically known that the optimal pulse depth at which the pulse is felt strongest is shallower at Chon than at Gwan or Cheok. Applying this clinical knowledge, the pulse pressure and the factors of 𝑓1, 𝑓4, and 𝑓5 are well correlated, which implies that the pulse pressure may be an appropriate quantity for the classification of the DEPs, which is one of our target variable to be discussed below.

We applied the coefficients in Table 4 to Fisher’s discriminant function and obtained the classification result in Table 5. The accuracy of the classification for the entire data set is 72.9% with the Matthews correlation coefficient of 0.46. We repeated the leave-one-out cross-validation test and obtained the accuracy of 61.4% (MCC=0.24). Reduced accuracy of about 11% in the cross-validation test indicates that the generated discriminant function overfits the training set, which opens the possibility for more efficient classification with less variables.

4.5. Fisher’s Discriminant Analysis with the Representative Pulse Quantities at Each Palpation Position

By factor analysis, we found that all the pulse amplitudes at various levels of applied pressures at the three palpation positions contributed with similar weight in the determination of the DEPs. OMDs rely mostly on the pulse force to determine the DEPs. A pulse may be considered forceful if either its maximum amplitude is large or the average amplitude over various applied pressures is large. The former is the pulse pressure and the latter is the mean pulse amplitude. In search of a simple form for the discriminant function with improved accuracy compared to the result using factor analysis, the two most relevant quantities are thought to be the pulse pressure PP𝑖=𝐻𝑖max and the mean pulse amplitude MPA𝑖=𝐻avg𝑖 at each (𝑖th) palpation position, whose detailed expressions in terms of 𝐻𝑖𝑗 are introduced in Section 3.5. With the pulse pressures and the mean pulse amplitudes at the three palpation positions, in the following, we apply Fisher’s discriminant analysis to determine an efficient classification model.

In the canonical correlation analysis, the coefficients of the standardized canonical discriminant function for 𝐻1max, 𝐻2max, and 𝐻3max were, respectively, given by 0.505, 0.282, and 0.555, which indicates that the three pulse pressures contribute to the discriminant function on equal orders of magnitude. The classification result by the linear discriminant function using the three pulse pressures is summarized in Table 6. The accuracy of the classification for the entire data set is 78.6%, and the Matthews correlation coefficient is 0.56. By a cross-validation test, we obtained the accuracy of 75.7% (MCC=0.51).

tab6
Table 6: Pulse classification result by 𝐻1𝑚𝑎𝑥, 𝐻2𝑚𝑎𝑥, and 𝐻3𝑚𝑎𝑥 into the DEPs.

By repeating the canonical correlation analysis using 𝐻avg1, 𝐻avg2, and 𝐻avg3, we again find that the 3 involved variables contribute to the discriminant function on equal orders of magnitude (the standardized canonical discriminant function coefficients are given by 0.248, 0.533, and 0.459 resp. for 𝐻avg1, 𝐻avg2, and 𝐻avg3). The accuracy of the classification for the entire data set was 75.7% (MCC=0.51), and the cross-validated classification accuracy was 68.6% (MCC=0.36). Using both the maxima and averaged variables together, we obtained a classification accuracy of 72.9% (MCC=0.45) for the entire data set and 65.7% (MCC=0.30) by the cross-validation test, which is worse than the three pulse pressures or the three averaged variables separately. In conclusion, the three pulse pressures 𝐻1max, 𝐻2max, and 𝐻3max were relevantly contributing to the pulse decisions for the DEPs, while additional or independent participation of the mean pulse amplitudes did not improve the concordance rate with OMDs’ pulse decisions.

4.6. Diagnostic Model with the Representative Pulse Quantities over All the Palpation Positions

Let us further reduce the number of pulse variables by taking the maximum or the average of pulse quantities 𝐻𝑖max and 𝐻avg𝑖 over the three palpation positions. This reduction procedure is based on the following reason. When a pulse is considered forceful over the three palpation positions, it may indicate that the maximum of a pulse quantity over the three palpation positions is large or that the average over the palpation positions is large.

Following this procedure, we obtain four pulse quantities, such as PP𝑖, PP𝑖max, MPA𝑖, and MPA𝑖max (see Section 3.5 for detailed relations with 𝐻𝑖𝑗’s), as competing variables for the pulse classification. By repeating Fisher’s discriminant analysis, as summarized in Table 7, we obtained good accuracy and Matthews correlation coefficient with the average pulse pressure (PP𝑖) or with the average of the MPAs (MPA𝑖), while less accurate classification was obtained with the maximum pulse pressure (PP𝑖max) or with the maximum of the MPAs. The combination of the four variables does not increase the accuracy of the discrimination, but rather reduces it due to overfitting.

tab7
Table 7: Fisher's discriminant analysis with 𝑃𝑃𝑖, 𝑃𝑃𝑖𝑚𝑎𝑥, 𝑀𝑃𝐴𝑖, and 𝑀𝑃𝐴𝑖𝑚𝑎𝑥.

We can improve the accuracy beyond the limit of the linear discriminant analysis by a mixed-variable classification model, where two competing variables participate sequentially in the decision process in complementary manner. To build an appropriate mixed-variable classification model, we review a simplest linear discriminant function with one participating variable; with PP𝑖, the rule is given by IfPP𝑖𝛼,thenitisanexcesspulsequality,whileifPP𝑖,<𝛽,thenitisadecientpulsequality(1) where parameters 𝛼 and 𝛽 are the criteria used for decision on the DEPs and they satisfy 𝛼𝛽. When 𝛼=𝛽, the accuracy and the MCC are as summarized in Table 7. In this simple discriminant analysis, there exists an intermediate regime between or near the criterion parameters 𝛼 and 𝛽, where the concordance of the pulse decisions with the OMDs’ decisions becomes abruptly poor; in the decision making for PP𝑖 far above 𝛼 as an excess pulse or for PP𝑖 far below 𝛽 as a deficient pulse, the concordance is usually very good.

An improved accuracy is possible by substituting the intermediate regime of discordance (the regime between 𝛼 and 𝛽) with another candidate variable such as PP𝑖max or MPA𝑖max. Despite the correlation between these variables is high (Pearson correlation coefficient ≥ 0.8), the intermediate regimes of discordance between themselves are slightly off-resonance. It gives an opportunity to improve accuracy by letting the secondary variable intervene in the intermediate regime of discordance of the primary variable. The standard classification procedure for this mixed-variable diagnostic model is outlined in Figure 7. It yielded an accuracy of 80.0% and MCC of 0.57 (after standardization, the criterion parameters were, resp., given by 𝛼=0.216, 𝛽=0.572, and 𝛾=𝛿=0.314), which was about 5.7% increase in the accuracy and 10% increase in the MCC compared to the original discriminant function with PP𝑖 alone. Applying MPA𝑖 and PP𝑖max as the two variables in the mixed-variable diagnostic model, it yielded similar improvement in performance (accuracy of 80%). The increased accuracy of the mixed-variable discriminant model implies that the average and maximum pulse amplitudes over the palpation positions are complementary to each other in OMDs’ pulse decisions.

fig7
Figure 7: A mixed-variable diagnostic method for the DEPs. (a) flowchart for the classification, and (b) classification regimes for the DEPs.

5. Conclusions

For the objectification and standardization of pulse diagnosis, reliable classification methods for the principal pulse qualities are urged to be developed. The deficient and excess pulse qualities (DEPs) are clinically important as they are the indicators representing the deficiency syndrome and the excess syndrome, respectively. In this work, we proposed a classification method for the DEPs. For this purpose, we conducted a clinical test and selected 70 subjects in their 20s either with the deficient pulse (26 samples) or with the excess pulse (44 samples), by concordant diagnoses between paired OMDs. By a Student’s 𝑡-test, among several candidate quantities such as the BMI or systolic/diastolic blood pressures measured at the brachial artery, we confirmed that the average pulse pressure defined over the three palpation positions was the most appropriate quantity in distinguishing the excess pulse group from the deficient pulse group. We continued to apply factor analysis and Fisher’s discriminant analysis and found that all the pulse amplitudes obtained at various applied pressures at the three palpation positions contributed relevantly to the diagnosis of the DEPs.

Next, we showed that either of the pulse pressure or the average pulse amplitude yielded as good accuracy as the original pulse amplitudes. It reflects that the diagnoses of the DEPs by OMDs rely mostly on the pulse force, as either of the two quantities is appropriate in representing the pulse force and good at describing the collective behavior of the original pulse amplitudes. Finally, we proposed a mixed-variable classification model, in which two complementary variables, for example, either two of the maximum or average of the pulse pressures, or the maximum or average of the mean pulse amplitudes, acted over the three palpation positions, were used sequentially to increase the classification accuracy in a reasonable degree. This study will contribute to the objectification and standardization of pulse diagnosis.

Acknowledgments

This work was supported by a grant from the Korea Institute of Oriental Medicine (KIOM) funded by the Korea government (K11070) and was partially supported by the Korea Ministry of Knowledge Economy (10028438). The authors thank Mr. Boncho Ku for his contribution to the statistical analysis and Dr. Young-Koo Lee for comments on the paper. The authors declare no commercial associations that may lead to conflicts of interests in connection with this work.

References

  1. Z. F. Fei, Contemporary Sphygmology in Traditional Chinese Medicine, People's Medical Publishing House, Beijing, China, 2003.
  2. S. Walsh and E. King, Pulse Diagnosis: A Clinical Guide, Elsevier, Sydney, Australia, 2008.
  3. L. Hammer, Chinese Pulse Diagnosis: A Contemporary Approach, Eastland Press, Seattle, Wash, USA, 2001.
  4. H. K. Huynh and G. M. Seifert, Pulse Diagnosis by Li Shi Zhen, Paradigm Publications, Sydney, Australia, 1985.
  5. B. Flaws, The Secret of Chinese Pulse Diagnosis, Blue Poppy Press, Denver, Colo, USA, 2006.
  6. H. H. Kim, J. Lee, K. W. Kim, and J. Y. Kim, “Proposal for pulse diagnosis positions(chon-kwan-chuk) for pulse analyzer based on literature review and anthropometry,” Korean Journal of Oriental Medicine, vol. 28, no. 3, pp. 13–22, 2007.
  7. Y. J. Lee, J. Lee, and J. Y. Kim, “Suggestion on an innovative pulse diagnosis system based on technical trend analysis,” Korean Journal of Oriental Physiology & Pathology, vol. 23, no. 1, pp. 174–179, 2009.
  8. J. Y. Kim, K. Y. Kim, and K. D. Ko, “A study on the problems and the method for improvement of pulse analyzers,” The Journal of the Korea Institute of Oriental Medical Diagnostics, vol. 3, no. 1, pp. 28–36, 1999.
  9. H. H. Ryu, S. W. Lee, J. Lee, Y. J. Lee, and J. Y. Kim, “Analysis of physical quantification of pulse types by pulse diagnosis literatures,” Korean Journal of Oriental Physiology & Pathology, vol. 21, no. 6, pp. 1381–1387, 2007.
  10. Y. J. Jeon, J. U. Kim, H. J. Lee, et al., “A clinical study of the pulse vave characteristics and the three pulse diagnosis positions of Chon, Gwan, and Cheok,” Evidence-Based Complementary and Alternative Medicine, vol. 2011, Article ID 904056, 9 pages, 2011. View at Publisher · View at Google Scholar · View at PubMed
  11. C.-M. Huang, C.-C. Wei, Y.-T. Liao, H.-C. Chang, S.-T. Kao, and T.-C. Li, “Developing the effective method of spectral harmonic energy ratio to analyze the arterial pulse spectrum,” Evidence-Based Complementary and Alternative Medicine, vol. 2011, Article ID 342462, 7 pages, 2011.
  12. C.-Y. Liu, C.-C. Wei, and P.-C. Lo, “Variation analysis of sphygmogram to assess cardiovascular system under meditation,” Evidence-Based Complementary and Alternative Medicine, vol. 6, no. 1, pp. 107–112, 2009. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  13. C.-M. Huang, H.-C. Chang, S.-T. Kao, et al., “Radial pressure pulse and heart rate variability in heat- and cold-stressed humans,” Evidence-Based Complementary and Alternative Medicine, vol. 2011, Article ID 751317, 9 pages, 2011. View at Publisher · View at Google Scholar · View at PubMed
  14. L. Xu, D. Zhang, and K. Wang, “Wavelet-based cascaded adaptive filter for removing baseline drift in pulse waveforms,” IEEE Transactions on Biomedical Engineering, vol. 52, no. 11, pp. 1973–1975, 2005. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  15. L. Xu, D. Zhang, K. Wang, and L. Wang, “Arrhythmic pulses detection using Lempel-Ziv complexity analysis,” EURASIP Journal on Applied Signal Processing, vol. 2006, Article ID 18268, pp. 1–12, 2006.
  16. P.-Y. Zhang and H.-Y. Wang, “A framework for automatic time-domain characteristics parameters extraction of human pulse signals,” EURASIP Journal on Applied Signal Processing, vol. 2008, Article ID 468390, 9 pages, 2008.
  17. Y. Chen, L. Zhang, D. Zhang, and D. Zhang, “Wrist pulse signal diagnosis using modified Gaussian models and fuzzy C-means classification,” Medical Engineering & Physics, vol. 31, no. 10, pp. 1283–1289, 2009. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  18. D.-Y. Zhang, W.-M. Zuo, D. Zhang, H.-Z. Zhang, and N.-M. Li, “Wrist blood flow signal-based computerized pulse diagnosis using spatial and spectrum features,” Journal of Biomedical Science and Engineering, vol. 3, no. 4, pp. 361–366, 2010.
  19. J.-J. Shu and Y. Sun, “Developing classification indices for Chinese pulse diagnosis,” Complementary Therapies in Medicine, vol. 15, no. 3, pp. 190–198, 2007. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  20. D. Zhang, W. Zuo, D. Zhang, H. Zhang, and N. Li, “Classification of pulse waveforms using edit distance with real penalty,” EURASIP Journal on Advances in Signal Processing, vol. 2010, Article ID 303140, 8 pages, 2010.
  21. J. U. Kim, Y. J. Jeon, Y. J. Lee, K. H. Kim, and J. Y. Kim, “Novel diagnostic algorithm for the floating and sunken pulse qualities and its clinical test,” Evidence-Based Complementary and Alternative Medicine, vol. 2011, Article ID 813427, 10 pages, 2011.
  22. S. H. Kim, J. U. Kim, Y. J. Lee, K. H. Kim, and J. Y. Kim, “New algorithm of determining the floating and sinking pulse with a pulse diagnosis instrument,” Korean Journal of Oriental Physiology & Pathology, vol. 3, no. 6, pp. 1221–1225, 2009.
  23. S. H. Kim, J. U. Kim, Y. J. Jeon, K. H. Kim, and J. Y. Kim, “Method for determining the deficient and solid pulse with a new pulse wave parameter,” Korean Journal of Oriental Physiology & Pathology, vol. 24, no. 1, pp. 42–47, 2010.
  24. S. S. Franklin, S. A. Khan, N. D. Wong, M. G. Larson, and D. Levy, “Is pulse pressure useful in predicting risk for coronary heart disease? The framingham heart study,” Circulation, vol. 100, no. 4, pp. 354–360, 1999. View at Scopus
  25. W. W. Nichols and M. F. O'Rourke, McDonald's Blood Flow in Arteries. Theoretical, Experimental and Clinical Principles, Hodder Arnold, London, UK, 5th edition, 2005.
  26. B. W. Matthews, “Comparison of the predicted and observed secondary structure of T4 phage lysozyme,” Biochimica et Biophysica Acta, vol. 405, no. 2, pp. 442–451, 1975. View at Scopus
  27. P. Baldi, S. Brunak, Y. Chauvin, C. A. F. Andersen, and H. Nielsen, “Assessing the accuracy of prediction algorithms for classification: an overview,” Bioinformatics, vol. 16, no. 5, pp. 412–424, 2000. View at Scopus
  28. The unit of pulse amplitudes Hijmax=100 (and the pulse pressure Himax) was predetermined by the manufacturer and was not known to the authors. The correspondence for, e.g., Hijmax=100, in millimeters of mercury estimated by us was 20 mmHg.
  29. Y. Lee, J. Lee, H. Lee, and J. Kim, “Study for correlation characteristics on radial artery and floating/sinking pulse with BMI,” Korean Journal of Oriental Medicine, vol. 14, no. 3, pp. 121–126, 2008.
  30. J. F. Hair, R. E. Anderson, R. L. Tatham, and W. C. Black, Multivariate Data Analysis, Prentice-Hall, New Jersey, NJ, USA, 5th edition, 1998.
  31. Each factor contains all original variables of H11 to H35, with their factor loadings ranging from −0.185 to +0.928. Since the behavior of a factor is dominated by variables with large factor loading, in Table 3, we presented only the leading components with factor loadings larger than 0.5.