3  Signal Typing

Note

Living draft. This chapter presents the classification that gates every measure in Parts II to IV. Later chapters refer to the types by number in their Caveats sections. All example signals are synthetic; the code is shown under each figure and is MIT licensed. Found an error? Use the Report an issue link in the sidebar.

3.1 Why signal typing comes first

F0 is the rate at which a cyclic waveform repeats. An aperiodic signal has no repeating waveform and therefore has no F0 (Baken and Orlikoff 2000, 152). Jitter, shimmer and harmonics-to-noise ratio all start from an identified F0 or period, so they inherit this limit. The acoustic signals of many voice disorders do not have a periodicity that allows automatic measurement of F0 or of cycle amplitudes (Ludlow et al. 2018, 378).

Analysis software does not protect the user from this. A computer-based system may calculate a putative F0 for any signal, even an aperiodic one, and the number it displays need not have any meaning (Baken and Orlikoff 2000, 152). Checking the signal is the user’s responsibility. It is often skipped: Ludlow, Kent and Gray report a finding by Behrman (2005) that 60 % of clinicians using acoustic analysis did not examine the voice signal first to see whether periodic analysis was possible (Ludlow et al. 2018, 378).

Signal typing is that check. It sorts voice signals by how periodic they are over the interval to be analysed, so that each measure is applied only where it is valid.

3.2 Origin

The classification comes from the Workshop on Acoustic Voice Analysis, held on 17 and 18 February 1994 in Denver, Colorado, whose summary statement was written by Titze (Titze 1995, 2). The statement presents the classification as central to all other considerations in acoustic voice analysis, and as following the general principles of the nonlinear dynamics of coupled oscillators (Titze 1995, 19). Its first recommendation is that it is useful to classify acoustic voice signals into three types (Titze 1995, 26).

Baken and Orlikoff adopted Titze’s three categories for clinical use and added subtypes within each (Baken and Orlikoff 2000, 152). Sprecher and colleagues (2010) proposed a fourth type, for signals predominated by stochastic noise, which are unsuitable for both perturbation and correlation dimension analysis (Sprecher et al. 2010, 3710). They classified 40 voice samples into the four types from narrowband spectrograms, and found that perturbation, perceptual and correlation dimension measures all worsened from Type 1 through Type 4 (Sprecher et al. 2010, 3710).

3.3 The types

The definitions for Types 1 to 3 are Titze’s (Titze 1995, 19). The subtypes are Baken and Orlikoff’s addition and are shown in italics (Baken and Orlikoff 2000, 152).

3.3.1 Type 1: nearly periodic

Nearly periodic (or, in the ideal case, periodic) waveforms that do not undergo qualitative changes during the interval analysed. Any modulating or subharmonic frequencies present have at least one order of magnitude less energy than the fundamental (Titze 1995, 19; Baken and Orlikoff 2000, 152).

  • Subtype A: minor random variation of F0 or waveshape.
  • Subtype B: minor monotonic modulation, with F0 or amplitude steadily increasing or decreasing.
  • Subtype C: minor cyclic variation (modulation) of F0 or amplitude.

(Baken and Orlikoff 2000, 152)

3.3.2 Type 2: organised, but no single F0

Signals with sudden qualitative changes (bifurcations) in the interval analysed, or with modulating or subharmonic frequencies whose energy is comparable to that of the fundamental. Type 2 signals are clearly ordered and organised, but no single F0 value validly characterises the whole segment (Titze 1995, 19; Baken and Orlikoff 2000, 153).

  • Subtype A: signal bifurcation(s) present.
  • Subtype B: F0 or amplitude discontinuity.
  • Subtype C: large-scale modulation of F0.

(Baken and Orlikoff 2000, 153)

A common Type 2 pattern is period doubling. Successive cycles alternate between two shapes, which puts a component at half the fundamental (F0/2) into the spectrum. The waveform then has two defensible periods: one per oscillation or one per pair (Baken and Orlikoff 2000, 151).

3.3.3 Type 3: no apparent periodicity

Signals without apparent periodicity, and hence not validly susceptible to F0 measurement. A Type 3 signal may nonetheless be chaotic, with a covert organising principle that the techniques of dynamical systems theory might reveal (Titze 1995, 19; Baken and Orlikoff 2000, 153).

  • Subtype A: no observable structure.
  • Subtype B: deep-level (chaotic) organisation.

(Baken and Orlikoff 2000, 153)

3.3.4 Type 4: dominated by noise

