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

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.

Luis M. T. Jesus 1, 2 and Joana Martinez 1 and Andreia Hall 1, 3 and Anibal Ferreira 4

Academic Editor:Haldun Oguz

1, Institute of Electronics and Informatics Engineering of Aveiro (IEETA), University of Aveiro, 3810-193 Aveiro, Portugal
2, School of Health Sciences (ESSUA), University of Aveiro, 3810-193 Aveiro, Portugal
3, Department of Mathematics (DMat), University of Aveiro, 3810-193 Aveiro, Portugal
4, Department of Electrical and Computer Engineering, University of Porto, 4200-465 Porto, Portugal

Received 22 January 2015; Accepted 24 April 2015; 18 October 2015

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.

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 (F 0) 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 (F 1) frequency and higher values of the second formant (F 2) 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 F 1 and F 2 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 F 1 and F 2 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 Antonio 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 F 0. 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 F 0; standard-deviation of F 0; 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 F 1 and F 2. 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 [figure omitted; refer to PDF] Coefficient. The Mann-Whitney [figure omitted; refer to PDF] test was used to analyse the GRBAS scale parameters. The acoustic parameters that had normal distribution (HNR, F 2[female], standard-deviation of F 0[male], F 1[male]) were statistically analysed using the [figure omitted; refer to PDF] -test and parameters that did not have normal distribution (Jitter (ppq5), Shimmer (apq11), F 0[female], standard-deviation of F 0[female], F 1[female], standard-deviation of F 1[female], standard-deviation of F 2[female], F 0[male], standard-deviation of F 1[male], F 2[male], and standard-deviation of F 2[male]) were analysed with the Mann-Whitney [figure omitted; refer to PDF] 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 Antonio 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 [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] '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 [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] test.

Scale parameter	[figure omitted; refer to PDF]	[figure omitted; refer to PDF] 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; [figure omitted; refer to PDF] : Kendall's [figure omitted; refer to PDF] .

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 [figure omitted; refer to PDF] 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	[figure omitted; refer to PDF]	[figure omitted; refer to PDF] value
	[figure omitted; refer to PDF]	Mean ± SD	[figure omitted; refer to PDF]	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; [figure omitted; refer to PDF] : number of cases; SD: standard deviation; [figure omitted; refer to PDF] : Mann-Whitney [figure omitted; refer to PDF] test.

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

[figure omitted; refer to PDF]

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	[figure omitted; refer to PDF] or [figure omitted; refer to PDF]	[figure omitted; refer to PDF] value
	[figure omitted; refer to PDF]	Mean ± SD	[figure omitted; refer to PDF]	Mean ± SD
Jitter ppq5 (%)	17	1.06 ± 1.02	85	0.26 ± 0.18	[figure omitted; refer to PDF]	<0.001
Shimmer apq11 (%)	17	10.16 ± 3.34	85	7.14 ± 3.22	[figure omitted; refer to PDF]	0.001
HNR (dB)	17	10.11 ± 4.94	85	14.75 ± 4.46	[figure omitted; refer to PDF]	<0.001

UVFP: Unilateral Vocal Fold Paralysis; SD: standard deviation; [figure omitted; refer to PDF] : [figure omitted; refer to PDF] -test; U: Mann-Whitney [figure omitted; refer to PDF] test.

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

[figure omitted; refer to PDF]

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

[figure omitted; refer to PDF]

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

