Evidence-Based Complementary and Alternative MedicineVolume 2011 (2011), Article ID 813427, 10 pagesdoi:10.1155/2011/813427

Research Article

Novel Diagnostic Algorithm for the Floating and Sunken Pulse Qualities and Its Clinical Test

Jaeuk U. Kim, Young Ju Jeon, Yu Jung Lee, Keun Ho Kim, and Jong Yeol Kim

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

Received 14 September 2010; Accepted 15 January 2011

Copyright © 2011 Jaeuk U. Kim et al. 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.

We propose a novel classification algorithm for the floating pulse and the sunken pulse using a newly defined coefficient (). To examine the validity of the proposed algorithm, we carried out a clinical test in which 12 oriental medical doctors made pairwise diagnoses on the pulses of volunteering subjects. 169 subjects were simultaneously diagnosed by paired doctors, and the diagnoses in 121 subjects were concordant, yielding an accuracy of 72% and a Matthews correlation coefficient of 0.42, which indicates reasonable agreement between doctors. Two sample T-tests showed that subjects in the sunken pulse group had significantly higher BMI and () than those in the floating pulse group. The pulse classification by the algorithm converged with the diagnoses of paired doctors with an accuracy up to 69%. With these results, we confirm the validity of the novel classification algorithm for the floating and sunken pulses.

1. Introduction

Pulse diagnosis is one of the four important diagnostic methods in oriental medicine, which are inspection, listening and smelling, inquiring, and palpation [1–4]. The primary site for pulse diagnosis is the radial artery placed on both wrists. To make a pulse diagnosis, oriental medical doctors (OMDs) divide the terminal region of the radial artery into three adjacent intervals of Chon, Gwan, and Cheok, and use index, middle, and ring fingers on each palpating position. By pressing with these fingers with varying pressure simultaneously or sequentially, doctors determine characteristic pulse qualities. By analyzing the pulse qualities, primarily, doctors obtain information on balance in health and illness or subtle changes in homeostasis [3].

Pulse diagnosis requires sophisticated palpation techniques and interpretation on the pulse qualities demands enduring efforts with many clinical trials and case studies. For this reason, experiential bias and subjective judgments remain difficult to avoid. To overcome technical difficulties, to answer doubts and skepticism on its clinical validity, and to go beyond its current state of the art, it is necessary to develop objective techniques for pulse-taking using standardized protocols and quantitative definitions of the pulse qualities in terms of a few measureable physical parameters [5–7].

To quantify and standardize pulse diagnosis, it is essential to develop palpation devices that provide objective pulse signals as a function of the hold-down pressure, and to develop classification algorithms of the pulse qualities by analyzing the measured pulses. At the same time, the validity and efficiency of pulse diagnosis need to be verified clinically by objective methods. Thus far, some progress has been made on the quantification and standardization of pulse diagnosis. For example, pulse devices equipped with arrays of piezoresistive pressure sensors have been developed. In addition, clinical data have been collected by pulse devices and it has been analyzed by statistical methods. Studies on the fundamental aspects such as the blood flow dynamics along the radial artery and classifications of pulse waveforms are currently under investigation [8–18].

Pulse-taking devices should be able to measure pulse waveforms for several different steps of hold-down pressure, and it should be able to record the spatial distribution of the pulse strength along the axis of the radial artery and in its cross-section [5]. With a pulse-taking device which supplies the above-described quantities with guaranteed intrarater reliability (repeatability) and interrater reliability (reproducibility), we can obtain some principal characteristics of the pulse such as the level of the depth, width, length, and strength of the pulse, as well as the speed or rate of the pulse. As there are some literature studies which reinterpreted the pulse qualities in the classic textbooks in terms of these measurable quantities [6, 7], a desirable direction to develop classification algorithms for the pulse qualities is to understand a pulse by simple physical quantities.

For developing a reliable classification algorithm for the pulse qualities, the algorithm needs to be validated by a clinical test with a relevant protocol. The pulse-taking device used should have guaranteed repeatability and reproducibility, and the pulse qualities determined by the algorithm should be tested for consistency with the OMDs’ diagnoses. As a first step towards developing reliable classification algorithms, in this work, we introduce a novel classification algorithm for the floating and sunken pulses. To study the validity of the suggested algorithm, we study clinically the correlation between algorithmic outcomes and the OMDs’ diagnoses. The floating and sunken pulses belong to four principal pulses together with the rapid and slow pulses [1]. Therefore, one cannot emphasize enough the importance of developing a reliable classification algorithm for the floating and sunken pulses.

