Next Article in Journal
In Situ Flow Cytometer Calibration and Single-Molecule Resolution via Quantum Measurement
Previous Article in Journal
A Tool for High-Resolution Volumetric Optical Coherence Tomography by Compounding Radial-and Linear Acquired B-Scans Using Registration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detecting Bulbar Involvement in Patients with Amyotrophic Lateral Sclerosis Based on Phonatory and Time-Frequency Features

1
CIMNE, Building C1, North Campus, UPC, Gran Capità, 08034 Barcelona, Spain
2
Department of Computer Science & INSPIRES, University of Lleida, Jaume II 69, 25001 Lleida, Spain
3
Neurology Department, Hospital Universitari de Bellvitge, L’Hospitalet de Llobregat, 08907 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(3), 1137; https://doi.org/10.3390/s22031137
Submission received: 27 October 2021 / Revised: 20 January 2022 / Accepted: 28 January 2022 / Published: 2 February 2022
(This article belongs to the Section Intelligent Sensors)

Abstract

:
The term “bulbar involvement” is employed in ALS to refer to deterioration of motor neurons within the corticobulbar area of the brainstem, which results in speech and swallowing dysfunctions. One of the primary symptoms is a deterioration of the voice. Early detection is crucial for improving the quality of life and lifespan of ALS patients suffering from bulbar involvement. The main objective, and the principal contribution, of this research, was to design a new methodology, based on the phonatory-subsystem and time-frequency characteristics for detecting bulbar involvement automatically. This study focused on providing a set of 50 phonatory-subsystem and time-frequency features to detect this deficiency in males and females through the utterance of the five Spanish vowels. Multivariant Analysis of Variance was then used to select the statistically significant features, and the most common supervised classifications models were analyzed. A set of statistically significant features was obtained for males and females to capture this dysfunction. To date, the accuracy obtained (98.01% for females and 96.10% for males employing a random forest) outperformed the models in the literature. Adding time-frequency features to more classical phonatory-subsystem features increases the prediction capabilities of the machine-learning models for detecting bulbar involvement. Studying men and women separately gives greater success. The proposed method can be deployed in any kind of recording device (i.e., smartphone).

1. Introduction

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease with an irregular and asymmetric progression, characterized by a progressive loss of both upper and lower motor neurons and that leads to muscular atrophy, paralysis and death, mainly from respiratory failure. The life expectancy of patients with ALS is between 3 and 5 years from the onset of symptoms.
ALS causes muscle weakness and movement, speech, eating and respiratory impediments, leaving the patient reliant on caretakers and relatives and causing considerable social costs. Currently, there is no cure for ALS, although early detection can lead to the use of more appropriate therapies that may slow progress [1].
When the disease starts in the arms and legs, it is called spinal ALS (limb or spinal onset; 80% of cases), and when it starts in the cranial nerve nuclei, it is called bulbar ALS (bulbar onset; 20%). The bulbar muscle is responsible for speech and swallowing, so patients with the later variant have a shorter life expectancy. However, dysarthria, or slurred or difficult speech articulation, affects 80% of all ALS patients [2]. In bulbar ALS, these symptoms usually appear at the onset of the disease, while in spinal ALS, they appear later. Early detection of bulbar involvement in those with ALS is crucial for better diagnosis and prognosis, and could be the key to effectively slowing the development of the disease.
The authors in [3,4,5] demonstrated that the deterioration of the bulbar muscle affected some phonatory-subsystem features. Among these were jitter, shimmer, harmonic-to-noise ratio (HNR), pitch, formant trajectories, correlations of formants with articulation patterns, fractal jitter, and Mel Frequency Cepstral Coefficients. Consequently, as suggested in previous works [6,7,8,9,10,11,12,13], imperceptible changes in speech and voice can be detected through objective measures.
Time-frequency representation (TFR), broadly applied to detecting several conditions [14,15,16,17,18], has been recently used to detect pathological changes in voice signals [19]. TFR enables the evolution of the periodicity and frequency components to be observed over time, allowing the analysis of non-stationary signals, such as voice signals [20]. The spectrogram is the most common TFR for the analysis of audio signals. This representation corresponds to Cohen’s class of time-frequency energy distributions in general. The depiction of a spectrogram is not optimal in terms of resolution quality. There are Cohen class representations that have greater resolution quality. They are all made by smoothing the Wigner distribution, which has the finest resolution but the most detrimental interference. The smoothing functions chosen strike a balance between resolution quality and the elimination of detrimental interference terms.
The authors in [21] used TFR representations from the Cohen class for the onset signal of the vocal fold to diagnose various phonation problems induced by pathological alterations. To assess the voice signal, Cohen class TFRs were combined with a cone kernel distribution to provide optimum smoothness across time. The authors demonstrated that even minor pathogenic alterations in the vocal folds can be seen in TFR, allowing for sensitive affection detection and diagnosis.
In ALS, voice and speech impairment can occur up to 3 years before a diagnosis [22], and when the bulbar muscle function is damaged, voice and speech deteriorates significantly as the disease advances [23]. Features obtained from Cohen class TFRs could aid in the identification of bulbar involvement even earlier than human hearing can.
Centering attention on the subject at hand, R. Norel et al. [24] developed machine-learning models that recognize the presence and severity of ALS based on a variety of frequency, spectral, and voice quality characteristics. An et al. [25] employed Convolutional Neural Networks to classify the intelligible speech of ALS patients compared to healthy people. Finally, Gutz et al. [26] combined SVM and feature filtering techniques.
Based on previous works, and starting from our recent studies [6,18], our paper suggests using phonatory-subsystem [6] and time-frequency [18] features jointly. This also is our hypothesis and main contribution. These features, extracted from a portion of the five Spanish vowels, could enhance the performance of the classification models for the early detection of bulbar involvement, for which the main goals (and contributions) of this research are:
1.
To design a new methodology for the automatic detection of bulbar involvement in males and females based on phonatory-subsystem and time-frequency features.
2.
To obtain a set of statistically significant features for diagnosing bulbar involvement efficiently.
3.
To analyze the performance of the most common supervised classification models to improve the diagnosis of bulbar involvement.

2. Methods

2.1. Participants

