BioMed Research International
Volume 2015 (2015), Article ID 704121, 9 pages
http://dx.doi.org/10.1155/2015/704121

Clinical Study

Acoustic Correlates of Compensatory Adjustments to the Glottic and Supraglottic Structures in Patients with Unilateral Vocal Fold Paralysis

Luis M. T. Jesus,1,2 Joana Martinez,¹ Andreia Hall,1,3 and Aníbal Ferreira⁴

¹Institute of Electronics and Informatics Engineering of Aveiro (IEETA), University of Aveiro, 3810-193 Aveiro, Portugal
²School of Health Sciences (ESSUA), University of Aveiro, 3810-193 Aveiro, Portugal
³Department of Mathematics (DMat), University of Aveiro, 3810-193 Aveiro, Portugal
⁴Department of Electrical and Computer Engineering, University of Porto, 4200-465 Porto, Portugal

Received 22 January 2015; Accepted 24 April 2015

Academic Editor: Haldun Oguz

Copyright © 2015 Luis M. T. Jesus et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The goal of this study was to analyse perceptually and acoustically the voices of patients with Unilateral Vocal Fold Paralysis (UVFP) and compare them to the voices of normal subjects. These voices were analysed perceptually with the GRBAS scale and acoustically using the following parameters: mean fundamental frequency (F0), standard-deviation of F0, jitter (ppq5), shimmer (apq11), mean harmonics-to-noise ratio (HNR), mean first (F1) and second (F2) formants frequency, and standard-deviation of F1 and F2 frequencies. Statistically significant differences were found in all of the perceptual parameters. Also the jitter, shimmer, HNR, standard-deviation of F0, and standard-deviation of the frequency of F2 were statistically different between groups, for both genders. In the male data differences were also found in F1 and F2 frequencies values and in the standard-deviation of the frequency of F1. This study allowed the documentation of the alterations resulting from UVFP and addressed the exploration of parameters with limited information for this pathology.

1. Introduction

A neural dysfunction of the larynx leads to alterations in voice, respiration, and airway protection. Usually, Unilateral Vocal Fold Paralysis (UVFP) is related to a set of well-documented perceptive alterations such as weak voice, breathiness, roughness, diminished voice intensity, vocal effort, low voice efficiency, voice breaks, diplophonia, and air loss [1–5]. Furthermore, vocal strain is a critical component in various vocal pathologies including UVFP. A neuronal dysphonia, such as UVFP, can alter the vibrational patterns of the Vocal Folds (VF) which leads to compensatory adjustments to the glottic and supraglottic structures that increase the vocal effort and vocal strain perception [6, 7]. In addition to the perceptive alterations, UVFP also results in higher values of jitter and shimmer and lower values of the harmonics-to-noise ratio (HNR) [1–4, 8]. Furthermore, values of standard-deviation of fundamental frequency (F0) are reported as higher than normal because of the diminished control of the vibrational pattern of the VF, causing greater variability [9–11]. According to Schwarz et al. [6], there is a need to describe and understand the UVFP patient’s larynx configuration for a better and more individualised vocal intervention, preventing compensatory adjustments. Formant frequencies provide acoustic cues about the vocal tract configuration [12–14]. According to Lee et al. [15] the formant’s values are relevant for discriminating normal from pathologic voices and the configuration of the vocal tract is different during phonation in people with vocal pathologies. The same authors [15] found slightly lower values of the first formant (F1) frequency and higher values of the second formant (F2) frequency in cases of UVFP. This indicates that UVFP subjects tend to have a more elevated and advanced tongue position during phonation [13, 14]. A breathy voice (common in UVFP) is reported to be associated with the same configuration referred to previously [16]. However, Titze [13] reports an approximation of the values of the frequency of F1 and F2 in cases of narrower vocal tract. These vocal tract modifications may result from the attempt to compensate the vocal alteration by patients exhibiting UVFP [2]. According to Lee et al. [15] the standard-deviations of the frequency of F1 and F2 have higher values in cases of UVFP indicating a higher instability of the vocal tract configuration during phonation.

The aim of this study was to compare perceptually and acoustically the voices of subjects with UVFP and the voices of subjects representing normal quality. Measures related to the vocal tract configuration, namely, formant frequencies, were also analysed and correlated with alterations caused by vocal pathology.

2. Materials and Methods

