A Wearable System for the Evaluation of the Human-Horse Interaction: A Preliminary Study

: This study reports on a preliminary estimation of the human-horse interaction through the analysis of the heart rate variability (HRV) in both human and animal by using the dynamic time warping (DTW) algorithm. Here, we present a wearable system for HRV monitoring in horses. Speciﬁcally, we ﬁrst present a validation of a wearable electrocardiographic (ECG) monitoring system for horses in terms of comfort and robustness, then we introduce a preliminary objective estimation of the human-horse interaction. The performance of the proposed wearable system for horses was compared with a standard system in terms of movement artifact (MA) percentage. Seven healthy horses were monitored without any movement constraints. As a result, the lower amount of MA% of the wearable system suggests that it could be proﬁtably used for reliable measurement of physiological parameters related to the autonomic nervous system (ANS) activity in horses, such as the HRV. Human-horse interaction estimation was achieved through the analysis of their HRV time series. Speciﬁcally, DTW was applied to estimate dynamic coupling between human and horse in a group of fourteen human subjects and one horse. Moreover, a support vector machine (SVM) classiﬁer was able to recognize the three classes of interaction with an accuracy greater than 78%. Preliminary signiﬁcant results showed the discrimination of three distinct real human-animal interaction levels. These results open the measurement and characterization of the already empirically-proven relationship between human and horse.