Of the 65 participants selected for this study, 14 of those with ALS had been diagnosed with bulbar involvement (11 females and 3 males; mean = 56.8 years, standard deviation = 12.3 years), 33 had ALS but had not been diagnosed with bulbar involvement (8 females and 25 males; mean = 57.6 years, standard deviation = 12.0 years) and 18 were healthy individuals (9 females and 9 males; mean = 45.2 years, standard deviation = 12.2 years). The main clinical records of the ALS participants are summarized in Table 1. It can be seen that the sample is well age-balanced.
The ALS patients’ voices were checked by a multidisciplinary clinical team and finally selected by a neurologist for this study.
The control subjects were recruited through personal advertisements in the hospital facilities by the researchers involved in this study. After contacting the volunteers, they received an information sheet explaining the procedure and goal of the study as well as the exclusion criteria. They were interviewed through a questionnaire and those who did not report any voice issue or relevant previous condition were selected for the study.
The control subjects were recruited through personal advertisements conducted in the hospital facilities by the researchers involved in this study. Most of them were companions of ALS patients. After contacting them, control subjects received an information sheet explaining the procedures and goals of the study as well as the exclusion criteria. Control subjects were informed that the study focused on voice analysis to distinguish bulbar involvement in ALS patients. They were interviewed through a questionnaire. Those who did not report any voice issue or relevant previous condition were selected for the study. When they were eligible and still willing to participate, they were invited to come to the hospital room where the voice samples were registered.

2.2. Vowel Recording

There are five vowel segments in the Spanish phonological system (a, e, i, o, u). These were obtained and analyzed from each ALS patient, all of whom were Spanish speakers.
Under medium vocal loudness conditions, each participant uttered a sustained sample of each Spanish vowel for 3–4 s. The recordings were made in a standard hospital room using a laptop and a USB EMITA Streaming GXT 252 microphone calibrated for dBSPL. It has a sensitivity of 35 dBSPL and a maximum sound pressure level of 135 dBSPL. The participants sat on a chair with the microphone positioned approximately 30 centimeters from their mouths. The voice signals were recorded using Audacity, an open-source application [27], at a sampling rate of 44.100 Hz and 32-bit quantization.
A visual inspection of the spectrograms of the voice signals was conducted similarly to the procedure in [28] to analyze the signal type of the participants’ voices. Their results suggested four voice types, of which only type 1 and type 2 were considered suitable for perturbation analysis.
In this study, all the control subjects presented type 1 voice signals, which were periodic without strong modulations or subharmonics. They showed multiple clearly and nearly straight defined harmonics.
Among the 14 ALS patients with bulbar involvement, 10 patients presented type 1 voice signals, which were nearly periodic and showed some clearly defined harmonics. However, a small amount of noise was observed in some voices (four of them). Four of the ALS patients with bulbar involvement presented type 2 voice signals. These had some strong modulations and subharmonics, yet still presented stable and periodic segments in their voices.
Among the 33 ALS patients without bulbar involvement, 29 presented type 1 voice signals, which were nearly periodic and showed multiple or at least some clearly defined harmonics. Instead, four of them presented type 2 voice signals with some strong modulations and subharmonics but still with stable and periodic segments.
It was observed that most of the information of the signal recordings was contained in the range from 0 to 4000 Hz. Therefore, it was decided to decimate all the recording signals sampled at 44.100 Hz using a decimated factor of 5. Signals re-sampled at 8820 Hz were obtained.
Then, each re-sampled signal was standardized by means of the z-score technique. The z-score measures the distance of a signal sample from the mean of the re-sampled signal in terms of the standard deviation. The resulting standardized signal had mean 0 and standard deviation 1, and retained the shape properties of the re-sampled signal. For the re-sampled signal with mean X ¯ and standard deviation S, the z-score of a signal sample x was computed as:
z = ( x X ¯ ) S
Finally, a segment of 150 ms of each re-sampled and standardized signal ( x ( t ) ) was chosen for analysis by tacking the midpoint at the center of the phonation.

2.3. Phonatory-Subsystem Features

A total of 15 features from the phonatory subsystem defined in [6,13] were used. They were computed by means of the standard methods used in Praat [29] and the setting details used were the same as in [6]. These features were:
  • Fundamental period cycle-to-cycle variation (Jitter(absolute), Equation (2)).
    J i t t e r ( a b s o l u t e ) = 1 N 1 i = 1 N 1 | T i T i 1 | ,
    where N is the number of cycles and T i the duration of the ith cycle.
  • Relative period (Jitter(relative), Equation (3)).
    J i t t e r ( r e l a t i v e ) = 1 N 1 i = 1 N 1 | T i T i 1 | 1 N i = 1 N T i × 100
  • Relative perturbation (Jitter(rap), Equation (4)).
    J i t t e r ( r a p ) = 1 N 1 i = 1 N 1 | T i 1 3 n = i 1 i + 1 T n | 1 N i = 1 N T i × 100
  • Five-point period perturbation quotient (Jitter(ppq5), Equation (5)).
    J i t t e r ( p p q 5 ) = 1 N 1 i = 2 N 2 | T i 1 5 n = i 2 i + 2 T n | 1 N i = 1 N T i × 100
  • Variability of the peak-to-peak amplitude (Shimmer(dB), Equation (6)).
    S h i m m e r ( d B ) = 1 N 1 i = 1 N 1 | 20 × l o g A i + 1 A i ,
    where A i is the extracted peak-to-peak amplitude data and N is the number of extracted fundamental periods.
  • Relative amplitudes of consecutive periods (Shimmer(relative), Equation (7)).
    S h i m m e r ( r e l a t i v e ) = 1 N 1 i = 1 N 1 | A i A i + 1 | 1 N i = 1 N A i × 100
  • Three-, five- and eleven-point amplitude perturbation (Shimmer(apqP), Equation (8)).
    S h i m m e r ( a p q P ) = 1 N 1 i = 1 N 1 | A i ( 1 P n = i 1 i + 1 A n ) | 1 N i = 1 N A i × 100 ,
    where P = {3, 5 and 11}.
  • Mean and standard deviation (HNR(mean) and HNR(SD)) of the harmonics-to-noise-ratio (HNR, Equation (9)).
    H N R = 10 × l o g 10 r ( T 0 ) 1 r ( T 0 ) ,
    where r ( T 0 ) is the second local maximum of the normalized auto-correlation function and T 0 is the period of the signal.
  • Mean, standard deviation, minimum and maximum value of the pitch (pitch(mean), pitch(SD), pitch(min) and pitch(max)). See [29] for more details about obtaining the pitch.

2.4. Time-Frequency Features