This is a quantitative, descriptive, and cross-sectional study [17–19]. The recordings were made in Hospital de Santo António and Hospital de São João, both in Porto, Portugal, and at the Speech, Language, and Hearing Laboratory (SLHlab) at the University of Aveiro, Portugal. This took place as part of the data collection process of the first representative European Portuguese pathological voice database [20]. Part of this data was divided into two groups: a group having vocal pathology (UVFP) and a group without vocal pathology. A group of 17 patients, evaluated with videolaryngoscopy and diagnosed with UVFP, formed the pathologic group. The inclusion criteria for this group were having diagnosis of UVFP, not having had speech and language therapy intervention, and being over 18 years old. The exclusion criteria were having other concomitant pathologies to UVFP and/or having been submitted to a surgical intervention to correct the vocal pathology. A group of 85 normal voice volunteers were included in the control group based on two distinct procedures: 43 subjects were evaluated with videolaryngoscopy and diagnosed as normal; 42 subjects were evaluated using a vocal anamnesis and summative evaluation (a similar procedure was used by Roark et al. [21]). The inclusion criteria for the control group were having normal voice quality and being over 18 years old. The exclusion criterion was having vocal or other pathologies that may interfere with normal voice production.

Each pathologic case was individually matched to five subjects of the control group in order to increase the power of statistical tests [17, 22]. The cases were matched according to gender and age. The first variable was gender because after puberty there is a set of different characteristics that differentiate male and female voices [23]. The second variable was age because with aging some functional and structural modifications occur at phonatory level [23, 24]. Taking into account the fact that there are notable voice changes if the subjects’ age difference is more than 10 years [25–29] the maximum allowed difference of age between the matched subjects was 5 years, in an attempt to reduce variability.

Four (4) subjects with UVFP were male (23.5%) and 13 subjects were female (76.5%). The youngest patient was 30 years old and the oldest 72. The mean age for the pathologic group was 56.7 years with a standard-deviation of 12.7 years. Nine (9) patients had left UVFP (52.9%) and 8 right UVFP (47.1%). In the control group 20 subjects were male (23.5%) and 65 were female (76.5%). The mean age of the control group was 56.1 years and the standard-deviation was 12.7 years.

The voice recordings were made in a clinical setting using Praat 5.3.56 (32-bit edition) [30]. A Behringer ECM8000 microphone and a Presonus AudioBox USB (16 bits and 48000 Hz) were used for all of the recordings. The subjects were seated and the microphone was aligned to the mouth at a distance of 30 cm [31, 32]. An informed consent was signed and the vowel [a] was recorded. A parcel of the vowel was then annotated according to criteria defined by Pinho et al. [3]: 200 ms after the onset of phonation and with approximately 100 cycles. This parcel was then manually analysed with Praat 5.3.56 (64-bit edition) with an autocorrelation method (used by default by the software) to estimate F0. There were some errors in the identification of the period, so a modification of the “octave cost” to a higher value was made (as suggested in Praat’s manual). The values of the parameters used to run the autocorrelation method are presented in Table 1.

Table 1: Values of the autocorrelation method used in Praat for the voice analysis.

Parameter
Value

Maximum number of candidates
15

Silence threshold
0.03

Voicing threshold
0.45

Octave cost
0.15

Octave-jump cost
0.35

Voiced/unvoiced cost
0.14

From the “voice report” Praat window the following values were extracted: mean F0; standard-deviation of F0; jitter (ppq5); shimmer (apq11); mean harmonics-to-noise ratio (HNR). The Burg [33] method (used by default by Praat) was used to track the formants. The “formant listing” for the same 100 cycles was obtained and the mean value and standard-deviation were calculated for the frequency of F1 and F2. The values were double-checked through the spectrogram of each segment.

Each voice was also perceptually assessed using the GRBAS scale [34]. For the pathologic voices a group of five speech and language therapists with expertise in voice assessment made the perceptive evaluation. For the normal voices one speech and language therapist made the perceptive assessment. For these procedures the experts used the following headphones connected to the internal soundcard of a laptop computer: Sennheiser HD 380 Pro; Sennheiser HD201; Sony MDR-CD270; Sony MDRZX100B; Sony MDR-ZX110NA. All of the assessments were made blindly regarding the group (patients or normal subjects).

For the statistical analysis IBM SPSS Statistics version 20 was used. The interrater consistency was analysed using the Kendall Coefficient. The Mann-Whitney test was used to analyse the GRBAS scale parameters. The acoustic parameters that had normal distribution (HNR, F2♀, standard-deviation of F0♂, F1♂) were statistically analysed using the -test and parameters that did not have normal distribution (Jitter (ppq5), Shimmer (apq11), F0♀, standard-deviation of F0♀, F1♀, standard-deviation of F1♀, standard-deviation of F2♀, F0♂, standard-deviation of F1♂, F2♂, and standard-deviation of F2♂) were analysed with the Mann-Whitney test. The normality was tested with the Shapiro-Wilk test. A level of significance of 0.05 was used for all statistical analyses.