Introduction
In the last few decades, the interest in decoding the human-horse relationship and interaction (HHRI) has increased dramatically.This was guided by the strong empirical evidence of the positive outcomes in equine assisted therapy (EAT) and horseback riding in therapeutic programs [1], as well as the positive impact of animal companionship on human quality of life [2], where the equine is an important element of these therapeutic practices, with its feelings and behavior.For this purpose, the investigation of the modality in which both human and horse can communicate might be crucial.Measuring and evaluating the impact of the interaction experience might be relevant [3].Some studies investigated the equine perception of humans in terms of positive, negative or neutral valence [4].A study on how human psychological and physiological state can be perceived by horses was performed in [5] via the study of the heart rate.A more relaxed equine behavior was observed when humans showed positive attitudes toward them [6,7], while an equine increased heart rate was observed when humans were engaged in negative thinking [8].A nervous mood can be transmitted from humans to horses under handling and riding conditions [9].Voice [10][11][12], posture [13,14], facial expression [15,16], autonomic signals [17][18][19], hormones [20][21][22][23] and pheromones [24] might be used to fruitfully describe and characterize the emotional content [25].Non-verbal communication between human and horse was also investigated in [26,27].Heart rate and behavior resulted in being sensitive and reliable indicators of fear or anxiety in horses [28,29].Horses that are in discomfort were observed to be more aggressive toward humans [30] or to be characterized by an increased heart rate, motor activity and vocalizations [31].The effect of the gender of the person interacting with the horse was discussed in [32].Although a recent review described a parallel behavior between the human multi-sensorial perception and the demonstrated equine cross-modal recognition [33,34], an interdisciplinary approach is mandatory to reach a deeper knowledge of human-horse interaction.In fact, these achieved experimental findings pointed out complex and multidimensional aspects of the interaction, which involve medical, bioengineering, physics and veterinary science [2,35,36].In [9], the heart rate of both human and horse were monitored simultaneously under different experimental handling or riding conditions.In this study, Keeling et al. asserted that the analysis of heart rate is an important tool to investigate horse-human interactions.Again, hormones, heart rate and some standard heart rate variability-related indices were investigated in [20] during both training and performance.Different feature trends were observed between human and horse when they were obtained from ECG (electrocardiographic) recordings related to public or private sessions.A body sensors network technology was used to real-time monitor the horse-rider dyad in [37].The aim of such a monitoring was evaluating the human-horse interaction.Specifically, based on a study concerning the equine emotional response during physical activity [38], the evaluation was based on the measurements of heart rate and physical activities via a mathematical model.Such a model was proposed to decompose the equine heart rate into two different components.The first one was concerning the physical component, while the second one contained information about equine emotional state [38].Finally, in [39], the link between horse and human was also investigated by studying their electroencephalograph signals (EEG), revealing a higher synchronicity in EEG waves at increasing interactions.A correlation analysis between human and equine hormone concentrations was performed in [20,21].
Therefore, the human-horse relationship appears to be a complex interaction affected by several psychological factors [40].The perception of humans by horses seems to be based on experience and repeated interactions, with horses that form a memory of humans that impacts their reactions in subsequent interactions [5].Hence, previous negative experiences with human contact could lead horses toward a negative emotional reaction [7] or, vice versa, previous positive experience could lead them toward positive feelings with humans [41].
In our hypothesis, human and horse are considered as complex systems that interact through a coupling process.In this frame, we hypothesize that coupling can be modulated by the type and time duration of the contact itself.Specifically, we analyzed the level of coupling by studying their heart rate variability (HRV) time series [42] through dynamic time warping (DTW).As is well known from the literature, HRV can be considered as a non-linear time series, in which complex oscillations are present [43].Therefore, we aimed at measuring this biological coupling over time [44].DTW is, generally, used for studying time series dynamics of non-stationary systems.It calculates the best possible warp alignment between two time series, by selecting the one with the minimum distortion.Specifically, DTW is also defined as a measure of similarity between two time series, and it is calculated as the minimum mapping distance between them.DTW was widely used in many contexts, including data mining [45], speech processing [46] and medicine [47].
In this study, DTW was adopted to evaluate how heart activities evolve in a similar or dissimilar way.For example, if during an experiment with increasing exciting levels we detect an increasing DTW between human and equine HRVs, it indicates that the distance of the two HRV dynamics is increasing, and therefore, the heart activities are following different patterns.In order to perform a continuous, comfortable and non-invasive monitoring of the interaction in a natural environment, we developed, and here present, a wearable monitoring system for horses.The amount of advantages that wearable systems have brought to physiological signal monitoring for humans is well known, occupying an even larger space in the research.Moreover, the continuous technological development and the increasing demand of smart systems push wearables as the most used and suitable systems for ubiquitous and pervasive investigations.In addition, their flexibility allows researchers and clinicians to face the large variability of biomedical signals and tasks in monitoring subjects during their daily activities [48][49][50][51][52].However, the biggest limitation in using wearable systems with humans and animals is due to motion artifacts (MA) [53][54][55][56][57][58], which are the major source of noise in biomedical signal acquisition, inducing the loss and alteration of informative content.For instance, the electrocardiographic (ECG) signal acquired in a naturalistic environment without movement constraints can be severely affected by important artifacts, and a great amount of data might be lost in contrast to the signal quality easily achievable in controlled environments and protocols [59].As a matter of fact, cardiac stress tests or simply respiration can generate a big amount of MAs that can alter the signal [59].Moreover, it is important to highlight that restraining horse is usually discouraged since it is unnatural and stressful and induces an increase of the sympathetic contribution of the heart control [60] that leads to misleading ECG interpretations [61].In this work, we present a textile-based system for the ECG monitoring in horses, where the electrodes are completely made of fabric (electro-textile or e-textile).Normally, textile materials are insulators, but for this application, conductive yarns are integrated into the fabric during the manufacturing process [62].In the human biomedical field, e-textiles are considered as higher value-added textiles and are prominently developed for being used in smart clothing.Smart clothing refers to a new garment that is able to acquire and process information, as well as actuate responses [63].The potentialities deriving from these kinds of textile sensors enable the application of wearable systems in a great variety of experimental settings.As a matter of fact, many human studies showed reliable recordings of biomedical signals [64][65][66][67][68]: for example, ECG in [69][70][71][72][73][74][75][76], respiration in [69,70,73,74,[77][78][79], electrodermal response (EDR) in [51,80] and, finally, PhotoPlethysmoGraphy (PPG)and blood pressure in [81,82].In this kind of application, physiological signals are monitored and recorded in order to evaluate or follow the health status of the person who is wearing the wearable device [48,83].
Similarly in the veterinary research field, some authors [84][85][86][87] proposed to use Holter devices in equine applications and subsequently systems with radiotelemetry.Here, we propose the use of wearable systems in both animals and humans simultaneously, in order to acquire their cardiac activity in a reliable and artifact-free way, as was generally demonstrated in [72,[88][89][90][91][92].The aim is to infer autonomic nervous system responses enabling the detection of uncontrolled responses of animals when elicited with emotional stimuli coming from humans and vice versa.The manuscript is organized as follows: Section 2 deals with the materials' and methods' description; specifically, it describes the wearable systems, the experimental protocol of the interaction and the signal processing chain.Section 3 reports the achieved results, and Section 4 is focused on the discussion and conclusion of the study highlighting the future perspectives of this pioneering work.

Materials and Methods
In this study, two wearable monitoring systems, i.e., one for humans and the other for horses, were employed.Each system was comprised of a smart garment and portable electronics.The human system (Figure 1) was designed as a sensorized t-shirt (Smartex, Pisa, IT, Italy), and it is exhaustively described in [48,79,[93][94][95].Differently, the equine system (Figure 2) was comprised of an elastic smart belt (Smartex, Pisa, IT, Italy) [89,91] fastened around the chest behind the shoulder area.In both systems, two textiles electrodes (Smartex, Pisa, IT, Italy) and a strain gauge sensor (Smartex, Pisa, IT, Italy) were integrated to acquire the ECG (with a sampling frequency of 250 Hz) and respiration activity (with a sampling frequency of 25 Hz).The strain gauge is carried out by textile sensors that monitor the cross-sectional variations of the rib cage.The respiration sensor along with electrodes are integrated in an elastic band through a one-step process in the fabric by means of a circular knitting machine [96].They are developed by Smartex s.r.l.; many details can be found in [92,96].Moreover, equine ECG was recorded by placing the electrodes according to the modified base-apex configuration [97].It is worthwhile to note that the use of a dry textile-based electrode provides several advantages.Firstly, the system is easy to use through an automatic placement of the sensors and allows high comfort.Secondly, electrodes are made of a special multilayer structure of textile material that increases the amount of sweat and reduces the rate of evaporation reaching very rapidly an electrochemical equilibrium between skin and electrode.This means that the signal quality [72] is remarkably improved and kept as constant as possible.These materials are knitted together and are fully integrated in the garment without any mechanical and physical discontinuity, creating areas with different functionalities (see Figure 1).For each system, the two ECG e-textile electrodes and the strain gauge sensor are finally connected to portable electronics through a simple plug that can be easily unplugged when necessary.Wearable monitoring system for humans [48].As is possible to note, from the box on the right, the e-textile electrodes for ECG (electrocardiographic) acquisition are knitted and completely integrated into the garment.Moreover, an inertial platform (triaxial accelerometer) (Smartex, Pisa, IT, Italy) integrated into the portable electronics that was positioned on the back of the horse allowed the monitoring of the horse's physical activity at a sampling frequency of 25 Hz.Finally, the wearable system was wirelessly connected to a smartphone where a dedicated app enabled checking the status of the sensors and remotely controlling the storing process of the physiological information in a secure digital (SD) card.

Equine Textile-Based ECG: Test and Validation
This section is focused on testing and validating the reliability of the equine ECG traces coming from textile electrodes.For this study, seven healthy standard bred mares in anestrus (mean age 8.4 ± 1.3 years) were enrolled.Equine subjects were socially housed in a paddock (75 × 75 m) and were provided ad libitum access to both hay and water.Horses were used as receivers in the embryo transfer program of the Department of Veterinary Sciences (University of Pisa, Italy) where this study was performed.Mares were in healthy condition and not pregnant at the time of this protocol.
The signal quality test of e-textile electrodes was performed by means of a comparative study in terms of the motion artifact (MA) between e-textiles and conventional Ag/AgCl electrodes.Each horse was simultaneously monitored by the wearable system, i.e., equine elastic smart belt and Ag/AgCl electrodes (see Figure 2).Two identical electronics were employed for both couples of electrodes Ag/AgCl and the textile ones.All data were acquired in a stall (4 × 4 m), where horses were left free to spontaneously move for 60 min.At the end of the recording session, all of the acquired ECGs were visually examined by one expert, and all of the ECG segments, corrupted by MAs, were marked.The goal was the estimation of the percentage of corrupted signal.Such a percentage was calculated as the sum of the time intervals in which an MA was observed over the whole time length of the recording.
A nonparametric Wilcoxon signed-rank test for paired samples was used to compare the percentage of corrupted signals between the signals acquired by the two kinds of electrodes, i.e., e-textile and Ag/AgCl.Specifically, the test was designed to compare the performance of the two kinds of electrodes while simultaneously recording the heart activity of each horse.A significant result of such a test would indicate a coherent different percentage of corrupted signal between the traces acquired by means of the two kinds of electrodes.

Protocol of Interaction
Here, the design of the interaction protocol is reported.Fourteen subjects (25 ± 3 years old, 4 males) were enrolled.The participants did not show any past or current experience of mental or personality disorder.In addition, one standard bred mare out of the seven previously described was enrolled (age: 8 years, weight: 560 kg, height: 160 cm).Informed consent was signed by the participants according to this specific protocol approved by the Ethical Committee of the University of Pisa.
During the whole experimental protocol, the autonomic nervous system (ANS) response of both human and horse were monitored.Specifically, two different systems were used to acquire the ECG signal.The fabric-based monitoring system [48,79,[93][94][95] was used to record the human ECG, while the elastic smart belt [89,91] was used to record the equine one.
The experimental protocol consisted of three phases, each one lasting 4 min.During the first phase, P 1 , the human and the mare were in different stalls (4 × 4 m).In this phase, considered as the resting phase, the subject sat on a chair, while the horse was free to move.Successively, the horse was moved from her stall to the human's.In this phase, P 2 , the subject was asked to keep himself/herself still on the chair, while the horse was free to move and explore the environment.This phase implies a visual and olfactory interaction.Finally, in the third phase, P 3 , all of the participants were asked to stand up and groom the horse.It is worthwhile to note that in order to keep the visual conditions of the horse as constant as possible, all of the subjects were asked to wear an azure plain t-shirt.

Signal Processing Chain for Human-Horse Interaction
This section deals with the description of the signal processing chain, which can be summarized in three steps.First, the HRV of both human and horse were estimated.Then, the DTW coupling index is computed.Finally, in the third step, the statistical analysis of DTW in the three sessions and the pattern recognition of the three phases is performed.

Heart Rate Variability Estimation
HRV is considered to be one of the most important ANS-related series.HRV describes the variations in the beat-to-beat intervals or in the instantaneous heart rate, since it reflects the regulatory mechanism of the cardiac activity by the ANS [98].ANS is divided into sympathetic and parasympathetic branches.According to a simple model, sympathetic activity is responsible for the increasing of the heart rate (HR) and of the decreasing of HRV, while parasympathetic activity is usually considered to be in charge of decreasing HR and increasing HRV [99].Many studies have focused on the significant relationship between ANS and HRV.Specially, frequency domain indexes, e.g., the LF/HF ratio (i.e., it is the ratio between the power in the low frequency bandwidth and the power in the high frequency bandwidth), have been deeply investigated [100], as well as temporal [19], spectral [101] and non-linear [102][103][104] indices.Several tools have been developed for its analysis [105], since its estimation is quite simple.In fact, HRV can be estimated as the series obtained from the interpolation of the beat-to-beat distances, i.e., RR distances (i.e. it means the distance between consecutive R-peaks of QRS complex).For this purpose, ECG signals were first of all digitally filtered with a zero-phase Butterworth infinite impulse response band-pass filter with cut-off frequencies between 0.5 and 40 Hz.Then, different R-peak detection algorithms were used to process human and equine ECG signals with the aim of estimating the HRV.More precisely, R-peaks related to the human ECG signals were detected by means of the well-known Pan-Tompkins method [106], while the method proposed in [91] was used to detect the R-peaks in the equine ECG signals.The Pan-Tompkins method [106] is an algorithm based on a pre-processing phase, including band-pass filtering, squaring of the data samples and moving average filtering, and on a decision rule phase, which includes an amplitude threshold to detect R-peaks.Differently, the algorithm proposed in [91] to detect R-peaks in equine ECG was based on the estimation of the energy of the second derivative of the ECG signal.R-peaks were detected by performing a thresholding of the obtained energy-signal.Finally, a cubic spline interpolation was applied to the irregular RR sampled series with a new sampling frequency equal to 10 Hz [105].To do that, the cubic spline interpolation was used to create an evenly-distributed sampled series.This procedure is necessary since RRseries are irregularly sampled, due to the physiological variability of the heart activity, and an evenly-distributed sampled series is a mandatory condition in several analyses.Then, since the RR interval sampling frequency is usually chosen among 2, 4, 6, 8 and 10 Hz [107], we decided to obtain a sampling frequency equal to 10 Hz as performed in [108,109].In fact, according to [107], it is important to highlight that the bandwidth within which the autonomic nervous system has a significant response is 0-1 Hz.

Feature Estimation
This section reports on the computation of the dynamic time warping [110].DTW is usually used for studying time series dynamics of non-stationary systems.It computes the best possible warp alignment between two time series, by selecting the one with the minimum distortion.DTW was widely used in many contexts including data mining [45], speech processing [46] and medicine [47].Specifically, DTW is a measure of the similarity between two temporal series.To estimate it, a temporal non-linear warping is performed to find the optimal match between the two sequences.In this frame, given two sequences, x and y of length N and M, respectively, it is possible to define a (N,M)-warping path as a sequence p = (p 1 , ..., p L ), with p l = (n l , m l ) ∈ [1 : N] × [1 : M] for l∈ [1 : L], able to align the elements of x and y.Such an alignment assigns the element x n l of x to the element y m l of y.Hence, the total cost c p (x, y) of a warping-path p between x and y can be defined, with respect to the local cost measure or local distance measure c, as reported in Equation (1).
Finally, the optimal warping path between x and y is the warping path p * corresponding to the minimal total cost.DTW distance is defined as the cost of p * , according to (Equation ( 2)): where p is a (N,M)-warping path.This means that a low DTW distance implies a short warping path, which results in a higher similarity between the two time series.In this study, DTW is used to estimate the degree of similarity/dissimilarity between the two subjects' heart activities (i.e., horse and human), as well as the degree of similarity/dissimilarity among all of the phases within each single subject (i.e., only human or horse).To this aim, the DTW method is applied to both human and horse HRV series.

Statistical Analysis and Pattern Recognition
A Wilcoxon signed rank test for paired data was used to compare the features estimated during the different phases of the experimental protocol.All of the possible combinations of the 2-class problems were taken into account.Specifically, the DTW, estimated between human and horse in each phase, was compared with the one obtained from the other phases, i.e., P 1 vs. P 2 , P 1 vs. P 3 and P 2 vs. P 3 , to evaluate possible statistically-significant differences between the human-horse interactions occurring during the different phases.Moreover, DTW was also estimated among the HRV series of each single enrolled subject between the different phases, i.e., only horse and only human.In this case, we aimed at comparing the degree of similarity/dissimilarity passing through the different phases within the same subject.A correction factor for multiple comparisons was applied according to the Benjamini and Hochberg method [111].Such a correction was necessary to cope with the rate of Type I errors in null hypothesis testing when conducting multiple comparisons.The Benjamini and Hochberg method, which is a false discovery rate (FDR)-controlling procedure, takes into account the expected proportion of rejected null hypotheses that were incorrect rejections ("false discoveries") [111].A significance level equal to 0.05 was used.
The classification process was performed by means of a supervised learning method, which aimed at performing the recognition of the experimental phases, i.e., the human-horse interaction, and evaluating the discriminant power of the DTW feature, as estimated between human and horse.Specifically, we implemented a support vector machine (SVM) classifier, with a radial-based kernel function (see Equation ( 3)): where γ is equal to 1 and x and x stand for the two samples.Moreover, a Leave-One-Subject-Out (LOSO) architecture was developed to apply our pattern recognition approach on the estimated DTW.Of note, the LOSO architecture is appositely designed to test the developed classifier on untreated and unknown data.Specifically, if N is the number of enrolled subjects (in this case, N = 14), the classifier is trained N − 1 times on the N − 1 subjects and, therefore, tested on the remaining 1 subject.To this aim, first of all, the data related to each subject were Z-transformed and then given as input to the LOSO architecture.

Results of the E-Textiles' Validation and Testing
The percentages of artifact signal recorded both by using Ag/AgCl electrodes and e-textile electrodes were estimated.Such an estimation was performed by manually placing the labels.A person with experience in the field of equine cardiology was in charge of this issue.The obtained results are displayed in Table 1.In addition, a non-parametric Wilcoxon signed-rank test for paired samples was implemented to analyze the statistical significance of the differences observed between the two kinds of electrodes in terms of the percentage of corrupted signal.A significant p-value equal to 0.0156 was achieved (Figure 3).This indicated a statistically-significant difference between the performance of the two kinds of electrodes.In fact, as is easily detectable in Figure 3 and in Table 1, the percentage of corrupted signal was significantly lower when ECG traces were recorded via e-textile electrodes.

Results of the Human-Horse Interaction
In Figure 4, an example of two HRVs estimated from both equine and human ECG traces is shown.There, it is possible to observe, within a smaller signal portion, how warping-paths are defined.The statistical analysis performed on the DTW estimated between the human and equine HRV traces of the three different phases is reported in Table 2 in terms of p-values.All of the pairwise comparisons between different phases were statistically significantly different after performing the p-value correction via the Benjamini and Hochberg method.Specifically, as highlighted in Table 2 by means of the rising arrows, a loss of similarity, i.e., an increased DTW distance value, can be observed as the human-horse interaction becomes progressively stronger.These trends can be also observed in Figure 5.In fact, in 12 out of 14 subjects, the DTW value can be observed to be lower or equal in P 1 compared to in P 2 , i.e., human and equine HRV are more similar in the P 1 phase than in P 2 .In 12 out of 14 subjects, the DTW value can be detected to be higher in the P 2 phase than in P 1 , but lower than in P 3 , while 12 out of 14 subjects show a higher DTW value in the P 3 phase with respect to the other ones.In Table 3, the statistical results obtained by comparing only the human DTW and only horse DTW, separately, between the different phases are reported.Observing the reported p-values and arrows, a significant difference is reported for horse DTW between P 1 -P 2 and P 2 -P 3 with a higher value for the P 2 -P 3 .For humans, a significant difference was found in the comparison of DTW between P 1 -P 2 and P 1 -P 3 , showing a higher value for the P 1 -P 3 , and the comparison of DTW between P 1 -P 3 and P 2 -P 3 , showing a higher value for the P 1 -P 3 .These results seem to show that the P 3 phase is the more responsive for both humans and horses.The results of the pattern recognition are returned as a confusion matrix (Table 4).In the main diagonal, each element represents the correct recognition percentage.Upper and lower triangular matrices report the percentage of incorrect recognition.The overall accuracy is about 85.71%.

Discussion and Conclusions
The objective of the study here presented was two-fold: the first one was to investigate the reliability of using wearable systems for monitoring physiological signal in horses; the second objective was to estimate and quantify the human horse interaction in a specific experimental protocol.As a matter of fact, concerning the first objective, the obtained results showed that the e-textile-based system outperformed the standard monitoring one, in terms of comfortability and reduction of the amount of movement artifacts.This latter one was demonstrated by the comparison of the blind visual inspection performed by an expert applied to both signals.Specifically, an amount of 7 h of equine ECG recording was analyzed, and the statistical analysis showed that ECG signals acquired with the textile electrodes had a lower percentage of corruption, revealing a greater robustness against movement artifacts.It is worthwhile to note that textile-based wearable monitoring systems offer interesting advantages in horse applications.In fact, it resulted in being very easy to use and very comfortable, and it could be automatically placed without any adhesive, glue and adhesive bandages.Furthermore, the multilayer structure of the textile electrode increased the amount of sweat and reduced the rate of evaporation, reaching an optimal electrochemical equilibrium between skin and electrode, rapidly.As a matter of fact, the signal quality was notably improved and kept as constant as possible.Moreover, wearable systems allow us to reliably monitor heart activity and parameters related to the involvement of ANS activity also without any restriction of movement.Indeed, it is possible to record an artifact-free ECG that enables non-linear signal processing techniques' application to horse's HRV signal.This allowed us to perform an objective measurement of the complex ANS responses for the investigation of human-horse interaction over a temporal dynamic evolution.Specifically, by analyzing the dynamic of the HRV of both human and horse, a similarity estimator, i.e., DTW, has been studied to describe different interaction phases.DTW comparisons for horse showed that the transition from P 2 to P 3 is more significant with respect to the transition from P 1 to P 2 , whilst for humans, we found that the difference between P 1 -P 3 is higher than both P 1 -P 2 and P 2 -P 3 .These results seem to show that the P 3 phase is the more responsive for both humans and horses.Moreover, the achieved results in the horse-human DTW estimation showed that it was able to significantly discriminate all of the three phases of the experimental protocol.Furthermore, a significantly continuous increment from P 1 to P 2 and then to P 3 was observed in human-horse DTW.This may indicate a monotonic variation of the coupling between the human and horse cardiovascular systems.Specifically, as highlighted in Table 2, a loss of similarity, i.e., an increasing DTW distance value, can be observed as the human-horse interaction changes.These trends can be also observed in Figure 5.In fact, almost all of the involved subjects showed a DTW value lower in P 1 than in P 2 and than in P 3 .A possible justification of this phenomenon could be found in the nature of the interaction.Probably, when human and horse are in physical contact, human and horse may suspiciously behave with a decrease of similarity and an increase of the distance between the two heart activities.The ability of the classification procedure of recognizing all of the three phases (i.e., accuracy greater than 85%) means that there exists a hyperplane able to differentiate these phases, whose distribution is not explained only by statistics, but it could be due to a non-linear dataset configuration in the feature space [112].As a matter of fact, the achieved results are not sufficient, at the moment, to generalize and confirm the described phenomena.Hence, an increased number of measures, including further variables, is needed to clarify the phenomena involved in the human-horse interaction.In conclusion, this work poses the basis for the application of novel high level signal processing techniques for stationary and non-stationary signals [113,114], already used in investigating human-to-human interaction [115][116][117], to animals, and particularly to horses, in order to objectively reveal interesting responses of both the central and autonomic nervous system (ANS) for a particular uncommon stimulation.Moreover, the achieved results lead us to conclude that a quantitative measure of the human-horse interaction is viable, and it could be very effective in many fields of application, for example in therapy assisted by equine (EAT) [2,118] or for controlling the effect of therapeutic horseback riding [119].In this field, our paradigm could permit analyzing the emotional interaction of a human patient with a "standardized horse" (a horse specifically trained and managed to create a controlled positive emotional background with a human).From a research point of view, these results can open a new scenario in which the emotional interaction [120] between human [121] and horse may be detectable and measurable, thus extending the current knowledge on comparative neuroscience.

Figure 1 .
Figure 1.Wearable monitoring system for humans[48].As is possible to note, from the box on the right, the e-textile electrodes for ECG (electrocardiographic) acquisition are knitted and completely integrated into the garment.

Figure 2 .
Figure 2. Systems placement on the horse.Elastic smart belt and Ag/AgCl electrodes placement.A scheme is presented on the left side and a real picture is on the right side.

Figure 3 .
Figure 3. Boxplot of the percentage of artifact signals observed in the two categories of data acquired.The central box represents the central 50% of the data.In fact, its lower and upper boundary lines are respectively at the 25% and 75% quantile of the data.A central red line indicates the median of the data.Finally, two vertical lines indicate the remaining data that are outside the central box and not considered as outliers.

Figure 4 .
Figure 4. Example of both human and equine heart rate variability (HRV) signals.Within the smaller portion of signal, it is shown how the warping-paths are defined.

Figure 5 .
Figure 5. DTW (dynamic time warping) values, as estimated between human and horse, obtained for each subject in each phase.

Table 1 .
Percentage of corrupted signal obtained by means of hand labeling of the two categories of signals: signals recorded by using Ag/AgCl electrodes and signals recorded by using textile electrodes.

Table 2 .
p-Values related to pairwise statistical analysis of the DTW (dynamic time warping) values, as estimated between human and horse, during the different phases (P 1 , P 2 , P 3 ).

Table 4 .
Confusion matrix obtained by applying the Support Vector Machine (SVM) classifier to the leave-one-subject-out architecture.: the percentage of the right recognised classes are marked in bold.