The methods employed to obtain the time-frequency features were inspired by the previous work, presented in [16,17], and implemented with MATLAB [30].
First, the Wigner distribution (WD) of the real signal x ( t ) of each voice segment was obtained and convoluted with the Choi-Williams exponential function. The resulting Choi-Williams distribution was normalized ( C W D N ( f , t ) ). For more details, see [18].
Then, the joint probability density distribution p D ( f , t ) (Equation (10)) was obtained.
p D ( f , t ) = m t ( t ) · m f ( f ) ,
where m t ( t ) , instantaneous power, and m f ( f ) , spectral energy density, are the marginal density functions of C W D N ( f , t ) .
According to Equation (10), p D can be only computed as the product of the marginal density functions m t ( t ) and m f ( f ) (of C W D N ) if they are statistically independent. To corroborate this assumption, we computed the joint time-frequency moments of the C W D N ( t n f m from n = 1 and m = 1 to n = 15 and m = 15 where n and m are the frequency and time moment orders) of the vowels of all the participants. All of these were 0 or very close to 0. This confirmed the statistical independence of m t ( t ) and m f ( f ) .
p D ( f , t ) is completely free of interference and negative values. Thus, it is very useful for extracting time-frequency features for classification.
Figure 1 shows the comparison of the p D ( f , t ) of the vowel “a” from three different patients. Non-undesirable effects were observed in the p D ( f , t ) . Figure 1a corresponds to a patient without bulbar involvement. The p D ( f , t ) shows a voice rich in harmonics. Figure 1b shows the p D ( f , t ) of the vowel “a” of a patient diagnosed with slight bulbar involvement. Significant differences can be observed. Voice harmonics appear attenuated. Figure 1c shows the p D ( f , t ) of an even more extreme case, diagnosed with severe bulbar involvement. It can be seen that its voice harmonics appear even more attenuated.The visual appraisal of these figures clearly shows the significant differences in the p D ( f , t ) between ALS patients with and without bulbar involvement.
From the p D ( f , t ) , a set of 30 features per vowel was obtained. Twenty-one features were computed by dividing the spectrum (0–4410 Hz) into 7 frequency bands. These were 1, 0–80 Hz; 2, 80–250 Hz; 3, 250–550 Hz; 4, 550–900 Hz; 5, 900–1500 Hz; 6, 1500–3000 Hz; 7, 3000–4410 Hz. These bands were selected to capture the differences observed in the time-frequency representations of the two groups of ALS patients by means of the visual appraisal of p D ( f , t ) in the range of these frequency bands. These features were:
  • Average instantaneous spectral energy (E(t), Equation (11)) for each frequency band (E_Bn1E_Bn7).
    E ( t ) = f 1 f 2 p D ( f , t ) d f ,
    where f 1 and f 2 are the lower and upper frequencies of each band.
  • Instantaneous frequency peak ( f _ C r e s ( t ) , Equation (12)) for each frequency band (f_Cres1f_Cres7).
    f _ C r e s ( t ) = 1 E ( t ) a r g m a x f f 1 f 2 f · p D ( f , t )
  • Average instantaneous frequency ( f m i ( t ) , Equation (13)) of the spectrum for each frequency band (f_Med1f_Med7).
    f m i ( t ) = f 1 f 2 1 E ( t ) f · p D ( f , t ) d f
10 additional features were added:
  • Instantaneous (H_t, Equation (14)) and spectral (H_f, Equation (15)) information entropies. Furthermore, the joint Shannon entropy (H_tf, Equation (16)) was also used.
    H _ t = l o g 2 ( m t N ( t ) ) · m t N ( t ) d t ,
    where m t N ( t ) is the quantified instantaneous marginal obtained from the m t ( t ) and m f N ( f ) is the quantified frequency marginal obtained from the m f ( f ) .
    H _ f = l o g 2 ( m f N ( f ) ) · m f N ( f ) d f
    H _ t f = H _ t + H _ f
  • Spectral information (IE(f), Equation (17)), for each frequency band (IE_Bn1IE_Bn7).
    I E ( f ) = l o g 2 ( m f N ( f ) )
  • Kurtosis (K, Equation (18)).
    K = m t ( t ) n m f ( f ) m ,
    where n = 4 and m = 0 .
  • Joint time-frequency moment ( t n f m ,18,31]) where n and m (n, m = 1, 7, 15) are the frequency and time moment orders, i.e., the following time-frequency moments were used: t 1 f 1 , t 7 f 7 and t 15 f 15 .

2.5. Feature Selection

From a total of 65 participants, 18 were labelled C (healthy group), 14 were labelled B (ALS patients with bulbar involvement) and 33 were labelled NB (ALS patients without bulbar involvement). Furthermore, every ALS participant was labelled A.
Accordingly, four classification problems were analyzed, males and females being studied separately, these being C vs. B, C vs. NB, B vs. NB and C vs. A.
The Multivariant Analysis of Variance (MANOVA), which uses the covariance between the features in testing the statistical significance of the mean differences, was performed in IBM SPSS Statistics [32] to select a subset of relevant features for use in constructing the classification model for these four cases. This procedure made it possible to contrast the null hypothesis in the features obtained.
To perform this statistical analysis, it was assumed that the features had a multivariable normal distribution, and no assumptions were made regarding the homogeneity of the variance or the correlation between the features. A significance value of p-value < 0.05 was considered sufficient to assume the existence of feature differences between the four groups analyzed.

2.6. Classification Models

Several supervised classification models were implemented in R [33] to measure the classification performance. These models were fitted with the features selected. These were standardized by subtracting the mean and centered at 0. Ten-fold cross-validation was implemented in R using the caret package to draw suitable conclusions. This consisted of dividing the dataset into 10 contiguous chunks, each containing approximately the same number of samples, and then performing 10 training-testing experiments as follows: for each chunk i { 1 , 2 , , 10 } , the current chunk was retained for testing the model and training was performed on the remaining 9 chunks, recording the results. The average performance of the 10 training-testing experiments was finally provided.
The upsampling technique with replacement was applied to the training data by making the group distributions equal to deal with the unbalanced dataset, which could bias the classification models [34].
The supervised models with classification thresholds of 50% were built in R [33]. In binary classification problems, the classification threshold is a value that converts the model prediction to positive or negative depending on whether the prediction is above or below the threshold.
The classification algorithms used were the most popular ones in ALS: Support Vector Machine (SVM), Neural Networks (NN), Linear Discriminant Analysis (LDA), Logistic Regression (LR) and Random Forest (RF). For more details, see [6].

2.7. Model Validation Metrics

There are various metrics for evaluating classification models [35]. The foremost among these, accuracy, sensitivity and specificity, were used to evaluate the performance of the classification models.

3. Results

First, the significant features from the four cases (C vs. B, C vs. NB, B vs. NB and C vs. A) were selected. Then, the performance of the classification models was evaluated.

3.1. Selecting the Significant Features

From the 50 features obtained, the MANOVA analysis was applied to select those that were statistically significant. Four comparisons were analyzed separately for males and females: C vs. B, C vs. NB, B vs. NB and C vs. A. Features not showing statistical significance (p-value ≥ 0.05) were discarded.
Table 2 shows the significant features obtained for males. In the C vs. B case, this was a set of 12 statistically (half phonatory) significant features (p-value < 0.05); in C vs. NB, there were 13 (10 of them phonatory); in B vs. NB, 9 (all time-frequency); and in C vs. A, 12 (10 of which were phonatory).
For females (Table 3), in the C vs. B case, a set of 20 statistically significant features (p-value < 0.05) was obtained (13 out of 20 were phonatory). In the C vs. NB case, a set of 10 statistically significant features was obtained (6 of them, phonatory). In the case B vs. NB, a set of 14 statistically significant features was obtained (12 of which were phonatory). In the C vs. A case, 20 statistically significant features were obtained (12 being phonatory).

3.2. Classification Models

The classification models were fitted with the significant features selected in Section 3.1. Table 4 and Table 5 show the classification performance for males and females, respectively. The results are presented for the accuracy, sensitivity and specificity of the models used for the four cases.
For males in C vs. B case, all the classifiers generally performed well. RF obtained the best accuracy, 96.1 % . For LDA and NN, accuracy was 95.0 % and for SVM and LR, 93.3 % and 91.9 % respectively. LR gave the best sensitivity ( 95.0 % ), and RF and LDA the best s p e c i f i c i t y = 97.5 % .
Similar performance was achieved in C vs. NB and C vs. A cases. In these, SVM was the best model (an accuracy of 93.1 % was reached for C vs. NB and 92.6 % for C vs. A).
Otherwise, the outcomes worsened in B vs. NB compared with the other cases. Despite RF obtained the best accuracy ( 91.8 % ), the sensitivity it achieved was the worst.
For females, in the C vs. B case, the results also indicate that the performance of all classifiers was excellent. RF gave the best accuracy, 98.1 % , sensitivity, 96.6 % , and specificity, 100 % .
Similar behavior was obtained in the C vs. NB and C vs. A cases. In these, RF was also the best model (obtaining accuracy of 94.1 % and 95.8 % for C vs. NB and C vs. A respectively). In both cases, LDA achieved the best specificity.
Meanwhile, the results were worse in B vs. NB compared with the other cases. Although RF obtained the best accuracy at 84.8 % , the outcomes obtained with it for specificity and especially sensitivity were very low.
In general, the best model was RF. Special attention should be paid to female outcomes. Poor results were obtained for both genders in the B vs. NB case.

4. Discussion

4.1. Principal Findings

The results obtained demonstrate that it is possible to diagnose bulbar involvement using supervised gender-specific models fitted to the significant phonatory and time-frequency features.
Table 2. Significant Features for males.
Table 2. Significant Features for males.
ComparisonFeaturep-Value
C vs. Bshimmer(dB)0.039
shimmer(apq11)<0.001
pitch(mean)0.001
pitch(SD)0.023
pitch(min)0.016
pitch(max)<0.001
f_Cres20.046
f_Cres60.046
f_Med2<0.001
f_Med60.008
K0.027
t 1 f 1 0.002
C vs. NBjitter(relative)0.008
shimmer(dB)0.001
shimmer(relative)0.008
shimmer(apq3)0.035
shimmer(apq11)<0.001
pitch(mean)0.001
pitch(SD)0.002
pitch(min)0.023
pitch(max)0.001
HNR(mean)0.037
IE_Bn10.045
H_tf0.015
H_f0.045
B vs. NBf_Cres10.044
f_Cres20.028
f_Med2<0.001
f_Med60.011
f_Med70.024
H_tf0.009
H_f0.009
K0.045
t 1 f 1 <0.001
C vs. Ajitter(relative)0.009
shimmer(dB)0.001
shimmer(relative)0.009
shimmer(apq3)0.044
shimmer(apq11)<0.001
pitch(mean)0.001
pitch(SD)0.002
pitch(min)0.015
pitch(max)<0.001
HNR(mean)0.046
H_tf0.048
t 1 f 1 0.034
Table 3. Significant Features for females.
Table 3. Significant Features for females.
ComparisonFeaturep-Value
C vs. Bjitter(relative)0.001
jitter(absolute)<0.001
jitter(rap)<0.001
jitter(ppq5)<0.001
shimmer(relative)<0.001
shimmer(dB)<0.001
shimmer(apq3)<0.001
shimmer(apq5)<0.001
shimmer(apq11)<0.001
pitch(mean)<0.001
pitch(SD)<0.001
pitch(max)<0.001
HNR(mean)<0.001
f_Cres20.004
f_Cres60.029
f_Cres70.020
E_Bn20.003
f_Med2<0.001
f_Med60.013
t 1 f 1 0.028
C vs. NBjitter(absolute)<0.001
shimmer(apq11)<0.001
pitch(mean)<0.001
pitch(SD)0.003
pitch(min)0.008
pitch(max)<0.001
f_Cres70.011
E_Bn20.015
f_Med10.014
t 7 f 7 0.022
B vs. NBjitter(relative)<0.001
jitter(absolute)<0.001
jitter(rap)<0.001
jitter(ppq5)<0.001
shimmer(relative)<0.001
shimmer(dB)<0.001
shimmer(apq3)<0.001
shimmer(apq5)<0.001
shimmer(apq11)<0.001
pitch(SD)<0.001
pitch(max)0.029
HNR(mean)<0.001
H_tf0.026
H_f0.048
C vs. Ajitter(relative)<0.001
jitter(rap)0.001
jitter(ppq5)0.004
shimmer(relative)<0.001
shimmer(dB)<0.001
shimmer(apq3)<0.001
shimmer(apq5)0.001
shimmer(apq11)<0.001
pitch(mean)<0.001
pitch(SD)<0.001
pitch(max)<0.001
HNR(mean)0.003
f_Cres20.006
f_Cres70.005
E_Bn20.003
f_Med10.049
f_Med20.001
f_Med70.049
H_t0.039
t 1 f 1 0.018
In the case of B vs. C, the a c c u r a c y achieved was up to 98.1% (RF) and 96.1% (RF) for females and males, respectively.
Lower performance was obtained in C vs. NB but this was still higher than expected. The voice performance in C or NB should be similar. This indicates that some participants in the NB group were probably incorrectly diagnosed. This is coherent with [6]. Similarly, the excellent performance achieved in C vs. A suggests that some of the members of A (14 out of 47) have bulbar involvement. Alternatively, although the most stable segments of the voice samples were selected for analysis, many co-articulatory effects could have influenced the results. Moreover, phonatory-subsystem features are subject to inherently large variability, even for Cs.
On the whole, huge uncertainty was observed in the evaluation concerning bulbar involvement among the participants in the NB group. The case of B vs. NB disclosed that the models did not differentiate between the B and NB subject groups as well as they did with the other groups. RF achieved the best overall performance (accuracy = 91.8%) in males. However, the model presented problems for spotting positive cases (sensitivity = 55.0%). In females, RF achieved an accuracy of 84.8%. These values are still far from the ones obtained in the C vs. B case. These outcomes additionally reinforce the idea that NB subjects were misdiagnosed.
The outcomes of each comparison between groups depend on the significant features chosen (between phonatory and time-frequency). In other words, the optimal results in each experiment are obtained with an ad-hoc set of features. This means the differentiation between the participants in different groups depends on different features. However, classifiers obtained very similar results for each experiment, showing a lesser influence.
The results obtained proved that combining phonatory-subsystem and time-frequency features improves the ability of the machine-learning models to detect bulbar involvement. In addition, detecting bulbar involvement also depends on the ad-hoc set of significant features found for such a case.
Table 4. Performance of male models. RF: Random Forest; LR: Logistic Regression; LDA: Linear Discriminant Analysis; NN: Neuronal Networks; SVM: Support Vector Machines.
Table 4. Performance of male models. RF: Random Forest; LR: Logistic Regression; LDA: Linear Discriminant Analysis; NN: Neuronal Networks; SVM: Support Vector Machines.
C vs. BC vs. NBB vs. NBC vs. A
RFAccuracy96.191.991.892.0
Sensitivity86.192.155.093.8
Specificity97.591.097.587.0
LRAccuracy91.989.288.591.3
Sensitivity95.090.375.090.7
Specificity92.086.989.594.0
LDAAccuracy95.091.181.392.0
Sensitivity85.088.690.090.7
Specificity97.598.080.596.0
NNAccuracy95.090.086.192.0
Sensitivity90.091.375.091.5
Specificity95.086.588.493.0
SVMAccuracy93.393.186.192.6
Sensitivity85.091.285.090.7
Specificity95.098.086.798.0