2. Quantification of the Floating and Sunken Pulses

2.1. Floating Pulse and Sunken Pulse

According to “Mai Jing”, the floating pulse is a pulse potent when felt with no pressure applied but impotent when felt with pressure applied, and the sunken pulse is a pulse impotent when felt with no pressure applied but potent when felt with pressure applied, while the hidden pulse is a pulse imperceptible until the fingers touch the bone with extremely heavy pressure [3]. Similar definitions can be found in other classics such as “Bin Hue Mai Xue” or “Yi Xue Ru Men”. A brief hypothetical explanation for the characteristic features of the floating and sunken pulses is as follows. The floating pulse is characterized as a Yang pulse and it is mostly seen in wind evil, external guest, and exterior patterns. The sunken pulse, on the other hand, is characterized as a Yin pulse and it is largely seen in those who have pathologic Qi hidden internally or interior patterns [2]. In general, superficial pulses are associated with acute diseases of exterior origin, while deep pulses are indicative of more severe and chronic illnesses related to the internal organs. For instance, the sunken pulse without other symptoms indicates a high potential for the impaired immune system. It is worth noting that the pulses of obese subjects are normally deeper than those of thin individuals [4].

2.2. Coefficient of the Floating and Sunken Pulses

There have been suggestions on how to quantify the pulse qualities defined only by the level of depth using pulse signals obtained by pulse-taking devices. As expounded by Fei [6], the floating/sunken/hidden pulses are distinguishable from the shape of the P-H curve (the curve of the pulse strength (H) as a function of the hold-down pressure (P)) as shown in Figure 1. Indeed, there are some pulse qualities other than the floating or sunken pulses that appear either at the superficial level or deep level, and they can be categorized as floating-type pulses or sunken-type pulses [1]. For instance, the scallion stalk pulse is a pulse exhibiting a floating characteristic combined with large but soft features (empty in the middle but solid at the sides when pressure is applied) [3]. In this aspect, strictly speaking, the P-H curve alone is insufficient to distinguish the floating or sunken pulses from more complex pulses. However, for the sake of simplicity, in this work, we do not distinguish the floating and sunken pulses from floating-type and sunken-type pulses and suppose them as the floating and sunken pulses. As shown in Figure 1, the hidden pulse can be considered the extreme of the sunken pulse and we will only focus on the floating and sunken pulses.

Figure 1: Floating versus sunken versus hidden pulses distinguished in the P-H curve [6].

Utilizing the idea in Figure 1, Lee et al. introduced the coefficient of the floating and sunken pulses () defined by [19] where () is the minimum (maximum) hold-down pressure in the P-H curve and is the optimal hold-down pressure at which the pulse strength becomes maximal (see Figure 2). The in (1) varies within [0,1], and the pulse is said to be floating-like if is close to 0 and sunken-like if is close to 1. In this scheme, the hidden pulse can be considered naturally the limiting case of the sunken pulse (). In a clinical test for determining the level of depth of the pulse by OMDs, Lee et al. showed that there was a significant difference in the averaged between subject groups diagnosed with floating pulse (for short, floating pulse group) and those with sunken pulse (for short, sunken pulse group), which verifies the method [19].

(a)

(b)

Figure 2: P-H curve with interpolation (solid curve) for a few discrete data points (circles) for the calculation of [19]. Examples for (a) a floating pulse, (b) an erroneous decision as a sunken pulse. Compare (b) with Figure 3(b).

The defined in (1) needs an interpolation technique since the pulse is recorded over a few discrete steps of the hold-down pressure. When the measured pulse data do not follow an ideal pattern of smooth down-hill after up-hill with increasing hold-down pressure, or when the pressure steps are too few that the data points are not sufficient for a proper interpolation, the calculated may yield an erroneous correspondence to the level of pulse depth. An example of such an erroneous decision by the method is illustrated in Figure 2(b), through which it is determined as a sunken-like pulse but an OMD would not diagnose it as a floating pulse or sunken pulse.

