International Journal of OtolaryngologyVolume 2012 (2012), Article ID 739523, 4 pagesdoi:10.1155/2012/739523

Clinical Study

Nonlinear Dynamic Analysis of Vowels in Cleft Palate Children with or without Hypernasality

Katsuaki Mishima,¹ Hiroyuki Nakano,¹ Tatsushi Matsumura,² Norifumi Moritani,² Seiji Iida,² and Yoshiya Ueyama¹

¹The Department of Oral and Maxillofacial Surgery, Graduate School of Medicine, Yamaguchi University, Minami-kogushi 1-1-1, Ube City, Yamaguchi 7558505, Japan²The Department of Oral and Maxillofacial Reconstructive Surgery, Graduate School of Medicine, Dentistry and Pharmaceutical Sciences, Okayama University, Shikata-cho 2-5-1, Okayama City, Okayama 7008525, Japan

Received 31 July 2011; Accepted 25 October 2011

Academic Editor: David R. White

Copyright © 2012 Katsuaki Mishima 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.

Objectives. To clarify the difference between Lyapunov exponents (LEs) for cleft palate (CP) patients with hypernasality versus without hypernasality and to investigate the relationship between their LEs and nasalance scores (NSs). Material and Methods. Six CP patients with severe hypernasality (mean age 9.2 years) and six CP patients without hypernasality (mean age 8.0 years) were enrolled. Five Japanese vowels were recorded at 44.1 KHz, and the NSs were measured simultaneously. The mean first LE (mLE₁) from all one-second intervals was computed. Results. The mLE₁ for /o/ in patients with hypernasality was significantly higher than that in patients without hypernasality. The correlation coefficients between the mLE₁ and NS for all vowels were not statistically different. Conclusion. The voice signal of /o/ for the patients with hypernasality was more instable than in those without hypernasality. The chaotic phenomenon was independent of nasal resonance in CP voice.

1. Introduction

Patients with cleft palate often exhibit nasality, which is a distinctive feature and an important target in speech therapy and rehabilitation. To evaluate velopharyngeal function, the aerodynamic and acoustic aspects of nasalization have been studied. An aerodynamic exam can diagnose the degree of velopharyngeal closure [1, 2], and acoustic measurements can categorize velopharyngeal insufficiency [3, 4]. The abnormal resonance generated by velopharyngeal insufficiency can be evaluated quantitatively using a nasometer [5].

On the other hand, the voice and speech of patients with cleft palates have been studied using many techniques including spectral analysis, perturbation analysis, and formant analysis. Zajac and Linville [6] and Lewis et al. [7] reported that cleft palate speakers have larger frequency perturbations (jitter) than normal controls. However, the methods used to calculate perturbations, jitter, and shimmer are only reliable for nearly periodic voice signals and cannot reliably analyze strongly aperiodic signals [8]. Recently, nonlinear dynamic methods have enabled the quantification of aperiodic and chaotic phenomena [9–11].

In our previous paper, we reported that the Lyapunov exponents (LEs) of the vowels /a/, /e/ and /o/ for adult cleft palate patients are higher than those for normal resonance adults and that there were no correlation coefficients between LEs and nasalance scores (NSs) [12]. These results suggested that vocal fold vibration may be less stable in adult cleft palate patients than in normal resonance subjects and that the LE may be a parameter independent of resonance. Subsequently, we investigated the nonlinear dynamic characteristics of cleft palate speech and voice. In the present paper, the purpose was to clarify the difference between the LEs for cleft palate patients with hypernasality versus without hypernasality and to investigate the relationship between their LEs and NSs.

2. Materials and Methods

2.1. Patients

Six repaired cleft palate patients with severe hypernasality (mean age 9.2 years; range 6 to 13, 2 boys and 4 girls) and six repaired cleft palate patients without hypernasality (mean age 8.0 years; range 6 to 13, 4 boys and 2 girls) were enrolled. The presence of hypernasality was perceptually judged by two speech therapists from Okayama University Hospital. The present study, which was approved by the Okayama University Institutional Ethical Board, was carried out after obtaining informed consents from the parents of all participants.