All of the procedures had the acceptance of the Ethical Commission of the Hospital de Santo António and Hospital de São João. An authorisation from the National Commission for Data Protection was also obtained.

3. Results and Discussion

3.1. Interrater Consistency

The consistency between the five judges that assessed the pathologic voices was analysed using Kendall’s test. Table 2 shows that there is consistency in all of the parameters of the GRBAS scale between judges. The fact that the judges presented consistency between them indicates that they have a similar internal understanding of the used instrument [35]. This consistency is likely to be related to the fact that the GRBAS scale is widely used, understood, and recommended worldwide by clinicians [36]. The ’s value, shown in Table 2, can vary between 0 (no general tendency of consistency between judges) and 1 (all judges responded equally) [37]. In Table 2 we can also see that the lowest value of was found for the R (Rough) parameter. This may be due to the fact that this parameter is a supraclass of perceptive parameters that can lead to various interpretations between different judges [38]. The fact that none of the parameters had a very good consistency was expected because the perceptive assessment is a very complex procedure that includes various subjective elements that are not totally understood [36, 39]. Despite the results varying from reasonable to good, perceptive evaluation is still a central procedure in the vocal assessment [40].

Table 2: Interrater consistency—Kendall’s test.

Scale parameter

value

G
0.263
0.001

R
0.160
0.033

B
0.381
<0.001

A
0.344
<0.001

S
0.438
<0.001

G: Grade; R: Rough; B: Breathy; A: Asthenic; S: Strained; : Kendall’s .

3.2. Comparison of GRBAS Scale Parameters between Normal and UVFP Voices

The results of the perceptive assessment of the voices of the normal and UVFP subjects were analysed using the Mann-Whitney test. Table 3 shows that all of the GRBAS parameters were statistically different between groups, being higher in the pathologic group as expected (see Figure 1). The control group had a mean score of zero, which was expected because the control group was intended to have a normal/nonaltered voice quality that would be associated to a 0 value (normal) of all parameters assessed in GRBAS. In the pathologic group we can see that the parameter with the highest values was G (Grade), which has been observed before by other authors [41, 42]. In this group of UVFP there were alterations in all of the GRBAS parameters, varying between a mild and moderate grade of perturbation. Grade (G), Rough (R), and Breathy (B) presented the highest mean scores, as previously observed by various authors [1–4, 43]. Another disturbance that is commonly found in subjects with UVFP is a weak voice [1, 4, 44] which is reflected in parameter A (Asthenic), also found in this sample. In addition to the previous parameters, according to Rosenthal et al. [7], it is usual to find vocal strain (parameter S) in these cases, which could also be observed in this study.

Table 3: Comparison of the results of GRBAS scale between UVFP and normal voice subjects.

UVFP
Normal

value

Mean ± SD

Mean ± SD

G
17
2.06 ± 0.827
85
0
0
<0.001

R
17
1.94 ± 0.899
85
0
0
<0.001

B
17
1.71 ± 0.772
85
0
0
<0.001

A
17
1.24 ± 0.437
85
0
0
<0.001

S
17
0.94 ± 0.556
85
0
25.5
<0.001

UVFP: Unilateral Vocal Fold Paralysis; G: Grade; R: Rough; B: Breathy; A: Asthenic; S: Strained; : number of cases; SD: standard deviation; : Mann-Whitney test.

Figure 1: Comparison between mean scores of the GRBAS scale for normal and UVFP subjects.

One of the major alterations caused by UVFP is the incomplete glottal closure that originates excess air during phonation that creates a breathy voice (parameter B is altered) [2, 4, 45]. This air leakage leads to a lower voice energy originating a weak voice (parameter A is altered) [2, 4, 45]. The irregularity of the VF cycles (parameter R reflects this) is due to the reduced mobility/immobility of the paralysed VF or to the fact that the unhealthy VF may present a passive vibration [4, 46]. In some cases, in an attempt to overcome the alterations caused by the UVFP, patients create compensations that can lead to strain in the supraglottic region, increasing the vocal effort and giving the voice a strained characteristic (parameter S) [4, 7]. Grade (G) is related with the other parameters and varies according to the severity of the overall voice perturbation [47].

3.3. Comparison of Acoustic Parameters between Normal and UVFP Voices