Recently, Kim et al. proposed a new method to determine the pulse qualities related to the level of pulse depth, by explicit comparison of pulse strengths at a low hold-down pressure () and at a high hold-down pressure () [20]. This method determines the pulse as a floating-like pulse if the pulse strength at low applied pressure is stronger than that at high hold-down pressure, and a sunken-like pulse if the pulse strength shows opposite behavior. This new proposal is more robust against anomalous data than the interpolated P-H curve method proposed by Lee et al., an example is shown in Figure 3(b).

(a)

(b)

Figure 3: New algorithm for [20]. It compares explicitly the pulse strengths at low hold-down pressure and high hold-down pressure. Example for (a) a floating pulse, (b) a middle-depth pulse (neither floating nor sunken). Compare (b) with Figure 2(b).

By updating the proposal by Kim et al., we introduce a normalization constant as the sum of the pulse strengths, and define the new as where () is the pulse strength at the pressure step of (). In Figure 3, it shows an appropriate definitions for and , given by and . Similar to the by Lee et al., the new in (2) varies within [0,1], and the pulse is floating-like if is close to 0 and sunken-like if is close to 1. For equal pulse strength of , the coefficient becomes . This method simulates OMDs who diagnose the pulse by comparing pulse strengths at two distinct pressure steps.

An appropriate assignment of pulse strengths and depends directly on the properties of the pressure-steps applied in the pulse-taking device. The pulse-taking device used in our study applied hold-down pressure in five different steps, which were , , , , and  mmHg on average (±standard deviation) over the subjects and over the palpation positions, for whom participated in this study. We found two appropriate choices for and which showed best concordance with OMD decisions for the floating and sunken pulses (detailed comparisons between different choices of and will not be addressed here); (1) () and (2) () = (, ). Let us distinguish these two scenarios by labeling their coefficients of the floating and sunken pulses as (1) for , and (2) for () = (, ). We will discuss the implications of these choices for () in Discussion and Conclusions (Section 5).

The method using the in (2) can be applied to simultaneous palpation by extending and to be averaged all over Chon, Gwan, and Cheok. For example, for (2), the pulse strengths can be extended by () = , and similarly for (1), where stands for the pulse strength at the th pressure-step at Chon for , at Gwan for , and at Cheok for . In our study, as conventionally done in clinics to diagnose the pulse, OMDs palpated simultaneously at Chon, Gwan, and Cheok of the left wrist of each subject. For this reason, we adopt the method with extended form of and over Chon, Gwan, and Cheok.

3. Subjects and Methods

3.1. Subjects

Subjects who were simultaneously pulse-diagnosed by paired OMDs were composed of 169 healthy volunteers in their 20s with no vascular deformity on the radial artery. The purpose and protocol of the study were explained to the subjects prior to study commencement. The study was carried out with the approval of the Institutional Review Board of the Korea Institute of oriental Medicine (IRB, I0903-01-02). The basic health conditions of each subject were examined by the OMDs, based on listening and smelling, inquiry, inspection, and palpation, with the help of questionnaire and blood pressure measurement. The basic physiological data of the subjects are summarized in Table 1.

Table 1: Subject information.

Characteristic (unit)
Value (mean ± SD)

Number (n)
169 (male = 88, female = 81)

Age (yr)
23.7 ± 2.4

Height (cm)
169.0 ± 8.2

Weight (kg)
63.6 ± 11.5

BMI (kg/m²)
22.1 ± 3.0

Systolic/diastolic blood pressure (mm Hg)
120.9/67.4 ± 16.2/11.2

3.2. Study Design

The overall flow of the clinical study is illustrated in Figure 4. We attempted to reduce experiential bias of individual OMD by allocating pairs of OMDs to each subject. Each pair of OMDs was cyclically chosen from an OMD pool comprised of 12 practitioners with more than 5 years of clinical experience. OMDs diagnosed the pulse by simultaneously palpating at three locations on the left wrist-Chon, Gwan and Cheok. OMDs examined the pulses only on the left wrist to reduce degrees of freedom in case that the left and right pulses differ in their characteristics. For the 121 subjects who were concordantly diagnosed as floating or sunken pulses by paired OMDs, we analyzed the measured pulses by the pulse-taking device.

Figure 4: Study design and the outcomes of OMDs' pulse diagnoses.

3.3. Pulse Measurement

The pulse was measured by 3-D MAC (Daeyomedi, Korea) approved by the Korea Food and Drug Administration. The 3-D MAC is equipped with an array-sensor with its elements arrayed crosswise within 10 × 10 mm². Each sensing element is a piezoresistive pressure sensor of about 2 × 3 mm².

A trained operator located the sensor at Gwan and initiated the measurement program. Before commencing actual pulse measurement, the device calibrated sensing performance and carried out fine adjustment to enable the central sensing element to locate at the center of the radial artery by automatically scanning nearby areas with varying hold-down pressures. It subsequently commenced measuring the pulse. The pulse was measured with applied hold-down pressures in five different steps. The applied pressure at each step was maintained constant for five seconds. For our study, the hold-down pressure at each pressure-step was set at 37  ±  6, 73  ±  10, 108  ±  13, 143  ±  17, and 181  ±  21 mmHg on average (±standard deviation). After measuring the pulse signal at Gwan, by an automated algorithm, the 3-D MAC shifted its sensor towards Cheok position and repeated measurement and subsequently it moved towards Chon position to complete the measurement. The 3-D MAC estimates that Cheok and Chon are about 10 mm away from Gwan which are comparable to the average palpation positions by OMDs for the pulse diagnosis (see Figure 5) [21]. We note that the way the 3-D MAC obtains pulse signals at Chon, Gwan, and Cheok is similar to the method of feeling the pulse with a single finger. The appearance of the 3-D MAC with its operational demo is shown in Figure 5(b). The performance of pulse measurement by the 3-D MAC is in the acceptable range of repeatability with the covariance for the maximum pulse strength at about 10%.

(a)

(b)

Figure 5: (a) palpation positions of Chon, Gwan, and Cheok, and (b) a demonstration of pulse-taking operation by the 3-D MAC.

3.4. Acquisition Process for the Pulse Strength for Each Pressure-Step

As shown in Figure 6 (top panel), raw data of the pulse waveform obtained from the measuring sensor contained noises due to breathing, uncontrolled movement of subject’s arm, and so forth. Baseline alignment was done by the 5th-order polynomial approximation, yielding the middle panel. By obtaining the average waveform for each pressure-step, we finally obtain the pulse strength ( to ) for each pressure-step ( to ) as shown in the bottom panel in Figure 6. By making appropriate choices for and , we proceeded to calculate (1) and (2).

Figure 6: Data acquisition process from raw data (top panel) to the pulse strength for each pressure-step (bottom panel).

4. Results

4.1. Diagnosis Result by Paired OMDs

Paired OMDs diagnosed 169 subjects simultaneously either as floating or sunken pulses. Among them, the diagnoses on 121 subjects (71.6%) were concordant, while the diagnoses of the remaining 48 subjects (28.4%) were divergent. Among the concordant cases, 49 subjects (40.5%) were diagnosed with floating pulses, and 72 subjects (59.5%) with sunken pulses.

Table 2 shows that the accuracy (alternatively, concordant diganoses) between OMD1 and OMD2 was 71.6% (121 agreements/169 simultaneous diagnoses), and the Matthews correlation coefficient (MCC) was 0.42. Note that Table 2 was obtained not by any fixed pairs of OMDs, but rather by cyclically paired OMDs among 12 OMD participants on each day of study. With these facts in mind, an accuracy above 70% and MCC above 0.4 are indicative of a good concordance between OMDs diagnoses.

Table 2: Consistency of the pulse diagnoses between paired OMDs.

OMD1

Floating
Sunken
Total

OMD2
Floating
49
22
71

Sunken
26
72
98

Total
75
94
169

4.2. Characteristics of the Floating Pulse Group and Sunken Pulse Group

We assumed diagnoses by paired OMDs returned reliable concordance. For the 121 subjects who were concordantly diagnosed with floating or sunken pulses by paired OMDs, we tabulate some physiological quantities stratified by pulse quality and gender in Table 3.

Table 3: Characteristics of 121 subjects stratified by pulse quality and gender.

