Skip to main content

Classification of Parkinson’s disease utilizing multi-edit nearest-neighbor and ensemble learning algorithms with speech samples

Abstract

Background

The use of speech based data in the classification of Parkinson disease (PD) has been shown to provide an effect, non-invasive mode of classification in recent years. Thus, there has been an increased interest in speech pattern analysis methods applicable to Parkinsonism for building predictive tele-diagnosis and tele-monitoring models. One of the obstacles in optimizing classifications is to reduce noise within the collected speech samples, thus ensuring better classification accuracy and stability. While the currently used methods are effect, the ability to invoke instance selection has been seldomly examined.

Methods

In this study, a PD classification algorithm was proposed and examined that combines a multi-edit-nearest-neighbor (MENN) algorithm and an ensemble learning algorithm. First, the MENN algorithm is applied for selecting optimal training speech samples iteratively, thereby obtaining samples with high separability. Next, an ensemble learning algorithm, random forest (RF) or decorrelated neural network ensembles (DNNE), is used to generate trained samples from the collected training samples. Lastly, the trained ensemble learning algorithms are applied to the test samples for PD classification. This proposed method was examined using a more recently deposited public datasets and compared against other currently used algorithms for validation.

Results

Experimental results showed that the proposed algorithm obtained the highest degree of improved classification accuracy (29.44%) compared with the other algorithm that was examined. Furthermore, the MENN algorithm alone was found to improve classification accuracy by as much as 45.72%. Moreover, the proposed algorithm was found to exhibit a higher stability, particularly when combining the MENN and RF algorithms.

Conclusions

This study showed that the proposed method could improve PD classification when using speech data and can be applied to future studies seeking to improve PD classification methods.

Background

Parkinson’s disease (PD) is a neurodegenerative disorder of the central nervous system that is characterized by a partial or full loss in motor reflexes, speech, behavior and other vital functions. It is generally observed in elderly people and causes disorders in patient speech and motor abilities (writing, balance, walk, etc.) [1]. PD is one of the most common neurodegenerative disorders, with an incidence rate of 20/100,000 and affecting approximately 5 million people worldwide, with half of PD cases found in China [2]. Moreover, these statistics are likely to be underestimating the incident rate due to difficulties in diagnosing PD. With populations growing, the number of diagnosed PD patients will continue to grow, thus increasing the damage of this disease in the future [2].

While identifying the causes of PD onset remain elusive, PD often referred to as an idiopathic disorder, genetic and environmental factors have been implicated [1]. Unfortunately, a reliable PD biomarker has yet to be identified, but the symptoms can often reflect PD occurrence and progression, such as tremor, rigidity, loss of muscle control, slowness in movement, poor balance and, especially, voice problems [16]. Therefore, diagnosing PD based on symptoms is reasonable and effective [714]. Among them, speech has been shown to be a useful signal for discriminating PWP (people with Parkinson’s) from healthy controls, with clinical evidence suggesting that the vast majority of PWP typically exhibit some form of vocal disorder [3, 7, 1416]. In fact, vocal impairment may be among the earliest prodromal PD symptoms and can be detectable up to five years prior to clinical diagnosis [2]. Thus, utilizing speech data can aid in the development of a noninvasive early PD diagnostic method, with speech alterations including reduced loudness, increased vocal tremor and breathiness (noise), while PD-associated vocal impairment is characterized by dysphonia (inability to produce normal vocal sounds) and dysarthria (difficulty in pronouncing words). Dysphonia can be measured by utilizing acoustic tools that detect voice abnormalities, such as aperiodic vibrations, non-Gaussian randomness, abnormality of vowel “a” phonations, etc. [2], with certain special words and sentences, such as Arabic numerals, special words, etc., used for detection.

While a reliable PD diagnosis is difficult, the effectiveness of speech based diagnostic approaches has inspired researchers to develop decision support tools able to extract dysphonic speech features and design classification algorithms to distinguish PD patients from healthy ones, with speech feature extraction, feature selection/transformation and classifier design [1724].

When examining feature extraction, the feature types commonly include pitch type, energy type, speed type and content type [1, 2, 1725]. As for feature transformation, the frequently used algorithms are PCA (principle component analysis) [26, 27, 31, 49]. As for feature selection, the frequently used algorithms are NN (neural network) based [2730, 32, 49], serial search based [2, 14, 29, 31], random based [32, 33, 48], p value based [2, 2734], relevance based [35, 36] or entropy based [37], discrimination algorithm (DA) based [47]. As for classifier design, the predominantly used classifiers include a support vector machine (SVM) [1, 2, 14, 29, 32, 35, 3841], KNN [1, 2, 26, 28, 40, 41, 47, 48, 49], random forest (RF) [2, 30, 36], Bayesian network [27, 28, 40, 42, 43, 48], discrimination algorithm (DA) [27, 29, 31, 37], probabilistic neural network (PNN) [27, 43] or decision tree [31, 40, 42, 4446]. Besides, several ensemble models were involved compared with single classifier [27, 47, 48, 50, 51].

While the methods mentioned above have been effective to some extent, obtaining optimal speech samples is difficult, with low quality samples being prone to misleading the classifiers and negatively impacting accuracy. Based on this potential pitfall, this study examined a speech sample selection algorithm to improve outcomes. First, a multi-edit-nearest-neighbor (MENN) algorithm was utilized to select the optimal speech samples iteratively, thereby enhancing the separability of the training samples. The MENN algorithm effectively optimizes training samples by removing samples that could mislead the classifier [50, 52, 54]. Next, an ensemble learning algorithm, random forest (RF) or decorrelated neural network ensembles (DNNE) was employed for classification on the selected training samples. These two ensembles were chosen for point of comparison, with RF being a classical ensemble learning algorithm with a proven stability [2, 34, 41], while DNNE is a newer ensemble learning algorithm able to effectively maximize differences between sub-learning algorithms [51]. Lastly, a classification was conducted on the test samples based on the trained ensemble learning algorithms.

Methods

Data descriptions

The data utilized in this study was obtained from the Parkinson speech dataset with multiple types of sound recordings [1], was deposited by Sakar et al. [1] and is available on the University of California, Irvine (UCI) machine learning dataset repository website. The deposited dataset (Table 1) included two datasets entitled “Training_Data” and “Test_Data”. The “Training_Data” set included 40 subjects, which were 20 PD (6 women, 14 men) and 20 healthy subjects (10 women, 10 men), with 26 speech samples collected per subject. The samples contained a variety of speech segments with varying pronunciations, to include continuous vowel sounds, digital pronunciation, word pronunciation and phrase sentence pronunciation. To set-up a feature vector, linear and nonlinear feature parameters were extracted for each speech sample in 26 dimensions. The “Test_Data” set included 28 subjects (all PD) with six speech samples collected per subject, with half containing pronunciation of the continuous vowel ‘a’, while the other half had pronunciation of the continuous vowel ‘o’. Feature vectors were then constructed by the 26 dimensions.

Table 1 Description of speech samples based on the same subject

The data collected in the context of this study belongs to 20 PWP (6 female, 14 male) and 20 healthy individuals (10 female, 10 male) who appealed at the Department of Neurology in Cerrahpasa Faculty of Medicine, Istanbul University. Test group consists of patients who are suffering from PD for 0 to 6 years. Individual ages vary between 43 and 77 (mean: 64.86, standard deviation: 8.97) along with 45 and 83 (mean: 62.55, standard deviation: 10.79) for test and control groups, respectively [1].

In addition to the dataset deposited by Sakar et al. [1], another broadly studied dataset also exists that was deposited by Little et al. [2, 14]. The major reasons the Sakar et al. dataset was chosen is as follows: (1) the Little et al. dataset has been more broadly studied, with a classification accuracy over 95% achievable, which is comparable to the close to 100% accuracy achieved herein; (2) the Sakar et al. dataset contains more samples, thus providing more statistically meaningful results; (3) only a few studies have examined the Sakar et al. dataset, thus making findings more insightful; and (4) corresponding classification accuracies for the Sakar et al. dataset are not promising at time. The aim of this study is to show that this type of data collection can lead to high classification accuracy of PD just by altering the classification method. Therefore, in the experimental section, most of the experiments are based on the PD data from Sakar et al. [1]. Since the PD data from Max Little was investigated frequently, the proposed algorithm is verified based on the data too.

Flowchart of the proposed algorithm (PD_MEdit_EL)

For simplification, the proposed algorithm was called PD_MEdit_EL and a flowchart can be seen below (Fig. 1). The PD_MEdit_EL algorithm includes two major parts, with one part optimizing speech samples via MENN while the other performs the classification using the ensemble learning algorithm, with either RF or DNNE used. Chronologically, training speech samples are optimized using the MENN algorithm, these training samples are then processed with the ensemble learning (EL) algorithm and a trained learning model is obtained. Finally, the trained learning algorithm is applied to aid in test sample classification.

Fig. 1
figure 1

Flowchart of the PD_MEdit_EL algorithm

Multi-edit-nearest-neighbor algorithm (MENN)

The MENN algorithm serves as a prototype selection algorithm when samples with different classes have overlapping distributions, the overlapping regions are more likely to be misclassified. Presumably, if the overlapping regions are accurately removed, the ‘noise’ or misleading data can be greatly suppressed. Thus, the remaining samples will better reflect the representative samples for different classes, thereby improving the classification accuracy [52].

The two-edit-nearest-neighbor algorithm mainly includes editing and classification. First, samples are divided into training and validation samples, with the validation samples classified and the misclassified samples removed. Generally, if the sample size is large enough, MENN should be applied. The main process of the algorithm is as follows:

Step 1: The original samples \(X = \left\{ {x_{i} = \left( {x_{i1} ,x_{i2} , \ldots ,x_{id} } \right)|i = 1,{ 2}, \ldots N} \right\}\) were randomly divided into s subsets, \(X_{ 1} ,X_{ 2} , \ldots ,X_{s} \left( {{\text{s}} > 3} \right)\), each containing \(M_{ 1} ,M_{ 2} , \ldots ,M_{s}\) samples. d is the dimension of the sample, N is the number of samples.

Step 2: Editing process. \(X_{i + 1}\) is used as the training set and KNN is used to classify the validation samples in set X i , with misclassified samples deleted. These set are repeated for i = 1 to i = s − 1. If i = s, X 1 is the training set.

Step 3: The samples retained after step 2 are merged to form a new sample set (X New ).

Step 4: Repeat the above steps until all misclassified samples are removed.

In Fig. 2, the results of the training samples being edited with the MENN algorithm can be seen, with the separability of on the training samples improved when comparing pre-MENN editing (Fig. 2a) and post-MENN editing (Fig. 2b).

Fig. 2
figure 2

Edited results of training samples: a prior to MENN and b with MENN.(1st dimension means the 1st feature; 2nd dimension means the 2nd feature of the PD data)

Previous studies have shown that the nearest neighbor algorithm can achieve a classification error rate close to the Bayesian error rate P* when a large enough sample size is utilized [53]. For the two-edit-nearest-neighbor algorithm, the progressive conditional false recognition rate is as follows:

$$P\,(e/x) = 1 - \frac{{p(\omega_{i} /x)^{2} }}{{\sum\nolimits_{i = 1}^{2} {p(\omega_{i} /x)^{2} } }} .$$
(1)

When the MENN algorithm is applied, the progressive conditional false recognition rate is as follows:

$$P_{M} (e/x) = 1 - \frac{{\sum\nolimits_{i = 1}^{2} {p(\omega_{i} /x)^{{2^{M + 1} }} } }}{{\sum\nolimits_{i = 1}^{2} {p(\omega_{i} /x)^{{2^{M} }} } }}$$
(2)

As seen in the formula, the progressive conditional false recognition rate will be greatly reduced when Mincreases. If \(M \to \infty ,{\text{ then lim}}_{M \to + \infty } P_{M} \left( {e/x} \right) = { \hbox{min} }\left[ {p\left( {\omega_{ 1} /x} \right),p\left( {\omega_{ 2} /x} \right)} \right] = P*\left( {e/x} \right)\).

Random forest

When the RF algorithm is utilized as the classifier, the corresponding training is conducted as follows [54]:

Step 1: Bootstrap method is used to re-sample and randomly generate T training sets: \(S_{1} ,\;S_{2} , \ldots S_{T}\).

Step 2: Corresponding decision trees are generated for each set: \(C_{1} ,\;C_{2} , \ldots C_{T}\) and mattributes are randomly selected from M attributes, with the determined optimal splitting used to split the node.

Step 3: Each tree is allowed to grow integrally, without being pruned.

Step 4: Based on each decision tree, the sample in the test set X is classified to generate the corresponding classes: \(C_{ 1} \left( X \right),\;C_{ 2} \left( X \right), \ldots C_{T} (X)\).

Step 5: The voting method is applied, output categories are determined based on the T decision trees and the category with the maximum votes becomes the assigned category.

Dnne

The main principle of the DNNE algorithm [51] is described as follows.

Initialization Training setD t ; base function G (option: sigmoid function); regularizing factor λ [0, 1]; the number of the sub-classifiers is M, dimension of sample is L.

Procedure Trained DNNE model.

1: Construct a M single-layer feedback neural network and an ensemble learning algorithm (model).

2: Randomly initialize weights w j and bias b j .

3: Input training samples D i and calculate the output g ij (x n ) of all the sub-classifiers,

$$g_{ij} (x_{n} ) = G_{ij} (w_{j} ,b_{j} ,x_{n} )$$
(3)

where the functional connection network (RVFL) for constructing random vector needs training samples D t with N dimensions, \(D_{t} = \left\{ {\left( {x_{1} ,y_{1} } \right),\left( {x_{2} ,y_{2} } \right), \ldots ,\left( {x_{N} ,y_{N} } \right)} \right\}\) and where \(\left( {x_{n} ,y_{n} } \right) \in R^{d} \times R\) belongs to pairs of observations. The G ij () means the ij sub-classifier.

4: Calculate constants C 1 , C 2 as Eq. (4):

$$\varphi \left( {i,j,k,l} \right) = \sum\limits_{n = 1}^{N} {g_{ij} \left( {x_{n} } \right)g_{kl} \left( {x_{n} } \right),\quad \varphi \left( {i,j} \right) = \mathop \sum \limits_{n = 1}^{N} g_{ij} \left( {x_{n} } \right)y_{n} }$$
(4)

\(\varphi \left( {i,j,k,l} \right)\) represents the correlation between the jth hidden neuron of ith independent RVFL network and the lth hidden neuron of kth independent RVFL network. \(\varphi \left( {i,j} \right)\) signifies the correlation between the jth hidden neuron of ith based network and object function ψ(x). Thus, C 1 and C 2 can be represented as Eqs. (5) and (6):

$$C_{1} = 1 - 2\lambda \frac{{(M - 1)^{2} }}{{M^{2} }},\quad C_{2} = 2\lambda \frac{M - 1}{{M^{2} }},$$
(5)
$$\mathop \sum \limits_{k = 1}^{L} C_{1} \varphi \left( {i,j,i,k} \right)\beta_{ik} + \mathop \sum \limits_{l \ne i}^{M} \mathop \sum \limits_{k = 1}^{L} C_{2} \varphi \left( {i,j,l,k} \right)\beta_{ik} = \varphi \left( {i,j} \right).$$
(6)

By following the formula above, the errore i , output weight β ij and M × L linear equation can be obtained.

So, the linear system can be calculated based on theβ ij and the RVFL ensemble model can be obtained. In order to simplify the calculation, the matrix form of the linear system is as follows Eq. (7):

$$H_{corr} B_{ens} = T_{h} .$$
(7)

H corr is hidden correlation matrix, B ens is the global output weights matrix and T h is the hidden-target matrix. H corr is defined as follows Eqs. (8) and (9):

$$H_{corr} = \left[ \begin{array}{ccccccc} C_{1} \varphi (1,1,1,1) & \cdots & C_{1} \varphi (1,1,1,L) & \cdots & C_{2} \varphi (1,1,M,1) & \cdots & C_{2} \varphi (1,1,M,L) \\ \vdots& & \vdots& & \vdots & & \vdots \\ C_{1} \varphi (1,L,1,1) & \cdots & C_{1} \varphi (1,L,1,L)& \cdots & C_{2} \varphi (1,L,M,1) & \cdots & C_{2} \varphi (1,L,M,L) \\ C_{2} \varphi (2,1,1,1) & \cdots & C_{2} \varphi (2,1,1,L)& \cdots& C_{2} \varphi (2,1,M,1) &\cdots & C_{2} \varphi (2,1,M,L) \\ \vdots & &\vdots & & \vdots & & \vdots \\ C_{2} \varphi (2,L,1,1)& \cdots & C_{2} \varphi (2,L,1,L) & \cdots & C_{2} \varphi (2,1,M,1) & \cdots & C_{2} \varphi (2,L,M,L)\\ \vdots & & \vdots & & \vdots & & \vdots \\ C_{2} \varphi (M,1,1,1) & \cdots &C_{2} \varphi (M,1,1,L) & \cdots & C_{1} \varphi (M,1,M,1) &\cdots & C_{1} \varphi (M,1,M,L) \\ \vdots & & \vdots & & \vdots & & \vdots \\ C_{2} \varphi (M,L,1,1) & \cdots& C_{2} \varphi (M,L,1,L) & \cdots & C_{1} \varphi (M,L,M,1)& \cdots & C_{1} \varphi (M,L,M,L) \end{array} \right]_{(ML \times ML) ,}$$
(8)
$$H_{corr(p,q)} = \left\{\begin{array}{l} C_{1} \varphi (m,n,k,l)\quad {\text{ if m}} = {\text{k;}} \\ C_{2} \varphi (m,n,k,l)\quad {\text{ otherwise;}} \\ \end{array}\right.$$
(9)

where p, \(q = 1, \ldots ,M \times L;\) \(m = \left\lceil {\frac{p}{L}} \right\rceil\), \(n = ((p - 1)\bmod \, L) + 1\); \(k = \left\lceil {\frac{q}{L}} \right\rceil\); \(l = ((q - 1)\bmod \, L) + 1\); and mod means modular operation. B ens and T h can be defined as follows Eq. (10):

$$B_{ens} = \left[ {\beta_{11} , \ldots ,\beta_{1L} ,\beta_{21} , \ldots ,\beta_{2L} ,\beta_{2L} , \ldots ,\beta_{M1} , \ldots ,\beta_{ML} } \right]_{ML \times 1}^{T} ,$$
(10)
$$T_{h} = \left[ {\varphi \left( {1, \, 1} \right), \ldots ,\varphi \left( {1,L} \right),\varphi \left( {2, \, 1} \right), \ldots ,\varphi \left( {2,L} \right), \ldots ,\varphi \left( {M, \, 1} \right), \ldots ,\varphi \left( {M,L} \right)} \right]_{ML \times 1}^{T} .$$

β ij can be modified using the optimization technique of gradient descent and is formulated as follows Eq. (11):

$$\hat{B}_{ens} = H_{corr}^{ - 1} T_{h} .$$
(11)

The DNNE algorithm can be programmed and realized. It’s pseudo code is as follows:

Results

Experimental condition

One thousand and forty training samples collected from 40 subjects within the “Train_Data” set, with each composed of 26 dimensional feature parameters, were evaluated with leave one out cross validation (LOO CV). If the majority of subject samples were patient, then the subject was deemed a patient; otherwise, the subject was deemed healthy. To reduce the effect of outlier samples from the same subject since not all samples can reflect speech characteristics equally, a cross validation method called leave-one-subject-out (LOSO) was employed. When completing the DNNE algorithm, the linear exhaustive searching algorithm was used to search for maximum parameter values, with the search range of M within [2, 15], the base function L of the RVFL network within [5, 50] and the penalty coefficient γ and boosting threshold value ϕ within [0, 1]. Their step sizes were 0.1 and 0.01, respectively. Additionally, 60% of the training samples are frequently used for bagging and boosting the ensemble. When performing RF, the statistically optimal number of decision trees was 500. The SVM was based on libsvm. When performing the MENN algorithm, the setup parameter (s) was based on statistical experiments. For MENN, the s = 4. For the weights of the classification model, they are set by supervised training.

In this paper, the experimental operating system platform was the Windows, version 7, 32-bit operating system, and the memory size was 4 GB. The data processing was completed in MATLAB, version 2014a. The SVM, RF and relief algorithms are from the toolboxs under MATLAB environment. The DNNE algorithm is from the official website of MATLAB. Original MENN algorithm is from the file exchange related websites (http://www.pudn.com), but it was modified and combined with the DNNE and RF by us. Other parts are designed by us.

Performance evaluation criteria

To verify the effectiveness of the proposed algorithm, the classification accuracy, sensitivity and specificity were used as an evaluation standard utilizing the terms TP: true positive, TN: true negative, FP: false positive and FN: false negative. Specific formulas are as follows:

$$Accuracy = \frac{TP + TN}{TP + FP + TN + FN}$$

Sensitivity, also called the true positive rate, determines the percentage of accurately identified disease subjects by the equation:

$$Sensitivity = \frac{TP}{TP + FN}$$

Specificity, also called the true negative rate, determines the percentage of accurately identified healthy subjects by the equation:

$$Specificity = \frac{TN}{TN + FP}$$

Classification performance of the PD_MEdit_EL algorithm

In the present study, PD_MEdit_EL algorithm, which utilized the ensemble learning algorithms RF and DNNE, was examined. Until now, only four studies have examined the Sakar et al. [1] dataset, with two only reporting on the classification accuracy of the validation set [42, 49]. Therefore, only two studies [1, 27] remain that are comparable to the PD_MEdit_EL algorithm examined herein. Currently, SVM is widely applied in the classification of the Parkinson’s disease, thus SVMs with linear and RBF kernel functions were compared to the proposed PD_MEdit_EL algorithm (Table 2). These results showed that the PD_MEdit_EL algorithm provided the most accurate classification in terms of classification accuracy (ACC), sensitivity (SEN) and specificity (SPE), regardless of LOO and LOSO. When examining all classification determinants, the RF (with MENN) was the best overall classification algorithm, with the highest mean SPE of LOSO seen with the DNNE (with MENN) algorithm. When comparing the LOO and LOSO methods, the classification performance of LOSO was higher than LOO, possibly due to LOSO effectively reducing outliers within a subject.

Table 2 Classification accuracy for the “Training_Data” set

Overall, these findings show that when MENN is combined with either ensemble learning algorithm, a notable improvement is seen over other classification methods. When examining ACC in the RF (with MENN) method, improved classification was seen when compared to the SVM (linear), increased 16.5%, SVM (RBF), increased 14%, method in [1], increased 29.44% and method in [27], increased 0.3%. When examining accuracy in the DNNE (with MENN), classification improved by 8.37% relative to SVM (linear), 5.87% relative to SVM (RBF) and 21.31% relative to the method in [1], but decreased by 2.5% relative to the method in [27]. As for the SEN for the RF (with MENN) method, a classification improvement of 27.5% was seen relative to the SVM (linear), 12.5% relative to the SVM (RBF), 37.58% relative to the method in [1] and 7% relative to the method in [27]. However, when examining the SEN for the DNNE (with MENN) method, a classification improvement was only seen relative to the SVM (linear), increased 6% and the method in [1], increased 16.08%. These results showed that the RF method was more robust than the DNNE algorithm in terms of classification performance.

Next, the classification accuracy of several classification algorithms was examined for the “Test_Data” set (Table 3). Since this set only contained patient subjects, the sensitivity and specificity could not be calculated. Furthermore, the method in Ref. [27] did not report the classification accuracy for the “Test_Data”, so it was omitted. The results when using this dataset further pointed out the strength of the proposed PD_MEdit_EL algorithm, with significantly higher classification accuracies seen in the RF (with MENN) and DNNE (with MENN). When comparing the DNNE and RF algorithms, no significant difference in terms of classification accuracy was noted. The results obtained in Tables 2 and 3 are also graphically displayed in Figs. 3 and 4.

Table 3 Classification accuracy of the “Test_Data” set
Fig. 3
figure 3

Classification performances between different classification algorithms using the “Training_Data” set. The y-axis: classification accuracy; DNNE_MENN and RF_MENN: reflect the proposed PD_MEdit_EL algorithm; SVM_Linear: SVM with the linear kernel function; SVM_RBF: SVM with radial basis function kernel; Method in [27]: classification algorithm from Ref. [27]

Fig. 4
figure 4

Classification performances between different classification algorithms using the “Test_Data” set. The y-axis: classification accuracy; DNNE_MENN and RF_MENN: reflect the proposed PD_MEdit_EL algorithm; SVM_Linear: SVM with the linear kernel function; SVM_RBF: SVM with radial basis function kernel; Method in [27]: classification algorithm from Ref. [27]

To further examine classification performance, differences between algorithms relative to run number were examined; both the ‘Training_Data’ and ‘Test_Data’ sets were examined 10 times each for each algorithm. These findings were consistent with the data presented in Tables 2 and 3. When examining the “Training_Data” set, both the RF (with MENN) and DNNE (with MENN) showed the highest classification accuracies. However, the RF (with MENN) showed a greater degree of stability compared to the DNNE (with MENN), possibly due to the DNNE being complex. While the DNNE (with MENN) method showed an overall improved classification accuracy, its low stability make it a less ideal candidate. While the method referenced in [1] was not included in this comparison due to the paper [1] did not include the accuracies during number of runs. Since it comprises both SVM_Linear and SVM_RBF components, an indirect comparison can still be drawn to show that the PD_MEdit_EL algorithm provides higher classification accuracy than the method in [1].

When examining the “Training_Data” set (Fig. 5), the RF (with MENN) again showed the highest degree of accuracy with a high stability, while both the RF (with MENN) and DNNE (with MENN) showed a higher classification accuracy than the SVM, which is currently used in the classification of Parkinson’s disease. When examining this dataset, it is worth noting that the RF becomes more stable, thus suggesting a higher compatibility with this set.

Fig. 5
figure 5

Classification accuracy of different algorithms for the “Training_Data” set. The y-axis: classification accuracy; the x-axis: the number of runs for each algorithm; DNNE_MENN and RF_MENN: reflect the proposed PD_MEdit_EL algorithm; SVM_Linear: SVM with the linear kernel function; SVM_RBF: SVM with radial basis function kernel

When examining the “Test_Data” set (Fig. 6), the RF (with MENN) again showed the highest degree of accuracy with a high stability, while both the RF (with MENN) and DNNE (with MENN) showed a higher classification accuracy than the SVM, which is currently used in the classification of Parkinson’s disease. When examining this dataset, it is worth noting that the DNNE becomes more stable, thus suggesting a higher compatibility with this set.

Fig. 6
figure 6

Classification accuracy of different algorithms for the “Test_Data” set. The y-axis: classification accuracy; the x-axis: the number of runs for each algorithm; DNNE_MENN and RF_MENN: reflect the proposed PD_MEdit_EL algorithm; SVM_Linear: SVM with the linear kernel function; SVM_RBF: SVM with radial basis function kernel

The proposed algorithm was verified in the PD data from [2]. The Max Little et al. introduced several machine learning algorithms into their dataset [2]. According to the results, the SVM with RBF kernel and relief algorithm was best. Therefore, it was realized here for comparison with the proposed algorithm. For fair comparison, the proposed algorithm was based on relief algorithm too. The CV method is tenfold CV as same as the paper [2] did.

Seen from the Table 4, the classification accuracy rates are better than those in the Table 2. For the proposed algorithm and LOSO, the classification accuracy is improved from 81.5 to 87.8%; the sensitivity is improved from 92.5 to 95.4% respectively. For the proposed algorithm and LOO, the classification accuracy is improved from 70.93 to 87.8%; the sensitivity is improved from 73.27 to 95.4% respectively. The possible reason is that the number of samples in the dataset is smaller than that in the dataset on the Table 2. Besides, it is worth noting that the sensitivity and the specificity are quite different. The possible reason is that the numbers of the patients and healthy people are different. In other words, the dataset is unbalanced.

Table 4 Classification accuracy for the PD data from Max Little et al.

It is worth noting that the 1040 samples are from the 40 subjects, therefore, the samples belonging to same subject are dependent each other to some extent. So, if the samples in training set and test set are independent, the classification accuracy rates possibly become worse. However, the verification section in the relevant papers did not consider this point except the paper the data originates from. In order to further study the problem, an experiment was conducted. In the experiment, when a sample is classified, the other 25 samples belonging to same subject are not used for building the classification model. It can guarantee the samples in training set and test set are independent each other. The Table 5 shows the differences between considering dependency and independency. Seen from the Table 5, the RF_MENN means the algorithm not considering the dependence of the samples; the RF_MENN_inDe means the algorithm considering the dependence of the samples; the SVM_RBF means the SVM with RBF kernel algorithm not considering the dependence of the samples; the RF_MENN_inDe means the SVM with RBF kernel algorithm considering the dependence of the samples.

Table 5 Classification result in terms of dependence

Seen from the Table 5, when considering the dependence of the samples, the classification accuracy rates become worse regardless of RF and SVM algorithm. The extent of becoming worse is large. However, the proposed algorithm considering dependence of samples is still better than the paper [1]. Compared with the proposed algorithm and the SVM with RBF kernel, the former is still better when considering dependence of samples. It means that the proposed algorithm is valuable.

Effect of the MENN algorithm

To examine the effects of the MENN algorithm, the “Training_Data” set, which included 20 PD (6 women, 14 men) and 20 healthy (10 women, 10 men) subjects with 26 speech samples per subject, were examined. In total, 1040 speech samples were utilized and examined pre- and post-MENN, with 1039 trained under LOO (Table 6). With the MENN algorithm, the total number of training samples was reduced from 1039 to 731, with the number of healthy subjects reduced from 519 to 364 and the number of patient subjects reduced from 520 to 367. Furthermore, after applying the MENN algorithm, the number of training subjects did not change despite reducing the sample number. Thus, the MENN algorithm can meet the requirement of the subsequent machine learning for number of training samples.

The observed effect of the MENN algorithm was further visualized (Fig. 7), with the three dimensions being the first three features in the datasets and used for coordinate determination. Prior to applying the MENN algorithm, the PWP mean patient samples and normal mean healthy samples were very mixed and difficult to separate (Fig. 7). After apply the MENN algorithm, the sample mixing was greatly improved, thus enabling better subsequent classifications.

Fig. 7
figure 7

Distribution of samples both (1) before and (2) after applying the MENN algorithm (1st dimension means the 1st feature of the PD data; 2nd dimension means the 2nd feature of the PD data; 3rd dimension means the 3rd feature of the PD data)

Next, an ensemble learning algorithm with and without the MENN algorithm was examined using the “Training_Data” set (Table 7), with the same abbreviations as in Tables 2 and 3 used. The MENN was shown to play an important role in improving classification performance with both DNNE and RF. For DNNE under LOO, the ACC improved 3.62% and the SPE improved 13.33%. For DNNE under LOSO, the ACC improved 4.12%, the SEN improved 2.5% and the SPE improved 5%. For RF under LOO, the ACC improved 3.54%, the SEN improved 1.65% and the SPE improved 5.23%. For DNNE under LOSO, the ACC improved from 3.25% and the SEN improved 7%. The highest improvement was 7% and the highest classification accuracy achieved was 88%.

The ensemble learning algorithms with and without the MENN algorithm were also compared using the “Test_Data” set (Table 8), with the same abbreviations as in Tables 2 and 3 used. These results were in accordance with those presented in Table 5. When examining the ACC, the DNNE for samples improved 14.88%, while the DNNE for subjects improved 10.72%. Also when examining ACC, the RF for samples improved 43.03%, while the RF for subjects improved 45.72%. The highest ACC that was obtained was 100% and the largest increase for the samples was 43.03% and for the subjects was 45.72%. This high degree of improvement may have been due to this dataset only containing patient samples, thus making the classification less difficult.

Table 6 Comparison before and after applying the MENN algorithm

Level of significance of PD_MEdit_EL algorithm

In an attempt to establish that a significant difference is present between the PD_MEdit_EL algorithm and the other examined algorithms, the p values of the ACCs, SENs and SPEs between the algorithms based on ten times experiments were calculated (Table 9). In terms of LOO and LOSO, the differences between the RF (with MENN) and DNNE (with MENN) when compared to the SVM (linear) or SVM (RBF) showed some highly significant differences. While these results show a significant difference between the PD_MEdit_EL algorithm and those examined, a significant difference was not noted when comparing the RF (with MENN) and DNNE (with MENN), thus suggesting that the two algorithms would perform in a comparable fashion.

Table 7 Improvement with the MENN algorithm for the “Training_Data” set

Statistically significant differences were also examined using the “Test_Data” set (Table 10). The same trend was seen with this dataset, with some highly significant differences noted between the RF (with MENN) and DNNE (with MENN) when compared to the SVM (linear) or SVM (RBF). However, when comparing the RF (with MENN) and DNNE (with MENN), no significant difference was noted for subjects, but the two were significantly statistically different for samples.

Table 8 Improvement with the MENN algorithm for the “Test_Data” set

Discussion

Herein, a classification algorithm for PD diagnosis, termed PD_MEdit_EL, was generated by combining a multi-edit-nearest neighbor (MENN) algorithm and ensemble learning algorithm. While research pertaining to PD classifications based on speech samples has been performed, few studies have utilized the Sakar et al. dataset [1]. Since this is a larger and more recently deposited dataset, it was chosen to examine the effectiveness of this new classification strategy. While present classification algorithms examine feature extraction, feature selection/transformation and classifier design, they are unable to consider sample optimization via selection. According to the theory of machine learning, the instance selection is very crucial to the quality of machine learning and in accurately obtaining a final classification, with outliers adding ‘noise’ that can mislead the classifier. In the present study, a sample selection algorithm, MENN, is utilized to reduce the impact of these outliers. Subsequent experimentation showed that the proposed algorithm can provide accurate classifications. While the RF method has been utilized in previous studies, a significant improvement was noted when it was combined with the MENN algorithm (Tables 9, 10).

Table 9 Establishing method significance using the “Training_Data” set
Table 10 Establishing method significance using the “Test_Data” set

Additionally, two ensemble learning algorithms, RF and DNNE, were utilized to enhance the classification accuracy and stability of classification for PD. While the RF algorithm has been more extensively examined, the DNNE is a new classification algorithm. The parameters for the two algorithms were based on prior knowledge and statistical experiments and the results showed that the two algorithms performed better than the algorithm from Ref. [1], even without the MENN algorithm. When combining the ensemble learning algorithms and the MENN algorithm, the classification performance was further improved in both the “Training_Data” (~30%) and “Test_Data” (~25%) sets. The two Sakar et al. [1] datasets were studied systematically in terms of ACC, SEN and SPE, to include examining samples and subjects. To establish statistical effectiveness of the proposed PD_MEdit_EL algorithm, significance differences between algorithms were examined, with some highly significant differences noted between the PD_MEdit_EL algorithm and the other examined algorithms.

Overall, PD_MEdit_EL algorithm achieved a higher classification performance when compared with examined algorithms currently utilized for classification. This could in part be attributed to the PD_MEdit_EL method processing the data, which the existing examined algorithms do not. Such processing techniques include feature selection, feature nonlinear transformation, multiple-classifier and so on. Thus the combining of these methods makes the PD_MEdit_EL algorithm effective method for classifying PD.

The main contributions and innovations of this paper can be described as follows:

  1. 1.

    Speech samples were optimized utilizing the MENN algorithm, thus improving PD classification accuracy.

  2. 2.

    The ensemble learning algorithms, RF and DNNE, were examined in conjunction with MENN and resulted in further optimized PD classifications.

  3. 3.

    The “Traning_Data” and “Test_Data” sets were studied and verified systematically, while the relevant studies did not involve the “Test_Data” set.

  4. 4.

    The statistical significance of the proposed PD_MEdit_EL algorithm was confirmed when compared to the other examined algorithms.

Conclusions

While speech based PD classifications have been shown to be effective, the existing methods lack the ability to optimize speech samples, which is crucial for improving PD classification performance. In this study, a new classification algorithm (PD_MEdit_EL) was proposed that combines the MENN algorithm and an ensemble learning algorithm. To evaluate the performance of the proposed algorithm, a public dataset provided by Sakar et al. was utilized, with the proposed algorithm compared to existing methods. Experimental results showed that the PD_MEdit_EL algorithm provides improved classification abilities when compared with the examined algorithms. Furthermore, the noted improvement was apparent irrespective of the LOO or LOSO or of the dataset utilized, with an improvement near 30% seen for the “Training_Data” set and 25% for the “Test_Data”. Overall, this study provides new insight that can be applied to subsequent research pertaining to PD classifications when utilizing speech data. In the near future, we will consider examining compressed speech feature data to further verify and possibly modify the PD_MEdit_EL algorithm.

Abbreviations

PD:

Parkinson disease

MENN:

multi-edit-nearest-neighbor

RF:

random forest

DNNE:

neural network ensembles

PWP:

people with Parkinson’s

PCA:

principle component analysis

NN:

neural network

SVM:

support vector machine

PNN:

probabilistic neural network

EL:

ensemble learning

RVFL:

functional connection network LOO

CV:

leave one out cross validation

LOSO:

leave-one-subject-out

References

  1. Sakar BE, Isenkul ME, Sakar CO, Sertbas A, Gurgen F, Delil S, Apaydin H, Kursun O. Collection and analysis of a Parkinson speech dataset with multiple types of sound recordings. IEEE J Biomed Health Inform. 2013;17(4):828–34.

    Article  Google Scholar 

  2. Tsanas A, Little MA, McSharry PE, Spielman J, Ramig LO. Novel speech signal processing algorithms for high-accuracy classification of Parkinson’s disease. IEEE Trans Biomed Eng. 2012;59(5):1264–71.

    Article  Google Scholar 

  3. de Rijk MC, Launer LJ, Berger K, Breteler MM, Dartigues JF, Baldereschi M, Fratiglioni L, Lobo A, Martinez-Lage J, Trenkwalder C, Hofman A. Prevalence of Parkinson’s disease in Europe: a collaborative study of population-based cohorts. Neurologic Diseases in the Elderly Research Group. Neurology. 2000;54(5):21–3.

    Article  Google Scholar 

  4. Van Den Eeden SK, Tanner CM, Bernstein AL, Fross RD, Leimpeter A, Bloch DA, Nelson LM. Incidence of Parkinson’s disease: variation by age, gender, and race/ethnicity. Am J Epidemiol. 2003;157(11):1015–22.

    Article  Google Scholar 

  5. Elbaz A, Bower JH, Maraganore DM, McDonnell SK, Peterson BJ, Ahlskog JE, Schaid DJ, Rocca WA. Risk tables for Parkinsonism and Parkinson’s disease. J Clin Epidemiol. 2002;55(1):25–31.

    Article  Google Scholar 

  6. O’Sullivan SB, Schmitz TJ. “Parkinson disease,” in physical rehabilitation. 5th ed. Philadelphia: F. A. Davis Company; 2007. p. 856–94.

    Google Scholar 

  7. Khan T. Parkinson’s disease assessment using speech anomalies: a review. Idt.mdh.se. 2014.

  8. Das R. A comparison of multiple classification methods for diagnosis of Parkinson disease. Expert Syst Appl. 2010;37(2):1568–72.

    Article  Google Scholar 

  9. Hirschauer TJ, Adeli H, Buford JA. Computer-aided diagnosis of Parkinson’s disease using enhanced probabilistic neural network. J Med Syst. 2015;39(11):1–12.

    Article  Google Scholar 

  10. Prashanth R, Roy SD, Mandal PK, Ghosh S. Parkinson’s disease detection using olfactory loss and REM sleep disorder features. IEEE Eng Med Biol Soc. Annual Conference; 2014. pp. 5764–67.

  11. Baghai-Ravary L, Beet SW. Automatic speech signal analysis for clinical diagnosis and assessment of speech disorders. Springer Br Electr Comput Eng. 2012;115(2):31–6.

    Google Scholar 

  12. Kim H, Lee HJ, Lee W, Kwon S, Kim SK, Jeon HS, Park H, Shin CW, Yi WJ, Jeon BS, Park KS. Unconstrained detection of freezing of Gait in Parkinson’s disease patients using smart phone. Annual International Conference of the IEEE engineering in medicine and biology society. IEEE engineering in medicine and biology society. Annual Conference; 2015. pp. 3751–54.

  13. Little MA, McSharry PE, Hunter EJ, Spielman J, Ramig LO. Suitability of Dysphonia measurements for telemonitoring of Parkinson’s disease. IEEE Trans Biomed Eng. 2009;56(4):1015–22.

    Article  Google Scholar 

  14. Hariharan M, Polat K, Sindhu R. A new hybrid intelligent system for accurate detection of Parkinson’s disease. Comput Methods Programs Biomed. 2014;113(3):904–13.

    Article  Google Scholar 

  15. Sapir S, Ramig LO, Spielman JL, Fox C. Formant centralization ratio: a proposal for a new acoustic measure of dysarthric speech. J Speech Lang Hear Res. 2010;53(1):114–25.

    Article  Google Scholar 

  16. Rusz J, Bonnet C, Klempíř J, Tykalová T, Baborová E, Novotný M, Ulseh A, Růžička E. Speech disorders reflect differing pathophysiology in Parkinson’s disease, progressive supra nuclear palsy and multiple system atrophy. J Neurol. 2010;262(2):992–1001.

    Google Scholar 

  17. Islam MS, Parvez I, Deng H, Goswami P. Performance comparison of heterogeneous classifiers for detection of Parkinson’s disease using voice disorder (dysphonia). 2014 International Conference on informatics, electronics and vision (ICIEV); 2014. pp. 1–7.

  18. Pérez CJ, Naranjo L, Martín J, Campos-Roca Y. A latent variable-based Bayesian regression to address recording replications in Parkinson’s disease. 22nd European signal processing conference; 2014. pp. 1447–51.

  19. Athanasios T, Little MA, Mcsharry PE, Ramig LO. Accurate telemonitoring of Parkinson’s disease progression by noninvasive speech tests. IEEE Trans Biomed Eng. 2010;57(4):884–93.

    Article  Google Scholar 

  20. Benba A, Jilbab A, Hammouch A, Sandabad S. Voiceprints analysis using MFCC and SVM for detecting patients with Parkinson’s disease. 2015 International Conference on electrical and information technologies (ICEIT), IEEE; 2015.

  21. Titze, Ingo R, Martin, Daniel W. Principles of voice production. 2nd ed. Iowa City: Natl. Center Voice Speech; 2000.

    Google Scholar 

  22. Schoentgen J, Guchteneere RD. Time series analysis of jitter. J Phon. 1995;23(23):189–201.

    Article  Google Scholar 

  23. Tsanas A, Little MA, McSharry PE, Ramig LO. Ramig. nonlinear speech analysis algorithms mapped to a standard metric achieve clinically useful quantification of average Parkinson’s disease symptom severity. J Royal Soc Interface. 2011;59(8):842–55.

    Article  Google Scholar 

  24. Yair E, Gath I. High resolution pole-zero analysis of parkinsonian speech. IEEE Trans Biomed Eng. 1991;38(2):161–7.

    Article  Google Scholar 

  25. Orozco-Arroyave JR, Hönig F, Arias-Londoño JD, Vargas-Bonilla JF, Daqrouq K, Skodda S, Rusz J, Nöth E. Automatic detection of Parkinson’s disease in running speech spoken in three different languages. J Acoust Soc Am. 2016;139(1):481–500.

    Article  Google Scholar 

  26. Ozkan H. A comparison of classification methods for telediagnosis of Parkinson’s disease. Entropy. 2016;18(4):115.

    Article  MathSciNet  Google Scholar 

  27. Mahnaz B, Ashkan S. A multiple-classifier framework for Parkinson’s disease detection based on various vocal tests. Int J Telemed Appl. 2016: 1–9.

  28. Froelich W, Wróbel K, Porwik P. Diagnosing Parkinson’s disease using the classification of speech signals. J Med Inform Technol. 2014;23:1642–6037.

    Google Scholar 

  29. Sharma RK, Gupta AK. Processing and analysis of human voice for assessment of Parkinson disease. J Med Imaging Health Inform. 2016;6(1):63–70.

    Article  Google Scholar 

  30. Ozcift A. SVM feature selection based rotation forest ensemble classifiers to improve computer-aided diagnosis of Parkinson disease. J Med Syst. 2012;36(2):2141–7.

    Article  Google Scholar 

  31. Yang S, Zheng F, Luo X, Cai S, Wu Y, Liu K, Wu M, Chen J, Krishnan S. Effective dysphonia detection using feature dimension reduction and kernel density estimation for patients with Parkinson’s disease. PLoS ONE. 2014;9(2):1–10.

    Google Scholar 

  32. Shahbakhti M, Taherifar D, Sorouri A. Linear and non-linear speech features for detection of Parkinson’s disease. The 2013 biomedical engineering international conference (BMEiCON-2013); 2013. pp. 1–3.

  33. Avci1 D, Dogantekin A. An Expert diagnosis system for Parkinson disease based on genetic algorithm-wavelet kernel-extreme learning machine. Parkinson’s disease; 2016: 1–9.

  34. Mekyska J, Rektorova I, Smekal Z. Selection of optimal parameters for automatic analysis of speech disorders in Parkinson’s disease. Telecommunications and signal processing (TSP), 2011 34th International Conference on 18–20 Aug; 2011. pp. 408–12.

  35. Sakar CO, Kursun O. Telediagnosis of Parkinson’s disease using measurements of dysphonia. J Med Syst. 2010;34(4):591–9.

    Article  Google Scholar 

  36. Galaz Z, Mekyska J, Mzourek Z, Smekal Z, Rektorova I, Eliasova I, Kostalova M, Mrackova M, Berankova D. Department prosodic analysis of neutral, stress-modified andrhymed speech in patients with Parkinson’sdisease. Comput Methods Programs Biomed. 2016;127:301–17.

    Article  Google Scholar 

  37. Su M, Chuang KS. Dynamic feature selection for detecting Parkinson’s disease through voice signal. IEEE Mtt-S 2015 International microwave workshop series on Rf and wireless technologies for biomedical and healthcare applications. IEEE; 2015.

  38. Khan T, Westin J, Dougherty M. Classification of speech intelligibility in Parkinson’s disease. Biocybern Biomed Eng. 2013;33(4):35–45.

    Google Scholar 

  39. Frid A, Safra EJ, Hazan H, Lokey LL,Hilu D, Manevitz L, Ramig LO, Sapir S. Computational diagnosis of Parkinson’s disease directly from natural speech using machine learning techniques. IEEE International Conference on software science, technology and engineering; 2014. pp. 50–3.

  40. Kaya E, Findik O, İsmail B, Arslan A. Effect of discretization method on the diagnosis of Parkinson’s disease. Int J Innov Comput Inform Control Ijicic. 2011;7(8):4669–78.

    Google Scholar 

  41. Benba A, Jilbab A, Hammouch A. Hybridization of best acoustic cues for detecting persons with Parkinson’s disease. World Conference on Complex Systems; 2014. pp. 622–5.

  42. Islam MS, Parvez I, Deng H, Goswami P. Performance comparison of heterogeneous classifiers for detection of Parkinson’s disease using voice disorder (dysphonia). International Conference on informatics, electronics and vision. IEEE; 2014. pp. 1–7.

  43. Naranjo L, Pérez CJ, Campos-Roca Y, Martín J. Addressing voice recording replications for Parkinson’s disease detection. Expert Syst Appl. 2016;46(C):286–92.

    Article  Google Scholar 

  44. Spadoto AA, Guido RC, Papa JP, Falcao AX. Parkinson’s disease identification through optimum-path forest. International Conference of the IEEE engineering in medicine and biology society. Conference Proceeding of IEEE engineering medicine biology society; 2010. pp. 6087–90.

  45. Froelich W, Wrobel K, Porwik P. Diagnosis of Parkinson’s disease using speech samples and threshold-based classification. J Med Imaging Health Inform. 2015;5(6):1358–63.

    Article  Google Scholar 

  46. Chung JW. Performance Comparison of algorithm through classification of Parkinson’s disease according to the speech feature. J Korea Multimed Soc. 2016;19(2):209–14.

    Article  Google Scholar 

  47. Gok M, Murat K. An ensemble of k-nearest neighbours algorithm for detection of Parkinson’s disease. Int J Syst Sci. 2015;46(6):1108–12.

    Article  MATH  Google Scholar 

  48. Mandal I, Sairam N. New machine-learning algorithms for prediction of Parkinson’s disease. Int J Syst Sci. 2014;45(3):647–66.

    Article  MathSciNet  MATH  Google Scholar 

  49. Chen H-L, et al. An efficient diagnosis system for detection of Parkinson’s disease using fuzzy k-nearest neighbor approach. Expert Syst Appl. 2013;40(1):263–71.

    Article  Google Scholar 

  50. Juan LI, Wang YPA. Fast Neighbor prototype selection algorithm based on local mean and class global information. Zidonghua Xuebao/acta Automatica Sinica. 2014;40(6):1116–25.

    Google Scholar 

  51. Alhamdoosh M, Wang D. Fast decorrelated neural network ensembles with random weights. Inf Sci. 2014;264(6):104–17.

    Article  MathSciNet  MATH  Google Scholar 

  52. Hattori K, Takahashi M. A new nearest-neighbor rule in the pattern classification problem. Pattern Recogn. 1999;32(3):425–32.

    Article  Google Scholar 

  53. Rico Juan JR, Iñesta JM. New rank methods for reducing the size of the training set using the nearest neighbor rule. Pattern Recogn Lett. 2012;33(5):654–60.

    Article  Google Scholar 

  54. Piyush R, Ramakrishnan S. Diffusion tensor based Alzheimer image analysis using region specific volume features and random forest classifier. Int Conf Biomed Eng. 2013;43:691–4.

    Article  Google Scholar 

Download references

Authors’ contributions

HZ and YL conceived of the whole study, and participated in its design and coordination and helped to draft the manuscript. LY, YL, XZ and MQ participated in the measurements of all subjects and drafted the complete manuscript. PW, JY and FY managed the trials and assisted with writing discussions in the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Authors thank the professor Yanling Zhang and Gui Li for their valuable suggests.

Competing interests

The authors declare that they have no competing interests.

Availability of data and supporting materials

The data utilized in this study was obtained from the Parkinson Speech Dataset with Multiple Types of Sound Recordings, was deposited by Sakar et al. and is available on the University of California, Irvine (UCI) machine learning dataset repository website. The URL is as follows:

http://archive.ics.uci.edu/ml/datasets/Parkinson+Speech+Dataset+with++Multiple+Types+of+Sound+Recordings.

Funding

This research is funded by National Natural Science Foundation of China NSFC (No: 61108086, 61171089, 91438104, 1304382), Basic and Advanced Research Project in Chongqing (cstc2016jcyjA0043, cstc2016jcyjA0134,cstc2016jcyjA0064), Chongqing Social Undertaking and People’s Livelihood Guarantee Science and Technology innovation Special Foundation (cstc2016shmszx40002, cstc2016shmszx130014, cstc2016shmszx0099), Fundamental Research Funds for the Central Universities (CDJZR155507), the China Postdoctoral Science Foundation (2013M532153), General fund in Sichuan province(14ZB0303), General fund of Dazhou science and technology bureau (KJJ201401), the Chongqing Postdoctoral Science Special Foundation of China, and The Ministry of education to return personnel research start fund; the Youth Training Project of Army Medical Technology (13QNP120).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yongming Li.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, HH., Yang, L., Liu, Y. et al. Classification of Parkinson’s disease utilizing multi-edit nearest-neighbor and ensemble learning algorithms with speech samples. BioMed Eng OnLine 15, 122 (2016). https://doi.org/10.1186/s12938-016-0242-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12938-016-0242-6

Keywords