Type 4 was added by Sprecher and colleagues for signals predominated by stochastic noise characteristics, split off from Type 3 (Sprecher et al. 2010, 3711). Such signals have an infinite correlation dimension, which is what makes them unsuitable for the nonlinear methods that can still describe a Type 3 signal (Sprecher et al. 2010, 3710).

On a narrowband spectrogram a Type 4 signal shows energy smeared across a broad range of frequencies, as with broadband white noise; an F0 may or may not be visible (Sprecher et al. 2010, 3711). A Type 3 signal, by contrast, shows smearing across the harmonics with the energy concentrated at lower frequencies, generally below about 1500 Hz, and often still shows F0 and one or two harmonics (Sprecher et al. 2010, 3711).

Adding Type 4 also narrows Type 3. Sprecher and colleagues define Type 3 as chaotic with a finite dimension, where Titze defined it as any signal with no apparent periodic structure (Sprecher et al. 2010, 3711; Titze 1995, 19). The four-type scheme is a subdivision of Titze’s Type 3 rather than a disagreement about the signals themselves.

3.4 What the types look like

Figure 3.1 shows one synthetic example of each type. Each is a pulse train passed through the same vowel filter; only the cycle-to-cycle pattern differs. The spectrograms use a long window, so harmonics show as horizontal lines when they exist.

Show the code for this figure
split_t = 0.4 * 0.6
fig, axes = plt.subplots(4, 2, figsize=(7.2, 7.6),
                         gridspec_kw={"width_ratios": [1, 1.35]})
for row, (name, x) in enumerate(examples.items()):
    aw, asg = axes[row]
    t0 = 0.40  # waveform excerpt start (s), after the Type 2 bifurcation
    seg = x[int(t0 * FS): int((t0 + 0.06) * FS)]
    aw.plot(np.arange(len(seg)) / FS * 1e3, seg, color=vs.INK, lw=0.8)
    aw.set_ylim(-1.05, 1.05)
    aw.set_yticks([])
    aw.set_title(name, color=vs.INK)
    vs.style_axes(aw, grid=False)
    aw.spines["left"].set_visible(False)

    vs.plot_spectrogram(asg, x, FS, win_s=0.050, fmax=1500, hop_s=0.002, dyn_db=50)
    asg.set_ylim(0, 1.5)
    asg.set_ylabel("kHz")
    if name == "Type 2":
        asg.axvline(split_t, color=vs.ORANGE, lw=1.2, ls="--")
        asg.text(split_t + 0.01, 1.38, "bifurcation", color=vs.INK, fontsize=8)
    vs.style_axes(asg, grid=False)
axes[-1, 0].set_xlabel("Time (ms)")
axes[-1, 1].set_xlabel("Time (s)")
fig.tight_layout()
plt.show()
Four rows, one per signal type. Each row has a short waveform on the left and a spectrogram on the right. Type 1 shows a regular waveform and parallel harmonic lines. Type 2 shows alternating big and small cycles and, after about 0.19 seconds, extra lines between the harmonics. Type 3 shows irregular cycles and a smeared spectrogram. Type 4 shows a noise-like waveform and a diffuse spectrogram.
Figure 3.1: Synthetic examples of the four signal types. Left: 60 ms of waveform. Right: narrowband spectrogram (50-ms Hamming window, 0–1.5 kHz) of the whole 0.6-s segment. Type 1: small random period and amplitude variation; harmonics are clean lines. Type 2: a bifurcation at 0.19 s (shaded) from period-1 to period-2 behaviour; new lines appear midway between the harmonics, at odd multiples of F0/2. Type 3: cycle lengths and amplitudes driven by a deterministic chaotic map; no stable harmonic lines. Type 4: weak irregular pulses buried in noise; energy is spread across frequency.

3.5 Applying the classification

The type belongs to an analysis interval, not to a speaker or a diagnosis. Both Type 1 and Type 2 are defined by what happens during the interval to be analysed (Baken and Orlikoff 2000, 152–53), so the same phonation can contain a Type 1 stretch and a Type 2 stretch. The practical question is always whether the segment you intend to measure is Type 1.

Ludlow, Kent and Gray describe the procedure. Select a segment of a prolonged vowel that is quasiperiodic, so that the periods of vibration can be detected accurately. Examine it with a narrowband spectrogram to see whether F0 and its harmonics are visible. If they cannot be identified, the signal is probably Type 3 or Type 4 and is not suitable for perturbation analysis (Ludlow et al. 2018, 378).

Figure 3.2 shows why the inspection cannot be delegated to the software. A simple autocorrelation F0 tracker without a voicing decision returns a value for every frame of every example. Only the Type 1 values describe the signal.