2.2. Methods

The voices were recorded through a microphone (Shure BG1.1, Niles, Ill) on a portable solid-state recorder (Marantz PMD 640, Itasca, Ill) with a nasometer II headset (model 6400, Kay Elemetrics Corp., Lincoln Park, NJ) in a quiet room designated for speech therapy in the Okayama University Dental Hospital. The nasalance scores (NSs) were measured simultaneously. The voice samples were recorded on the compact flush medium of the recorder at a sampling rate of 44.1 KHz, at 16 bits, in a *.wav file format. The Japanese vowels /a/; [a], /i/; [i], /u/; [ɯ], /e/; [e], and /o/; [o] were used as voice samples. Each vowel was naturally phonated during approximately one-second three times.

The voice data were processed on a personal computer (NEC mate MA30Y, Tokyo) with a modified Chaos Analyzing Program (ver. 1.0.4, CCI Corporation, Fukuoka), which used the algorithm from Sano and Sawada [13]. The first Lyapunov exponent (LE₁) was computed for each one second interval, while the interval was being shifted by 100 msec. The mean first Lyapunov exponents for all intervals (mLE₁) were then calculated.

The delay time was determined as follows. The first zero-crossing points of autocorrelation were calculated for each vowel. As a result, the delay time was estimated at 15, 32, 27, 22, and 21 points (1 point = 1/44.1 msec.) for the vowels /a/, /i/, /u/, /e/, and /o/, respectively (Table 1). The embedding dimension was determined as follows. The fractal dimensions were computed using the Grassberger-Procaccia algorithm [14], and convergent diagrams, in which the embedding dimensions were assumed, were then constructed for each vowel. Thus, the embedding dimensions were estimated at 5 for all vowels (Table 1).

Table 1: First zero-crossing points of autocorrelation and estimated embedding dimensions.

Hypernasality (+)
6

Hypernasality (−)
6

P value

0.429
1.000
0.931
0.931
0.662

0.063
0.190
0.286
0.063
0.111

Total
12

Upper column: first zero-crossing points of autocorrelation (unit point; 1 point = 1/44.1 msec.),

lower column: estimated embedding dimensions.

The differences of the first zero-crossing points of autocorrelation, estimated the embedding dimensions, and the mLE₁s between the two groups with versus without hypernasality were analyzed statistically using the Mann-Whitney test. The correlation coefficients between the mLE₁ and NS were also calculated for each vowel. The statistical package SPSS (ver.16.0) was utilized, and differences with values of less than 0.05 were considered to be statistically significant.

3. Results

There were no significant differences between the first zero-crossing points and the estimated embedding dimensions of those patients with or without hypernasality for all vowels (Table 1). The mLE₁ for /o/ in the patients with hypernasality was significantly higher than in patients without hypernasality () (Table 2). The NSs for /i/, /u/, /e/, and /o/ in the patients with hypernasality were significantly higher than in patients without hypernasality (Table 3). The correlation coefficients between the mLE₁ and NS for all vowels were not statistically different (Table 4).

Table 2: Lyapunov exponents (mLE₁).

Hypernasality (+)
6

Hypernasality (−)
6

P value

0.093
0.310
0.589
0.699
0.015

mLE₁; mean of the first Lyapunov exponents.

Table 3: Nasalance scores (NSs).

Hypernasality (+)
6

Hypernasality (−)
6

P value

0.240
0.026
0.026
0.026
0.041

Table 4: Correlation coefficients between the mLE₁ (mean of the first Lyapunov exponents) and NS (nasalance score).

mLE₁ of /a/
mLE₁ of /i/
mLE₁ of /u/
mLE₁ of /e/
mLE₁of /o/

NS of /a/
0.177

0.582

NS of /i/

0.430

0.163

NS of /u/

0.420

0.174

NS of /e/

0.131

0.685

NS of /o/

0.424

0.170

Upper column: correlation coefficients, lower column: P value.

4. Discussion

Although nasality can be evaluated using a spectral analysis of speech signals, voice acoustic measures of nasality are not universally used in clinical or empirical work because of ambiguity in the literature regarding the appropriate acoustic methodology, the amount of labor involved as compared with the Nasometer, and so forth [15]. However, Vogel et al. [15] demonstrated the potential for the wider application of acoustic investigation into nasality.

Several authors have described laryngeal disorders, including organic and functional disorders, in cleft palate speakers [16, 17]. Zajac and Linville [6] and Lewis et al. [7] reported that cleft palate speakers have higher frequency perturbations (jitter) than normal controls. Nicollas et al. [18] demonstrated that neither jitters nor shimmers significantly differed with age or gender. Van Lierde et al. [19] reported a multiparameter approach to vocal quality but stated that the nature of the vocal quality and the voice range measurement differences cannot be explained from their study. Therefore, we concluded that future studies on the voice of cleft palate subjects using nonlinear analysis may be beneficial in gaining further insight into the mechanics of phonation.

To our knowledge, there have been no reports on the application of nonlinear dynamic analysis to cleft palate speech. Our previous study demonstrated that the mLE₁ for /a/, in both males and females with CP, is significantly higher than in normal resonance individuals and that the mLE₁ for /e/ in males with CP and for /o/ in females with CP are significantly higher than in normal resonance individuals [12]. Since the mLE1 is a measure of the instability of the voice signal, these results suggest that the vocal fold vibration is less stable in CP speakers than in normal resonance subjects. In addition, the correlation coefficients between the mLE₁ and NS for all vowels were not statistically different in both normal and CP subjects. Therefore, the Lyapunov exponents may be independent of resonance [12]. Subsequently, the nonlinear dynamic characteristics were investigated in the present study. The mLE₁ for /o/ in the patients with hypernasality was significantly higher than in patients without hypernasality; in other words, the voice signal of /o/ for the patients with hypernasality was more instable than in those without hypernasality. This may contribute to the instability of the vocal fold. On the other hand, the correlation coefficients between the mLE₁ and NS for all vowels were not statistically different in patients both with versus without hypernasality. This supported the independence of chaotic phenomenon and nasal resonance in cleft palate speech and voice, which was demonstrated in our previous paper [12].

Nicollas et al. [18] reported that the large LE seems to decrease with age from their studies of children between 6 and 12 years of age. It was also suggested that the large LE is lower in boys than in girls overall but varies for each age [18]. In our present study, the boys and girls were not separated because of the small sample size. A further investigation is necessary in this respect.

5. Conclusions

The voice signal of /o/ for the patients with hypernasality was more instable than in those without hypernasality. The chaotic phenomenon was independent of nasal resonance in cleft palate speech and voice.

Conflict of Interests

The authors declare that they have no conflict of interests.

Acknowledgment

The authors thank Mr. Shozo Hirose from Ergovision Co. Ltd. (Osaka) for his advice.

References

D. L. Jones, “The relationship between temporal aspects of oral-nasal balance and classification of velopharyngeal status in speakers with cleft palate,” Cleft Palate-Craniofacial Journal, vol. 37, no. 4, pp. 363–369, 2000.

D. L. Jones, H. L. Morris, and D. R. Van Demark, “A comparison of oral-nasal balance patterns in speakers who are categorized as “almost but not quite” and ‘sometimes but not always’,” Cleft Palate-Craniofacial Journal, vol. 41, no. 5, pp. 526–534, 2004.

D. W. Warren, R. M. Dalston, W. C. Trier, and M. B. Holder, “A pressure-flow technique for quantifying temporal patterns of palatopharyngeal closure,” Cleft Palate Journal, vol. 22, no. 1, pp. 11–19, 1985.

D. W. Warren, R. M. Dalston, K. E. Morr, W. M. Hairfield, and L. R. Smith, “The speech regulating system: temporal and aerodynamic responses to velopharyngeal inadequacy,” Journal of Speech and Hearing Research, vol. 32, no. 3, pp. 566–575, 1989.

S. G. Fletcher and M. E. Bishop, “Measurement of nasality with tonar,” Cleft Palate Journal, vol. 7, pp. 610–621, 1970.

D. J. Zajac and R. N. Linville, “Voice perturbations of children with perceived nasality and hoarseness,” Cleft Palate Journal, vol. 26, no. 3, pp. 226–232, 1989.

J. R. Lewis, M. L. Andreassen, H. A. Leeper, D. L. Macrae, and J. Thomas, “Vocal characteristics of children with cleft lip/palate and associated velopharyngeal incompetence,” Journal of Otolaryngology, vol. 22, no. 2, pp. 113–117, 1993.

Y. Zhang and J. J. Jiang, “Acoustic analyses of sustained and running voices from patients with laryngeal pathologies,” Journal of Voice, vol. 22, no. 1, pp. 1–9, 2008.

A. Giovanni, M. Ouaknine, and J. M. Triglia, “Determination of largest Lyapunov exponents of vocal signal: application to unilateral laryngeal paralysis,” Journal of Voice, vol. 13, no. 3, pp. 341–354, 1999.

Y. Zhang, C. McGilligan, L. Zhou, M. Vig, and J. J. Jiang, “Nonlinear dynamic analysis of voices before and after surgical excision of vocal polyps,” Journal of the Acoustical Society of America, vol. 115, no. 5, pp. 2270–2277, 2004.

J. J. Jiang, Y. Zhang, and C. McGilligan, “Chaos in voice, from modeling to measurement,” Journal of Voice, vol. 20, no. 1, pp. 2–17, 2006.

K. Mishima, A. Sugii, T. Yamada, T. Matsumura, and N. Moritani, “Lyapunov exponents for Japanese vowels in normal adult and cleft palate speakers,” Folia Phoniatrica et Logopaedica, vol. 63, no. 3, pp. 129–133, 2010.

M. Sano and Y. Sawada, “Measurement of the lyapunov spectrum from a chaotic time series,” Physical Review Letters, vol. 55, no. 10, pp. 1082–1085, 1985.

P. Grassberger and I. Procaccia, “Characterization of strange attractors,” Physical Review Letters, vol. 50, no. 5, pp. 346–349, 1983.

A. P. Vogel, H. M. Ibrahim, S. Reilly, and N. Kilpatrick, “A comparative study of two acoustic measures of hypernasality,” Journal of Speech, Language, and Hearing Research, vol. 52, no. 6, pp. 1640–1651, 2009.

S. B. Leder and J. W. Lerman, “Some acoustic evidence for vocal abuse in adult speakers with repaired cleft palate,” Laryngoscope, vol. 95, no. 7, pp. 837–840, 1985.

L. L. D'Antonio, H. R. Muntz, M. A. Province, and J. L. Marsh, “Laryngeal/voice findings in patients with velopharyngeal dysfunction,” Laryngoscope, vol. 98, no. 4, pp. 432–438, 1988.

R. Nicollas, R. Garrel, M. Ouaknine, A. Giovanni, B. Nazarian, and J. M. Triglia, “Normal voice in children between 6 and 12 years of age: database and nonlinear analysis,” Journal of Voice, vol. 22, no. 6, pp. 671–675, 2008.

K. M. Van Lierde, S. Claeys, M. De Bodt, and P. Van Cauwenberge, “Vocal quality characteristics in children with cleft palate: a multiparameter approach,” Journal of Voice, vol. 18, no. 3, pp. 354–362, 2004.