Although perceptive assessment is the most used technique for vocal assessment, it is a subjective process that leads to some variability issues [8]. Contrary to this, acoustic data allows objective and noninvasive measures about the behaviour of the VF [8, 15, 48–50]. Table 4 shows statistically different values of jitter (ppq5), shimmer (apq11), and HNR between the normal and pathologic voices. Jitter, which is related to the absolute difference between the durations of consecutive cycles [43], is higher in UVFP subjects (see Figure 2). These results were also obtained by other authors [2, 3, 8]. These higher values may be due to the asymmetry at the VF level, caused by the UVFP that leads to vibration irregularities in frequency altering the jitter values [2]. Similarly shimmer, which is related to the absolute difference between the amplitudes of consecutive cycles [43], is also higher in UVFP cases (see Figure 3). These results were also obtained by other authors [2, 3, 8]. The asymmetry caused by UVFP leads to vibration irregularities in amplitude altering shimmer values [2]. This parameter is also increased by a poor and inconsistent contact between VF, which is very common in UVFP [51]. Thus, these UVFP subjects present more cyclic irregularity at frequency and amplitude level compared to the normal voice subjects. It should be noted that we also had higher than normal values of shimmer in the normal sample. This may be due to the fact that the recordings were made in a clinical setting that is not entirely noise-free and this could have interfered with the data calculation of this parameter. Regarding the HNR, which is obtained from the ratio between the harmonic and noise components of the signal [43], the results indicate a lower value in the pathologic group (see Figure 4). These results were consistent with the literature [2, 8]. The alterations in periodicity caused by the UVFP originate a lower ratio between the two components, diminishing the HNR values in the pathologic cases [2]. These results indicate that patients with UFVP have higher relative noise amplitude during phonation (than the normal subjects) lowering the HNR value.

Table 4: Comparison of jitter, shimmer, and HNR between normal and UVFP subjects.

UVFP
Normal
or
value

Mean ± SD

Mean ± SD

Jitter ppq5 (%)
17
1.06 ± 1.02
85
0.26 ± 0.18

<0.001

Shimmer apq11 (%)
17
10.16 ± 3.34
85
7.14 ± 3.22

0.001

HNR (dB)
17
10.11 ± 4.94
85
14.75 ± 4.46

<0.001

UVFP: Unilateral Vocal Fold Paralysis; SD: standard deviation; : -test;  U: Mann-Whitney test.

Figure 2: Jitter (%) values for UVFP and normal subjects.

Figure 3: Shimmer (%) values for UVFP and normal subjects.

Figure 4: HNR (dB) values for UVFP and normal subjects.

The parameters presented in Tables 5 and 6 were divided by gender because females and males have different inherent vocal tract and VF characteristics, especially in terms of size and mass [52]. For F0 (see Figures 5 and 6) we can see that there are no significant statistical differences between pathologic and normal voices in both genders. This fact was also previously described by Oguz et al. [8]. Fundamental frequency is directly related to and dependent of length, tension, mass, rigidity, and the interaction with the subglottic pressure [53]. The fact that there are no differences between the two groups indicates that, in this sample, the modifications at VF level caused by UVFP are not sufficient to create real alterations in F0. Also according to Woo et al. [54] the majority of UVFP subjects present F0 values close to normal.

Table 5: Fundamental frequency () and first and second formant frequencies ( and ) and their standard-deviations, for normal and UFVP female participants.

♀
UVFP
Normal
or
value

Mean ± SD

Mean ± SD

(Hz)
13
218.38 ± 72.36
65
195.36 ± 33.02

0.335

SD (Hz)
13
6.65 ± 12.28
65
2.62 ± 2.15

0.018

(Hz)
13
826.16 ± 171.73
65
819.03 ± 164.75

0.495

SD (Hz)
13
117.18 ± 99.91
65
71.28 ± 53.54

0.116

(Hz)
13
1522.69 ± 96.00
65
1453.51 ± 139.20
= 1.58
0.059

SD (Hz)
13
156.31 ± 146.20
65
62.26 ± 48.67

0.002

UVFP: Unilateral Vocal Fold Paralysis; : number of cases; SD: standard deviation; : -test;  : Mann-Whitney test; : fundamental frequency; SD : standard-deviation of the fundamental frequency; : first formant frequency; SD : standard-deviation of first formant frequency; : second formant frequency; SD : standard-deviation of second formant frequency.

Table 6: Fundamental frequency () and first and second formant frequencies ( and ) and their standard-deviations, for normal and UFVP male participants.

♂
UVFP
Normal
or
value

Mean ± SD

Mean ± SD

(Hz)
4
121.43 ± 12.70
20
128.27 ± 23.85

0.485

SD (Hz)
4
3.41 ± 1.25
20
1.36 ± 0.58

<0.001

(Hz)
4
821.48 ± 331.80
20
677.33 ± 84.95

0.043

SD (Hz)
4
191.91 ± 105.62
20
30.58 ± 18.84

<0.001

(Hz)
4
1629.09 ± 474.06
20
1282.95 ± 104.35

0.011

SD (Hz)
4
263.68 ± 144.45
20
36.89 ± 27.53

0.001