4.2. Comparison with Prior Work

This study is consistent with [6,7,8,36] which demonstrated that such phonatory-subsystem features as jitter, shimmer, pitch and HNR were sensitive indicators for describing pathological voices in ALS. It is also consistent with [6] where great uncertainty was found in the diagnosis of NBs participants.
Table 5. Performance of female models. RF: Random Forest; LR: Logistic Regression; LDA: Linear Discriminant Analysis; NN: Neuronal Networks; SVM: Support Vector Machines.
Table 5. Performance of female models. RF: Random Forest; LR: Logistic Regression; LDA: Linear Discriminant Analysis; NN: Neuronal Networks; SVM: Support Vector Machines.
C vs. BC vs. NBB vs. NBC vs. A
RFAccuracy98.194.184.895.8
Sensitivity96.692.592.395.8
Specificity10095.575.096.0
LRAccuracy91.493.074.791.3
Sensitivity91.390.075.093.4
Specificity91.595.575.087.0
LDAAccuracy93.190.472.190.7
Sensitivity87.682.570.087.3
Specificity86.697.575.098.0
NNAccuracy93.286.971.190.6
Sensitivity93.385.072.393.6
Specificity94.089.070.084.5
SVMAccuracy95.191.674.293.6
Sensitivity93.390.073.694.7
Specificity97.593.075.091.5
Besides the 15 phonatory-subsystem features obtained in [6], this study also provides 35 time-frequency features. The combination of phonatory-subsystem and time-frequency features, after performing MANOVA for feature selection, enhanced the outcomes of [6], which achieved the best results to date for detecting bulbar involvement in ALS using only acoustic features, ahead of [8,13,24].
Accuracies of up to 98.1% (RF) and 96.1% (RF) for females and males respectively were achieved when comparing the bulbar and control participants (case B vs. C). This accuracy exceeded the one obtained in [24] with SVM (79.0%) by 17.1% for males and 15.1% for females. The other studies found did not distinguish the classification problems by gender. In [6], SVM obtained an accuracy of 95.8%. In [13], NN based on Mel Frequency Cepstral Coefficients (coefficients for speech representation based on human auditory perception) obtained 90.7%. In [8], NN based on phonatory-subsystem features obtained 91.7% and adding motion sensors for both lip and tongue data increased the accuracy to 96.5% at the expense of including more invasive measurements. For females, our results outperformed those from the aforementioned studies by 2.3%, 7.4% and 6.4% respectively. For males, ours were 0.3% above those obtained in [6] and 5.4% and 4.4% above those obtained in [8,13].
When comparing ALS patients diagnosed with bulbar involvement with those patients in whom bulbar involvement has yet to be detected (B vs. NB), the outcomes outperformed the ones obtained in [6]. The respective accuracy for males and females increased by 16.3% and 9.3% with the same classifier (RF) (91.8% and 84.8% as against 75.5%). This is an important outcome which indicates that the use of time-frequency features increases the identification of bulbar involvement among patients with ALS.
The outcomes obtained in the C vs. NB and C vs. A cases were very similar to those in [6], reinforcing the idea that some NBs could have bulbar involvement.
The most important gains were obtained when comparing B and NB. The selection of the significant features for this comparison improved the outcomes. Thus, involvements (i.e., bulbar) could be detected through a separate, and more closely adjusted, set of features. Consequently, by increasing the identification of particular features, treatment could be better customized for each ALS patient.
In addition, only studies showing C vs. B have been presented in the literature (except in [6]). No attempts to distinguish other subjects have been made to date. We highlight this differentiating issue, and the importance of future research into it.

4.3. Limitations