Show the code for this figure
fig, axes = plt.subplots(1, 4, figsize=(7.2, 2.6), sharey=True)
for ax, (name, x) in zip(axes, examples.items()):
    tt, f0 = vs.ac_f0_track(x, FS)
    ax.axhline(125, color=vs.MUTED, lw=0.8, ls="--")
    ax.plot(tt, f0, "o", color=vs.BLUE, ms=3)
    ax.set_title(name, color=vs.INK)
    ax.text(0.04, 0.03, f"mean {np.mean(f0):.0f} Hz\nSD {np.std(f0):.1f} Hz",
            transform=ax.transAxes, va="bottom", fontsize=8, color=vs.INK)
    ax.set_xlabel("Time (s)")
    vs.style_axes(ax)
axes[0].set_ylabel("Reported F0 (Hz)")
axes[0].set_ylim(0, 420)
fig.tight_layout()
plt.show()
Four small panels of F0 against time. Type 1 is a flat line at 125 Hz. Type 2 is flat at 125 Hz, then steps up to about 137 Hz. Type 3 scatters between about 95 and 150 Hz. Type 4 jumps between widely different values up to 400 Hz.
Figure 3.2: F0 returned every 10 ms by a naive autocorrelation tracker (40-ms frames, search range 60–400 Hz, no voicing decision) for the four synthetic signals of Figure 3.1. The dashed line is the 125-Hz rate of the pulse generator. Every frame gets a number. Type 1: the numbers describe the signal. Type 2: after the bifurcation the tracker settles on about 137 Hz, the rate of the shorter cycle of each pair; it describes neither the alternating pattern nor the 62.5-Hz pair rate. Types 3 and 4: the numbers are produced by the algorithm, not by any period in the signal.

Figure 3.3 summarises the decision.

Show the code for this figure
flowchart TD
    A["Select a quasiperiodic segment<br/>of the sustained vowel"] --> B["Inspect the waveform and a<br/>narrowband spectrogram"]
    B --> C{"F0 and harmonics<br/>visible?"}
    C -- "No" --> G["Type 3 or 4:<br/>no F0-based measures"]
    C -- "Yes" --> D{"Qualitative change, subharmonics,<br/>or strong modulation<br/>in the interval?"}
    D -- "No" --> E["Type 1:<br/>F0, perturbation and<br/>noise measures apply"]
    D -- "Yes" --> F{"Can a shorter Type 1<br/>segment be isolated?"}
    F -- "Yes" --> A
    F -- "No" --> H["Type 2:<br/>report perturbation as<br/>not applicable"]
    G --> I["Use measures that do not<br/>need F0, such as CPPS;<br/>describe the signal"]
    H --> I
flowchart TD
    A["Select a quasiperiodic segment<br/>of the sustained vowel"] --> B["Inspect the waveform and a<br/>narrowband spectrogram"]
    B --> C{"F0 and harmonics<br/>visible?"}
    C -- "No" --> G["Type 3 or 4:<br/>no F0-based measures"]
    C -- "Yes" --> D{"Qualitative change, subharmonics,<br/>or strong modulation<br/>in the interval?"}
    D -- "No" --> E["Type 1:<br/>F0, perturbation and<br/>noise measures apply"]
    D -- "Yes" --> F{"Can a shorter Type 1<br/>segment be isolated?"}
    F -- "Yes" --> A
    F -- "No" --> H["Type 2:<br/>report perturbation as<br/>not applicable"]
    G --> I["Use measures that do not<br/>need F0, such as CPPS;<br/>describe the signal"]
    H --> I
Figure 3.3: Signal-typing workflow before measurement. Branch labels follow the definitions of Baken and Orlikoff (2000, pp. 152–153) and the inspection procedure described by Ludlow, Kent and Gray (2018, p. 378).

3.6 Which measures apply to which types

Table 3.1: What each signal type supports.
Type F0 Perturbation (jitter, shimmer)
1 The definition of F0 is unambiguous only for Type 1 signals (Titze 1995, 26). Evaluation of F0 presupposes Type 1, subtype A, within reasonable limits (Baken and Orlikoff 2000, 153). Considerable utility and reliability (Titze 1995, 26).
2 No obvious single F0 throughout the segment (Titze 1995, 19). Per-cycle F0 and per-cycle amplitude, and their perturbation, cannot be reliably extracted. Perturbation measures by themselves are unreliable and contain little pattern information (Titze 1995, 26).
3 Not validly susceptible to F0 measurement (Titze 1995, 26; Baken and Orlikoff 2000, 153). Cannot be reliably extracted (Titze 1995, 26).
4 May or may not show an F0 on the spectrogram (Sprecher et al. 2010, 3711). Not reliable, and the signal is unsuitable for correlation dimension analysis as well (Sprecher et al. 2010, 3710).

Four points of detail behind the table.

  • How much perturbation is still reliable. Titze’s practical guideline for Type 1 signals is that perturbation measures below about 5% have been found reliable (Titze 1995, 26).
  • Where the contour breaks. Average cyclic parameter contours can be extracted from both Type 1 and Type 2 signals, but where a bifurcation occurs in a Type 2 signal, arbitrary decisions by the extraction algorithm are likely to affect the contour in a non-unique way (Titze 1995, 27). Event-based contours are highly susceptible to error in Types 2 and 3. They may be used for visual display, but are not recommended for perturbation measures on Type 3 signals (Titze 1995, 27).
  • An empirical cutoff. Sprecher and colleagues used the TF32 unreliability measures trk and err, taking an err above 10 to mark a sample ill-suited to perturbation analysis (Sprecher et al. 2010, 3711). Percent jitter and percent shimmer rose and SNR fell from Type 1 to Type 4, and both trk and err grouped Types 1 and 2 together against Types 3 and 4, with no significant difference between Types 3 and 4 (Sprecher et al. 2010, 3713).
  • What to use instead. For Type 2, visual displays are the most useful, and assessment is best based on the entire display rather than a single measure (Titze 1995, 26). For Type 3, perceptual ratings of roughness and of other auditory manifestations of aperiodicity are likely the best clinical measures, with fractal dimension, attractor dimension or the Lyapunov exponent as possible acoustic complements, and phase portraits as visual confirmation of high dimensionality (Titze 1995, 26). For Type 4, correlation dimension could not be calculated at all, the curves failing to converge, and spectrographic and perceptual analysis remain valid but yield no objective quantitative data (Sprecher et al. 2010, 3714; 2010, 3715).

Other measure families relate to the types as follows.

  • Harmonics-to-noise ratio. HNR needs a filter that separates F0 and its harmonics from the rest of the signal. Like jitter and shimmer, it can only be applied to a steady phonation such as a prolonged vowel (Ludlow et al. 2018, 378). It therefore shares their dependence on a Type 1 segment.
  • Cepstral peak prominence. CPP does not require direct computation of F0, and it can be extracted from connected speech as well as sustained vowels (Murton et al. 2020, 1596). Because no pitch detection or tracking is used, the cepstral measures are robust in dysphonic voices, particularly breathy ones (Barsties v. Latoszek et al. 2018, 316). The 2018 ASHA guidance recommended CPP as a general measure of dysphonia in place of jitter, shimmer and HNR (Murton et al. 2020, 1596). See the CPPS chapter.
  • Composite indices. AVQI includes shimmer local, shimmer local dB, and HNR among its components (Barsties and Maryn 2016, 574). ABI includes jitter local (Delgado Hernández et al. 2018, 320). These components carry the restrictions in Table 3.1 into the composite. See the AVQI and ABI chapters.
  • Dynamical-systems methods. For Type 3 signals, any covert organisation might be discoverable with the techniques of dynamical systems theory (Baken and Orlikoff 2000, 153).

3.7 Caveats

  • Typing is a judgement. The types are defined by inspection, and the boundary between Type 1 and Type 2 is stated as an energy ratio (an order of magnitude below the fundamental versus comparable to it) (Titze 1995, 19; Baken and Orlikoff 2000, 152–53). None of the ingested sources reports inter-rater agreement for signal typing itself. In Sprecher and colleagues the typing was done by one author with consultation from the others, and the Spearman correlation of 0.867 they report is for their three listeners’ perceptual ratings, which are a separate judgement (Sprecher et al. 2010, 3711; 2010, 3713). They also note that the cutoff for strong modulations is somewhat subjective, and excluded wavy-harmonic samples for that reason (Sprecher et al. 2010, 3711).
  • Perturbation can move the wrong way with treatment. When a voice is very irregular before treatment, few periods can be tracked, and jitter and shimmer may come out low. As the voice improves, more periods are identified and the measured values may rise (Ludlow et al. 2018, 378). Type the signal at every session before comparing perturbation values over time.
  • Occasional discontinuities. Occasional discontinuities may be less damaging to perturbation measurement than was once believed, at least for some software (Baken and Orlikoff 2000, 153). The requirement for a Type 1 segment still stands.
  • Subtypes are an extension. The A/B/C subtypes are Baken and Orlikoff’s addition to Titze’s three categories (Baken and Orlikoff 2000, 152). Cite a subtype only where the analysis depends on it.
  • The samples behind Type 4. Sprecher and colleagues’ 40 samples, ten per type, were selected as the most representative of each type from about 150 sustained phonations in a disordered-voice database (Sprecher et al. 2010, 3711). They were chosen to present the concept of a Type 4 signal, not as a consecutive clinical series, and the authors defer the association between pathology and type to future work (Sprecher et al. 2010, 3715).