[figure omitted; refer to PDF]

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 F 0 (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 F 0. Also according to Woo et al. [54] the majority of UVFP subjects present F 0 values close to normal.

Table 5: Fundamental frequency ( [figure omitted; refer to PDF] ) and first and second formant frequencies ( [figure omitted; refer to PDF] and [figure omitted; refer to PDF] ) and their standard-deviations, for normal and UFVP female participants.

[female]	UVFP	Normal	[figure omitted; refer to PDF] or [figure omitted; refer to PDF]	[figure omitted; refer to PDF] value
[figure omitted; refer to PDF]	Mean ± SD	[figure omitted; refer to PDF]	Mean ± SD
[figure omitted; refer to PDF] (Hz)	13	218.38 ± 72.36	65	195.36 ± 33.02	[figure omitted; refer to PDF]	0.335
SD [figure omitted; refer to PDF] (Hz)	13	6.65 ± 12.28	65	2.62 ± 2.15	[figure omitted; refer to PDF]	0.018
[figure omitted; refer to PDF] (Hz)	13	826.16 ± 171.73	65	819.03 ± 164.75	[figure omitted; refer to PDF]	0.495
SD [figure omitted; refer to PDF] (Hz)	13	117.18 ± 99.91	65	71.28 ± 53.54	[figure omitted; refer to PDF]	0.116
[figure omitted; refer to PDF] (Hz)	13	1522.69 ± 96.00	65	1453.51 ± 139.20	[figure omitted; refer to PDF] = 1.58	0.059
SD [figure omitted; refer to PDF] (Hz)	13	156.31 ± 146.20	65	62.26 ± 48.67	[figure omitted; refer to PDF]	0.002

UVFP: Unilateral Vocal Fold Paralysis; [figure omitted; refer to PDF] : number of cases; SD: standard deviation; [figure omitted; refer to PDF] : [figure omitted; refer to PDF] -test; [figure omitted; refer to PDF] : Mann-Whitney [figure omitted; refer to PDF] test; [figure omitted; refer to PDF] : fundamental frequency; SD [figure omitted; refer to PDF] : standard-deviation of the fundamental frequency; [figure omitted; refer to PDF] : first formant frequency; SD [figure omitted; refer to PDF] : standard-deviation of first formant frequency; [figure omitted; refer to PDF] : second formant frequency; SD [figure omitted; refer to PDF] : standard-deviation of second formant frequency.

Table 6: Fundamental frequency ( [figure omitted; refer to PDF] ) and first and second formant frequencies ( [figure omitted; refer to PDF] and [figure omitted; refer to PDF] ) and their standard-deviations, for normal and UFVP male participants.

[male]	UVFP	Normal	[figure omitted; refer to PDF] or [figure omitted; refer to PDF]	[figure omitted; refer to PDF] value
[figure omitted; refer to PDF]	Mean ± SD	[figure omitted; refer to PDF]	Mean ± SD
[figure omitted; refer to PDF] (Hz)	4	121.43 ± 12.70	20	128.27 ± 23.85	[figure omitted; refer to PDF]	0.485
SD [figure omitted; refer to PDF] (Hz)	4	3.41 ± 1.25	20	1.36 ± 0.58	[figure omitted; refer to PDF]	<0.001
[figure omitted; refer to PDF] (Hz)	4	821.48 ± 331.80	20	677.33 ± 84.95	[figure omitted; refer to PDF]	0.043
SD [figure omitted; refer to PDF] (Hz)	4	191.91 ± 105.62	20	30.58 ± 18.84	[figure omitted; refer to PDF]	<0.001
[figure omitted; refer to PDF] (Hz)	4	1629.09 ± 474.06	20	1282.95 ± 104.35	[figure omitted; refer to PDF]	0.011
SD [figure omitted; refer to PDF] (Hz)	4	263.68 ± 144.45	20	36.89 ± 27.53	[figure omitted; refer to PDF]	0.001

UVFP: Unilateral Vocal Fold Paralysis; [figure omitted; refer to PDF] : number of cases; SD: standard deviation; [figure omitted; refer to PDF] : [figure omitted; refer to PDF] -test; [figure omitted; refer to PDF] : Mann-Whitney [figure omitted; refer to PDF] test; [figure omitted; refer to PDF] : fundamental frequency; SD [figure omitted; refer to PDF] : standard-deviation of the fundamental frequency; [figure omitted; refer to PDF] : first formant frequency; SD [figure omitted; refer to PDF] : standard-deviation of first formant frequency; [figure omitted; refer to PDF] : second formant frequency; SD [figure omitted; refer to PDF] : standard-deviation of second formant frequency.

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

[figure omitted; refer to PDF]

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

[figure omitted; refer to PDF]

The standard-deviation of F 0 (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 F 0 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 F 0 (Hz) values for female UVFP and normal subjects.

[figure omitted; refer to PDF]

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

[figure omitted; refer to PDF]

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 F 1 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 F 1 frequencies in UVFP cases were expected (based on data reported previously [15]); however, Table 6 clearly shows that the F 1 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 F 1 frequency values for voices with similar characteristics to UVFP patients. Since the frequency of F 1 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: F 1 and F 2 frequency (Hz) values for female UVFP and normal subjects.

[figure omitted; refer to PDF]

Figure 10: F 1 and F 2 frequency (Hz) values for male UVFP and normal subjects.

[figure omitted; refer to PDF]

The second formant (F 2) 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 [figure omitted; refer to PDF] 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 F 2 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 F 1 shown in Table 5, there were significant differences between the two groups being SD of F 1 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 F 2 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 F 2, 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 F 1 and F 2 frequency (Hz) for female UVFP and normal subjects.

[figure omitted; refer to PDF]

Figure 12: SD of F 1 and F 2 frequency (Hz) for male UVFP and normal subjects.

[figure omitted; refer to PDF]

Overall results related to the vocal tract configuration (F 1 and F 2) 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 F 0 is higher, and this increases the difficulty in F 1 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 F 0 values between normal and UVFP voices in this sample. Jitter (ppq5), shimmer (apq11), HNR, and SD of F 0 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.

Acknowledgments

The authors would like to thank the Otorhinolaryngology Team from Hospital de Santo Antonio 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.

Conflict of Interests

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

References

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

[2] D. M. Hartl, S. Hans, J. Vaissière, M. Riquet, 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.

[3] C. M. R. Pinho, L. M. T. Jesus, A. Barney, "Aerodynamic measures of speech in unilateral vocal fold paralysis (UVFP) patients," Logopedics Phoniatrics Vocology , vol. 38, no. 1, pp. 19-34, 2013.

[4] L. Sulica, A. Blitzer Vocal Fold Paralysis , Springer, New York, NY, USA, 2006.

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

[6] K. Schwarz, C. A. Cielo, N. Steffen, J. Becker, 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.

[7] A. L. Rosenthal, S. Y. Lowell, R. H. Colton, "Aerodynamic and acoustic features of vocal effort," Journal of Voice , vol. 28, no. 2, pp. 144-153, 2014.

[8] H. Oguz, M. Demirci, M. A. Safak, N. Arslan, A. Islam, 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.

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

[10] C. Madill, P. McCabe, "Acoustic analysis using freeware: praat," Handbook of Voice Assessments , Singular, San Diego, Calif, USA, 2011.

[11] A. Vogel, "Multidimensional analysis of voice: computerized speech lab," Handbook of Voice Assessments , Singular, San Diego, Calif, USA, 2011.

[12] 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.

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

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

[15] J.-W. Lee, H.-G. Kang, J.-Y. Choi, Y.-I. Son, "An investigation of vocal tract characteristics for acoustic discrimination of pathological voices," BioMed Research International , vol. 2013, 2013.

[16] A. J. Lotto, L. L. Holt, K. R. Kluender, "Effect of voice quality on perceived height of English vowels," Phonetica , vol. 54, no. 2, pp. 76-93, 1997.

[17] L. Gordis Epidemiology , Elsevier, Philadelphia, Pa, USA, 2004., 3rd.

[18] G. Breakwell, J. Smith, D. Wright Research Methods in Psychology , SAGE Publications, London, UK, 2012., 4th.

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

[20] L. Jesus, "University of Aveiro's advanced voice function assessment databases (AVFAD)," Revista de Saúde Pública , vol. 48, pp. 291, 2014.

[21] R. M. Roark, B. C. Watson, R. J. Baken, D. J. Brown, J. M. Thomas, "Measures of vocal attack time for healthy young adults," Journal of Voice , vol. 26, no. 1, pp. 12-17, 2012.

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

[23] J. Beck, W. Hardcastle, J. Laver, F. Gibbon, "Organic variation of the vocal apparatus," The Handbook of Phonetic Sciences , pp. 155-201, Blackwell Publishing, Oxford, UK, 2010., 2nd.

[24] S. Linville, R. Kent, M. Ball, "The aging voice," Voice Quality Measurements , pp. 359-376, Singular, San Diego, Calif, USA, 2000.

[25] I. Chatterjee, H. Halder, S. Bari, S. Kumar, A. Roychoudhury, 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.

[26] H. B. Fisher, S. E. Linville, "Acoustic characteristics of women's voices with advancing age," Journals of Gerontology , vol. 40, no. 3, pp. 324-330, 1985.

[27] E. Perrin, C. Berger-Vachon, 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.

[28] L. A. Ramig, 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.

[29] T. Shipp, H. Hollien, "Perception of the aging male voice," Journal of Speech and Hearing Research , vol. 12, no. 4, pp. 703-710, 1969.

[30] P. Boersma, "Praat, a system for doing phonetics by computer," Glot International , vol. 5, pp. 341-345, 2001.

[31] B. Boyanov, 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.

[32] J. G. Svec, 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.

[33] J. Burg, "Maximum entropy spectral analysis," in Proceedings of the 37th Meeting of the Society of Exploration Geophysicists, 1967.

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

[35] R. Artstein, M. Poesio, "Inter-coder agreement for computational linguistics," Computational Linguistics , vol. 34, no. 4, pp. 555-596, 2008.

[36] F. Jalalinajafabadi, C. Gadepalli, F. Ascott, J. Homer, M. Lujan, 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.

[37] M. G. Kendall, B. B. Smith, "The problem of m rankings," Annals of Mathematical Statistics , vol. 10, no. 3, pp. 275-287, 1939.

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

[39] C. Sellarsa, A. Stantona, A. McConnachiea, C. Dunneta, L. Chapmana, C. Bucknalla, K. MacKenziea, "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.

[40] M. P. Karnell, S. D. Melton, J. M. Childes, T. C. Coleman, S. A. Dailey, 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.

[41] M. Hirano, S. Hibi, R. Terasawa, M. Fujiu, "Relationship between aerodynamic, vibratory, acoustic and psychoacoustic correlates in dysphonia," Journal of Phonetics , vol. 14, pp. 445-456, 1986.

[42] P. Yu, R. Garrel, R. Nicollas, M. Ouaknine, 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.

[43] M. A. Little, D. A. E. Costello, 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.

[44] S. Bielamowicz, 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.

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

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

[47] H. Takahashi, Japan Society of Logopedics and Phoniatrics, "Assessment of auditory impression of dysphonia," Voice Examination Methods , Interna, Tokyo, Japan, 1979.

[48] P. H. Dejonckere, P. Bradley, P. Clemente, G. Cornut, L. Crevier-Buchman, G. Friedrich, P. Van De Heyning, M. Remacle, V. Woisard, "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.

[49] N. Yan, L. Wang, 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.

[50] V. Teles, 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.

[51] P. Reijonen, S. Lehikoinen-Söderlund, 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.

[52] D. G. Childers, 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.

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

[54] P. Woo, R. Colton, D. Brewer, J. Casper, "Functional staging for vocal cord paralysis," Otolaryngology-Head and Neck Surgery , vol. 105, no. 3, pp. 440-448, 1991.

[55] R. Kent, C. Read The Acoustic Analysis of Speech , Singular, San Diego, Calif, USA, 2002., 2nd.

[56] D. H. Klatt, 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.

[57] M. Higashikawa, K. Nakai, A. Sakakura, 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.

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

[59] I. R. Titze, 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.

[60] H. M. Hanson, "Glottal characteristics of female speakers: acoustic correlates," Journal of the Acoustical Society of America , vol. 101, no. 1, pp. 466, 1997.