The use of classification models with small datasets hinders the full assessment of the importance of the findings. The size of the dataset was, in part, determined by the low prevalence of ALS, which is considered a rare disease. The small number of samples in the B group was influenced by the heterogeneity of the ALS disease in which patients’ symptomatology is very diverse.
Furthermore, hand editing the segments of the voice recordings is inherently subjective and may introduce subtle and unintended selection biases. Although automatic instruments have been created, these methods are currently insufficiently accurate and require manual correction.

5. Conclusions and Future Work

This research directly addresses a recent statement released by the NEALS bulbar subcommittee regarding the need for methodologies based on objective measurements [37]. The outcomes achieved reinforce the idea that machine learning can be a suitable tool for helping with the diagnosis of ALS with bulbar involvement using common recording or mobile (i.e., smartphone) devices.
We demonstrate the usefulness of assessing bulbar involvement properly using phonatory-subsystem and time-frequency features from a study of the Spanish vowels that outperformed previous works, specifically [6,8,13,24]. It was also demonstrated that group identification depends on the significant features found for such an experiment.
The main contribution is the differentiation of diagnosis by gender. This outperformed all the results in the literature.
The next steps of this work will consist of improving the corpus for diagnosing bulbar dysfunction. It is planned to increase the sample size and enhance the annotation of the ALS patients without bulbar involvement. Novel methods based on the creation of vowel patterns and semi-supervised classification models will be developed to provide hints for distinguishing those ALS patients without bulbar involvement who may have been misdiagnosed.
Vowel patterns could be generated from the quasi-periodic components of a short stable segment of the five Spanish vowels. Principal and independent component analysis of these patterns is also envisioned.
Moreover, additional research is required to develop this concept properly. Longitudinal research studies are conceived in which patients’ diagnoses are obtained at multiple follow-ups. Several repetitions of the sustained phonations will be required to minimize sampling variability even for the control subjects.

Author Contributions

Conceptualization, A.T., F.C. and F.S.; methodology, A.T. and F.C.; software, A.T. and F.C.; validation, A.T., F.C., F.S. and M.P.; formal analysis, A.T. and F.C.; investigation, A.T., F.C. and F.S.; resources, A.T., F.C., F.S. and M.P.; data curation, A.T. and F.C.; writing—original draft preparation, A.T. and F.S.; writing—review and editing, A.T. and F.S.; visualization, A.T., F.C. and F.S.; supervision, F.S.; project administration, F.S.; funding acquisition, F.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Ministerio de Economía, Industria y Competitividad and the Ministerio de Ciencia e Innovación under contracts TIN2017-84553-C2-2-R and PID2020-113614RB-C22, respectively. AT is a member of CIMNE, a Severo Ochoa Centre of Excellence (2019-2023) under grant CEX2018-000797-S, funded by MCIN/AEI/10.13039/501100011033.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Research Ethics Committee for Biomedical Research Projects (CEIm) at the Bellvitge University Hospital in Barcelona. (ref.: 065/19). Date: 29 May 2019.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. Written informed consent has been obtained from the patient(s) to publish this paper.

Acknowledgments

Neurology Department of the Bellvitge University Hospital in Barcelona allowed the recording of the voices of the participants in its facilities. Carlos Augusto Salazar Talavera reported the clinical records of the participants. Marta Fulla and Maria Carmen Majos Bellmunt gave advise about the process of eliciting the sounds.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Source and Synthetic Data Repository. All the models implemented in R are available in a GitHub repository in ELA source files (https://github.com/atenad/ALS), accessed on 27 January 2022. A synthetic dataset of the features obtained is also provided. It was built in R from the original values using the synthpop package [38].

Abbreviations

The following abbreviations are used in this manuscript:
TFRTime-frequency representation
Jitter (absolute)inter-cycle variation of the fundamental period
Jitter(relative)relative period
Jitter(rap)relative perturbation
Jitter(ppq5)five-point period perturbation quotient
Shimmer(dB)Variability of the peak-to-peak amplitude
Shimmer(relative)relative amplitudes of consecutive periods
Shimmer(apqP)three, five and eleven-point amplitude perturbation
HNRharmonics-to-noise ratio
HNR(mean)mean HNR
HNR(SD)standard deviation of HNR
WDWigner distribution
CWDChoi-Williams exponential function
pDjoint probability density distribution
E(t)average instantaneous spectral energy
f m i ( t ) average instantaneous frequency
f _ C r e s ( t ) instantaneous frequency peak
H_tinstantaneous information entropy
H_fspectral information entropy
H_tfjoint Shannon entropy
IE(f)spectral information
Kkurtosis
t n f m joint time-frequency moment
SVMSupport Vector Machine
NNNeural Networks
LDALinear Discriminant Analysis
LRLinear Logistic Regression
RFRandom Forest
Ccontrol group
Bgroup of ALS participants diagnosed with bulbar involvement
NBgroup of ALS participants not diagnosed with bulbar involvement
Agroup of ALS participants with or without bulbar involvement

References

  1. Kiernan, M.C.; Vucic, S.; Cheah, B.C.; Turner, M.R.; Eisen, A.; Hardiman, O.; Burrell, J.R.; Zoing, M.C. Amyotrophic lateral sclerosis. Lancet 2011, 377, 942–955. [Google Scholar] [CrossRef] [Green Version]
  2. Tomik, B.; Guiloff, R. Dysarthria in amyotrophic lateral sclerosis: A review. Amyotroph. Lateral Scler. Off. Publ. World Fed. Neurol. Res. Group Mot. Neuron Dis. 2010, 11, 4–15. [Google Scholar] [CrossRef] [PubMed]
  3. Shellikeri, S.; Green, J.R.; Kulkarni, M.; Rong, P.; Martino, R.; Zinman, L.; Yunusova, Y. Speech Movement Measures as Markers of Bulbar Disease in Amyotrophic Lateral Sclerosis. J. Speech Lang. Heart Res. JSLHR 2016, 59, 887–899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Lee, J.; Dickey, E.; Simmons, Z. Vowel-Specific Intelligibility and Acoustic Patterns in Individuals With Dysarthria Secondary to Amyotrophic Lateral Sclerosis. J. Speech Lang. Heart Res. 2019, 62, 1–26. [Google Scholar] [CrossRef] [PubMed]
  5. Carpenter, R.; McDonald, T.; Howard, F. The Otolaryngologic Presentation of Amyotrophic Lateral Sclerosis. Otolaryngology 1978, 86, ORL-479–ORL-484. [Google Scholar] [CrossRef]
  6. Tena, A.; Clarià, F.; Solsona, F.; Meister, E.; Povedano, M. Detection of Bulbar Involvement in Patients With Amyotrophic Lateral Sclerosis by Machine Learning Voice Analysis: Diagnostic Decision Support Development Study. JMIR Med. Inform. 2021, 9, e21331. [Google Scholar] [CrossRef]
  7. Silbergleit, A.K.; Johnson, A.F.; Jacobson, B.H. Acoustic analysis of voice in individuals with amyotrophic lateral sclerosis and perceptually normal vocal quality. J. Voice 1997, 11, 222–231. [Google Scholar] [CrossRef]
  8. Wang, J.; Kothalkar, P.V.; Heitzman, D. Towards Automatic Detection of Amyotrophic Lateral Sclerosis from Speech Acoustic and Articulatory Samples. In Proceedings of the InterSpeech, San Francisco, CA, USA, 8–12 September 2016. [Google Scholar]
  9. Chiaramonte, R.; Luciano, C.; Chiaramonte, I.; Serra, A.; Bonfiglio, M. Multi-disciplinary clinical protocol for the diagnosis of bulbar amyotrophic lateral sclerosis. Acta Otorrinolaringol. 2019, 70, 25–31. [Google Scholar] [CrossRef]
  10. Tomik, J.; Tomik, B.; Wiatr, M.; Skladzien, J.; Strek, P.; Szczudlik, A. The Evaluation of Abnormal Voice Qualities in Patients with Amyotrophic Lateral Sclerosis. Neuro-Degener. Dis. 2015, 15, 225–232. [Google Scholar] [CrossRef]
  11. Horwitz-Martin, R.L.; Horwitz-Martin, R.L.; Quatieri, T.F.; Lammert, A.C.; Williamson, J.R.; Yunusova, Y.; Godoy, E.; Mehta, D.D.; Green, J.R. Relation of automatically extracted formant trajectories with intelligibility loss and speaking rate decline in amyotrophic lateral sclerosis. In Proceedings of the InterSpeech, San Francisco, CA, USA, 8–12 September 2016; pp. 1205–1209. [Google Scholar]
  12. Spangler, T.; Vinodchandran, N.V.; Samal, A.; Green, J.R. Fractal features for automatic detection of dysarthria. In Proceedings of the 2017 IEEE EMBS International Conference on Biomedical Health Informatics (BHI), Orland, FL, USA, 16–19 February 2017; pp. 437–440. [Google Scholar] [CrossRef]
  13. Suhas, B.; Patel, D.; Rao, N.; Belur, Y.; Reddy, P.; Atchayaram, N.; Yadav, R.; Gope, D.; Ghosh, P.K. Comparison of Speech Tasks and Recording Devices for Voice Based Automatic Classification of Healthy Subjects and Patients with Amyotrophic Lateral Sclerosis. In Proceedings of the Interspeech, Graz, Austria, 15–19 September 2019; pp. 4564–4568. [Google Scholar] [CrossRef] [Green Version]
  14. Melia, U.; Vallverdú, M.; Jospin, M.; Jensen, E.W.; Valencia, J.F.; Clariá, F.; Gambus, P.L.; Caminal, P. Prediction of nociceptive responses during sedation by time-frequency representation. In Proceedings of the 2013 35th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), Osaka, Japan, 3–7 July 2013; pp. 2547–2550. [Google Scholar] [CrossRef]
  15. Melia, U.; Claria, F.; Vallverdu, M.; Caminal, P. Measuring Instantaneous and Spectral Information Entropies by Shannon Entropy of Choi-Williams Distribution in the Context of Electroencephalography. Entropy 2014, 16, 2530–2548. [Google Scholar] [CrossRef] [Green Version]
  16. Claria, F.; Vallverdú, M.; Baranowski, R.; Chojnowska, L.; Caminal, P. Heart rate variability analysis based on time-frequency representation and entropies in hypertrophic cardiomyopathy patients. Physiol. Meas. 2008, 29, 401–416. [Google Scholar] [CrossRef]
  17. Clariá, F.; Vallverdú, M.; Riba, J.; Romero, S.; Barbanoj, M.J.; Caminal, P. Characterization of the cerebral activity by time-frequency representation of evoked EEG potentials. Physiol. Meas. 2011, 32, 1327–1346. [Google Scholar] [CrossRef] [PubMed]
  18. Tena, A.; Clarià, F.; Solsona, F. Automated detection of COVID-19 cough. Biomed. Signal Process. Control 2022, 71, 103175. [Google Scholar] [CrossRef] [PubMed]
  19. Mika, D.; Józwik, J. Advanced Time-Frequency Representation in Voice Signal Analysis. Adv. Sci. Technol. Res. J. 2018, 12, 251–259. [Google Scholar] [CrossRef]
  20. Yang, Y.; Peng, Z.; Zhang, W.; Meng, G. Parameterised time-frequency analysis methods and their engineering applications: A review of recent advances. Mech. Syst. Signal Process. 2019, 119, 182–221. [Google Scholar] [CrossRef]
  21. Davies, M.; Daudet, L. Sparse audio representations using the MCLT. Signal Process. 2006, 86, 457–470. [Google Scholar] [CrossRef]
  22. Makkonen, T.; Ruottinen, H.; Puhto, R.; Helminen, M.; Palmiol, J. Speech deterioration in amyotrophic lateral sclerosis (ALS) after manifestation of bulbar symptoms: Speech deterioration in ALS. Int. J. Lang. Commun. Disord. 2017, 53. [Google Scholar] [CrossRef] [Green Version]
  23. Tomik, B.; Krupinski, J.; Glodzik-Sobanska, L.; Bala-Slodowska, M.; Wszolek, W.; Kusiak, M.; Lechwacka, A. Acoustic analysis of dysarthria profile in ALS patients. J. Neurol. Sci. 1999, 169, 35–42. [Google Scholar] [CrossRef]
  24. Norel, R.; Pietrowicz, M.; Agurto, C.; Rishoni, S.; Cecchi, G. Detection of Amyotrophic Lateral Sclerosis (ALS) via Acoustic Analysis. bioRxiv 2018. [Google Scholar] [CrossRef]
  25. An, K.; Kim, M.; Teplansky, K.; Green, J.; Campbell, T.; Yunusova, Y.; Heitzman, D.; Wang, J. Automatic Early Detection of Amyotrophic Lateral Sclerosis from Intelligible Speech Using Convolutional Neural Networks. In Proceedings of the Interspeech, Hyderabad, India, 2–6 September 2018; pp. 1913–1917. [Google Scholar] [CrossRef] [Green Version]
  26. Gutz, S.E.; Wang, J.; Yunusova, Y.; Green, J.R. Early Identification of Speech Changes Due to Amyotrophic Lateral Sclerosis Using Machine Classification. In Proceedings of the Interspeech, Graz, Austria, 15–19 September 2019; pp. 604–608. [Google Scholar] [CrossRef] [Green Version]
  27. Audacity Manual. Available online: https://manual.audacityteam.org (accessed on 27 January 2022).
  28. Sprecher, A.; Olszewski, A.; Jiang, J.J.; Zhang, Y. Updating signal typing in voice: Addition of type 4 signals. J. Acoust. Soc. Am. 2010, 127, 3710–3716. [Google Scholar] [CrossRef] [Green Version]
  29. Boersma, P.; Weenink, D. Praat: Doing Phonetics by Computer, V6.1.01; Technical Report; University of Amsterdam: Amsterdam, The Netherlands, 2019. [Google Scholar]
  30. MATLAB. Version 9.9.0.1495850 (R2020b); The MathWorks Inc.: Natick, MA, USA, 2020. [Google Scholar]
  31. Loughlin, P. What are the time-frequency moments of a signal? Proc. SPIE Int. Soc. Opt. Eng. 2001, 4474, 35–44. [Google Scholar] [CrossRef]
  32. IBM Corp. IBM SPSS Statistics for Windows; IBM Corp.: Armonk, NY, USA, 2016. [Google Scholar]
  33. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  34. Kuhn, M.; Johnson, K. Applied Predictive Modeling; Springer: Berlin/Heidelberg, Germany, 2013. [Google Scholar] [CrossRef]
  35. Tharwat, A. Classification assessment methods. Appl. Comput. Inform. 2021, 17, 168–192. [Google Scholar] [CrossRef]
  36. Vashkevich, M.; Petrovsky, A.; Rushkevich, Y. Bulbar ALS Detection Based on Analysis of Voice Perturbation and Vibrato. In Proceedings of the 2019 Signal Processing: Algorithms, Architectures, Arrangements and Applications (SPA), Poznan, Poland, 18–20 September 2019; pp. 267–272. [Google Scholar]
  37. Plowman, E.K.; Tabor, L.C.; Wymer, J.; Pattee, G. The evaluation of bulbar dysfunction in amyotrophic lateral sclerosis: Survey of clinical practice patterns in the United States. Amyotroph. Lateral Scler. Front. Degener. 2017, 18, 351–357. [Google Scholar] [CrossRef] [PubMed]
  38. Nowok, B.; Raab, G.M.; Dibben, C. synthpop: Bespoke Creation of Synthetic Data in R. J. Stat. Softw. 2016, 74, 1–26. [Google Scholar] [CrossRef] [Green Version]