3.8 Points of disagreement

  • Perturbation analysis of Type 2 signals. Titze recommends against it: for Type 2 signals, perturbation measures by themselves are unreliable and contain little pattern information (Titze 1995, 26). Sprecher and colleagues found that most of their Type 2 samples fell below the err cutoff of 10, and concluded that Type 1 and some Type 2 voices, those without bifurcations, could be analysed with perturbation (Sprecher et al. 2010, 3714). They state that this is inconsistent with Titze’s recommendations and may be attributable to their own typing procedure, and advise that future work consider pitch-tracking measures before including Type 2 voices (Sprecher et al. 2010, 3715). This chapter follows Titze for Type 2, and Table 3.1 reflects that.
  • Date of the classification. The workshop was held on 17 and 18 February 1994 (Titze 1995, 2), and the summary statement carries no publication year of its own. Baken and Orlikoff cite the scheme as Titze
    1. (Baken and Orlikoff 2000, 152); Ludlow, Kent and Gray cite it as Titze
    2. (Ludlow et al. 2018, 378); Sprecher and colleagues cite it as Titze (1995) (Sprecher et al. 2010, 3710). This reference uses 1995, the more common form, and gives the workshop date where it matters.
  • Character of Type 2 and Type 3 in the secondary accounts. Baken and Orlikoff define Type 2 by bifurcations or strong modulating or subharmonic components, and Type 3 by the absence of apparent periodicity, possibly with chaotic organisation (Baken and Orlikoff 2000, 153). Ludlow, Kent and Gray describe Type 2 as highly irregular periodicity with subharmonics and modulations such as tremor, and Type 3 as chaotic and usually perceived as rough (Ludlow et al. 2018, 378). Both are compatible with Titze’s wording (Titze 1995, 19), which this chapter now uses directly.
Baken, Ronald J., and Robert F. Orlikoff. 2000. Clinical Measurement of Speech and Voice. 2nd ed. Singular Thomson Learning.
Barsties, Ben, and Youri Maryn. 2016. “External Validation of the Acoustic Voice Quality Index Version 03.01 with Extended Representativity.” Annals of Otology, Rhinology & Laryngology 125 (7): 571–83. https://doi.org/10.1177/0003489416636131.
Barsties v. Latoszek, Ben, Youri Maryn, Ellen Gerrits, and Marc De Bodt. 2018. “A Meta-Analysis: Acoustic Measurement of Roughness and Breathiness.” Journal of Speech, Language, and Hearing Research 61 (2): 298–323. https://doi.org/10.1044/2017_JSLHR-S-16-0188.
Delgado Hernández, Jonathan, Nieves M. León Gómez, Alejandra Jiménez, Laura M. Izquierdo, and Ben Barsties v. Latoszek. 2018. “Validation of the Acoustic Voice Quality Index Version 03.01 and the Acoustic Breathiness Index in the Spanish Language.” Annals of Otology, Rhinology & Laryngology 127 (5): 317–26. https://doi.org/10.1177/0003489418761096.
Ludlow, Christy L., Raymond D. Kent, and Lincoln C. Gray. 2018. Measuring Voice, Speech, and Swallowing in the Clinic and Laboratory. Plural Publishing.
Murton, Olivia, Robert Hillman, and Daryush Mehta. 2020. “Cepstral Peak Prominence Values for Clinical Voice Evaluation.” American Journal of Speech-Language Pathology 29 (3): 1596–607. https://doi.org/10.1044/2020_AJSLP-20-00001.
Sprecher, Alicia, Aleksandra Olszewski, Jack J. Jiang, and Yu Zhang. 2010. “Updating Signal Typing in Voice: Addition of Type 4 Signals.” The Journal of the Acoustical Society of America 127 (6): 3710–16. https://doi.org/10.1121/1.3397477.
Titze, Ingo R. 1995. Workshop on Acoustic Voice Analysis: Summary Statement. National Center for Voice; Speech.