SEX
Pulse (subject#)
BMI (kg/m²)
(mmHg)
〈Bp〉 (mmHg)
Heart rate (beats/min)

(1)
(2)

Male
Floating (26)
21.3 ± 1.6
126.0 ± 16.0
85.1 ± 11.6
80.4 ± 14.5
170.2 ± 40.1
0.46 ± 0.13
0.55 ± 0.13

Sunken (32)
24.2 ± 3.2
129.4 ± 14.1
91.2 ± 10.8
77.1 ± 12.2
160.5 ± 30.6
0.58 ± 0.10
0.66 ± 0.09

-test (-value)
1.2E-4**
0.38
0.044*
0.35
0.30
5.1E-4**
3.2E-4**

Female
Floating (23)
19.4 ± 2.0
109.0 ± 12.2
79.0 ± 7.2
83.0 ± 12.4
143.4 ± 33.5
0.43 ± 0.14
0.52 ± 0.16

Sunken (40)
22.3 ± 2.8
115.3 ± 13.0
82.9 ± 9.9
79.9 ± 11.0
136.0 ± 35.3
0.49 ± 0.16
0.61 ± 0.16

-test (-value)
6.2E-5**
0.063
0.10
0.30
0.42
0.12
0.029*

All
Floating (49)
20.4 ± 2.0
118.0 ± 16.6
82.2 ± 10.2
81.6 ± 13.5
157.6 ± 39.2
0.45 ± 0.14
0.54 ± 0.14

Sunken (72)
23.1 ± 3.1
121.6 ± 15.2
86.6 ± 11.1
78.6 ± 11.5
146.9 ± 35.3
0.53 ± 0.14
0.63 ± 0.14

-test (-value)
4.3E-7**
0.22
0.031*
0.20
0.12
1.8E-3**
1.9E-4**

Data presented are the mean ± SD except for the -test result. *: , **: . Abbreviated: “” = systolic blood pressure, “〈Bp〉” = average blood pressure, and “” stands for the maximum pulse strength among all pressure steps, whose unit is device-specific. “-test” = the two-sample -test between the floating and sunken pulse groups.

Table 3 shows that BMI and (1) and (2) were significantly different between the two pulse groups and in both genders at the significance level of , while blood pressures, heart rate, or remained not different significantly. Here is the maximum value among the measured pulse strengths ranging from to , where is the pulse strength at the th pressure-step averaged over the three palpation locations. Other parameters such as diastolic blood pressure, difference between systolic and diastolic blood pressure were not significantly different between the floating pulse group and sunken pulse group.

The Pearson correlation coefficient (PCC) between (1) and (2) for the 121 concordant subjects is 0.96. Despite strong correlation between (1) and (2), Table 3 shows that (2) is more meaningfully different than (1) between the floating and sunken pulse groups in both genders. It implies that (2) may give improved performance over (1) in predicting the floating and sunken pulses. Based on this observation, we examined the diagnostic performance of pulse qualities using (1) and (2), with reference to the OMDs’ diagnoses. The performance of the algorithm with reference to the OMDs’ diagnoses is slightly better in male volunteers than in female volunteers. Since the difference between two gender groups is minor, hereafter we do not distinguish genders and present results only for the entire subject group.

4.3. Determination of Pulse Qualities by (1) and (2)

We attempt to classify the measured pulses into the floating or sunken pulses using (1) and (2). A simple decision rule can be made by where is the discriminant for the floating and sunken pulses. Applying this rule, we obtain Figure 7.

(a)

(b)

Figure 7: The number of decisions by (3) that are concordant with the diagnoses by paired OMDs, for the floating (asterisked curve) and sunken pulses (squared curve), as a function of the discriminant . Result by using (1) is shown in (a) and (2) in (b), where the plain curve corresponds to the total concordant decisions as the floating or sunken pulses, between the algorithm and OMDs. The boxed area in (a) is the regime showing both the maximal accuracy and MCC, while in (b) the accuracy is shown maximal in the solid box and the MCC is shown maximal in the dotted box.

Figure 7 shows the number of concordant decisions between the classification by the rule in (3) and diagnosis by paired OMDs as a function of the discriminant . Results by using (1) is shown in (a) and (2) in (b). As plotted by the asterisked curve, the number of concordant decisions for the floating pulse increase with increasing , while it was the opposite for the sunken pulse (squared curve). The overall maximum in the concordant decisions for the floating and sunken pulses lies in the central regime of . Specifically, the accuracy (ACC), which is the ratio between the concordant decisions and total decisions, is maximally given by ACC = 82/121 = 67.8% at with (1) (boxed region in (a)), and % at with (2) (solid box in (b)). On the other hand, the Matthews correlation coefficient (MCC) is maximally given by MCC = 0.31 at with (1) (boxed region in (a)), and MCC = 0.40 at with (1) (dotted box in (b)).

Decisions by both (1) and (2) show similar patterns in accuracy with maximum values given similarly, that is, ACC = 67.8% with (1) and ACC = 69.4% with (2). However, compared to (1), (2) shows improved correlations with decisions by OMDs. Noting that the MCC between pairs of OMDs is 0.42, the MCC between diagnoses by (2) and OMDs is reasonably good. Both the Matthews correlation coefficient and the two-sample -test in Table 3 indicate that (2) yields better concordance with the pulse determinations by OMDs than (1). It implies that an adequate pair of the low and high applied pressures to discriminate the floating pulse from the sunken pulse is given by (, ) = (37 mmHg, 143 mmHg). We will discuss this result further in the section of Discussion and Conclusions.

By setting the discriminant for the floating and sunken pulses independently, we can improve the accuracy with the regime between the two discriminants left undetermined. In this case, the Matthews correlation coefficient is poorly defined and we do not take this correlation into account. We denote the discriminant for the floating pulse by and that for the sunken pulse by . Then, the classification rule (3) is modified as where the determinants satisfy , and it reduces to (3) in the case of . The result is shown in Figure 8 for (2) (we skip the result for (1)). Since the rule in (4) may have the regime of left undetermined, we introduce the selection rate (SR) to count the determined pulses either as floating or sunken pulses with respect to the total number of subjects to be determined, that is,

(a)

(b)

Figure 8: Decision by (4) with (2) for (a) the floating pulse and (b) sunken pulse as a function of the discriminant (a) and (b) . The asterisked curve corresponds to the total number of decisions, and the squared curve the concordant decisions with the OMDs (scale on the left vertical axis), while the plain curve is the accuracy (scale on the right vertical axis).

Figure 8 shows classifications into (a) the floating pulse and (b) the sunken pulse by the rule in (4) with (2). Both the total number of decisions (asterisked curve) and the concordant decisions with OMDs’ (squared curve) increase with increasing for the floating pulse, while they increase with decreasing for the sunken pulse. On the other hand, the accuracy (plain curve) tends to increase with decreasing number of decisions made. The sudden drop in accuracy for the sunken pulse close to 1 () is unimportant anomaly that may occur with small sample size. An optimal pair of two discriminants turns out to be (, ) = (0.58, 0.68) (boxed areas) with respective accuracy given by ACC = 29/47 = 61.7% for the floating pulse, and ACC = 35/40 = 87.5% for the sunken pulse. The overall accuracy for the 87 determined cases is ACC = 64/87 = 73.6%, which is improvement of about 4.2% compared to the result by (3), while 34 subjects (28.1%) remained undetermined as the opportunity cost. A similar result can be obtained with (1). For instance, for an optimal pair of (, ) = (0.45, 0.52), ACC = 24/38 = 63.3% for the floating pulse and ACC = 43/57 = 75.4% for the sunken pulse. The overall accuracy for the 95 determined cases is ACC = 67/95 = 70.5%, which is improvement of about 2.7% in accuracy with 26 subjects (21.5%) left undetermined.

Depending on one’s preference for accuracy versus selection rate, one can have different combinations for an optimal pair of discriminants (, ). For instance, choosing much larger than is expected to increase the accuracy but decrease the selection rate, while choosing an equal value for and will maximize the selection rate (SR = 1) but at the expense of the accuracy. To examine gain and loss between the selection rate and accuracy as a function of the discriminants (, ), we calculated the number of concordant decisions and the number of decisions made on two-dimensional space spanned by and .

As a result, in Figure 9, we present the accuracy for the maximum concordant decisions as a function of the selection rate, where the asterisked curve presents results by (1) and the squared curve presents results by (2). There exists an optimal pair of discriminants (, ) with which the concordant decisions become maximal for a given number of decisions made as the floating or sunken pulses (refer Figure 8). Information on the optimal pairs of discriminants is not shown in the figure. It shows better accuracy with (2) than with (1) in the entire range of the selection rate. Especially, for (2), a meaningful improvement in accuracy is observed at the expense of selection rate.

Figure 9: Decision by (4) with (1) (asterisked curve) and (b) (2) (squared curve). The graph shows the accuracy for the maximum number of coincidences with OMDs’ decisions as a function of the decision rate, that is, the number of decisions made into the floating or sunken pulses.

Indeed, the undetermined regime of satisfying can be classified into the group of middle-depth pulse. Then (4) can be modified to classify the pulse into the three groups in the level of depth by Equation (6) classifies a pulse into a floating, a sunken, or a middle-depth pulse, which simulates OMDs clinical diagnosis in practice. Applying this rule to Figure 9, the number of the middle-depth pulse becomes “the number of subjects − number of decision” (reading x-axis). It shows a clear tendency for (2) (squared curve) that increasing the number of the decisions into the middle-depth pulse group increases the accuracy in the decisions into the floating and the sunken pulse groups.

5. Discussion and Conclusions

We proposed a novel classification algorithm for the floating pulse and sunken pulse and carried out a clinical test to examine its validity. First, we introduced the coefficient of the floating and sunken pulses () which is bounded between 0 and 1. When is close to 0, it is classified as a floating pulse, When is close to 1, it is classified as a sunken pulse. Using the classification algorithm based on the new method, we performed a clinical study to determine correlation between diagnoses made by the OMDs and by the algorithm.

For the pulse decisions by OMDs, out of 12 OMD participants, we paired two OMDs in different combination on each day of the study. Among 169 simultaneous diagnoses made by paired OMDs, concordant diagnoses were made on 121 subjects (72%), and the rest (28%) showed divergent diagnoses, yielding a Matthews correlation coefficient of 0.42. It lies within a reasonable range of agreement between the OMDs. Out of the 121 concordant diagnoses, 72 were diagnosed as the sunken pulses (sunken pulse group), and 49 were the floating pulses (floating pulse group).

Two-sample -tests between the floating pulse group and the sunken pulse group showed that BMI and were significantly different () in both genders. This result is in accordance with an experiential knowledge that obese individuals are likely to have sunken pulses [4]. We can assume that there is a positive correlation between BMI and skin thickness. This leads to the conclusion that a major factor affecting the pulse depth may be tissue thickness between the radial artery and skin surface. A complementary experiment, for example, measuring skin thickness by ultrasonography, will help us verify the assumed correlation. On the other hand, average blood pressure between the floating and sunken pulse groups was different at the significance level of , while it did not show differences in heart rate, maximum pulse strength, and so forth It indicates that average blood pressure is another factor which determines the pulse depth. By these observations, we suggest that two major factors determining the pulse depth are the thickness of skin tissue above the radial artery and the average blood pressure.

For the 121 subjects diagnosed concordantly by paired OMDs, we tested for consistency with the novel classification algorithm. For the entire study cohort, the maximum accuracy reached was about using coefficient (1) or (2), while the Matthews correlation coefficient reached was up to 0.40 for (1) and 0.31 for (1). It suggests that the classification algorithm using either coefficient yields an acceptable concordance with OMDs’ diagnoses, and (2) yields an optimal performance in its concordance. It suggests that an adequate pair of low and high applied pressures for our algorithm is given by () = (37 mmHg, 143 mmHg). Since our algorithm simulates diagnosis by OMDs who compare pulse strengths at low and high hold-down pressures to determine the pulse depth, the above result implies that OMDs are likely to determine the pulse depth by comparing pulse strengths at a weak applied pressure of  mmHg and at a strong pressure of  mmHg, one of which lies much below the normal range of the blood pressure and the other one lies above.

When the pulse strength is larger at than at , the maximal pulse strength will be positioned well below the average blood pressure ( mmHg) and it shows a typical floating characteristic, and when the pulse strength at is larger, it shows a sunken characteristic. On the other hand, if pulse strengths at and are not so different, the maximum pulse will be observed around , yielding a middle-depth pulse. A more accurate choice for the low and high pressures to be used in the algorithm can be made with a pulse-taking device that offers more narrowly spaced hold-down pressures. A pulse-taking device with 9 steps of the hold-down pressure that are spaced equally up to 180 mmHg is an appropriate candidate for such purpose.

Finally, we suggested an updated version of the above-tested algorithm by allowing a regime of the middle-depth pulse group for the pulses satisfying , which is more useful since OMDs in clinics diagnose a pulse either as floating, sunken, or neither of them corresponding to the middle-depth pulse. The implication of the middle-depth pulse is a sign of balance between Yang/function (Qi) and Yin/form (blood) [3].

In conclusion, we proposed a novel algorithm to classify pulse into floating, sunken, or middle-depth categories using a coefficient called . By a clinical test, primarily we confirmed the validity of the novel classification algorithm and secondly we found appropriate numerical values of the pressure steps which are apt for determining the level of pulse depth. To our knowledge, this study is the first report on developing classification algorithm for the floating and sunken pulses based on a massive clinical data with human subjects of the order of hundred. This work may motivate further research activities towards developing objective classification algorithms for other pulse qualities.

Acknowledgments

This work was supported by a grant from the Korea Institute of oriental Medicine (KIOM) funded by the Korea government (K10070) and was partially supported by the Korea Ministry of Knowledge Economy (10028438). The authors thank Dr. Hee Jeong Jin, Boncho Ku, Jaechul Lee, and Dr. Jun-Hyeong Do for discussions. The authors declare no commercial associations that may lead to conflicts of interests in connection with the submitted paper.

References

H. K. Huynh and G. M. Seifert, Pulse Diagnosis by Li Shi Zhen, Paradigm Publications, Sydney, Australia, 1985.

B. Flaws, The Secret of Chinese Pulse Diagnosis, Blue Poppy Press, Denver, Colo, USA, 2006.

S. Walsh and E. King, Pulse Diagnosis: A Clinical Guide, Elsevier, Sydney, Australia, 2008.

L. I. Hammer, Chinese Pulse Diagnosis: A contemporary approach, Eastland Press, Seattle, Wash, USA, 2001.

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.

Z. F. Fei, Contemporary Sphygmology in Traditional Chinese Medicine, People’s Medical Publishing House, Beijing, China, 2003.

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.

J. Lee, Y. J. Jeon, H. J. Lee, and J. Y. Kim, “Variations in non-invasive radial artery pulse waveform at different measuring positions,” in Proceedings of the 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '08), pp. 1415–1418, August 2008.

Y.-J. Lee, J. Lee, and J.-Y. Kim, “A Study on Characteristics of Radial Arteries through Ultrasonic Waves,” in Proceedings of the 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBS '08), pp. 2453–2456, August 2008.

Y. J. Lee, J. Lee, H. J. Lee, and J. Y. Kim, “A study on correlation between BMI and oriental medical pulse diagnosis using ultrasonic,” in Proceedings of the 13th International Conference on Biomedical Engineering, pp. 2052–2055, 2008.

Y. J. Jeon, J. U. Kim, H. J. Lee, et al., “A clinical study of the pulse wave characteristics and the three pulse diagnosis positions of Chon, Gwan, and Cheok,” Evidence-Based Complementary and Alternative Medicine. In press.

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.

S. Tan, K. Tillisch, and E. Mayer, “Functional somatic syndromes: emerging biomedical models and traditional Chinese medicine,” Evidence-Based Complementary and Alternative Medicine, vol. 1, no. 1, pp. 35–40, 2004.

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.

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.

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.

Y. Chen, L. Zhang, D. Zhang, and D. Zhang, “Wrist pulse signal diagnosis using modified Gaussian models and Fuzzy C-Means classification,” Medical Engineering and Physics, vol. 31, no. 10, pp. 1283–1289, 2009.

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, pp. 361–366, 2010.

S. W. Lee, Y. J. Lee, H. J. Lee, H. J. Kang, and J. Y. Kim, “Clinical study on the floating and sunken pulse detection with piezoresistive sensors and contact pressure control robot,” Korean Journal of Oriental Physiology < Pathology, vol. 19, no. 6, pp. 1673–1675, 2005.

S. H. Kim, J. U. Kim, Y. J. Lee, K. H. Kim, and J. Y. Kim, “New algorithm of determining the floating and sunken pulse with a pulse diagnosis instrument,” Korean Journal of Oriental Physiology & Pathology, vol. 23, no. 6, pp. 1221–1225, 2009.

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,” Journal of Korean Oriental Medicine, vol. 28, no. 3, pp. 13–22, 2007.