Figure 1. p D ( f , t ) of vowel “a” for 3 different patients with bulbar involvement. The marked difference in the graphic representation of the time-frequency between the subjects can be clearly appreciated. (a) Patient pD without bulbar involvement. (b) Patient pD with slight bulbar involvement. (c) Patient pD with severe bulbar involvement.
Figure 1. p D ( f , t ) of vowel “a” for 3 different patients with bulbar involvement. The marked difference in the graphic representation of the time-frequency between the subjects can be clearly appreciated. (a) Patient pD without bulbar involvement. (b) Patient pD with slight bulbar involvement. (c) Patient pD with severe bulbar involvement.
Sensors 22 01137 g001
Table 1. ALS participants clinical records. Notation: Age (in years). ALSFR-R (Rating Scale-Revised): scores (0–48) the severity of ALS; Bulbar: Bulbar involvement; NA: Data not available.
Table 1. ALS participants clinical records. Notation: Age (in years). ALSFR-R (Rating Scale-Revised): scores (0–48) the severity of ALS; Bulbar: Bulbar involvement; NA: Data not available.
AgeSexALSFR-RBulbarBulbar Onset Symptoms
37F37NONo Symptoms
38M6YESNA
39M43NONo Symptoms
41M34NONo Symptoms
41M34NONo Symptoms
43F21YESDysphagia
44F19NONo Symptoms
48F36NONo Symptoms
48F29YESDysphagia
48M31NONo Symptoms
48M45NONo Symptoms
49MNANONo Symptoms
50M39NONo Symptoms
52M43NONo Symptoms
52F27YESDysphagia
52M33NONo Symptoms
53F29YESDysphagia/Dysarthria
55M26NONo Symptoms
55M24NONo Symptoms
56M35NONo Symptoms
56M27NONo Symptoms
58F46YESDysarthria
58M28YESNA
59F33YESNA
60M46YESNA
63M22NONo Symptoms
63M42NONo Symptoms
63MNANONo Symptoms
65M24NONo Symptoms
66F41NONo Symptoms
67MNANONo Symptoms
67F33YESDyspnoea
68MNANONo Symptoms
68F21NONo Symptoms
69M37NONo Symptoms
70F28YESDysphagia
70F17NONo Symptoms
70M46NONo Symptoms
70M27NONo Symptoms
70F23YESDysphagia/Dysarthria
71M39NONo Symptoms
71F32YESDysphagia
72M30NONo Symptoms
72F38NONo Symptoms
76F30NONo Symptoms
81M36NONo Symptoms
81M28NONo Symptoms
84F30YESNA
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tena, A.; Clarià, F.; Solsona, F.; Povedano, M. Detecting Bulbar Involvement in Patients with Amyotrophic Lateral Sclerosis Based on Phonatory and Time-Frequency Features. Sensors 2022, 22, 1137. https://doi.org/10.3390/s22031137

AMA Style

Tena A, Clarià F, Solsona F, Povedano M. Detecting Bulbar Involvement in Patients with Amyotrophic Lateral Sclerosis Based on Phonatory and Time-Frequency Features. Sensors. 2022; 22(3):1137. https://doi.org/10.3390/s22031137

Chicago/Turabian Style

Tena, Alberto, Francesc Clarià, Francesc Solsona, and Mònica Povedano. 2022. "Detecting Bulbar Involvement in Patients with Amyotrophic Lateral Sclerosis Based on Phonatory and Time-Frequency Features" Sensors 22, no. 3: 1137. https://doi.org/10.3390/s22031137

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop