2.1. Study 1: Impact of familiar companion
The median value of companion’s reactivity score in the novel object test was 0.5 (Quartile 1 (Q1) = 0; Quartile 3 (Q3) = 1) for habituated companions versus 1 (Q1 = 0; Q3 = 2.5) for naïve companions. These values were not significantly different (Mann-Whitney test: W=23.5, n1=8, n2=8, P = 0.372). In the umbrella test, the median value of companion’s reactivity score was significantly lower at 1 (Q1 = 1; Q3 = 1.5) for habituated companions compared to 2 (Q1 = 2; Q3 = 3) for naïve companions (Mann-Whitney test: W=4, n1=8, n2=7, P = 0.004).
Only during the novel object test, did any subjects investigated the ball, and these numbers were very low (2 in the alone condition and 2 others in the condition with the companion).
2.1.1. Influence of context
Reactivity scores were analysed for 16 pairs for the novel object test. In the novel object test, subjects in pairs had significantly lower reactivity scores than those alone (Cumulative link mixed models (CLMM): estimate ± S.E.: −2.12 ± 0.88, P = 0.016) (see Figure 1a).
The mean value of the subject’s pre-test heart rate measurement before the experiment was 42.29 (SD = 4.58) and significantly increased to 48.33 (SD = 12.54) during the novel object test (Wilcoxon tests: V = 6, n = 13, p = 0.011).
The analyses of heart rate recovery were run on 12 pairs in the novel object test. The difference in time to recover obtained from the subjects in pairs and when alone was not significantly different (Linear mixed effect models (LMM): t-value, df: −1.336, 11, P = 0.208) (see Figure 1b).
The companion’s reactivity score did not have an impact on the subject’s reactivity score in the novel object test (Cumulative link model (CLM): estimate ± S.E.: 0.31 ± 0.51, P = 0.535).
Reactivity scores were analysed for 14 pairs for the umbrella test. In this test, considering only subjects who showed a reaction, the reactivity score for the subjects in pairs was not significantly lower than that obtained from those alone (CLMM: estimate ± S.E.: −1.14 ± 0.78, P = 0.145) (see Figure 2a).
The mean value of the subject’s pre-test heart rate measurement before the experiment was 41.50 (SD = 3.66) and significantly increased to 50.68 (SD = 8.66) during the umbrella test (Wilcoxon tests: V = 7, n = 13, p = 0.004).
The data on heart rate recovery were analysed for 14 pairs for the umbrella test. The presence of a companion was associated with quicker heart rate recovery (LMM: t-value, df: −2.319, 13, P = 0.037) (see Figure 2b).
The reactivity score of the companion did not have an impact on the subject’s reactivity score in the umbrella test by pair (CLM: estimate ± S.E.: −0.37 ± 0.77, P = 0.623).
2.1.2. Impact of companion habituation
The difference in reactivity scores and heart rate recovery were not significantly different between subjects with habituated companions and subjects with naïve companions regardless of the test (Reactivity score (RS), CLM: Novel object test (NO), estimate ± S.E.: 0.36 ± 0.95, P = 0.705; Umbrella test (UM), estimate ± S.E.: 0.20 ± 1.36, P = 0.881; Heart rate recovery (HRR), Linear effect model (LM): NO, t-value, df: −1.027, 10, P = 0.329; UM, t-value, df: 0.216, 12, P = 0.833).
2.2. Study 2: Impact of unfamiliar companion
The median value of companion’s reactivity score in the novel object test was 0 (Quartile 1 (Q1) = 0; Quartile 3 (Q3) = 0) for habituated companions versus 0 (Q1 = 0; Q3 = 1) for naïve companions. These values were not significantly different (Mann-Whitney test: W=25, n1=8, n2=8, P = 0.369).
In the umbrella test, the median value of companion’s reactivity score was significantly lower at 1 (Q1 = 0.5; Q3 = 3) for habituated companions versus 4 (Q1 = 2.5; Q3 = 4) for naïve companions (Mann-Whitney test: W=11, n1=8, n2=7, P = 0.047).
Only one subject investigated the umbrella during the umbrella test with a companion. Also, it has been observed that during the tests 11 pairs had exchanged agonistic behaviour such as head threat when they ate, i.e. one of the horses turned its heard toward the other with the ears pinned back.
2.2.1. Influence of context
For the novel object test, the analysis of reactivity scores was run on 12 pairs. The subject’s reactivity score in the novel object test for paired subjects was significantly lower than when alone (CLMM: estimate ± S.E.: −2.87 ± 1.19, P = 0.015) (see Figure 3a).
The mean value of the subject’s pre-test heart rate measurement before the experiment was 48.34 (SD = 4.73) versus 50.53 (SD = 10.52) for the novel object test. These values are not significantly different (Wilcoxon tests: V = 56, n = 16, p = 0.552).
The data on heart rate recovery were analysed for 11 pairs in the novel object test. The time to recover for the subjects in pairs was not significantly different to those alone (LMM: t-value, df: −1.209, 10, P = 0.254) (see Figure 3b).
In the novel object test pairs, the companion’s reactivity score did not have an impact on the subject’s reactivity score (CLM: estimate ± S.E.: 0.07 ± 0.52, P = 0.885).
In the umbrella test, the reactivity scores were analysed for 15 pairs. The reactivity score for the subjects in pairs was not significantly lower than that obtained from those alone (CLMM: estimate ± S.E.: −1.08 ± 0.96, P = 0.259) (see Figure 4a).
The mean value of the subject’s heart rate pre-test measurement before the experiment was significantly lower at 48.34 (SD = 4.73) versus 71.54 (SD = 33.04) during the umbrella test (Wilcoxon tests: V = 28.5, n = 16, p = 0.044).
The data on heart rate recovery were analysed for 15 pairs for the umbrella test. Subjects in pairs had significantly lower recover times compared than those alone (LMM: t-value, df: −2.229, 14, P = 0.043) (see Figure 4b).
The subject’s reactivity score was not impacted by the reactivity score of the companion in the umbrella test (CLM: estimate ± S.E.: −0.68 ± 0.46, P = 0.137).
2.2.2. Impact of companion habituation
Regardless of the test, there were no significant differences in reactivity scores and heart rate recovery, between subjects with a habituated companion and subjects with a naïve companion (RS, CLM: NO, estimate ± S.E.: 0.37 ± 1.08, P = 0.731; UM, estimate ± S.E.: 1.92 ± 1.42, P = 0.174; HRR, LM: NO, t-value, df: −1.477, 9, P = 0.173; UM, t-value, df: 1.287, 13, P = 0.221).