Next Article in Journal
Determination of Prolactin in Canine Saliva: Is It Possible to Use a Commercial ELISA Kit?
Previous Article in Journal
Motivations for Industry Stakeholders in China, Vietnam, Thailand and Malaysia to Improve Livestock Welfare
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Talking to Dogs: Companion Animal-Directed Speech in a Stress Test

1
Department of Cognitive Biology, University of Vienna, Althanstrasse 14, 1090 Vienna, Austria
2
Department of Behavioral Biology, Konrad Lorenz Research Station, University of Vienna, Althanstrasse 14, 1090 Vienna, Austria
3
Department of Special Education, Institut für Sonderpädagogische Entwicklungsförderung und Rehabilitation, University of Rostock, August-Bebel-Str. 28, 18055 Rostock, Germany
*
Authors to whom correspondence should be addressed.
Animals 2019, 9(7), 417; https://doi.org/10.3390/ani9070417
Submission received: 3 June 2019 / Revised: 25 June 2019 / Accepted: 27 June 2019 / Published: 4 July 2019
(This article belongs to the Section Companion Animals)

Abstract

:

Simple Summary

Companion animal-directed speech is a current topic of research, interesting due to its similarity to infant-directed speech. Dog owners seem to almost subconsciously use this high-pitched and repetitive way of speaking, slightly adapted, for dogs. The aim of this study was to investigate dog-directed speech in different contexts and examine whether owner personality and relationship quality affect it. We found that owners’ personality and gender affect their dog-directed speech. The majority of the modifications of dog-directed speech could be explained by a differential use of voice pitch and range. Our study supports the idea that voice pitch was used to communicate affect, whereas pitch range was used as an attention-getting strategy. Based on our results, we conclude that dog-directed speech is adjusted depending on context, gender, and personality. Societal value in this study consists of its contribution to basic knowledge of how we talk to animals, which may help in preventing accidents (e.g., dog bites) as well as improving animal training.

Abstract

Companion animal-directed speech (CADS) has previously been investigated in comparison to infant-directed speech and adult-directed speech. To investigate the influence of owner caregiving, attachment pattern, and personality on CADS, we used the Ainsworth strange situation procedure. It allowed us to assess voice source parameters of CADS across different contexts. We extracted speech parameters (voicing duration, voice pitch, pitch range, and jitter) from 53 dog owners recorded during the procedure. We found that owner personality and gender but not caregiving/attachment behavior affect their voice’s pitch, range, and jitter during CADS. Further, we found a differential and context-specific modification of pitch and range, consistent with the idea that pitch communicates affect, whereas range is more of an attention-getting device. This differential usage, and the increased pitch, emphasize and support the parallels described between CADS and infant-directed speech. For the first time, we also show the effect of personality on CADS and lay the basis for including jitter as a potentially useful measure in CADS.

1. Introduction

It has long been known that the human voice carries substantial cues to the emotional state of the owner, completely independent of any verbal or linguistic content [1]. With the exception of a few familiar words (e.g., their name, “good”, “bad”, etc.), animals listening to human speech presumably rely primarily upon such nonverbal information to interpret speaker state and intentions [2]. However, the specific forms of information available in companion animal-directed speech (CADS) remain understudied. In this study, we analyzed dog-directed speech in different situations to gain insight into this question. We focused attention on source characteristics, particularly voice fundamental frequency (f0, often colloquially termed “voice pitch”), which has reliably been found to vary with emotional state, though often in individual-specific manners [3,4]. Based on this previous research, we hypothesized that voice pitch variables would increase in increasingly arousing situations. In addition to summary statistics on f0 itself (mean, range, duration), we also analyzed voice perturbation using jitter, which quantifies local period-to-period deviation in the length of successive f0 periods. Measures of voice pitch irregularities (jitter) have a long history in voice emotion analysis, and high jitter has been suggested to correlate with high stress or arousal [3,5]. However, inconsistencies across studies have been argued to reflect individual personality differences [1] and may vary with coping style (e.g., for inhibited vs. outgoing children [6]). We therefore had only a weak a priori prediction that jitter should increase with arousal, perhaps in individual- or personality-specific ways.
Research trying to identify information encoded in CADS has used infant-directed speech (IDS) as a comparison. When compared to ordinary speech, CADS has a higher fundamental frequency, higher pitch range, and is more repetitive; these characteristics have also been found in IDS [7,8,9]. CADS may be similar to IDS in quality but differs from IDS in an apparent lack of hyper-articulation and lack of a language learning context [7,10,11]. These similarities and subtle differences make the comparison all the more intriguing and fascinating. Not only do people use CADS with dogs, but dogs also seem to pay attention to it. CADS has been evaluated for its attention-getting quality with adult dogs and seems to draw adult dogs’ attention more than adult-directed speech [12]. It appears that we adapt CADS (e.g., omitting hyper articulation), knowing that our four-legged friends will never be able to speak, but we still use it in navigating our relationship with them. One of the most important functions within relationships is giving and receiving emotional social support. How we give and receive this emotional and social support depends on and varies with arousal level and attachment representations [13,14]. Humans biophilic nature combined with dogs’ adaptation to humans (domestication) allows both sides to use and respond to the aforementioned human offspring–caregiver attachment behavior [15,16].
Any attachment relationship is characterized by maintaining proximity, separation distress, and using the caregiver as a secure base to explore from, and as a safe haven to return to [17,18]. Mary Ainsworth developed a classification tool to diagnose the quality of an infant’s attachment [19,20], the Ainsworth strange situation procedure (ASSP). The ASSP is designed to assess attachment patterns by increasing the subjects’ mental stress load and arousal over consecutive episodes [21,22]. It is created to increasingly activate attachment behavior on the side of the child or dog; this behavior is usually reciprocated by the caregiver in a complementary manner [23], which means the experience of mental stress has the subject turn towards support and comfort, which is found in the form of social support [19]. This is why we want to investigate CADS in this setup: Social communication has a highly relevant function in providing emotional and social support in any relationship. Both attachment theory and the ASSP have previously been applied to human–dog dyads, with promising results [24,25,26,27]. Since literature suggests that CADS can be influenced by owner gender and context [28,29], the ASSP is an excellent tool for testing CADS in varying contexts. To date, analyses of CADS did not take into account important factors in social communication, such as attachment pattern, caregiving behavior or personality. Here, we want to investigate their potential impact on the voice source during CADS.

