Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Context-Related Acoustic Variation in Male Fallow Deer (Dama dama) Groans

Abstract

While social and behavioural contexts are known to affect the acoustic structure of vocal signals in several mammal species, few studies have investigated context-related acoustic variation during inter-sexual advertisement and/or intra-sexual competition. Here we recorded male fallow deer groans during the breeding season and investigated how key acoustic parameters (fundamental frequency and formant frequencies) vary as a function of the social context in which they are produced. We found that in the presence of females, male fallow deer produced groans with higher mean fundamental frequency when vocal males were also present than they did when no vocal males were in close vicinity. We attribute this to the increased arousal state typically associated with this context. In addition, groan minimum formant frequency spacing was slightly, but significantly lower (indicating marginally more extended vocal tracts) when males were alone than when potential mates and/or competitors were nearby. This indicates that, contrary to our predictions, male fallow deer do not exaggerate the acoustic impression of their body size by further lowering their formant frequencies in the presence of potential mating partners and competitors. Furthermore, since the magnitude of the variation in groan minimum formant frequency spacing remains small compared to documented inter-individual differences, our findings are consistent with the hypothesis that formants are reliable static cues to body size during intra- and inter-sexual advertisement that do not concurrently encode dynamic motivation-related information.

Introduction

Recent studies of mammal communication have highlighted the multi-component nature of vocal signals (for review see [1]), and shown how the inter and intra-individual variation of the spectral components that compose these signals can simultaneously convey information on both static (long-term) and dynamic (short-term) attributes of callers. Vocalisations have been found to encode information on static attributes of callers such as their sex [2], [3], body size [3][9] or identity [10][15], as well as information on dynamic attributes of callers such as their emotional or motivational state [10], [16], [17]. Dynamic attributes are typically affected by a combination of internal factors such as the physiological state of the caller [18], [19] and external factors such as the social contexts of emission [20][28]. While several detailed investigations of the acoustic structure of mammal signals have identified the simultaneous presence of static and dynamic variation, there is an obvious potential conflict between the two types of information: components that encode information on static attributes should be characterised by limited variability, while components which encode information on dynamic attributes should be free to co-vary with the traits they express [29].

In the contexts of agonistic and sexual interactions, the “size” or “frequency” code hypothesis offers an evolutionary scenario whereby specific features of vocal signals known to encode honest information on static attributes of callers, such as their size, can simultaneously encode short-term, dynamic, motivational information [30]. This theory predicts that callers may evolve the ability to make ritualised use of relatively small variations of otherwise constrained parameters during social interactions. For example, in a species where spectral components correlate negatively with body size, callers may lower these components by a small amount in order to exaggerate the acoustic expression of their body size in order to sound more aggressive, or raise these components in order to minimise the acoustic expression of their body size and therefore sound more submissive [30]. However, the possible ritualised use of acoustic cues to static attributes has not been investigated in the context of nonhuman mammal vocal interactions.

The sexually-selected vocal displays of male fallow deer (Dama dama) constitute a useful model for investigating short-term variation of vocal behaviour. Early studies of vocal behaviour in fallow bucks focused on the function of calling rates: males who initiated vocal activity and remained vocal on most days during the rut achieved higher numbers of copulations [31] and average groaning rates were found to be related to mating success [32]. Investigations of short-term variation in calling rates found that males groaned at higher rates when females were present, and that males with females increased their groaning rates when vocal males were in close vicinity [33]. More recent studies of groaning in fallow deer have extensively described inter-individual variation in the acoustic structure of groans [12] as well as medium-term (during the rut) and long-term (lifespan) intra-individual variation [29], [34].

The acoustic structure of the fallow deer groan is best described in terms of its source- (fundamental frequency) and filter- (formant frequency) related components [6], [12], [29], [34], [35]. Male fallow deer groans are short, stereotypical vocalisations with an extremely low and sexually dimorphic fundamental frequency (F0) [11], [35], [36]. Although the basis for the inter-individual F0 variation in fallow deer groans remains unknown, in humans a lower F0 is associated with higher androgen levels [37], and correlational investigations have shown that higher-ranking, more dominant males produce groans with lower minimum F0 [6]. In addition, vocal tract resonances (hereafter termed formants) are reliable cues to body size in male fallow deer groans [6], as they are in the calls of several other male mammals [3][5], [7][9].

Here, we recorded free-ranging male fallow deer during the breeding season in order to determine whether variation in social context affects the acoustic structure of fallow deer ‘common groans’. More specifically, following McElligott & Hayden's [33] investigation of the effect of context on short-term rates of groaning, we investigated the effect of key proximate social factors on the acoustic structure of male fallow deer groans: the presence or absence of female(s) and vocal male(s). Along the “frequency code” theory [30], and previous findings in red deer (Cervus elaphus) [38], [39] we predict that fallow bucks should maximally extend their vocal tracts and therefore fully lower their formant frequencies when other vocal males are nearby and/or when females are in close proximity, as a means of maximising their apparent body size to competitors and potential mates. We have no strong a priori prediction for the affect of social context on groan F0. However, since lower minimum F0 is an established indicator of dominance in fallow bucks [6], we may expect mean F0 to be lower when males are in close proximity.

Results

Groan mean F0

The social context of recording had a significant effect on groan F0 (F3, 246 = 3.637, p = 0.013). The pair-wise comparisons showed that males with females produced groans with higher mean F0 when another vocal male was present (mean ± S.E. = 32.1±0.93 Hz) than they did when no other vocal male was close by (mean ± S.E. = 29.2±0.97 Hz; P = 0.032) or when they were alone (mean ± S.E. = 29.4±0.92 Hz; P = 0.030), but not significantly higher than when only a vocal male was present (mean ± S.E. = 30.2±1.11 Hz; P = 0.359) (see Figure 1A). No significant differences were detected in the mean F0 of male groans between the other social contexts (see Figure 1A).

thumbnail
Figure 1. Estimated marginal means ± SE for the effect of social context (alone, with a vocal male, with a female, with a female and a vocal male) on mean F0 (A) and minimum ΔF (B).

Mean responses sharing the same letter are not significantly different.

https://doi.org/10.1371/journal.pone.0021066.g001

Groan minimum ΔF

The social context also had a significant effect on groan minimum ΔF (F3, 238 = 8.640, p<0.001): groan minimum ΔF was lower (indicating more fully extended vocal tracts) when males were alone (mean ± S.E. = 300.8±2.05 Hz) than it was when vocal males (mean ± S.E. = 316.5±2.78 Hz; P<0.001), females (mean ± S.E. = 307.1±1.99 Hz; P = 0.041) or vocal males and females (mean ± S.E. = 308.6±1.89 Hz; P = 0.003) were in close attendance (see Figure 1B). In addition, males produced groans with higher minimum ΔF when they were with other vocal males (mean ± S.E. = 316.5±2.78 Hz) than they did in any other context (see Figure 1B).

Discussion

The results of this study indicate that the presence of mating partners and competitors has a significant effect on the vocal behaviour of male fallow deer, with frequency components of male groans varying between different social contexts. Specifically, we have shown that the mean F0 of male fallow deer groans increases when males are with females and other vocal males are in close proximity. In addition, fallow bucks that were alone produced groans with lower minimum ΔF than they did when vocal males and/or females were in close attendance, and higher minimum ΔF when vocal males were in close vicinity than they did in any other social context.

The higher mean F0 of groans produced by males with females when other vocal males were nearby could plausibly reflect greater arousal levels around potentially dangerous rivals. Indeed, vocal males are considered to be actively rutting [40] and, therefore, likely to constitute a threat. Moreover, previous work on this species showed that the rate of male groans increased when males were with females and other vocal males were close by [33], and higher rates of vocal activity are likely to signal increased motivation [41]. In addition, work on nonhuman mammals has shown that increased mean F0 is associated with increased arousal [10] and, from a vocal production perspective, high arousal levels could lead to increased sub-glottal pressure and dynamic strain on the vocal folds, leading to higher mean F0 [42]. Interestingly, the mean F0 of male fallow deer groans did not appear to change according to whether a female was present or not (see Figure 1), consistent with previous studies suggesting that dynamic variation in groaning behaviour is primarily directed at male competitors [33].

Contrary to our predictions, we found that bucks extended their vocal tracts less fully in the presence of potential mating partners and competitors than they did in the absence of both and thus, do not appear to make ritualistic use of size-related formant information to maximise the acoustic impression of their body size in these contexts. An explanation for this could be the very high groaning rates that male fallow deer are known to produce when there are females or vocal males nearby [circa 55 groans per minute: 33], which may make it more difficult for them to always lower the larynx to the maximum [unlike in the closely related red deer which only roars up to 8 times per minute: [43]]. However, the total magnitude of the observed variation (∼5%) remains small compared to the natural variation in the minimum formant spacing of adult male fallow deer groans of around 16% [minimum = 281 Hz, maximum = 326 Hz: 12] and, therefore, may not reflect a meaningful difference between a large, medium, and small male representative of the population i.e. males would probably still be categorised as small, medium or large by receivers.

It is also important to note that we cannot automatically assume that all the context-related differences in minimum ΔF and F0 that we report here are perceptible by fallow deer. Studies of human voice perception show that while listeners can discriminate differences in F0 as low as 2%, this rises to about 9% when F0 approaches 40 Hz [44], below which the mean F0 of most fallow deer groans falls [12], [34]. Similarly, human listeners cannot perceive shifts in formant spacing or apparent vocal tract length of less than 4% [44], [45]. While we cannot generalise human perceptual abilities to fallow deer, whether the variation reported above is perceptually relevant for this species remains to be demonstrated. We suggest that future studies investigate the perceptual abilities of male and female fallow deer, using habituation–discrimination tests and re-synthesised groan stimuli with different levels of F0 and formant spacing manipulation [46].

In conclusion, the results of the current study show that the mean F0 of male fallow deer groans increases when other vocal males are in close proximity, and do not support the hypothesis that males maximise the acoustic impression of their body size during close-range intra-sexual assessment and inter-sexual contexts. While the observed variation in F0 probably reflects variation in arousal, the higher minimum ΔF of groans produced by males when other vocal males are nearby may reflect a trade off between groaning rate and laryngeal lowering in these contexts. Nevertheless, the observed variation in minimum ΔF remains small compared to inter-individual variation reported for this species, and its perceptual and functional significance remain to be investigated. Furthermore, our findings are consistent with the hypothesis that formants are mostly static cues to body size during intra- and inter-sexual advertisement [6] that do not concurrently encode dynamic motivation-related information. Indeed, recent work on red deer has shown that formants are likely to function as static indices of body size in female mate choice contexts [39], [47] and when harem holders assess rivals [38]. Accordingly, we suggest that future studies on fallow deer use resynthesis techniques and playback experiments to further investigate the relative importance of size-related formant information during intra-sexual assessment and female mate choice contexts.

Materials and Methods

Ethical statement

This work follows the Association for the study of Animal Behaviour/Animal Behaviour Society guidelines for the use of animals in research, and was approved by the Ethics Committee of the Department of Psychology at the University of Sussex.

Study site and animals

A free-ranging population of fallow deer were observed during the 2004 breeding season (7th to 26th October) at Petworth Park, West Sussex, in the U.K. During the study period there were approximately 1000 deer in the park. Because the males in our study were not individually marked, unique antler shape and coat markings were used to identify individuals from multiple digital photographs taken at the time of each recording using a Nikon Coolpix 8700 (8 megapixel, 10× optical zoom). Only recordings from individuals that could be identified with absolute certainty were retained in the analysis.

Recordings and contextual settings

Groans were recorded from each of 18 male fallow deer at distances of 10–30 metres using a Marantz PMD 670 CF-card recorder and a Sennheiser MKH 416 P48 directional microphone (sampling rate of 44.1 kHz, 16 bits amplitude resolution). The sound files were transferred to an Apple Macintosh Macbook computer for visual inspection using narrow band spectrograms (see Figure 2: FFT method, window length = 0.03 s, time steps = 250, frequency steps = 1000, Gaussian window shape, dynamic range = 50 dB) and recordings with high levels of background noise were discarded. In total we used 269 groans from 11 males recorded on nine separate days. For each of the recording sessions we noted whether a vocally active male was within approximately 30 metres of the focal male. Vocal males are likely to constitute a greater threat because they are actively rutting [40]. We also noted whether any females were within 30 metres of the focal male (following [33]), and if the focal male was a harem holder.

thumbnail
Figure 2. Spectrogram (A) and waveform (B) of a male fallow deer groan.

Male fallow deer groans are low-pitched vocalisations characterised by a rapid downward shift in formant frequencies across the call (spectrogram settings: FFT method, window length = 0.05 s, time step = 0.01 s, Gaussian window shape, dynamic range = 30 dB).

https://doi.org/10.1371/journal.pone.0021066.g002

Acoustic analysis

We used Praat DSP package version 5.0.29 (www.praat.org) for the acoustic analysis. To measure F0 parameters we extracted the F0 contour of groans using the ‘To pitch (ac) command’ (time step = 0.01 s; minimum and maximum F0 = 20 Hz and 50 Hz, respectively). Time-varying numerical representations of the F0 contour were compared with the F0 contour as visualized on a spectrogram and checked to see if the autocorrelation algorithm was correctly tracking the F0. In cases where a harmonic or a sub-harmonic were tracked instead of the fundamental frequency (octave jumps), numerical representations of the F0 contour were manually adjusted using the ‘Edit’ window in Praat, and the resulting F0 contour was played back as a sine wave for comparison with the original recording.

As stated in the introduction, male fallow deer retract their larynx towards the sternum during vocalisation, resulting in a lowering of formants and formant spacing (ΔF) across the call [35]. Consequently, to quantify the minimum ΔF achieved when the larynx is fully retracted, we took the last pulses of each groan (∼0.05 s in duration) for our analysis. The frequency values of the first six formants were measured using Linear Predictive Coding (LPC; ‘To Formants (Burg)’ command in Praat) and the following analysis parameters: time step: 0.01 s; window analysis: 0.05 s; maximum formant value: 2300 Hz; maximum number of formants: 7; pre-emphasis: 6000 Hz. Since the first two formants are often poorly represented and vocal tract elongation has a minimal effect on them [35], the upper formant values (F3, F4, F5 & F6) were used to estimate ΔF during each groan using the linear regression method described by Reby and McComb [9].

Statistical analysis

We used linear mixed-effect models (LMM's) fitted with maximum likelihood estimation to examine relationships between acoustic features of male groans and the different contextual settings. Our acoustic measures (mean F0 and minimum ΔF) were entered as dependent variables and the social context was entered as a fixed factor categorical independent variable. Q-Q plots and scatter plots of residuals were inspected to ensure the dependant variables were normally distributed. The different contextual settings were as follows: no vocal adult male(s) or female(s) within 30 metres of the focal male; vocal adult male within 30 metres of the focal male; female(s) within 30 metres of the focal male; female(s) and vocal male(s) within 30 metres of the focal male. Pair-wise comparisons with Bonferroni adjustments were used to examine differences between the different social contexts. In addition, we entered subject identity as a random factor in each of our LMM's to control for repeated measures taken from the same individual. We also included the recording date and if the focal male was a harem holder as covariates, and used a model selection criteria based on the Akaike's Information Criteria (AIC), in which the model having the lowest AIC value is chosen [48]. Non-significant interaction terms were always removed from the model [49]. For the effect of social context on mean F0 the model with the lowest AIC value did not include the recording date or if the focal male was a harem holder as covariates. For groan minimum ΔF, the best model included the recording date, whether the focal male was a harem holder or not, and the interaction terms harem*recording date and social context*recording date. SPSS version 16 for Mac OS 10.5 was used to run the linear mixed effect models, significance levels were set at 0.05, and two-tailed probability values are quoted.

Acknowledgments

We thank Karen McComb and Paul Myles for help collecting the recordings, and the gamekeeper at Petworth Park, David Whitby, for giving us access to the animals. We would also like to thank Martha Tate for helping to confirm the identity of fallow bucks using photographs taken at the time of the recordings.

Author Contributions

Conceived and designed the experiments: DR BC. Performed the experiments: DR BC. Analyzed the data: DR BC. Contributed reagents/materials/analysis tools: DR BC. Wrote the paper: DR BC.

References

  1. 1. Taylor A, Reby D (2010) The contribution of source-filter theory to mammal vocal communication research. Journal of Zoology 280: 221–236.
  2. 2. Rendall D, Owren MJ, Weerts E, Hienz RD (2004) Sex differences in the acoustic structure of vowel-like grunt vocalizations in baboons and their perceptual discrimination by baboon listeners. Journal of the Acoustical Society of America 115: 411–421.
  3. 3. Charlton BD, Zhihe Z, Snyder RJ (2009) The information content of giant panda, Ailuropoda melanoleuca, bleats: acoustic cues to sex, age and size. Animal Behaviour 78: 893–898.
  4. 4. Fitch WT (1997) Vocal tract length and formant frequency dispersion correlate with body size in rhesus macaques. Journal of the Acoustical Society of America 102: 1213–1222.
  5. 5. Riede T, Fitch WT (1999) Vocal tract length and acoustics of vocalization in the domestic dog (Canis familiaris). Journal of Experimental Biology 202: 2859–2867.
  6. 6. Vannoni E, McElligott AG (2008) Low frequency groans indicate larger and more dominant fallow deer (Dama dama) males. PLoS ONE 3: e3113.
  7. 7. Sanvito S, Galimberti F, Miller EH (2007) Vocal signalling in male southern elephant seals is honest but imprecise. Animal Behaviour 73: 287–299.
  8. 8. Harris TR, Fitch WT, Goldstein LM, Fashing PJ (2006) Black and white colobus monkey (Colobus guereza) roars as a source of both honest and exaggerated information about body mass. Ethology 112: 911–920.
  9. 9. Reby D, McComb K (2003) Anatomical constraints generate honesty: acoustic cues to age and weight in the roars of red deer stags. Animal Behaviour 65: 519–530.
  10. 10. Rendall D (2003) Acoustic correlates of caller identity and affect intensity in the vowel-like grunt vocalizations of baboons. Journal of the Acoustical Society of America 113: 3390–3402.
  11. 11. Reby D, Joachim J, Lauga J, Lek S, Aulagnier S (1998) Individuality in the groans of fallow deer (Dama dama) bucks. Journal of Zoology 245: 79–84.
  12. 12. Vannoni E, McElligott AG (2007) Individual acoustic variation in fallow deer (Dama dama) common and harsh groans: a source-filter theory perspective. Ethology 113: 223–234.
  13. 13. McComb K, Reby D, Baker L, Moss C, Sayialel S (2003) Long-distance communication of acoustic cues to social identity in African elephants. Animal Behaviour 66: 317–329.
  14. 14. Blumstein DT, Armitage KB (1997) Alarm calling in yellow-bellied marmots: I. The meaning of situationally variable alarm calls. Animal Behaviour 53: 143–171.
  15. 15. Charlton BD, Zhihe Z, Snyder RJ (2009) Vocal cues to identity and relatedness in giant pandas (Ailuropoda melanoleuca). Journal of the Acoustical Society of America 126: 2721–2732.
  16. 16. Theis KR, Greene KM, Benson-Amram SR, Holekamp KE (2007) Sources of variation in the long-distance vocalizations of spotted hyenas1–29.
  17. 17. Manser MB, Seyfarth RM, Cheney DL (2002) Suricate alarm calls signal predator class and urgency. Trends in Cognitive Sciences 6: 55–57.
  18. 18. Charlton BD, Keating JL, Kersey D, Rengui L, Huang Y, et al. (2011) Vocal cues to male androgen levels in giant pandas. Biology Letters 7: 71–74.
  19. 19. Charlton BD, Keating JL, Huang Y, Rengui L, Swaisgood RR (2010) Female giant panda (Ailuropoda melanoleuca) chirps advertise the caller's fertile phase. Proceedings of the Royal Society of London, Series B: Biological Sciences 277: 1101–1106.
  20. 20. Weilgart L, Whitehead H (1990) Vocalizations of the North Atlantic pilot whale (Globicephala melas) as related to behavioral contexts. Behavioral Ecology and Sociobiology 26: 399–402.
  21. 21. Yin S, McCowan B (2004) Barking in domestic dogs: context specificity and individual identification. Animal Behaviour 68: 343–355.
  22. 22. Fischer J, Hammerschmidt K, Cheney DL, Seyfarth RM (2002) Acoustic features of male baboon loud calls: Influences of context, age, and individuality. Journal of the Acoustical Society of America 111: 1465–1474.
  23. 23. Rendall D, Seyfarth RM, Cheney DL, Owren MJ (1999) The meaning and function of grunt variants in baboons. Animal Behaviour 57: 583–592.
  24. 24. Notman H, Rendall D (2005) Contextual variation in chimpanzee pant hoots and its implications for referential communication. Animal Behaviour 70: 177–190.
  25. 25. Sugiura H (2007) Effects of proximity and behavioral context on acoustic variation in the coo calls of Japanese macaques. American Journal of Primatology 69: 1412–1424.
  26. 26. Taylor A, Reby D, McComb K (2009) Context-related variation in the vocal growling behaviour of the domestic dog (Canis familiaris). Ethology 115: 905–915.
  27. 27. Wich SA, Koski S, Vries H, Schaik CP (2003) Individual and contextual variation in Thomas langur male loud calls. Ethology.
  28. 28. Poole J, Payne K, Langbauer W, Moss C (1988) The social contexts of some very low frequency calls of African elephants. Behavioral Ecology and Sociobiology 22: 385–392.
  29. 29. Briefer E, Vannoni E, McElligott AG (2010) Quality prevails over identity in the sexually selected vocalisations of an ageing mammal. BMC Biology 8:
  30. 30. Ohala J, J. (1984) An ethological perspective on common cross-language utilization of F0 of voice. Phonetica 41: 1–16.
  31. 31. McElligott AG, O'Neill KP, Hayden TJ (1999) Cumulative long-term investment in vocalization and mating success of fallow bucks, Dama dama. Animal Behaviour 57: 1159–1167.
  32. 32. Clutton-Brock TH, Green D, Hiraiwahasegawa M, Albon SD (1988) Passing the buck - resource defense, lek breeding and mate choice in fallow deer. Behavioral Ecology and Sociobiology 23: 281–296.
  33. 33. McElligott AG, Hayden TJ (1999) Context-related vocalization rates of fallow bucks, Dama dama. Animal Behaviour 58: 1095–1104.
  34. 34. Vannoni E, McElligott AG (2009) Fallow bucks get hoarse: vocal fatigue as a possible signal to conspecifics. Animal Behaviour 78: 3–10.
  35. 35. McElligott AG, Birrer M, Vannoni E (2006) Retraction of the mobile descended larynx during groaning enables fallow bucks (Dama dama) to lower their formant frequencies. Journal of Zoology 270: 340–345.
  36. 36. Vannoni E, McElligott AG (2006) Individual Acoustic Variation in Fallow Deer (Dama dama) Common and Harsh Groans: A Source-Filter Theory Perspective. Ethology 113: 1–12.
  37. 37. Dabbs JM, Mallinger A (1999) High testosterone levels predict low voice pitch among men. Personality and individual differences 27: 801–804.
  38. 38. Reby D, McComb K, Cargnelutti B, Darwin C, Fitch WT, et al. (2005) Red deer stags use formants as assessment cues during intrasexual agonistic interactions. Proceedings of the Royal Society of London, Series B: Biological Sciences 272: 941–947.
  39. 39. Charlton BD, Reby D, McComb K (2007) Female red deer prefer the roars of larger males. Biology Letters 3: 382–385.
  40. 40. Moore N, Kelly P, Cahill J, Hayden T (1995) Mating strategies and mating success of fallow (Dama dama) bucks in a non-lekking population. Behavioral Ecology and Sociobiology 36: 91–100.
  41. 41. Berger J, Cunningham C (1991) Bellows, copulations, and sexual selection in bison (Bison bison). Behavioral Ecology 2: 1.
  42. 42. Titze IR (1994) Principles of Voice Production. Englewood Cliffs, New Jersey: Prentice Hall.
  43. 43. Clutton-Brock TH, Albon S, Guinness FE (1982) Red Deer: Behavior and Ecology of Two Sexes. Chicago: University of Chicago Press.
  44. 44. Smith DRR, Patterson RD, Turner R (2005) The processing and perception of size information in speech sounds. Journal of the Acoustical Society of America 117: 305–318.
  45. 45. Puts DA, Hodges CR, Cardenas RA, Gaulin SJC (2007) Men's voices as dominance signals: vocal fundamental and formant frequencies influence dominance attributions among men. Evolution and Human Behavior 28: 340–344.
  46. 46. Charlton BD, Reby D, McComb K (2007) Female perception of size-related formants in red deer, Cervus elaphus. Animal Behaviour 74: 707–714.
  47. 47. Charlton BD, Reby D, McComb K (2008) Effect of combined source (F0) and filter (formant) variation on red deer hind responses to male roars. Journal of the Acoustical Society of America 123: 2936–2943.
  48. 48. Burnham KD, Anderson DR (1998) Model selection and multimodal inference: a practical information-theoretic approach. New York: Springer-Verlag.
  49. 49. Engqvist L (2005) The mistreatment of covariate interaction terms in linear model analyses of behavioural and evolutionary ecology studies. Animal Behaviour 70: 967–971.