UVFP: Unilateral Vocal Fold Paralysis; : number of cases; SD: standard deviation; : -test; : Mann-Whitney test; : fundamental frequency; SD : standard-deviation of the fundamental frequency; : first formant frequency; SD : standard-deviation of first formant frequency; : second formant frequency; SD : standard-deviation of second formant frequency.

Figure 5: F0 (Hz) values for female UVFP and normal subjects.

Figure 6: F0 (Hz) values for male UVFP and normal subjects.

The standard-deviation of F0 (see Figures 7 and 8), which is related to the variations in vibration and muscular control of the VF, is higher in the pathologic group indicating important alterations in the described aspects [53]. Thus, subjects with UVFP present more F0 variability indicating a poorer muscular control and lower vibrational stability of VF. These results are supported by other authors [10, 11, 46, 53].

Figure 7: SD of F0 (Hz) values for female UVFP and normal subjects.

Figure 8: SD of F0 (Hz) values for male UVFP and normal subjects.

The vocal tract configuration interacts with VF oscillation; that is, vocal tract configuration constrains VF functioning during phonation [15, 55]. After the onset of UVFP patients usually develop some compensatory adjustments at glottic and supraglottic level altering voice and vocal tract configuration [6]. The description of vocal tract configurations in subjects with UVFP could guide treatments and help prevent negative compensations [6, 7].

Regarding F1 frequency, Table 5 shows that for females there are no statistically significant differences between pathologic and normal subjects (see Figures 9 and 10). A similar result was obtained by Lee et al. [15]. Formant frequency values shown in Table 6 reveal that, for males, differences between groups are statistically significant. Lower values of F1 frequencies in UVFP cases were expected (based on data reported previously [15]); however, Table 6 clearly shows that the F1 frequency values were higher in the pathologic group. However, authors such as Hartl et al. [2] and D. H. Klatt and L. C. Klatt [56] have also reported higher F1 frequency values for voices with similar characteristics to UVFP patients. Since the frequency of F1 is inversely related to the vertical movement of the tongue, higher values of this formant (in UVFP subjects) indicate a lower tongue position during phonation for the pathologic subjects. This result is in line with what was found by Higashikawa et al. [57] for whispered voices.

Figure 9: F1 and F2 frequency (Hz) values for female UVFP and normal subjects.

Figure 10: F1 and F2 frequency (Hz) values for male UVFP and normal subjects.

The second formant (F2) frequency, which is related to the horizontal tongue movement, is higher in UVFP male subjects (see Figures 9 and 10). This result was also obtained in other studies [2, 15]. For females, although the value is very close to the significance level, there are no statistically significant differences between normal and UVFP subjects. However, we can see a slightly higher value of F2 frequency in the female pathologic group compared to normal females. Therefore results indicate that there could be a tendency to a more advanced tongue position during phonation in cases of UVFP. This is consistent with the results presented by Lotto et al. [16] who studied breathy voices (typical of UVFP).

As for the SD of the frequency of F1 shown in Table 5, there were significant differences between the two groups being SD of F1 frequency higher in the patients, for male participants. There were no significant differences between groups for females (see Table 6). As for the SD of the frequency of F2 there were statistically significant higher values in the UVFP group for both genders (see Figures 11 and 12). Therefore these parameters, especially the SD of the frequency of F2, may have an important role in discriminating normal and UVFP voices. Pathologic voices showed higher values of formant frequency SD. These results were also obtained by Lee et al. [15]. This indicates a greater instability of the vocal tract configuration in UVFP during phonation.

Figure 11: SD of F1 and F2 frequency (Hz) for female UVFP and normal subjects.

Figure 12: SD of F1 and F2 frequency (Hz) for male UVFP and normal subjects.

Overall results related to the vocal tract configuration (F1 and F2) show great potential to discriminate between normal and UVFP voices (especially for males) in spite of the localisation of the lesion being at the VF level. This is in agreement with the literature which clearly indicates that the behaviour of the VF is not entirely independent of the vocal tract [55, 58, 59]. Thus, these parameters can add useful information to the assessment procedure and may be used as a complement to the more traditional VF behavioural assessment.

It should be noted that the overall results obtained for females distance themselves from what was initially expected. These differences between genders may be due to a greater technical difficulty in analysing female voices [56, 60]. To a large extent, these difficulties are associated with the identification of formants, due to the fact that F0 is higher, and this increases the difficulty in F1 estimation [56].

4. Conclusions

In this study various ways of assessing the UVFP voice were combined. Since vocal therapy is one of the first noninvasive treatment options with potential to help the client to reacquire a functional voice, it is fundamental to know in detail the alterations created by the pathology at VF and vocal tract level to better guide the treatment. Perceptual differences between normal and UVFP voices were found. The perceptual parameters that better characterised this data of UVFP subjects were Rough (R) and Breathy (B), but altered values of Asthenic (A) and Strained (S) were also found. As far as acoustic parameters are concerned there were no differences in F0 values between normal and UVFP voices in this sample. Jitter (ppq5), shimmer (apq11), HNR, and SD of F0 had an important role in discriminating normal and UVFP voices. Measures related to the vocal tract configuration were also indicative of alterations at VF level; therefore the analysis of formant frequencies values and their SD may have an important role in a clinical setting contributing to a better knowledge of the alterations caused by the vocal pathology. Future work should continue to explore formants and their relation to vocal pathology.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors would like to thank the Otorhinolaryngology Team from Hospital de Santo António and Hospital de São João. This work was partially funded by National Funds through FCT (Foundation for Science and Technology), in the context of the projects UID/CEC/00127/2013 and Incentivo/EEI/UI0127/2014. The Advanced Voice Function Assessment Databases (AVFAD) project is supported by the School of Health Sciences (ESSUA), University of Aveiro, Portugal.

References

A. Blitzen, M. Brin, and L. Ramig, Neurologic Disorders of the Larynx, vol. 101, Thieme, New York, NY, USA, 2nd edition, 2009.

D. M. Hartl, S. Hans, J. Vaissière, M. Riquet, and D. F. Brasnu, “Objective voice quality analysis before and after onset of unilateral vocal fold paralysis,” Journal of Voice, vol. 15, no. 3, pp. 351–361, 2001.

C. M. R. Pinho, L. M. T. Jesus, and A. Barney, “Aerodynamic measures of speech in unilateral vocal fold paralysis (UVFP) patients,” Logopedics Phoniatrics Vocology, vol. 38, no. 1, pp. 19–34, 2013.

L. Sulica and A. Blitzer, Vocal Fold Paralysis, Springer, New York, NY, USA, 2006.

K. Verdolini, C. Rosen, and R. Branski, Classification Manual for Voice Disorders-I, Lawrence Erlbaum Associates, Mahwah, NJ, USA, 2006.

K. Schwarz, C. A. Cielo, N. Steffen, J. Becker, and G. P. Jotz, “Voice and laryngeal configuration of men with unilateral vocal fold paralysis before and after medialization,” Journal of Voice, vol. 25, no. 5, pp. 611–618, 2011.

A. L. Rosenthal, S. Y. Lowell, and R. H. Colton, “Aerodynamic and acoustic features of vocal effort,” Journal of Voice, vol. 28, no. 2, pp. 144–153, 2014.

H. Oguz, M. Demirci, M. A. Safak, N. Arslan, A. Islam, and S. Kargin, “Effects of unilateral vocal cord paralysis on objective voice measures obtained by Praat,” European Archives of Oto-Rhino-Laryngology, vol. 264, no. 3, pp. 257–261, 2007.

D. K. Chhetri, J. Neubauer, J. L. Bergeron, E. Sofer, K. A. Peng, and N. Jamal, “Effects of asymmetric superior laryngeal nerve stimulation on glottic posture, acoustics, vibration,” Laryngoscope, vol. 123, no. 12, pp. 3110–3116, 2013.

C. Madill and P. McCabe, “Acoustic analysis using freeware: praat,” in Handbook of Voice Assessments, Singular, San Diego, Calif, USA, 2011.

A. Vogel, “Multidimensional analysis of voice: computerized speech lab,” in Handbook of Voice Assessments, Singular, San Diego, Calif, USA, 2011.

P. Alku, “Glottal inverse filtering analysis of human voice production—a review of estimation and parameterization methods of the glottal excitation and their applications,” Sadhana, vol. 36, no. 5, pp. 623–650, 2011.

I. Titze, Principles of Voice Production, National Center for Voice and Speech, Iowa City, Iowa, USA, 2nd edition, 2000.

G. Fant, Acoustic Theory of Speech Production—With Calculations Based on X-Ray Studies of Russian Articulations, Mouton, The Hague, The Netherlands, 2nd edition, 1970.

J.-W. Lee, H.-G. Kang, J.-Y. Choi, and Y.-I. Son, “An investigation of vocal tract characteristics for acoustic discrimination of pathological voices,” BioMed Research International, vol. 2013, Article ID 758731, 11 pages, 2013.

A. J. Lotto, L. L. Holt, and K. R. Kluender, “Effect of voice quality on perceived height of English vowels,” Phonetica, vol. 54, no. 2, pp. 76–93, 1997.

L. Gordis, Epidemiology, Elsevier, Philadelphia, Pa, USA, 3rd edition, 2004.

G. Breakwell, J. Smith, and D. Wright, Eds., Research Methods in Psychology, SAGE Publications, London, UK, 4th edition, 2012.

F. Kerlinger and H. Lee, Foundations of Behavioral Research, Wadsworth and Thomson, Orlando, Fla, USA, 4th edition, 2000.

L. Jesus, “University of Aveiro's advanced voice function assessment databases (AVFAD),” Revista de Saúde Pública, vol. 48, p. 291, 2014.

R. M. Roark, B. C. Watson, R. J. Baken, D. J. Brown, and J. M. Thomas, “Measures of vocal attack time for healthy young adults,” Journal of Voice, vol. 26, no. 1, pp. 12–17, 2012.

N. Breslow and N. Day, Statistical Methods in Cancer Research: Vol.1—The Analysis of Case-Control Studies, IARC, Lyon, France, 1980.

J. Beck, “Organic variation of the vocal apparatus,” in The Handbook of Phonetic Sciences, W. Hardcastle, J. Laver, and F. Gibbon, Eds., pp. 155–201, Blackwell Publishing, Oxford, UK, 2nd edition, 2010.

S. Linville, “The aging voice,” in Voice Quality Measurements, R. Kent and M. Ball, Eds., pp. 359–376, Singular, San Diego, Calif, USA, 2000.

I. Chatterjee, H. Halder, S. Bari, S. Kumar, A. Roychoudhury, and P. Murthy, “An analytical study of age and gender effects on voice range profile in bengali adult speakers using phonetogram,” International Journal of Phonosurgery and Laryngology, vol. 1, no. 2, pp. 65–70, 2011.

H. B. Fisher and S. E. Linville, “Acoustic characteristics of women's voices with advancing age,” Journals of Gerontology, vol. 40, no. 3, pp. 324–330, 1985.

E. Perrin, C. Berger-Vachon, and L. Collet, “Acoustical recognition of laryngeal pathology: a comparison of two strategies based on sets of features,” Medical and Biological Engineering and Computing, vol. 37, no. 5, pp. 652–658, 1999.

L. A. Ramig and R. L. Ringel, “Effects of physiological aging on selected acoustic characteristics of voice,” Journal of Speech and Hearing Research, vol. 26, no. 1, pp. 22–30, 1983.

T. Shipp and H. Hollien, “Perception of the aging male voice,” Journal of Speech and Hearing Research, vol. 12, no. 4, pp. 703–710, 1969.

P. Boersma, “Praat, a system for doing phonetics by computer,” Glot International, vol. 5, pp. 341–345, 2001.

B. Boyanov and S. Hadjitodorov, “Acoustic analysis of pathological voices: a voice analysis system for the screening and laryngeal diseases,” IEEE Engineering in Medicine and Biology Magazine, vol. 16, no. 4, pp. 74–82, 1997.

J. G. Švec and S. Granqvist, “Guidelines for selecting microphones for human voice production research,” American Journal of Speech-Language Pathology, vol. 19, no. 4, pp. 356–368, 2010.

J. Burg, “Maximum entropy spectral analysis,” in Proceedings of the 37th Meeting of the Society of Exploration Geophysicists, 1967.

M. Hirano, Clinical Examination of Voice, Springer, New York, NY, USA, 1981.

R. Artstein and M. Poesio, “Inter-coder agreement for computational linguistics,” Computational Linguistics, vol. 34, no. 4, pp. 555–596, 2008.

F. Jalalinajafabadi, C. Gadepalli, F. Ascott, J. Homer, M. Lujan, and B. Cheetham, “Perceptual evaluation of voice quality and its correlation with acoustic measurement,” in Proceedings of the European Modelling Symposium (EMS ’13), pp. 283–286, IEEE, Manchester, UK, November 2013.

M. G. Kendall and B. B. Smith, “The problem of m rankings,” Annals of Mathematical Statistics, vol. 10, no. 3, pp. 275–287, 1939.

C. Moers, B. Möbius, F. Rosanowski, E. Nöth, U. Eysholdt, and T. Haderlein, “Vowel- and text-based cepstral analysis of chronic hoarseness,” Journal of Voice, vol. 26, no. 4, pp. 416–424, 2012.

C. Sellarsa, A. Stantona, A. McConnachiea et al., “Reliability of perceptions of voice quality: evidence from a problem asthma clinic population,” The Journal of Laryngology & Otology, vol. 123, no. 7, pp. 755–763, 2009.

M. P. Karnell, S. D. Melton, J. M. Childes, T. C. Coleman, S. A. Dailey, and H. T. Hoffman, “Reliability of clinician-based (GRBAS and CAPE-V) and patient-based (V-RQOL and IPVI) documentation of voice disorders,” Journal of Voice, vol. 21, no. 5, pp. 576–590, 2007.

M. Hirano, S. Hibi, R. Terasawa, and M. Fujiu, “Relationship between aerodynamic, vibratory, acoustic and psychoacoustic correlates in dysphonia,” Journal of Phonetics, vol. 14, pp. 445–456, 1986.

P. Yu, R. Garrel, R. Nicollas, M. Ouaknine, and A. Giovanni, “Objective voice analysis in dysphonic patients: new data including nonlinear measurements,” Folia Phoniatrica et Logopaedica, vol. 59, no. 1, pp. 20–30, 2007.

M. A. Little, D. A. E. Costello, and M. L. Harries, “Objective dysphonia quantification in vocal fold paralysis: comparing nonlinear with classical measures,” Journal of Voice, vol. 25, no. 1, pp. 21–31, 2011.

S. Bielamowicz and S. V. Stager, “Diagnosis of unilateral recurrent laryngeal nerve paralysis: laryngeal electromyography, subjective rating scales, acoustic and aerodynamic measures,” Laryngoscope, vol. 116, no. 3, pp. 359–364, 2006.

C. Baylor, K. Yorkston, E. Strand, T. Eadie, and J. Duffy, “Measurment of treatment outcome in unilateral vocal fold paralysis: a systematic review,” UVFP Technical Report 5, Academy of Neurologic Communication Disorders and Sciences, Washington, DC, USA, 2005.

S. Pinho, D. Tsuji, and S. Bohadana, Fundamentos em Laringologia e Voz, Guanabara Koogan, Rio de Janeiro, Brasil, 2006.

H. Takahashi, “Assessment of auditory impression of dysphonia,” in Voice Examination Methods, Japan Society of Logopedics and Phoniatrics, Ed., Interna, Tokyo, Japan, 1979.

P. H. Dejonckere, P. Bradley, P. Clemente et al., “A basic protocol for functional assessment of voice pathology, especially for investigating the efficacy of (phonosurgical) treatments and evaluating new assessment techniques: guideline elaborated by the Committee on Phoniatrics of the European Laryngological Society (ELS),” European Archives of Oto-Rhino-Laryngology, vol. 258, no. 2, pp. 77–82, 2001.

N. Yan, L. Wang, and M. L. Ng, “Acoustical analysis of voices produced by Cantonese patients of unilateral vocal fold paralysis acoustical analysis of voices by Cantonese UVFP,” in Proceedings of the IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC '13), pp. 1–5, IEEE, Kunming, China, August 2013.

V. Teles and A. Rosinha, “Acoustic analysis of formants and measures of the sonorous signal disturbance in non-smoker and non-alcoholic women without vocal complaints,” International Archives of Otorhinolaryngology, vol. 12, no. 4, pp. 523–530, 2008.

P. Reijonen, S. Lehikoinen-Söderlund, and H. Rihkanen, “Results of fascial augmentation in unilateral vocal fold paralysis,” Annals of Otology, Rhinology and Laryngology, vol. 111, no. 6, pp. 523–529, 2002.

D. G. Childers and K. Wu, “Gender recognition from speech. Part II: fine analysis,” Journal of the Acoustical Society of America, vol. 90, no. 4, pp. 1841–1856, 1991.

M. Behlau, Voz: O Livro do Especialista, Revinte, Rio de Janeiro, Brazil, 2001.

P. Woo, R. Colton, D. Brewer, and J. Casper, “Functional staging for vocal cord paralysis,” Otolaryngology—Head and Neck Surgery, vol. 105, no. 3, pp. 440–448, 1991.

R. Kent and C. Read, The Acoustic Analysis of Speech, Singular, San Diego, Calif, USA, 2nd edition, 2002.

D. H. Klatt and L. C. Klatt, “Analysis, synthesis, and perception of voice quality variations among female and male talkers,” Journal of the Acoustical Society of America, vol. 87, no. 2, pp. 820–857, 1990.

M. Higashikawa, K. Nakai, A. Sakakura, and H. Takahashi, “Perceived pitch of whispered vowels-relationship with formant frequencies: a preliminary study,” Journal of Voice, vol. 10, no. 2, pp. 155–158, 1996.

M. Rothenberg, “Source-tract acoustic interaction in breathy voice,” in Vocal Fold Fisiology—Biomechanics, Acoustic and Phonatory Control, pp. 465–481, The Denver Center for the Performing Arts, Denver, Colo, USA, 1980.

I. R. Titze and B. H. Story, “Acoustic interactions of the voice source with the lower vocal tract,” Journal of the Acoustical Society of America, vol. 101, no. 4, pp. 2234–2243, 1997.

H. M. Hanson, “Glottal characteristics of female speakers: acoustic correlates,” Journal of the Acoustical Society of America, vol. 101, no. 1, p. 466, 1997.