2. Materials and Methods

2.1. Subjects

A total of 59 human–dog teams participated in this study. Six out of the 59 teams had to be excluded from the analysis due to recording equipment malfunctions. The remaining 53 human–dog teams (28 women, 25 men, 28 female dogs, 25 male dogs) in balanced combinations of male–male, female–female, male–female, and female–male were a subset of 132 human–dog dyads who had already participated in an experimental study of interaction styles and human–dog relationships [30,31]. The recruitment for our study was based on voluntary participation of teams from the pool of 132 dyads. Our subset of dyads represents those who responded positively to our request for participation in this study. All dogs lived with the owner who was also the main attachment figure from puppyhood onward. All dogs were intact, i.e., were neither spayed nor neutered. Their mean age was 4 years ± 1.5 SD and their mean weight was 30 kg ± 13.2 SD. Mean owner age was 46.2 years ± 10.2 SD. All dyads were recruited from Vienna and surrounding areas in Austria.

2.2. Owner Personality Axis

Dog owners were asked to fill in the German version of the NEO Five-Factory Inventory evaluating their own personality [32,33], a 60-item psychometric instrument designed to evaluate nonclinical adult personality structures along five major dimensions: neuroticism, extraversion, openness, agreeableness, and conscientiousness.

2.3. Ainsworth Strange Situation Procedure

The Ainsworth strange situation procedure [20] was adapted to assess the human–dog relationship. Within this assessment, dog attachment behavior is activated via increasing stress during the procedure, which in turn potentially activates caregiving behavior in the owner. Before testing, the room was prepared by placing two color-tagged chairs next to each other in the middle of the room, closing the windows and shades, and depositing several toys between and in front of the chairs. The dog owners were guided through the procedure by the experimenter in the following fixed order:
0. Controls (~10 s): Two control settings (reading and speaking/adult-directed speech) are recorded in the waiting room, adjacent the experimental room: First, the owner is asked to read a predefined text in presence of the dog out loud. (“…der beste Freund des Schäferhundes Rex ist die Ente Oskar…”—“… the German shepherd Rex’s best friend is the duck Oscar…”). For the second control, the experimenter engages the owner in small talk about their dog.
1. Introduction (~20 s): The experimenter leads the owner–dog dyad into the room, introduces them to the surroundings, asks the owner to unleash the dog for the duration of the procedure, and leaves the room again.
2. Exploration (3 min): Both dog and owner can move around freely and explore the room and the toys. The owner is free to interact and talk with the dog normally.
3. Encounter (1 + 2 min): A person previously unknown to the dog (‘the stranger’) enters the room and asks the owner to take their designated seat if they were not sitting yet. Other than this, there is no further interaction initially: For one minute, the stranger sits down on the other chair, remaining motionless. In the second minute, the stranger initiates small talk with the dog owner. In the third and last minute of this first encounter, the stranger tries to interact with and elicit play from the dog. At the end of minute three, the stranger asks the dog owner to leave the room.
4. Separation (3 min): The dog stays alone in the experimental room with the stranger. The stranger tries to engage the dog by playing and interacting.
5. Call (~5 s): The experimenter instructs the dog owner to stand in front of the experimental room door and to loudly call the dog’s name.
6. Reunion (3 min): The owner re-enters the room and reunites with the dog. The owner is free to interact and to play with the dog. The stranger quietly leaves the room. At the end of minute three, a vibrating phone indicates that the owner should leave the experimental room.
7. Separation (3 min): The dog remains alone in the experimental room for three minutes.
8. Encounter (3 min): The stranger from stage 3 re-enters and the dog remains in the experimental room with the stranger. The stranger tries to comfort the dog by playing and interacting.
9. Call (~5 s): The experimenter again instructs the dog owner to stand in front of the experimental room door and loudly call the dog’s name.
10. Reunion (3 min): The owner enters the room, reunites with the dog, and provides at least one physical contact. The owner is free to interact and play with the dog. The stranger quietly leaves the room. At the end of minute three, the experimenter enters the room and ends the procedure.
The procedure is illustrated graphically in Figure 1, marking all episodes (0, 2, 3, 5, 6, 9, 10) where the owner’s utterances were analyzed with the owner icon.

2.4. Video Recording, Surveillance and Animal Welfare

The test room was equipped with a wide-angle lens video surveillance system (Canon Inc., Canon Austria GmbH, 1100 Vienna, Austria) to record and monitor the proceedings throughout from outside the room. The experimenter monitored the situation directly outside of the test room to ensure the dog and dyad’s safety and to intervene if needed. The dog owners could terminate the experiment at any point without giving any reason and were permitted to monitor the dog after exiting the room from the video monitor outside, with the experimenter. In cases where any of the dogs experienced extensive stress, the experiment would have been stopped by the experimenter independently of the dog owner; however, this was never the case. As judged by behavioral parameters, none of the dogs experienced stress levels outside of ordinary levels, so in no case did the owner or experimenter intervene to terminate the test situation. All participants were asked to sign consent forms and were informed about their voluntary participation and their right to stop the experiment at any time without providing a reason (see Section 2.11).

2.5. Dog Attachment Classification

The dog’s attachment classification, based on the adapted ASSP, was analyzed according to methods detailed in Schöberl et al. and Solomon et al. [30,34]. Based on the video recordings of the ASSP, coauthors AB and JS, two psychologists trained in attachment categorization of human toddlers, together classified the dogs (with 89% interrater reliability) into five categories: secure, insecure avoidant, insecure ambivalent, insecure disorganized, and unclassifiable. KK and IS assisted with their expertise in dog behavior. Sample sizes and criteria were as follows:
A total of 27 out of 53 dogs were classified as ‘securely attached’: They eagerly approached their owners during the reunion and actively searched for and tolerated physical contact, while also showing interest in exploring the room. Three dogs were classified as ‘insecure-avoidant’, showing little tendency to approach their owner during the reunion but spending much time exploring the room and the toys. Six dogs were classified as ‘insecure-ambivalent’ based on their tendency to not explore their environment after the reunion at all but instead remaining in close proximity to their owner. Nine dogs were categorized as ‘insecure-disorganized’ who showed odd behavioral elements not being part of the normal attachment repertoire, such as freezing, staring or evident stereotypies. For eight dogs, no consensus could be reached; they were therefore labelled ‘unclassifiable’.

2.6. Owner Caregiving Rating

The owner caregiving rating was developed by Solomon et al. [34] in the context of previous studies [30,31] and is based on Ainsworth’s maternal sensitivity scale [20] and the ‘supportive presence’ scale [35]. The scale was designed to capture the caregiver’s responsiveness and sensitivity to the dog’s needs in a threatening situation [34]. During the threat task, an unfamiliar person entered the test room wearing a black coat, hat, and ski mask (with only the eyes visible). The unfamiliar person took three steps (in an interval of three seconds) towards the tethered dog while staring at its face. This process was repeated twice. After the second encounter, the unfamiliar person de-escalated the threat by stepping back, taking off the disguise, talking to the dog in a calming matter and offering cheese. The dog owner was present in only one encounter with the unfamiliar threat. The order of the owner’s presence or absence during the encounters was randomized for each dyad. The experimenter (and in one test scenario, also the owner) observed both conditions on monitors outside the room. The dyad participated in the mild threat task four weeks to a year prior to entering the ASS procedure. Caregiving behavior by owners was rated based on the threat task on a seven-point caregiving scale for all dyads. The highest score of seven was given if the owner showed a consistent, quick, and flexible response in their caregiving behavior. The minimum score of one was given if the dog owner did not respond at all to the dog, or if they responded in a negative or punishing way. The ethics of the caregiving scales development was reviewed by the ‘Faculty of Life Sciences’ at the University of Vienna (case number: 2014-015).

2.7. Audio Recording

The dog owners’ vocalizations were recorded during the entire ASSP with an H4N recorder connected to a small Sennheiser (ew 100 G3) microphone attached to the clothing in the chest area. The sampling frequency was 48 kHz with 16-bit quantization, and the sensitivity was adjusted prior to recording to prevent clipping. The recording was stopped by the experimenter shortly after entering the room after episode ten.

2.8. Audio Treatment and Analysis

All 53 recordings had an adequate signal-to-noise ratio for audio analysis. The audio files of all dyads were prepared for semi-automated, acoustic analysis by hand editing with Audacity (Version 2.1.1; www.audacityteam.org), removing ASSP episodes and episode parts where either a talking second person was in the test room with the owner, or where the dog owner quietly waited outside of the test room (Figure 1). Following this procedure, episodes 1, 4, 7, and 8 as well as minutes two to three of episode 3 were excluded with Audacity (1: experimenter talking in room; 3: strange person talking in room; 4, 7, and 8: owner waiting outside). The remaining episodes 0, 2, 3, 5, 6, 9, and 10 were run through a semi-automated analysis pipeline in Praat (Version 6.0.23; www.fon.hum.uva.nl/praat/). All analysis was conducted using five custom written Praat scripts.
I. A: Normalization and prefiltering: Each audio file was adjusted to a maximum amplitude peak of 0.99 and treated with a spectral subtraction background noise filter using Praat commands ‘Scale peak’ and ‘Remove noise’.
I. B: Speech segmentation: Praat’s ‘To TextGrid (silences)’ was used to label those portions of the audio containing the owners’ speech. The thresholds for the speech stream segmentation were based on silence in between intervals with a threshold of −35 dB silence, a minimum duration of 0.5 s for silent intervals, and a minimum duration of 0.07 s for spoken intervals. This process labelled the relevant spoken intervals automatically by marking them in the Praat TextGrid. The analyst then verified these automatic labels and made manual accuracy adjustments if necessary, in case of bursts of sound created by the dog vocalizing or playing with the toys.
II: Pitch contour tracking: All spoken intervals were saved, and their pitch contours were tracked and extracted using the Praat command ‘Extract visible pitch contour’. The Praat internal commands ‘Get mean’, ‘Get minimum’, and ‘Get maximum’ were used to measure the acoustic parameters mean, minimum, and maximum of each individual pitch contour.
III: Episode matching: A TextGrid tier with episode identifiers was added to match each spoken interval with the corresponding episode. The start and end time of all spoken intervals were extracted.
IV: Jitter measurements: The command ‘To PointProcess’ was used to create a point process object out of the pitch contour extracted in script II. Those ‘PointProcesses’ and the command ‘Get jitter (local)’ were used to measure voice’s jitter. The default settings were used except for the ‘Longest period (s)’, where 0.033 (minimum frequency measured) was input.

2.9. Analyzed Parameters

Four summary parameters of the f0 (‘voice pitch’) track were used to describe the owner’s speech in spoken intervals throughout the test procedure (Figure 2): voiced duration, mean, range, and jitter. The semiautomated system measured each speech parameter below for every spoken interval. Intervals labelled as spoken but without measurable f0 were automatically excluded from further analysis. For readability, the fundamental frequency (f0) is referred to as ‘pitch’ throughout the rest of the paper.
Duration: Voiced utterance duration measured in seconds for each interval and calculated by subtracting the start time from the end time.
Mean Pitch: Mean f0 of the spoken interval in hertz.
Pitch Range: The minimum f0 measured within the interval, subtracted from the f0 maximum.
Jitter: Jitter was measured using Praat’s algorithm ‘To Jitter (local)’. This parameter is used as a measure of voice quality and is ‘the average absolute difference between consecutive periods, divided by the average period’ (Version 6.0.23; http://www.fon.hum.uva.nl/praat/manual/PointProcess__Get_jitter__local____.html).

2.10. Statistical Analysis

Statistical analyses were done using R (Version 3.3.3; www.r-project.org/) and R-Studio (Version 1.0.13; www.rstudio.com/) and the packages ggthemes, ggplot2, psych, data.table, visreg, dplyr, doBy, xlsx, usdm, tidyverse, G.Gally, stats, lme4, car, nlme, lme4, cowplot, multcomp, MuMIn, piecewiseSEM, sjstats, and MASS. The data set for this study is available in the Supplementary Materials (Spreadsheet S1). Of all 53 dyads included in this study, over 4100 measurements were analyzed for each of the four response variables. Prior to running detailed analysis, all fixed factors tested negatively for multicollinearity. After visual inspection of the residuals of each response variable, the basic assumption of linear mixed models that the residuals follow a normal distribution could not be confirmed for jitter, mean pitch, pitch range, and voiced duration. Therefore, generalized linear mixed models (glmm) fit using maximum likelihood were calculated for all these four variables. The best distribution fit for each individual model’s response variable was established by investigating the residuals’ histogram and plotting the full models fitted residuals over the estimated residuals. Since the four response variables jitter, mean pitch, pitch range, and voiced duration were continuous and positively skewed towards the right, a Gamma distribution with a log link was used. Dyad was included as random effect to control for individual variances in each model, and the centered and scaled NEO FFI character traits (agreeableness, conscientiousness, extraversion, openness, and neuroticism), centered and scaled caregiving rating, attachment pattern, and owner gender were added as fixed factors and episode as a covariate of the full model. Because voice pitch varies by roughly an octave between adult men and women, gender was highly important to include as a fixed factor for each model. All attachment categories, except for the secure group, were collapsed into the nonsecure attachment group. This was a necessary simplification to reduce the number of levels within categories and facilitate model convergence. Null models included the covariate episode and the random effect dyad. All full models, except the duration model, tested to be a significantly better fit than the null model. Based on the results of the full-null model testing, duration was thus excluded from further analysis. Up to this point, standard practice of fitting linear mixed models was followed. Further decisions for the statistical analysis require some explanation.
A widely-used traditional but criticized approach for finding the best fitting model is a method called stepwise model reduction [36]. This process excludes one variable at a time from the full model based on p-values or AICc scores until reaching the null model. Apart from the problems arising from the usage of p-values, another issue arises in this context [37]: This stepwise reduction prevents an overview of all possible combinations of fixed factors explaining the response variable. To circumvent the problems and limitations arising through stepwise reduction, we used the function ‘dredge’ of the MuMIn package. This function creates every single model possible out of the fixed factors of the full model with (in our case) AICc scores. The model with the lowest AICc score is picked to provide a baseline, and within a delta of 2 upwards, all models are considered mathematically equally good fits [38]. The models within this range are put through model averaging to create averaged coefficients. This process allows evaluation of the influence all fixed factors have on the response variable without the restriction of p-values. For compatibility with the described approach of model averaging based on AICc values, a confidence interval of 85% is used to judge each factor’s impact on the speech parameters [39]. Factors with a confidence interval not including 0 were chosen to be of importance in predicting voice parameters. Relative importance, a measure describing each factor’s relative importance compared to the most valuable factor within the averaged coefficients, is used as a second measure and confirmation of the confidence intervals.
This approach of comparing all models within the delta 2 of the lowest AICc through model averaging was used to gain insight into the complex framework of CADS and its interactions with human personality and the human–dog attachment/caregiving system, without the unnecessary restriction of eliminating valuable comparisons from step one. For a review on null hypotheses significance testing and information theory based approaches and their possible combinations as used here, see Mundry 2011 [40].

2.11. Ethics

All participants were asked to sign consent forms, were informed about the procedure, and could terminate participation at any time. The ethics regarding human participation in the ASS procedure was reviewed and approved by the German Society for Psychology (Deutsche Gesellschaft für Psycholgie, AB 07_2011). All human/animal data collection and analysis was done in accordance with the declaration of Helsinki and the EU Directive 2010/63/EU for animal experiments. The ethics for this study was reviewed and approved by the ‘Faculty of Life Sciences’ animal welfare committee at the University of Vienna (case number: 2014-015).

3. Results

3.1. Mean Pitch

We found owner gender to affect pitch across episodes (Table 1). CADS was consistently higher in voice pitch than read speech (male mean: 121 ± 20 Hz SD; female mean: 203 ± 33 Hz SD) or conversational speech (male mean: 119 ± 19 Hz SD; female mean: 204 ± 38 Hz SD; Figure 3). Both male and female owners’ highest median pitch was recorded during the call episodes (male mean: 176 ± 32 Hz SD; female mean: 300 ± 57 Hz SD). In both genders, the reunions, exploration, and encounter had a similar mean pitch. No effect of personality traits on CADS mean pitch was observed for male and female owners.

3.2. Pitch Range

We found gender and openness to affect pitch range across episodes (Table 2). Pitch range in CADS and conversation was reduced relative to reading. Both male and female owners showed the highest median pitch range in the reading control condition (Figure 4A). The median pitch range in the reunions, exploration, and encounter episodes in both men and women was lower than the speaking (male mean: 73 ± 48 Hz SD; female mean: 124 ± 84 Hz SD; Figure 4A) and reading controls (male mean: 84 ± 42 Hz SD; female mean: 146 ± 54 Hz SD; Figure 4A). Female but not male owners high in openness showed an increased frequency range (Figure 4B).

3.3. Jitter

Jitter results in the CADS condition were highly variable. We found owner gender and openness to influence jitter across episodes (Table 3). The median jitter was lowest in the encounter (male mean: 0.014 ± 0. 009 % SD; female mean: 0.01 ± 0.007 % SD) and the call episodes (male mean: 0.009 ± 0.004 % SD; female mean: 0.007 ± 0.003 % SD) in both male and female owners. The speaking control condition had the highest median percentage of jitter in both men (male mean: 0.018 ± 0.008 % SD) and women (female mean: 0.012 ± 0.005 % SD). Male owners’ reunion and exploration episodes had a lower percentage of voice jitter than their speaking control condition; women’s percentage of jitter in the same episodes had similar levels to the speaking control (Figure 5A). Male owners’ openness values scaled positively with the percentage of jitter in their vocalizations; this effect was comparatively weaker in female owners (Figure 5B).

4. Discussion

We analyzed vocal parameters in dog-directed speech during a series of controlled encounters in dog–human dyads and compared these to normal speech (a read passage, or between-human conversation). The staged encounters were designed to elicit arousal in the dogs and caregiving from their owners. In general, we found that CADS was higher in f0 (“voice pitch”) than normal speech but showed a narrower pitch range. Results concerning pitch perturbation (jitter) were quite variable and showed no clear effect of arousal; the most pronounced effect was a considerable decrease in jitter associated with the owner calling the dog’s name in the two “call” episodes.
Compared to earlier work on CADS, our results support the findings of an increased pitch but fail to support the previously described broader pitch range. Jeannin et al. [29] partially used an approach similar to the ASSP and found voice pitch in dog-directed speech to increase in reunion episodes compared to adult-directed speech. Our results showed the exact same pattern, with a strong increase in voice pitch during the reunion episodes. Gergely et al. [28] reported a broader frequency range in CADS in mothers compared to fathers. We found the same effect in women to using a broader frequency range during CADS than men did. A well-described phenomenon in IDS and CADS is a differential use of voice pitch and frequency range. Pitch is said to communicate affect, whereas the range is used more in attention getting [10,28,41,42,43]. This different function of pitch and range might also be the explanation as to why we did not find a broader pitch range. Pitch range might simply be adjusted strongly in accordance with context. Namely, if the owner wants to draw the dogs’ attention away from something, a broad pitch range could be used. The opposite would be true if the owner already has the dog’s attention and is trying to soothe and calm it [41]. This would explain why the owners used a broader pitch range and a high pitch while calling their dogs and therefore hoping to draw their attention while also communicating positive intentions. The pitch in exploration, encounter, and the reunion episodes was lower in comparison to the call episodes but still increased compared to the read and adult-directed speech. This elevated pitch coincided with the narrowest pitch ranges in the same episodes. Exploration, encounter, and both reunions share in common that most dogs had already directed their attention towards their owner (in a 2013 study, 85% of the dogs followed their owners around in the reunion and encounter episodes [44]); no broad pitch range was needed to keep their focus. The owners continued to use an elevated pitch to communicate positive affect and a nonthreatening situation.
The results regarding jitter were less clear and more variable. Both males’ and females’ adult-directed speech had the highest percentage of jitter compared to the lowest in the call episodes. Female owners’ median jitter slowly returns to the levels of the reading control throughout the exploration, encounter, and reunions. Male owners’ jitter stays at an elevated level during the exploration and reunion episodes. Interestingly, jitter was highest in the adult-directed speech and behaved contrary to what we would have expected. An explanation for this might be found in the hypothesis that an increase in the speakers’ stress leads to a decreased jitter due to higher tension on the vocal folds [1,4]. This might be a parabola-like phenomenon with relaxed and extremely stressed speakers producing the highest jitter values. This hypothesis would explain why we would see jitter to be lowest when the speaker is stressed but not overly so; enough to cause tension (in the body and the vocal folds) yet mild enough to be dealt with. This hypothesis fits best with our results, but further research must be done to empirically assess this claim.
Our work also illustrates the importance of owner gender and owner personality on CADS. Not only does our study support the idea that CADS is gender-specific and variable enough among dyads to constitute an important component of human–dog caregiving and attachment strategy, we also found openness to drive CADS modulation within this ASSP setup. The influence of personality on CADS seems to be gender-specific and may be considered part of gender-specific performance. Men’s openness was positively correlated with utterance jitter, while women showed almost no correlation. Higher scores in openness correlated with a higher pitch range in women but showed no correlation in men. The limited literature regarding personality in the context of CADS restricts us in doing more than speculating as to why openness might influence jitter and pitch range. There might be one intriguing hypothesis explaining the influence of this one personality axis. The ASSP is designed to evoke arousal and a stress response (and therefore cortisol excretion) in the dogs, which in turn causes a similar stress response in owners. Research suggests a link between cortisol stress response, personality, and gender [45,46]. We propose this link to be reflected in the vocal parameters of jitter and pitch range. This possible interaction of personality, cortisol, and CADS might be a potentially interesting field of further research.

5. Conclusions

To summarize, our data partially support and partially contradict our initial hypothesis of an increase in voice source variables with increasingly aroused situations. We did find voice pitch to increase in the ASSP setup in comparison to normal speech, but we found the opposite to be true for pitch range. The decrease in pitch range might be caused by the differential use and function of mean pitch and range. With the narrower range, the owners tried to calm the dogs in this stressful setting. Due to limited reports in literature and close to no precedent, our predictions and hypothesis on jitter were less clear. The idea of voice perturbation to increase with arousal was not supported by our data, but it fit the hypothesis of being (at least partially) personality-dependent.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2615/9/7/417/s1, Spreadsheet S1: data set.

Author Contributions

Conceptualization, R.L.; Data curation, I.S.; Formal analysis, R.L.; Funding acquisition, K.K. and W.T.F.; Investigation, I.S.; Methodology, R.L., K.K., I.S., A.B., J.S., and W.T.F.; Project administration, I.S.; Resources, K.K. and W.T.F.; Software, R.L.; Supervision, K.K. and W.T.F.; Validation, I.S., A.B., and J.S.; Visualization, R.L.; Writing—original draft, R.L.; Writing—review and editing, K.K. and W.T.F.

Funding

The study was funded by the Austrian Science Fund (FWF), Project P23345, by Waltham grant (FA 566001) (both to K.K.), by the PhD program Cognition and Communication 1 and 2 from the Austrian Science Fund (FWF): Project W1234-G17 (to K.K./I.S.), Project W1262-B29 (to R.L./W.T.F.) and by a J. William Fulbright award (to J.S.). The funding agencies were not involved in the collection, analysis, interpretation of data, writing the report, and decision to submit the article.

Acknowledgments

We would like to thank all owners and their dogs for their participation in our study. We thank Manuela Wedl for her important impact during the developmental stage of this project and all students, who were involved: Evi Myska, Sarah Prettner, Philipp Stöger, Sigrid Amon, and Verena Ziemen. We thank Dan C. Mann for support with Praat and Matthias Loretto for input on statistical analysis. Further, we thank Arik Kershenbaum and one anonymous reviewer for comments and suggestions of improvement on a previous version of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The sponsors had no role in the design, execution, interpretation, or writing of the study.

References

  1. Protopapas, A.; Lieberman, P. Fundamental frequency of phonation and perceived emotional stress. J. Acoust. Soc. Am. 1997, 101, 2267–2277. [Google Scholar] [CrossRef] [PubMed]
  2. Siniscalchi, M.; D’Ingeo, S.; Fornelli, S.; Quaranta, A. Lateralized behavior and cardiac activity of dogs in response to human emotional vocalizations. Sci. Rep. 2018, 8, 77. [Google Scholar] [CrossRef] [PubMed]
  3. Williams, C.E.; Stevens, K.N. Emotions and Speech: Some Acoustical Correlates. J. Acoust. Soc. Am. 1972, 52, 1238–1250. [Google Scholar] [CrossRef] [PubMed]
  4. Scherer, K. Vocal affect expression: A review and model for further research. Psychol. Bull. 1986, 99, 143–165. [Google Scholar] [CrossRef] [PubMed]
  5. Smith, G.A. Voice analysis for the measurement of anxiety. Br. J. Med. Psychol. 1977, 50, 367–373. [Google Scholar] [CrossRef] [PubMed]
  6. Kagan, J.; Reznick, J.S.; Snidman, N. Biological bases of childhood shyness. Science 1988, 240, 167–171. [Google Scholar] [CrossRef]
  7. Burnham, D.; Kitamura, C.; Vollmer-Conna, U. What’s new, pussycat? On talking to babies and animals. Science 2002, 296, 1435. [Google Scholar] [CrossRef]
  8. Burnham, D.; Francis, E. Are you my little pussy-cat? Acoustic, phonetic and affective qualities of infant-and pet-directed speech. In Proceedings of the 5th International Conference on Spoken Language Processing (ICSLP 98), Sydney, Australia, 30 November–4 December 1998. [Google Scholar]
  9. Ben-Aderet, T.; Gallego-abenza, M.; Reby, D.; Mathevon, N. Dog-directed speech: Why do we use it and do dogs pay attention to it? Proc. R. Soc. B Biol. Sci. 2017, 284, 20162429. [Google Scholar] [CrossRef]
  10. Xu, N.; Burnham, D.; Kitamura, C.; Vollmer-Conna, U. Vowel hyperarticulation in parrot-, dog- and infant-directed speech. Anthrozoos 2013, 26, 373–380. [Google Scholar] [CrossRef]
  11. Hoff, E. Language Development; Jon-David Hague: Wadsworth, OH, USA, 2009. [Google Scholar]
  12. Jeannin, S.; Gilbert, C.; Amy, M.; Leboucher, G. Pet-directed speech draws adult dogs’ attention more efficiently than Adult-directed speech. Sci. Rep. 2017, 7, 4980. [Google Scholar] [CrossRef]
  13. Julius, H.; Beetz, A.; Kotrschal, K.; Turner, D.C.; Uvnäs-Moberg, K. Attachment to Pets; Hogrefe: Göttingen, Germany, 2012. [Google Scholar]
  14. Siniscalchi, M.; d’Ingeo, S.; Quaranta, A. Orienting asymmetries and physiological reactivity in dogs’ response to human emotional faces. Learn. Behav. 2018, 46, 574–585. [Google Scholar] [CrossRef] [PubMed]
  15. Archer, J. Why do people love their pets? Evol. Hum. Behav. 1997, 18, 237–259. [Google Scholar] [CrossRef]
  16. Kotrschal, K. Hund-Mensch. Das Geheimnis der Seelenverwandtschaft; Brandstätter: Wien, Austria, 2016. [Google Scholar]
  17. Cassidy, J.; Shaver, P.R. The Nature of the Child’s Ties. In Handbook of Attachment: Theory, Research, and Clinical Applications; The Guilford Press: New York, NY, USA, 1999; Volume XVII. [Google Scholar]
  18. Bowlby, J. Attachment and Loss. Key Concepts Fam. Stud. 1980, 54, 9–12. [Google Scholar]
  19. Ainsworth, M. Attachments Beyond Infancy. Am. Psychol. 1989, 44, 709–716. [Google Scholar] [CrossRef] [PubMed]
  20. Ainsworth, M.D.S.; Blehar, M.C.; Waters, E.; Wall, S. Patterns of Attachment: A Psychological Study of the Strange Situation; Lawrence Erlbaum: Oxford, UK, 1978. [Google Scholar]
  21. Palestrini, C.; Prato Previde, E.; Custance, D.M.; Spiezio, C.; Sabatini, F.; Verga, M. Heart rate and behavioural responses in dogs (Canis familiaris) in the Ainsworth’s Strange situation stress: A pilot study. In Proceedings of the 3th International Congress on Veterinary Behavioural Medicine, Vancouver, BC, Canada, 7–8 August 2001; Volume 94, pp. 89–91. [Google Scholar]
  22. Main, M.; Solomon, J. Procedures for idenifying infants as disorganaized/disoriented during the Ainsworth Strange Situation. In Attachment in the Preschool Years: Theory, Research, and Intervention; University of Chicago Press: Chicago, IL, USA, 1990; pp. 161–182. [Google Scholar]
  23. Ainsworth, M.D.S.; Wittig, B.A. Attachment and the exploratory behavior of one-year olds in a strange situation. Determ. Infant Behav. 1969, 4, 111–136. [Google Scholar]
  24. Gácsi, M.; Maros, K.; Sernkvist, S.; Miklósi, Á. Does the owner provide a secure base? Behavioral and heart rate response to a threatening stranger and to separation in dogs. J. Vet. Behav. Clin. Appl. Res. 2009, 4, 90–91. [Google Scholar] [CrossRef]
  25. Palmer, R.; Custance, D. A counterbalanced version of Ainsworth’s Strange Situation Procedure reveals secure-base effects in dog-human relationships. Appl. Anim. Behav. Sci. 2008, 109, 306–319. [Google Scholar] [CrossRef]
  26. Prato-Previde, E.; Custance, D.M.; Spiezio, C.; Sabatini, F. Is the dog-human relationship an attachment bond? An observational study using Ainsworth’s strange situation. Behaviour 2003, 140, 225–254. [Google Scholar] [CrossRef]
  27. Topál, J.; Miklósi, Á.; Csányi, V.; Dóka, A. Attachment Behavior in Dogs (Canis familiaris): A New Application of Ainsworth’s (1969) Strange Situation Test. J. Comp. Psychol. 1998, 112, 219–229. [Google Scholar] [CrossRef]
  28. Gergely, A.; Faragó, T.; Galambos, Á.; Topál, J. Differential effects of speech situations on mothers’ and fathers’ infant-directed and dog-directed speech: An acoustic analysis. Sci. Rep. 2017, 7, 13739. [Google Scholar] [CrossRef]
  29. Jeannin, S.; Gilbert, C.; Leboucher, G. Effect of interaction type on the characteristics of pet-directed speech in female dog owners. Anim. Cogn. 2017, 20, 499–509. [Google Scholar] [CrossRef] [PubMed]
  30. Schöberl, I.; Beetz, A.; Solomon, J.; Wedl, M.; Gee, N.; Kotrschal, K. Social factors influencing cortisol modulation in dogs during a strange situation procedure. J. Vet. Behav. Clin. Appl. Res. 2016, 11, 77–85. [Google Scholar] [CrossRef]
  31. Schöberl, I.; Wedl, M.; Beetz, A.; Kotrschal, K. Psychobiological factors affecting cortisol variability in human-dog dyads. PLoS ONE 2017, 12, e0170707. [Google Scholar] [CrossRef] [PubMed]
  32. Costa, P.T.; McCrae, R.R. Revised NEO Personality Inventory (NEO-PI-R) and NEO Five-Factor; The Sage Handbook of Personality Theory and Assessment; Psychological Assessment Resources: Odessa, Ukraine, 1992. [Google Scholar]
  33. Borkenau, P.; Ostendorf, F. NEO-Fünf-Faktoren Inventar:(NEO-FFI); nach Costa und McCrae; Hogrefe: Göttingen, Germany, 1993. [Google Scholar]
  34. Solomon, J.; Beetz, A.; Schöberl, I.; Gee, N.; Kotrschal, K. Attachment security in companion dogs: Adaptation of Ainsworth’s strange situation and classification procedures to dogs and their human caregivers. Attach. Hum. Dev. 2018, 21, 389–417. [Google Scholar] [CrossRef] [PubMed]
  35. Matas, L.; Arend, R.A.; Sroufe, L.A. Continuity of Adaptation in the Second Year: The Relationship between Quality of Attachment and Later Competence. Child Dev. 1978, 49, 547–556. [Google Scholar] [CrossRef]
  36. Grueber, C.E.; Nakagawa, S.; Laws, R.J.; Jamieson, I.G. Multimodel inference in ecology and evolution: Challenges and solutions. J. Evol. Biol. 2011, 24, 699–711. [Google Scholar] [CrossRef] [PubMed]
  37. Amrhein, V.; Korner-Nievergelt, F.; Roth, T. The earth is flat (p > 0.05): Significance thresholds and the crisis of unreplicable research. PeerJ 2017, 5, e3544. [Google Scholar] [CrossRef] [PubMed]
  38. Jørgensen, S.E. Model Selection and Multimodel Inference. Ecol. Model. 2004, 172, 96–97. [Google Scholar] [CrossRef]
  39. Arnold, T.W. Uninformative Parameters and Model Selection Using Akaike’s Information Criterion. J. Wildl. Manag. 2010, 74, 1175–1178. [Google Scholar] [CrossRef]
  40. Mundry, R. Issues in information theory-based statistical inference—A commentary from a frequentist’s perspective. Behav. Ecol. Sociobiol. 2011, 65, 57–68. [Google Scholar] [CrossRef]
  41. Trainor, L.J.; Austin, C.M.; Desjardins, N. Is infanct-directed speech prosody a result of the vocal expression of emotion? Psychol. Sci. 2000, 11, 188–195. [Google Scholar] [CrossRef] [PubMed]
  42. Singh, L.; Morgan, J.L.; Best, C.T. Infants’ Listening Preferences: Baby Talk or Happy Talk? Infancy 2002, 3, 365–394. [Google Scholar] [CrossRef]
  43. Kitamura, C.; Burnham, D. Pitch and communicative intent in mother’s speech: Adjustments for age and sex in the first year. Infancy 2009, 4, 85–110. [Google Scholar] [CrossRef]
  44. Siniscalchi, M.; Stipo, C.; Quaranta, A. “Like owner, like dog”: Correlation between the owner’s attachment profile and the owner-dog bond. PLoS ONE 2013, 8, e78455. [Google Scholar] [CrossRef] [PubMed]
  45. Oswald, L.M.; Zandi, P.; Nestadt, G.; Potash, J.B.; Kalaydjian, A.E.; Wand, G.S. Relationship between cortisol responses to stress and personality. Neuropsychopharmacology 2006, 31, 1583–1591. [Google Scholar] [CrossRef] [PubMed]
  46. Schöberl, I.; Wedl, M.; Bauer, B.; Day, J.; Möstl, E.; Kotrschal, K. Effects of owner-dog relationship and owner personality on cortisol modulation in human-dog dyads. Anthrozoos 2012, 25, 199–214. [Google Scholar] [CrossRef]
Figure 1. The Ainsworth strange situation (ASS) procedure visualized in the form of a timeline. The icons represent all parties of this study. The dog was present at all times and is therefore marked as such in each episode. The episodes marked with the owner icon were used for audio analysis.
Figure 1. The Ainsworth strange situation (ASS) procedure visualized in the form of a timeline. The icons represent all parties of this study. The dog was present at all times and is therefore marked as such in each episode. The episodes marked with the owner icon were used for audio analysis.
Animals 09 00417 g001
Figure 2. Spectrogram of a female owner talking to her dog during a reunion of the ASS procedure with illustrations of the spoken intervals and the measured variables mean f0 (average f0), voiced duration, and f0 range.
Figure 2. Spectrogram of a female owner talking to her dog during a reunion of the ASS procedure with illustrations of the spoken intervals and the measured variables mean f0 (average f0), voiced duration, and f0 range.
Animals 09 00417 g002
Figure 3. Change in mean pitch by gender over the ASS procedure episodes. The male owners are represented in light blue and the female owners in dark blue. Due to physiological differences, pitch (f0) is separated by about an octave between men and women. Therefore, results are presented in separate graphs (mean pitch plotted by gender over ASS procedure (ASSP) episodes, with the dashed line representing the reading controls median. The two calls and reunions are visualized combined as “calls” and “reunions”.) n = 53; number of observations = 4121.
Figure 3. Change in mean pitch by gender over the ASS procedure episodes. The male owners are represented in light blue and the female owners in dark blue. Due to physiological differences, pitch (f0) is separated by about an octave between men and women. Therefore, results are presented in separate graphs (mean pitch plotted by gender over ASS procedure (ASSP) episodes, with the dashed line representing the reading controls median. The two calls and reunions are visualized combined as “calls” and “reunions”.) n = 53; number of observations = 4121.
Animals 09 00417 g003
Figure 4. Pitch range by gender over the NEO FFI character trait openness and throughout the ASSP. The light blue line represents male and the dark blue line represents female owner. (A: Pitch range plotted by gender over ASSP episodes with the dashed line representing the reading controls median. The two calls and reunions are visualized combined as “calls” and “reunions”. B: Pitch range plotted by gender over the NEO FFI openness score.) n = 53, number of observations = 4121.
Figure 4. Pitch range by gender over the NEO FFI character trait openness and throughout the ASSP. The light blue line represents male and the dark blue line represents female owner. (A: Pitch range plotted by gender over ASSP episodes with the dashed line representing the reading controls median. The two calls and reunions are visualized combined as “calls” and “reunions”. B: Pitch range plotted by gender over the NEO FFI openness score.) n = 53, number of observations = 4121.
Animals 09 00417 g004
Figure 5. Percentage of jitter by gender over the dog owners’ NEO FFI character trait openness and throughout the ASSP. The light blue line represents male and the dark blue line represents female owner. (A: Jitter plotted by gender over ASSP episodes with the dashed line representing the reading controls median. The two calls and reunions are visualized combined as “calls” and “reunions”. B: Jitter plotted by gender over the NEO FFI openness score.) n = 53; number of observations = 4121.
Figure 5. Percentage of jitter by gender over the dog owners’ NEO FFI character trait openness and throughout the ASSP. The light blue line represents male and the dark blue line represents female owner. (A: Jitter plotted by gender over ASSP episodes with the dashed line representing the reading controls median. The two calls and reunions are visualized combined as “calls” and “reunions”. B: Jitter plotted by gender over the NEO FFI openness score.) n = 53; number of observations = 4121.
Animals 09 00417 g005
Table 1. Model averaged coefficients for mean pitch models. Only variables with a minimum relative importance of 1 and a confidence interval not ranging over 0 are taken into consideration for model interpretation.
Table 1. Model averaged coefficients for mean pitch models. Only variables with a minimum relative importance of 1 and a confidence interval not ranging over 0 are taken into consideration for model interpretation.
ParameterEstimateStd. ErrorConfidence Interval (85%)Relative Importance
0.0750.925
(Intercept)4.76580.04294.70414.8276
speaking control−0.0050.0249−0.04090.03091
exploration0.19040.01960.16210.2187´´
stranger0.18590.02630.1480.2237´´
call (1)0.38630.03530.33550.4371´´
reunion (1)0.18220.01940.15440.2101´´
call (2)0.37560.03480.32550.4257´´
reunion (2)0.14530.01970.1170.1736´´
gender (female)0.55130.05370.47390.62861
caregiving−0.00270.0121−0.05960.01970.13
extraversion0.0020.0101−0.01860.05060.13
attachment (non-secure)−0.00350.0093−0.0301−0.02590.13
neuroticism0.00150.0106−0.02810.05450.12
openness0.00090.0099−0.0330.05020.11
agreeableness0.00080.0089−0.02990.04510.11
Intercept includes reading control, secure attachment, and owner gender male.
Table 2. Model averaged coefficients for frequency range models. Only variables with a minimum relative importance of >0.8 and a confidence interval not ranging over 0 are taken into consideration for model interpretation.
Table 2. Model averaged coefficients for frequency range models. Only variables with a minimum relative importance of >0.8 and a confidence interval not ranging over 0 are taken into consideration for model interpretation.
ParameterEstimateStd. ErrorConfidence Interval (85%)Relative Importance
0.0750.925
(Intercept)4.34260.09114.21144.4738
speaking control−0.16560.1013−0.3114−0.01991
exploration−0.60980.0799−0.7248−0.4948´´
stranger−0.75180.1072−0.9062−0.5974´´
call (1)−0.29620.1436−0.5028−0.0895´´
reunion (1)−0.59740.0787−0.7107−0.4841´´
call (2)−0.33510.1415−0.5388−0.1315´´
reunion (2)−0.6180.079−0.7317−0.5042´´
gender (female)0.7880.07820.67550.90061
openness0.06130.04360.01790.1270.85
extraversion0.01570.0303−0.00160.09970.32
agreeableness0.00380.0159−0.01580.08580.11
neuroticism0.00590.0217−0.02590.09450.17
attachment (non-secure)0.00530.0281−0.04890.16680.09
conscientiousness0.00180.012−0.02930.07510.08
caregiving−0.00060.0099−0.06240.04350.07
Intercept includes reading control, secure attachment, and owner gender male.
Table 3. Model averaged coefficients for jitter models. Only variables with a minimum relative importance of >0.8 and a confidence interval not ranging over 0 are taken into consideration for model interpretation.
Table 3. Model averaged coefficients for jitter models. Only variables with a minimum relative importance of >0.8 and a confidence interval not ranging over 0 are taken into consideration for model interpretation.
ParameterEstimateStd. ErrorConfidence Interval (85%)Relative Importance
0.0750.925
(Intercept)−4.30290.0599−4.389−4.2167
speaking control0.21930.05320.14270.29591
exploration0.17950.04180.11930.2397´´
stranger0.03090.0563−0.05010.1119´´
call (1)−0.41040.0754−0.5189−0.3019´´
reunion (1)0.1460.04130.08650.2055´´
call (2)−0.29830.0746−0.4057−0.191´´
reunion (2)0.16790.04150.10820.2276´´
gender (female)−0.30870.0652−0.4025−0.21481
openness0.05310.03580.01560.10630.87
caregiving0.00370.0145−0.01810.06860.15
extraversion−0.00230.0115−0.05750.02120.13
agreeableness0.0010.01−0.03310.05090.11
attachment (non-secure)−0.00130.0207−0.10250.07840.1
neuroticism−0.00040.0103−0.04960.04220.1
Intercept includes reading control, secure attachment, and owner gender male.

Share and Cite

MDPI and ACS Style

Lesch, R.; Kotrschal, K.; Schöberl, I.; Beetz, A.; Solomon, J.; Fitch, W.T. Talking to Dogs: Companion Animal-Directed Speech in a Stress Test. Animals 2019, 9, 417. https://doi.org/10.3390/ani9070417

AMA Style

Lesch R, Kotrschal K, Schöberl I, Beetz A, Solomon J, Fitch WT. Talking to Dogs: Companion Animal-Directed Speech in a Stress Test. Animals. 2019; 9(7):417. https://doi.org/10.3390/ani9070417

Chicago/Turabian Style

Lesch, Raffaela, Kurt Kotrschal, Iris Schöberl, Andrea Beetz, Judith Solomon, and W. Tecumseh Fitch. 2019. "Talking to Dogs: Companion Animal-Directed Speech in a Stress Test" Animals 9, no. 7: 417. https://doi.org/10.3390/ani9070417

APA Style

Lesch, R., Kotrschal, K., Schöberl, I., Beetz, A., Solomon, J., & Fitch, W. T. (2019). Talking to Dogs: Companion Animal-Directed Speech in a Stress Test. Animals, 9(7), 417. https://doi.org/10.3390/ani9070417

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop