Next Article in Journal
Kinetic Sensors for Ligament Balance and Kinematic Evaluation in Anatomic Bi-Cruciate Stabilized Total Knee Arthroplasty
Next Article in Special Issue
Improving Animal Monitoring Using Small Unmanned Aircraft Systems (sUAS) and Deep Learning Networks
Previous Article in Journal
Statistical Characterization of Temperature and Pressure Vertical Profiles for the Analysis of Laser Heterodyne Radiometry Data
Previous Article in Special Issue
An Evaluation of the Factors Affecting ‘Poacher’ Detection with Drones and the Efficacy of Machine-Learning for Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Using Machine Learning for Remote Behaviour Classification—Verifying Acceleration Data to Infer Feeding Events in Free-Ranging Cheetahs

Leibniz-Institute for Zoo- and Wildlife Research, Alfred-Kowalke-Straße 17, 10315 Berlin, Germany
*
Author to whom correspondence should be addressed.
Contributed equally to supervising this work.
Sensors 2021, 21(16), 5426; https://doi.org/10.3390/s21165426
Submission received: 16 June 2021 / Revised: 1 August 2021 / Accepted: 5 August 2021 / Published: 11 August 2021
(This article belongs to the Special Issue Sensors and Artificial Intelligence for Wildlife Conservation)

Abstract

:
Behavioural studies of elusive wildlife species are challenging but important when they are threatened and involved in human-wildlife conflicts. Accelerometers (ACCs) and supervised machine learning algorithms (MLAs) are valuable tools to remotely determine behaviours. Here we used five captive cheetahs in Namibia to test the applicability of ACC data in identifying six behaviours by using six MLAs on data we ground-truthed by direct observations. We included two ensemble learning approaches and a probability threshold to improve prediction accuracy. We used the model to then identify the behaviours in four free-ranging cheetah males. Feeding behaviours identified by the model and matched with corresponding GPS clusters were verified with previously identified kill sites in the field. The MLAs and the two ensemble learning approaches in the captive cheetahs achieved precision (recall) ranging from 80.1% to 100.0% (87.3% to 99.2%) for resting, walking and trotting/running behaviour, from 74.4% to 81.6% (54.8% and 82.4%) for feeding behaviour and from 0.0% to 97.1% (0.0% and 56.2%) for drinking and grooming behaviour. The model application to the ACC data of the free-ranging cheetahs successfully identified all nine kill sites and 17 of the 18 feeding events of the two brother groups. We demonstrated that our behavioural model reliably detects feeding events of free-ranging cheetahs. This has useful applications for the determination of cheetah kill sites and helping to mitigate human-cheetah conflicts.

1. Introduction

One of the most complex aspects to study in animals is their behaviour [1]. Understanding species-specific behaviour is essential for their management, welfare and conservation [2,3,4,5,6]. Studying the behaviour of free-ranging animals often poses several challenges. Traditional ethology is based on direct observation, but the presence of an observer might alter the natural behaviour of the study animals if they are not fully habituated to the observer [7]. Furthermore, the animals might be difficult to observe when they are nocturnal, shy and cryptic [8], or are from an aquatic or migratory species [9].
Using biologgers has fueled research on remote monitoring of animal behaviour and unlocked the former limitations of direct observations [10,11,12]. Collecting data on fine-scale movement via accelerometers (ACCs) is an effective way to continuously monitor free-ranging animals [13,14]. ACCs measure the change in velocity of the sensors attached to the body and allow translation of the quantified movement patterns into distinct behavioural categories or activities [13]. This technique has been successfully used for a variety of bird species [15,16,17,18,19], marine animals [20,21,22,23], and terrestrial mammals [24,25,26,27,28], including cheetahs, Acinonyx jubatus [29,30,31,32].
The training of machine learning algorithms for pattern recognition and data classification is conducted with observational verified data. For this, a ground-truthed behavioural data set is used to train the algorithm and another data set is used to infer behaviour from it. Campbell and colleagues [30] used ground-truthed data of a tame surrogate dog (Canis lupus familiaris) to train their support vector machine (SVM) algorithm to predict walking, running, sitting, standing and lying in ACC data collected from different species, including the cheetah. Using the cheetah dataset generated by Campbell et al. [30], Bidder and colleagues [31] successfully predicted sitting and standing in cheetahs using another algorithm, the k-nearest neighbour (KNN) algorithm. While these two studies had strict validation purposes, Grünewälder et al. [29] observed six free-ranging cheetahs and recorded stationary, mobile and feeding behaviour during an average time of 31 h. A SVM and a hidden Markov model were trained with this ground-truthed dataset of each cheetah and used to predict the three behaviour categories for the remaining unobserved data collected during an average time of 332 days of the six cheetahs [29]. To our knowledge, there is currently no study that transferred a behaviour classification model trained on one group of cheetahs to another group of cheetahs. If the inferred behaviour classification in the latter group is verified, the classification model has a wide application for the species.
Cheetahs occur at low densities and cover large home ranges, making it difficult to observe them in the wild [33,34,35]. Southern Africa hosts the largest free-ranging cheetah population in the world, where cheetahs predominantly occur on farmland, i.e., outside of protected areas such as national parks or game reserves [36,37,38]. On farmland, cheetahs might come into conflict with landowners when they predate on their livestock and/or valuable game species [35,39]. Because cheetahs are not protected on farmland, landowners regularly kill cheetahs to prevent such losses [35,39]. The cheetah is listed as vulnerable by the International Union for Conservation of Nature and Natural Resources (IUCN) [36], and human–cheetah conflicts are amongst the major threats to its existence in the wild.
To investigate the extent of the farmer–cheetah conflict, it is important to understand the feeding ecology of the cheetahs, thus detecting feeding behaviours in the wild. Cheetahs feed very rarely on carcasses, thus feeding events can be interpreted as cheetahs having killed the prey animals [40]. Grünewälder et al. [29] used individual-based behaviour predictions from ACC data to investigate feeding-related behaviours. Other studies successfully used clusters of GPS positions to detect potential kill sites of carnivores, e.g., of leopards (Panthera pardus) [41,42] and jaguars (Panthera onca) [43]. Such clusters represent locations to which the animal returned repeatedly and/or spent an extended amount of time. From these spatial data alone, it is not possible to determine whether the animal was indeed feeding or rather resting. It is therefore required that GPS positions are retrieved and visualised in regular intervals and detected clusters visited within a time frame that allows finding of prey remains in the field [41,42]. Additional spatial and temporal information can be helpful in discriminating between potential kill sites and non-kill sites [43]. Combining both ACC and GPS data, Wang et al. [44] were able to identify and verify five predation events in free-ranging cougars (Puma concolor) by associating periods of high acceleration movements with subsequent GPS clusters. This was a useful approach because their behavioural model was trained with ground-truthed data of captive pumas and was weak at predicting feeding behaviour. Also, the combination of ACC and GPS data was successfully used to differentiate feeding from non-feeding events in a single leopard [45] and non-active (resting) from active (potentially feeding) stationary events in cougars [46].
The goal of this study was to detect six behaviours, including feeding behaviour, from acceleration data recorded in free-ranging cheetahs by using behaviour classification models trained on ACC data of captive cheetahs. For this, we first tested the potential of six common supervised machine learning algorithms (MLA) in inferring the behaviours from ACC data from captive cheetahs. We trained the MLAs on a dataset which we ground-truthed by direct observation. We then used the validated dataset to identify behaviour in ACC data from free-ranging cheetahs and, by adding GPS cluster analysis, detect feeding events. We verified feeding events identified by the model with kill sites that were previously confirmed by visiting GPS clusters in the field and identifying the prey remains.

2. Materials and Methods

2.1. Study Area and Animals

For the data validation, we used five captive cheetahs held in two enclosures on a private game farm in east-central Namibia (−22.5803° S, 18.1875° E). The enclosures were set within their natural environment. One enclosure (1.0 ha) contained three brothers of approximately 8 years of age and the other enclosure (1.3 ha) a brother and sister of approximately 4 years of age. All five animals were born in the wild but came into captivity at an early age. They were accustomed to human presence and allowed observers in their vicinity (Figure 1). The animals had ad libitum access to clean water in an artificial water trough and were fed once a day in the late afternoon. On 3 October 2017, we fitted GPS collars with ACC sensors (e-obs GmbH, Grünwald, Germany [47]) on the three brothers during immobilisations for medical check-ups and on the two other cheetahs without immobilisation since they accepted this handling by one of the authors (DB). The collars were colour-coded to help with the identification of the animals during the behavioural observations. After data collection, the collars were removed again with the same procedures.
To test the applicability of the behavioural model, we used data from four free-ranging territorial males in two coalition groups. They were equipped with the same collar types as the five captive cheetahs on 18 November 2014 and 7 December 2014 (coalition 1) and 21 May 2014 (coalition 2). Their GPS data were used to visually detect kill sites (by GPS clusters) and verify them by checking for prey remains in the field. Their ACC and GPS data were then used to detect feeding events correlating to the confirmed kill sites.
All experimental procedures described were approved by the Internal Ethics Committee of the Leibniz Institute for Zoo and Wildlife Research (IZW, Berlin, Germany) (permit number: 1 April 2002) and the Ministry of Environment, Forestry and Tourism of Namibia (permit numbers: 1514/2011, 1689/2012, 1813/2013), 1914/2014, 2067/2015, 2194/2016, 2208/2017).

2.2. Data Collection

2.2.1. Data Collection on Captive Cheetahs for Model Validation

Observations of the three brothers were conducted during 36 full or half days between the 6 October and the 10 December 2017, whereas the brother and sister were observed during 20 full or half days between the 20 November and the 19 December 2017. The e-obs collars were set to record ACC data from 4:00 until 17:00 Universal Time Coordinated (UTC), corresponding to the hours of daylight at the time of data collection to ensure good visibility for the behavioural observations. The ACC data were recorded in ‘bursts’ occurring every 30 s. The collars were equipped with two different tag firmware and the settings for the burst lengths and sampling frequencies were chosen from the configuration options provided by the manufacturer. We programmed the collars in such a way that the dataset collected with the two firmwares were as similar as possible and fitted best with the data collected from the free-ranging animals (Table 1). The sampling frequency for the three captive brothers was temporarily increased from 10 Hz to 33.3 Hz during 7 to 15 days to assess differences in model performance due to different sampling frequencies in another study. To also make use of the 33.3 Hz data in the present study, we resampled the data using the resample_poly function from the signal package within the scipy library (Version 1.4.1 [48]). The resampling was done in python 3.7.10 [49]. This function includes a low pass filter that is applied before down sampling. The collars did not have an analog anti-aliasing filter, but these low sampling frequencies were nevertheless chosen (see Section 4). We focused on the heave-axis (z-axis, Figure 1) because the collars of the free-ranging cheetahs only recorded ACC data from this axis (see Section 2.2.2).
All observations were carried out by two observers simultaneously. An ultra-high frequency (UHF) pinger with a unique frequency for each collars was running during the on-time period, emitting a rhythmical beeping signal which was picked up by a UHF receiver (AR8200, AOR, Tokyo, Japan). A change in the rhythm of the beeping signal indicated the start and the end of each ACC burst. During the burst, the behaviours of the animals were observed and recorded on a sheet. Using a digital radio clock, each recorded burst was linked to a unique timestamp. Since the collars of the cheetahs sharing the same enclosure were synchronised in their sampling time of ACC data, we recorded the behaviour of all individuals in sight during any given burst.
We recorded six behavioural categories during the observations: drinking (D), feeding (F), grooming (G), resting (R, including all non-motion behaviours), trotting/running (T) and walking (W). To simulate feeding events as naturally as possible, we used either big pieces of meat, bones and organs or entire carcasses obtained from animals hunted on the game farm. We aimed to record at least 30 bursts per behaviour and animal to ensure good trainability of the MLAs. The ACC data were downloaded manually from the collars with a UHF connection using an e-obs basestation (BaseStation b5) [47] at the end of each observation day.
For each behaviour, one ACC burst was randomly chosen to create sample plots of the raw data (Appendix A Figure A1).

2.2.2. Data Collection on Free-Ranging Cheetahs to Test Model Application

For the model application part, GPS and ACC data from four collared free-ranging cheetah males in two coalitions were used over 91 days from 6 February 2015 to 7 May 2015. The collars recorded ACC data of the heave-axis (z-axis, Figure 1) in bursts of 3.6 s every two minutes at a frequency of 10 Hz (Table 1). Due to battery lifetime constraints and fast, field-applicable data download, we focused on one axis. We decided on the heave-axis, because our main interest was to detect feeding behaviour and the ripping movements of the head during feeding are likely to be well display in this axis. In each coalition, GPS data were recorded every 15 min from both coalition members, while for one coalition member the schedule was increased temporarily to every three minutes. The GPS collars with the higher resolution, recorded fixes from the 22 January 2015 to the 16 April 2015 for coalition 1 and from the 6 February 2015 to the 25 March 2015 for coalition 2, respectively. This higher frequency of GPS data recording was used to facilitate a higher accuracy in determining a feeding event and the duration thereof. To save battery lifetime we used an ACC informed GPS recording scheme, such that GPS data were recorded only every 6 h when a predetermined threshold of ACC data variance between consecutive GPS fixes were not reached, implying the animal was inactive. As soon as the animal became active again, the GPS interval raised again to 15 min or 3 min, respectively. Due to the low battery lifetime of one member of coalition 2, we set the GPS interval from 3 min to 6 h for the last 12 days.
The GPS data were downloaded on a weekly basis during aerial tracking flights as described in Melzheimer et al. [34]. When initially analysing the GPS positions on a map using geographical information systems, distinct clusters were identified visually at places where the cheetahs spent an extended amount of time at the same geographical location. Due to the ACC informed GPS measurement, the collars only recorded consecutive GPS fixes at the same position if the animal was active. The dynamic GPS scheduling was therefore used as an additional clue to distinguish resting and feeding spots. Once such potential feeding spots were identified, the sites were visited within a week of the suspected kill and examined for prey remains to confirm that a prey animal was killed and eaten there. The duration of the feeding event was set based on the first and last GPS position taken at the specific location where the prey was consumed. In cheetahs, members of the same coalition are always together [34]. Thus, we assigned feeding events of members of the same coalition to a common kill site.

2.3. Data Analyses

2.3.1. Algorithms Used and Data Validation

We worked with supervised machine learning algorithms (MLAs) to analyse the ACC data [50]. Supervised MLAs are common tools in pattern recognition and are characterised by the two processes of training and testing. In short, we used six supervised MLAs, namely linear discriminant analysis (LDA) [51], quadratic discriminant analysis (QDA) [16], the KNN algorithm [31,52], the classification and regression tree (CART) [16,53] algorithm, SVM algorithm [54,55], and the random forest (RF) [56,57] algorithm. All six MLAs were applied on the same set of six predictors: mean of z-axis (mnz), standard deviation of z-axis (sdz), weighted mean of the Fourier transformation of the z-axis (wmz = weighted mean of the periodogram ordinates for all Fourier frequencies w = 2πj/n, j = 1,..., q with q = n/2 (n even) or q = n/2 − 1 (n odd), with n being the number of data points in the sample), inverse coefficient of variance of z-axis (ICVz = mnz / sdz), kurtosis of z-axis (kz) and skewness of z-axis (sz). They were calculated for each burst of raw ACC data. We also implemented two ensemble learning approaches, the mean and the majority voting, that are based on the predictions derived from the six algorithms to improve the classification performance of our behavioural model. More details are provided in the Appendix A Text A1.
The validation of ACC data was done using two different approaches. Firstly, we estimated the classification performance of our behavioural model with a classic leave-one-out cross-validation (LOOCV) using the complete dataset of the five captive cheetahs, i.e., the ground-truthed dataset. Secondly, we conducted a per-animal cross-validation (PACV) for each of the five cheetahs by training the behavioural model with the complete datasets of four animals to predict behaviour in the complete dataset of the fifth animal.
To further improve the performance of the model, we added a probability threshold. For this the model labelled a burst as a particular behaviour only when an algorithm classified this behaviour with the defined or a higher probability. Any burst below the probability threshold was labelled “not conclusive”. To determine the probability threshold with the highest gain in predictive performance in relation to the percentage of bursts that did not exceed the probability threshold, we performed the LOOCV with thresholds ranging from 0.3 to 0.8. For the PACV we only used the probability threshold that performed best. For the LOOCV, we assessed the performance of each algorithm and of the two ensemble approaches by calculating overall precision, precision per behaviour with PR = TP/(TP+FP) and recall per behaviour using the respective confusion matrices, with RE = TP/(TP+FN) (TP = true positive, FP = false positive, FN = false negative). For the PACV, we did the same but assessed the performance only for the two ensemble approaches since they performed best in the LOOCV (Table 2).

2.3.2. Model Application and Verification

To test the applicability of the behavioural model, we used ACC and GPS data of the four free-ranging cheetahs, covering one full day before and one full day after the date of each confirmed kill site (see Section 2.2.2) (that were previously identified by visually checking for GPS clusters and visiting these clusters in the field to check for prey remains). This time-frame was based on the average kill rates of free-ranging cheetahs, which is less than, or equal to, one kill every two days in southern Africa [29,58] and never exceeded a kill per day in the Serengeti National Park in Tanzania, East Africa [33]. Using this time frame, we limited the chances of including data from unconfirmed and unverifiable feeding events.
After calculating the six predictors for all raw ACC bursts of the free-ranging cheetahs, we trained the behavioural model with the previously established complete ground-truthed dataset to predict behaviour in this unknown data using the six algorithms, the two ensemble learning approaches and the probability threshold that performed best. The resulting file contained all eight behaviour predictions (from the six MLAs and the two ensemble learning approaches) for each burst. The unique timestamp of each burst allowed us to ‘read’ the predicted behaviours in chronological order.
Next, we searched for clusters of feeding bursts in the chronological sequence of all predicted behaviours, in the following termed as ‘feeding cluster’. Since both the ensemble learning approaches with probability threshold and the SVM algorithm with probability threshold performed best in predicting feeding behaviour in the LOOCV (Table 2), all bursts that were identified to be ‘feeding’ by either one of the ensemble learning approaches or the SVM algorithm were defined to be ‘true feeding’. We used a sliding window approach to find feeding cluster with at least a third of bursts per sliding window being ‘true feeding’. We determined start and end time of these feeding clusters allowing a minimum length of 30 min per cluster to cover medium to large kills following Grünewälder et al. [29]. They determined for cheetahs a feeding time of 10 to 15 min for small prey animals and one of up to 3 h for large prey animals. We allowed for a maximum gap of 120 min of ‘not feeding’ within a cluster.
We then screened the corresponding GPS data for clusters. GPS clusters were determined based on the condition that consecutive coordinates did not exceed a distance of 50 m to each other and lasted at least 30 min. The centroid of each GPS cluster was determined by calculating the mean longitude and latitude of all coordinates of the cluster.
The GPS clusters and feeding clusters were then checked for overlaps in time. If multiple GPS clusters fell into the time frame of one feeding cluster, we either merged these GPS clusters, if their centroids were less than 50 m apart from each other, or we determined the GPS cluster with the highest time overlap with the feeding cluster to match a feeding event. If this process was inconclusive, we examined the data manually and assigned GPS coordinates to the ‘unmatched’ feeding cluster or dropped the feeding cluster as ‘false positive’. If feeding clusters started/ended more than 3 or 15 min, respectively (i.e., the interval between consecutive GPS fixes) before/after the animal was stationary according to a GPS cluster, start and end time of the feeding cluster were corrected accordingly.
Since free-ranging cheetahs kill their prey themselves and only rarely scavenge [40], they typically engage in hunting activity prior to a feeding event. We therefore also scanned the behaviour predictions one hour prior to the feeding events for bursts identified as ‘high ACC variation’, i.e., trotting and running, indicating a potential hunt.
All feeding events determined this way were then descriptively compared to the previously confirmed nine kill sites.
All analyses were conducted using the statistical computing software R Version 3.5.1 [59].

3. Results

3.1. Data Validation

For the five captive cheetahs, a total of 7760 bursts were recorded, of which 6,641 were single-behaviour bursts. The remaining 1119 bursts were discarded as they contained mixed behaviours, i.e., a change from one behaviour to another within a burst. The most sampled behaviour was resting with 2717 bursts (35.0%), followed by walking with 1673 bursts (21.6%), feeding with 1319 bursts (16.9%), and grooming with 487 bursts (6.3%). Trotting/running and drinking were the least sampled behaviours with 260 (3.4%) and 185 (2.4%) bursts, respectively.
The LOOCV was run with thresholds from 0.3 to 0.8 to find the most effective one. The probability threshold of 0.5 had the best ratio in increasing the precision and at the same time decreasing the recall (Appendix A Table A1). Confusion matrices for all algorithms with the 0.5 probability threshold are in the appendix (Appendix A Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8 and Table A9). Adding this threshold to the LOOCV increased the performance of the model, particularly in predicting feeding behaviour by the SVM algorithm (performance increase 2.7%) by simultaneously only missing 5.4% of true feeding bursts by labeling them as “not conclusive” (Appendix A Table A1). Thus, the SVM algorithm played an essential role in the determination of ‘true feeding’ in the model application.
With the threshold of 0.5, the LOCCV achieved overall precisions ranging from 57.5% to 91.4%, depending on the algorithm and ensemble learning method used (Table 2). The algorithms achieved high accuracies for resting, walking and trotting/running behaviour, but performed less well for drinking and grooming behaviour (Table 2). The combination of precision and recall for feeding behaviour was best at (78.8%/82.4% (precision/recall)) with the SVM algorithm (Table 2). The two ensemble learning approaches, both performed better than the SVM in respect to the precision but worst in respect to the recall (Table 2). The highest precisions for behaviours were 99.0% and 100% for resting and trotting/running at a recall of 90.7% and 99.2%, respectively. A comparison of all values in Table 2 with the results for the LOOCV without adding a probability threshold is presented in Appendix A Table A10.
The PACV with a probability threshold of 0.5 also reached high overall precision, ranging from 69.9% to 92.5% for the two ensemble learning approaches, while the overall recall, ranging between 57.5% and 70.0%, did not reach such high values (Appendix A Table A11). The precision for feeding behaviour ranged from 65.6% to 91.4% and the recall from 47.0% to 81.1%.

3.2. Model Application

The visual inspection of the GPS data of the two free-ranging male coalitions for clusters and the checking of these spots in the field resulted in nine GPS clusters that we confirmed to be kill sites of the cheetah males because we found in all cases prey remains at the sites.
We successfully identified in the unknown ACC data set of the two coalitions all nine kill sites and 17 of the 18 feeding events (Table 3). The time ranges of these feeding events fitted into the timespans of the feeding events deduced from the confirmed kill sites. Differences in the durations of GPS clusters between the feeding events detected visually and the model application was based on the different approaches of determining the GPS clusters, i.e., visual vs. automated approach (Table 3). We identified bursts with high ACC variation in the hour before the respective feeding event in eight out of 17 cases (47.1%) and in seven out of nine kill sites (77.8%, Table 3). We also identified one feeding event that was not confirmed in the field with prey remains, although a GPS cluster was identified (Table 3, Figure 2).
Figure 2 illustrates nine days of GPS coordinates of two males of coalition 2 with GPS intervals of 3 min and 15 min, respectively. Both males revealed almost an identical movement path, and the same feeding and resting clusters were detected in the corresponding ACC data. While two feeding events (marked red in Figure 2) were previously confirmed in the field with prey remains, the third one (marked orange in Figure 2) was not confirmed in the field, but was identified in both males. Walking events were only detected in the male with the 3-min GPS schedule. All walking clusters corresponded with lacking stationarity, i.e., no GPS clusters, indicating that the animal was indeed on the move. The calculation of resting and walking clusters is described in Appendix A Text A2.

4. Discussion

In this study, we validated a method to remotely and automatically detect with MLAs feeding events/kill sites of free-ranging cheetahs by using ACC and GPS data (Figure 3). We used a two-step process by demonstrating in a first step the applicability of ACC data to predict basic behaviours in five captive cheetahs by using six supervised MLAs. To improve prediction accuracy, we introduced two ensemble learning methods and a probability threshold. In a second step, we successfully used the trained model to predict behaviour in free-ranging cheetahs with the specific intention to detect feeding behaviour and, by adding GPS cluster analysis, identify feeding events. Such feeding events were verified by previously confirmed kill sites in the field (Figure 3).

4.1. Data Validation

Observing captive animals inherently goes along with a reduced range of behaviour. In our study, the missing behaviours were mainly hunting, killing, fighting and mating. Nevertheless, our ground-truthed dataset mirrored the natural activity budget of cheetahs, with resting (and observing), walking (moving) and feeding (hunting and eating) making up by far the largest proportion of their daily activities [33,60].
While the six MLAs (LDA, QDA, KNN, CART, SVM, RF) differed in their success in predicting each behaviour, we found the two ensemble learning approaches to be the most promising ones to overall improve precision and recall and a good compromise for the successful identification of the six behaviours. On the other hand, the performance differences of the MLAs allow for a specific use of the best performing algorithm when concentrating on a particular behaviour, as we did in the model application: We found the SVM algorithm to work best in identifying feeding behaviour and used it to determine instances of ‘true feeding’. To account for unverified and mixed behaviour bursts, we implemented a probability threshold of 0.5 [31,61]. This further improved precision and recall.
The results of the LOOCV and the PACV with overall precisions ranging from 57.5% to 91.4% and 69.9% to 92.5%, respectively, are similar to other studies that examined cheetah behaviour predictions: Campbell et al. [30] achieved with the SVM algorithm precision and recall values of >90% for correctly classifying sitting and standing behaviour in tri-axial ACC data of one cheetah. Using the cheetah dataset provided by Campbell et al. [30], Bidder et al. [31] reached precision and recall values of 90% and 97% using the KNN algorithm to predict sitting and standing behaviour. Grünewälder et al. [29] predicted stationary, mobile and feeding behaviour in bi-axial ACC data from six free-ranging cheetahs, for which both training and testing data originated from the same individuals. They achieved overall accuracy between 83.9% and 94.0%, while feeding behaviour had a larger range and was correctly identified in 22.6% to 100.0% of cases [29]. Our precision for feeding behaviour with the uni-axial ACC data using the LOOCV were also high and ranged between 74.4% and 81.6% at a recall between 54.8% and 82.4%. The PACV, i.e., a model trained with n-1 animals and applied to the remaining animal for which the model was not trained, generally leads to a reduction in the prediction accuracy [62,63]. Nevertheless, our precision for the prediction of feeding behaviour were between 65.6% and 91.4% at a recall between 47.0% and 81.1%. Thus, performed in the upper range compared with the results of Grünewälder et al. [29].
Behaviours that were not very variable in their movements and lasted relatively long such as resting, walking and trotting (Appendix A Figure A1) reached highest values for precision and recall, whereas more variable and shorter lasting behaviours such as grooming and drinking (Appendix A Figure A1) performed least. The latter two behaviours are highly flexible concerning body postures and duration, thus inter-individual differences are expected to be most pronounced for these behaviours [61]. Feeding behaviour reached intermediate to high precision and recall. Previous studies on other species also had low predictive performance for behaviours that had either structural similarities in movements and/or a high variance within each behaviour (cattle [54], cougar [44], African elephant (Loxodonta Africana) [25], roe deer (Capreolus capreolus) [28]).
A reduction in the prediction accuracy of inter-individual models such as LOOCV and PACV might also be caused by different sensor tags, variations in individual-specific behaviours, sex or changing environmental conditions [28]. The female cheetah had lower prediction accuracy for walking behaviour and a lower overall prediction accuracy compared to the males. This might be due to having predicted her behaviours by a model trained with only male data or having worn a tag firmware with slightly different ACC setting (10.54 Hz frequency, 3.8 s burst length) from the ones worn by the males (10.0 Hz frequency, 4.0 s burst length). In addition, this female had a higher body mass index and shorter legs than the males (pers. observation), which might have affected gait and/or head movements and amplitudes.

4.2. Model Application

We applied our behavioural model to data from free-ranging cheetah. Similar to other studies combining GPS cluster analysis with ACC data [44,45,46], we used the GPS data corresponding to particular behaviour in the ACC data as additional information to support the model predictions. We used the occurrence of GPS clusters to confirm the predicted feeding events because cheetahs that feed on a typical prey animal of 14 to 56 kg [64] are stationary for approximately one to three hours, resulting in GPS clusters [40].
We successfully identified all nine kill sites and 17 of 18 feeding events of the two male coalitions, which we previously verified in the field by using GPS clusters and checking for prey remains. In addition, we detected a feeding event in both members of one coalition that was not confirmed in the field (Figure 2). Perhaps there was a feeding event and kill site in the field, but the prey remains were carried away by other sympatric carnivores such as leopards, brown hyenas (Parahyaena brunnea) or black-backed jackals (Canis mesomelas).
There might be several reasons for a reduced model performance on free-living conspecifics such as behaviours that do not occur in captivity and thus are not in the training data set or recordings with mixed behaviours within one burst. Nevertheless, we assess our model using the SVM algorithm and the two ensemble learning approaches as sufficiently successful to detect feeding events in free-ranging cheetahs.
The detection of the chasing part of a hunt before a feeding event would increase the reliability of the identified feeding event [44]. Such chasing behaviour, which in our model was represented by the trotting/running behaviour, would be indicated by bursts with high ACC data variation. The chances to identify bursts during chasing behaviour are limited since bursts only cover 3.0% of the time (3.6 s per 120 s), and the time for a chase only lasts on average 37.9 s [65]. Nevertheless, we detected bursts of high ACC variation in 47.1% of the feeding events in the hour before the feeding event started, which represented 77.8% of the kill sites. To avoid missing bursts during chasing time periods, continuous ACC data recordings or the implementation of additional sensors might be useful. For example, Hetem et al. [66] used a combination of body temperature and activity patterns to identify hunting in cheetahs.
Since there are trade-offs between ACC sampling rate, behaviour predictability, memory usage, battery capacity, and data download duration, these parameters need to be considered in the specific context the study is performed in. Setting a low sampling frequency might result in aliasing effects when patterns of higher frequencies occur. Our results on the data validation and the model application using ACC data at approximately 10 Hz sampling frequency at only one axis, however, are particularly promising and important for future applications on free-ranging cheetahs when low sampling frequencies are necessary to ensure long battery and storage lifetime. Similarly, behaviour model accuracies in free-ranging cougars were not affected significantly when decreasing sampling frequency from 65 Hz to 16 Hz, and only a 2% reduction in model accuracy was noted when down-sampling further to 8 Hz [44]. A higher sampling frequency would have disadvantageous practical implications in the field in that the increased data volume would reduce the time period to reach the memory capacity. It would also increase the time needed to download the data from the collars, which is typically done during costly aerial tracking flights. Additionally, our study demonstrated that the use of only one ACC axis was sufficient to reliably identify feeding, resting, trotting/running and walking.

4.3. Method Application

Understanding the behaviour of carnivore species is pivotal to their management [67]. Thus, including accelerometry in ecological studies is providing novel insights into carnivore behaviour, with conservation implications. For example, a study on African lions (Panthera leo) in a human-dominated landscape in Kenya demonstrated that lions adopted flexible feeding patterns to avoid conflicts with humans, enabling co-existence [68]. Another study showed how human housing density negatively affected the movement and energetics of cougars in the USA, with cougars increasing their kill rate to compensate for anthropogenic disturbance [69]. In the context of human–wildlife conflict, collaring and tracking of conflict species already provided successful conflict mitigation solutions such as adapting livestock management to avoid livestock grazing in high activity areas of carnivores [35] and using geofencing for early warning systems that facilitate quick responses of livestock herders or owners to an approaching carnivore [70]. With human–carnivore conflicts posing one of the greatest challenges in conservation [71], it is essential to develop tools that help improve our understanding of carnivore behaviour and can be applied in conservation to safeguard their future.

5. Conclusions

In this study, we demonstrated that MLAs are a useful approach that provide detailed information on the behaviour, particularly the feeding behaviour of cheetahs. We demonstrated the strong potential and practical applicability of ACC data and MLAs for continuous, automated, and high-resolution behaviour monitoring of cheetahs and showed that their feeding behaviour is reliably detectable. Such information can be used to give new insights into the human–wildlife conflict in Namibia and elsewhere in the cheetah range.

Author Contributions

Conceptualization, J.M., B.W. (Bettina Wachter) and A.B.; methodology, L.G., D.B., J.M., B.W. (Bettina Wachter) and A.B.; software, L.G., W.R. and A.B.; validation, L.G. and D.B.; formal analysis, L.G., W.R. and D.B.; investigation, L.G., D.B. and B.W. (Bernd Wasiolka); resources, B.W. (Bettina Wachter); data curation, L.G., D.B. and B.W. (Bernd Wasiolka); writing—original draft preparation, L.G., B.W. (Bettina Wachter) and D.B.; writing—review and editing, L.G., B.W. (Bettina Wachter), D.B., B.W. (Bernd Wasiolka), W.R., J.M. and A.B.; visualization, L.G. and D.B.; supervision, J.M., B.W. (Bettina Wachter) and A.B.; project administration, B.W. (Bettina Wachter); funding acquisition, B.W. (Bettina Wachter). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. The publication has been partly funded by the Leibniz Open Access Publishing.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of the Leibniz Institute for Zoo and Wildlife Research (IZW, Berlin, Germany) (permit number: 2002-04-01) and the Ministry of Environment, Forestry and Tourism of Namibia (permit numbers: 1514/2011, 1689/2012, 1813/2013), 1914/2014, 2067/2015, 2194/2016, 2208/2017).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to conservation status of the species.

Acknowledgments

We thank the Ministry of Environment, Forestry and Tourism (MEFT) in Namibia for permission and support to conduct this study and Christian Schmitt and family from the Okambara Elephant Lodge for allowing us to conduct the study on the captive cheetahs on their farm. We also thank Heribert Hofer for stimulating discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Sample raw acceleration (ACC) plots of one burst for the six behaviours (a) drinking (D), (b) feeding (F), (c) grooming (G), (d) resting (R), (e) trotting/running (T), and (f) walking (W) with corresponding pictograms.
Figure A1. Sample raw acceleration (ACC) plots of one burst for the six behaviours (a) drinking (D), (b) feeding (F), (c) grooming (G), (d) resting (R), (e) trotting/running (T), and (f) walking (W) with corresponding pictograms.
Sensors 21 05426 g0a1
Table A1. Output of the leave-one-out cross-validations (LOOCV) with probability thresholds for behaviour classification ranging from 0.3 to 0.8. Percentage points of performance change were determined by calculating the difference between precision and recall reached without and with the threshold. Percentage points of bursts below the probability threshold for each behaviour were calculated by dividing the number of bursts below the threshold by the total number of observed counts. Performance change and bursts below the threshold are presented as the range over all six algorithms for all ACC data (i.e., linear discriminant analysis, quadratic discriminant analysis, the k-nearest neighbour algorithm, the classification and regression tree algorithm, support vector machine algorithm and the Support Vector Machine (SVM), the range over all six algorithms for only feeding (F) data and the result for the SVM algorithm for only feeding data. The probability threshold with the best trade-off between precision increase and recall decrease and the percentage of bursts that did not exceed the probability threshold was used for the analyses. This was assessed to be the probability threshold 0.5 (bold).
Table A1. Output of the leave-one-out cross-validations (LOOCV) with probability thresholds for behaviour classification ranging from 0.3 to 0.8. Percentage points of performance change were determined by calculating the difference between precision and recall reached without and with the threshold. Percentage points of bursts below the probability threshold for each behaviour were calculated by dividing the number of bursts below the threshold by the total number of observed counts. Performance change and bursts below the threshold are presented as the range over all six algorithms for all ACC data (i.e., linear discriminant analysis, quadratic discriminant analysis, the k-nearest neighbour algorithm, the classification and regression tree algorithm, support vector machine algorithm and the Support Vector Machine (SVM), the range over all six algorithms for only feeding (F) data and the result for the SVM algorithm for only feeding data. The probability threshold with the best trade-off between precision increase and recall decrease and the percentage of bursts that did not exceed the probability threshold was used for the analyses. This was assessed to be the probability threshold 0.5 (bold).
Probability ThresholdMetricRange All Algorithms, All ACC DataRange All Algorithms, only F DataSVM Algorithm, only F Data
0.3Precision0.1–0.8−0.2–0.90.1
Recall−0.3–0.1−0.2–0.80.0
Below threshold bursts 0.0–1.10.0–0.80.1
0.4Precision0.0–4.40.0–6.10.5
Recall−2.2–0.0−3.3–0.0−0.7
Below threshold bursts 0.0–6.30.0–5.51.1
0.5Precision0.0–4.00.0–13.02.7
Recall−14.7–−0.2−14.2–0.0−2.3
Below threshold bursts2.4–15.44.2–24.55.4
0.6Precision0.0–7.70.0–19.85.1
Recall−21.7–−0.4−27.7–0.0−9.5
Below threshold bursts 5.8–24.56.4–45.616.8
0.7Precision0.0–11.50.0–24.514.4
Recall−27.3–−0.8−40–0−32.6
Below threshold bursts 5.8–33.76.4–64.328.0
0.8Precision0.0–7.7−74.5–28.117.7
Recall−43.8–−1.2−75.4–−25.3−43.5
Below threshold bursts 17.7–47.237.4–81.844.3
By adding the probability threshold ranging from 0.3 to 0.8 to the model, precision increased for most behaviours. Drinking behaviour was never classified correctly by the LDA and CART which is why calculating changes in all algorithms was not possible resulting in 0.0 change. There was also an improvement for precision for feeding behaviour, except for the ones with a threshold of 0.3 and 0.8. The SVM algorithm performed best in correctly identifying feeding behaviour and therefore played an essential role in the determination of ‘true feeding’ in the model application. The probability threshold of 0.5 best trade-off between precision increase and recall decrease and the percentage of bursts that did not exceed the probability threshold, particularly in predicting feeding behaviour by the SVM algorithm (performance increase 2.7%) by simultaneously only missing 5.4% of true feeding bursts by labeling them as “not conclusive”. While higher thresholds showed higher precision improvements, we assessed the lower recall and percentage of feeding bursts below the probability threshold to be too high to compensate for the improvement. As a useful compromise, we decided to conduct all analyses with a probability threshold of 0.5.
Table A2. Confusion matrix for the LDA algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
Table A2. Confusion matrix for the LDA algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
LDAExpected
Predicted RestingWalkingTrottingDrinkingFeedingGrooming
resting237210135385
walking121583287321719
trotting00232000
drinking000000
feeding133007723113
grooming7300333
unknown31356092323237
Table A3. Confusion matrix for the QDA algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
Table A3. Confusion matrix for the QDA algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
QDAExpected
Predicted RestingWalkingTrottingDrinkingFeedingGrooming
resting2464001169
walking1158802720523
trotting00258060
drinking79401042537
feeding2766117941150
grooming12461780217
unknown2290195651
Table A4. Confusion matrix for the KNN algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
Table A4. Confusion matrix for the KNN algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
KNNExpected
Predicted RestingWalkingTrottingDrinkingFeedingGrooming
resting257540254068
walking3156705011310
trotting00258000
drinking4002310
feeding1466016981127
grooming32112397
unknown11834170161185
Table A5. Confusion matrix for the CART algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
Table A5. Confusion matrix for the CART algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
CARTExpected
Predicted RestingWalkingTrottingDrinkingFeedingGrooming
resting253310775189
walking72153817718734
trotting01254022
drinking000000
feeding25116320995177
grooming100000
unknown861721184185
Table A6. Confusion matrix for the RF algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
Table A6. Confusion matrix for the RF algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
RFExpected
Predicted RestingWalkingTrottingDrinkingFeedingGrooming
resting263000271524
walking11527014782
trotting01256020
drinking5307025
feeding7990151027129
grooming1250963216
unknown6238450132111
Table A7. Confusion matrix for the SVM algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
Table A7. Confusion matrix for the SVM algorithm with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
SVMExpected
Predicted RestingWalkingTrottingDrinkingFeedingGrooming
resting265300322635
walking11542022924
trotting00258000
drinking33068210
feeding8960161087173
grooming1431741188
unknown38291407177
Table A8. Confusion matrix for the ensemble learning method “mean” with a probabillity threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
Table A8. Confusion matrix for the ensemble learning method “mean” with a probabillity threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
MeanExpected
Predicted RestingWalkingTrottingDrinkingFeedingGrooming
resting263100211125
walking115640291034
trotting00255010
drinking1003300
feeding760013988144
grooming650019125
unknown7144589197189
Table A9. Confusion matrix for the ensemble learning method “majority” with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
Table A9. Confusion matrix for the ensemble learning method “majority” with a probability threshold of 0.5. The diagonal shows all correctly predicted behaviour classes. All predictions that were below the probability threshold were labeled as unknown.
MajorityExpected
Predicted RestingWalkingTrottingDrinkingFeedingGrooming
resting265000302839
walking215720391249
trotting01258000
drinking4104913
feeding9791171047166
grooming841226141
unknown441604893129
Table A10. Results of the leave-one-out cross-validations (LOOCV) without (no) and with (yes) a probability threshold (PT) of 0.5 for behaviour classification for each of the six machine learning algorithms (MLA), i.e., linear discriminant analysis (LDA), quadratic discriminant analysis (QDA), the the k-nearest neighbour (KNN) algorithm, the classification and regression tree (CART) algorithm, support vector machine (SVM) algorithm and the random forest (RF) algorithm, and the two ensemble-learning approaches, i.e., the majority and mean voting. Precision (PR) across all acceleration data and for the behaviours drinking (D), feeding (F), grooming (G), resting (R), trotting/running (T) and walking (W), and recall (RE) for all behaviours are shown. PR and RE were calculated using the respective confusion matrices, with PR = TP/(TP+FP) and RE = TP/(TP+FN). TP = true positive, FP = false positive, FN = false negative.
Table A10. Results of the leave-one-out cross-validations (LOOCV) without (no) and with (yes) a probability threshold (PT) of 0.5 for behaviour classification for each of the six machine learning algorithms (MLA), i.e., linear discriminant analysis (LDA), quadratic discriminant analysis (QDA), the the k-nearest neighbour (KNN) algorithm, the classification and regression tree (CART) algorithm, support vector machine (SVM) algorithm and the random forest (RF) algorithm, and the two ensemble-learning approaches, i.e., the majority and mean voting. Precision (PR) across all acceleration data and for the behaviours drinking (D), feeding (F), grooming (G), resting (R), trotting/running (T) and walking (W), and recall (RE) for all behaviours are shown. PR and RE were calculated using the respective confusion matrices, with PR = TP/(TP+FP) and RE = TP/(TP+FN). TP = true positive, FP = false positive, FN = false negative.
MLAPTOverallDFGRTW
PRREPRREPRREPRREPRREPRREPRRE
LDAno63.160.8<0.010.068.669.043.618.992.291.7100.089.274.296.2
yes71.555.5 <0.10.0 81.654.8 71.76.8 94.087.3 10089.2 82.094.6
QDAno74.978.240.460.577.173.349.749.998.990.997.799.285.695.1
yes75.576.2 41.856.2 78.271.3 50.044.6 99.090.7 97.899.2 86.194.9
KNNno81.872.168.126.576.280.367.234.592.597.1100.099.287.094.7
yes87.565.7 12.4 12.4 78.3 74.4 76.419.9 94.994.8 10099.2 95.093.7
CARTno65.566.0<0.010.074.575.447.938.092.1 93.298.197.780.691.9
yes57.559.7 <0.10.074.475.4 0.00.092.193.2 98.097.7 80.191.9
RFno83.079.169.446.577.383.062.855.996.497.898.599.293.492.1
yes87.474.4 82.437.8 80.477.9 70.844.4 97.696.8 98.898.5 94.191.3
SVMno83.977.572.244.976.084.767.745.296.098.2100.099.291.792.8
yes86.974.5 79.136.8 78.882.4 74.038.6 96.697.6 10099.2 92.892.2
Meanno84.175.573.136.876.882.371.241.195.7 98.599.299.288.595.1
yes91.467.8 97.117.8 81.574.9 80.625.7 97.996.8 99.698.1 91.993.5
Majorityno84.575.676.136.277.182.370.842.995.798.599.299.288.394.7
yes87.970.9 84.426.5 79.479.4 77.529.0 96.597.5 99.699.2 90.094.0
Table A11. Results of the per-animal cross-validations (PACV) with a probability threshold of 0.5 for the five cheetahs (M1, M2, M3, M4 and F1). Precision (PR) of the two ensemble-learning approaches (ELA), mean and majority voting, across all acceleration data and for the behaviours drinking (D), feeding (F), grooming (G), resting (R), trotting/running (T) and walking (W), and recall (RE) for all behaviours. PR and RE were calculated using the respective confusion matrices, with PR = TP/(TP + FP) and RE = TP/(TP + FN). TP = true positive, FP = false positive, FN = false negative.
Table A11. Results of the per-animal cross-validations (PACV) with a probability threshold of 0.5 for the five cheetahs (M1, M2, M3, M4 and F1). Precision (PR) of the two ensemble-learning approaches (ELA), mean and majority voting, across all acceleration data and for the behaviours drinking (D), feeding (F), grooming (G), resting (R), trotting/running (T) and walking (W), and recall (RE) for all behaviours. PR and RE were calculated using the respective confusion matrices, with PR = TP/(TP + FP) and RE = TP/(TP + FN). TP = true positive, FP = false positive, FN = false negative.
Test DatasetELAOverallDFGRTW
PRREPRREPRREPRREPRREPRREPRRE
M1Mean80.359.433.32.586.347.069.439.197.895.597.486.097.486.5
Majority77.766.652.922.580.349.749.149.190.798.597.590.795.789.2
M2Mean74.160.5<0.10.082.772.470.06.399.794.3100.097.892.492.6
Majority82.664.575.06.080.575.557.114.497.697.897.897.887.395.3
M3Mean80.361.5<0.10.0 91.465.6100.017.299.493.4100.0100.090.892.5
Majority92.570.0100.01.684.481.191.244.897.898.2100.0100.081.594.1
M4Mean76.161.5<0.10.072.778.4100.05.797.797.4100.0100.086.087.8
Majority73.065.2<0.10.072.179.994.720.593.699.0100.0100.077.691.8
F1Mean75.257.7<0.10.069.079.5100.01.6 98.780.8100.0100.083.684.3
Majority69.962.40.00.065.680.877.811.398.191.8100.0100.078.290.4
Discussion of Table A11: The PACV with a probability threshold of 0.5 reached, similar to POCCV, high overall precision but lower overall recall. They ranged from 69.9% to 92.5% and from 57.5% and 70.0% for precision and recall respectively for the two ensemble learning approaches. Resting and trotting/running behaviours were predicted with very high precision and recall of up to 100%. Feeding behaviour reached precision and recall ranging from 65.6% to 91.4% and from 47.0% to 81.1% respectively. Precision for feeding was lowest in the female F1 while the M1 showed the lowest recall for feeding behaviour. Drinking and grooming showed low recall scores in all cheetahs. High precision values result from a very low rate of these behaviours being assigned to a burst at all.
Text A1: Details on data analyses
Supervised machine learning algorithms (MLAs) use feature characteristics of predictor variables calculated from a ground-truthed dataset (training dataset) with different states/behaviours, to construct a model (=training process). The model is then used for predicting the different states on another ground-truthed dataset (=testing or validating process) or on data with unknown states/behaviours (=application). References and R packages for the six supervised MLAs were as follows: linear discriminant analysis (LDA) [51], quadratic discriminant analysis (QDA) [16], the KNN algorithm [31] using the R package kknn v. 1.3.1 [52], the classification and regression tree (CART) algorithm [16] using the package rpart v. 4.1-13 [53], SVM algorithm [54] using the package e1071 v. 1.6-8 [55], and the random forest (RF) algorithm [56] using package randomForest v. 4.6-14 [57].
The two ensemble-learning approaches were used to improve performance classification since the six MLAs showed differences in prediction accuracy per behaviour. For the mean voting, the burst gets assigned the behaviour that has the highest mean probability of all six MLAs. For the majority voting, each burst gets assigned the behaviour which is classified most by the six MLAs. In case of no majority, i.e., a draw, the behaviour gets assigned randomly to the most classified ones.
Text A2: Calculation of resting and walking clusters
The detection of resting and walking clusters was conducted very similarly to the procedure described for the detection of feeding clusters, again using the chronological order of behaviour predictions paired with the calculated GPS clusters (see Section 2.3). We started by determining all bursts that were identified as resting/walking by either one or both of the ensemble approaches to be bursts of ‘true resting’/’true walking’. We then again used a sliding window approach to find clusters of ‘true resting’/’true walking’ instances on the condition that the majority of bursts (at least 12 of 15) per sliding window were voted to be ‘true resting’/’true walking’. We then determined start and end time of these resting/walking clusters. Detected clusters that were at least 30 min long were then paired with their respective GPS coordinates via the timestamp. In case of the resting clusters we again assume overlaps with the already determined GPS clusters, while walking clusters are expected to be associated with missing GPS clusters, e.g., at least 2-3 coordinates that are more than 50 m apart from each other.

References

  1. Wilson, E.O. Foreword. The Behavior Guide to African Mammals, 20th Anniversary ed.; University of California Press: Berkeley, CA, USA, 2012. [Google Scholar]
  2. Sutherland, W.J. The importance of behavioral studies in conservation biology. Anim. Behav. 1998, 56, 801–809. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Blumstein, D.T.; Fernandez-Juricic, E. The Emergence of Conservation Behavior. Conserv. Biol. 2004, 18, 1175–1177. [Google Scholar] [CrossRef]
  4. Buchholz, R. Behavioral biology: An effective and relevant conservation tool. Trends Ecol. Evol. 2007, 22, 401–407. [Google Scholar] [CrossRef] [PubMed]
  5. Swaisgood, R.R. Current status and future directions of applied behavioral research for animal welfare and conservation. Appl. Anim. Behav. Sci. 2007, 102, 139–162. [Google Scholar] [CrossRef]
  6. Joubert, S.C.J. Animal behaviour. In Game Ranch Management, 5th ed.; Bothma, J.d.P., du Toit, J.G., Eds.; Van Schaik Publishers: Pretoria, South Africa, 2014; pp. 384–390. [Google Scholar]
  7. Crofoot, M.C.; Lambert, T.D.; Kays, R.; Wikelski, M.C. Does watching a monkey change its behaviour? Quantifying observer effects in habituated wild primates using automated radiotelemetry. Anim. Behav. 2010, 80, 475–480. [Google Scholar] [CrossRef] [Green Version]
  8. Sliwa, A. Home range size and social organisation of black-footed cats (Felis nigripes). Mammal. Biol. 2004, 69, 96–107. [Google Scholar] [CrossRef]
  9. Block, B.A.; Jonsen, I.D.; Jorgensen, S.J.; Winship, A.J.; Shaffer, S.A.; Bograd, S.J.; Hazen, E.L.; Foley, D.G.; Breed, G.A.; Harrison, A.-L.; et al. Tracking apex marine predator movements in a dynamic ocean. Nature 2011, 475, 86. [Google Scholar] [CrossRef] [PubMed]
  10. Cagnacci, F.; Boitani, L.; Powell, R.A.; Boyce, M.S. Animal ecology meets GPS-based radiotelemetry: A perfect storm of opportunities and challenges. Phil. Trans. R. Soc. B 2010, 365, 2157–2162. [Google Scholar] [CrossRef] [Green Version]
  11. Wilson, A.D.M.; Wikelski, M.; Wilson, R.P.; Cooke, S.J. Utility of biological sensor tags in animal conservation. Conserv. Biol. 2015, 29, 1065–1075. [Google Scholar] [CrossRef] [Green Version]
  12. Kays, R.; Crofoot, M.C.; Jetz, W.; Wikelski, M. Terrestrial animal tracking as an eye on life and planet. Science 2015, 348, aaa2478. [Google Scholar] [CrossRef] [Green Version]
  13. Brown, D.D.; Kays, R.; Wikelski, M.; Wilson, R.; Klimley, A.P. Observing the unwatchable through acceleration logging of animal behaviour. Anim. Biotelem. 2013, 1, 20. [Google Scholar] [CrossRef] [Green Version]
  14. Shepard, E.L.C.; Wilson, R.P.; Quintana, F.; Gómez-Laich, A.; Liebsch, N.; Albareda, D.A.; Halsey, L.G.; Gleiss, A.; Morgan, D.T.; Myers, A.E.; et al. Identification of animal movement patterns using tri-axial accelerometry. Endanger. Species Res. 2010, 10, 47–60. [Google Scholar] [CrossRef] [Green Version]
  15. Sakamoto, K.Q.; Sato, K.; Ishizuka, M.; Watanuki, Y.; Takahashi, A.; Daunt, F.; Wanless, S. Can Ethograms Be Automatically Generated Using Body Acceleration Data from Free-Ranging Birds? PLoS ONE 2009, 4, e5379. [Google Scholar] [CrossRef] [Green Version]
  16. Nathan, R.; Spiegel, O.; Fortmann-Roe, S.; Harel, R.; Wikelski, M.; Getz, W.M. Using tri-axial acceleration data to identify behavioral modes of free-ranging animals: General concepts and tools illustrated for griffon vultures. J. Exp. Biol. 2012, 215, 986–996. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Williams, H.J.; Shepard, E.L.C.; Duriez, O.; Lambertucci, S.A. Can accelerometry be used to distinguish between flight types in soaring birds? Anim. Biotelem. 2015, 3, 45. [Google Scholar] [CrossRef] [Green Version]
  18. Korner, P.; Sauter, A.; Fiedler, W.; Jenni, L. Variable allocation of activity to daylight and night in the mallard. Anim. Behav. 2016, 115, 69–79. [Google Scholar] [CrossRef]
  19. Sur, M.; Suffredini, T.; Wessells, S.M.; Bloom, P.H.; Lanzone, M.; Blackshire, S.; Sridhar, S.; Katzner, T. Improved supervised classification of accelerometry data to distinguish behaviors of soaring birds. PLoS ONE 2017, 12, e0174785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Viviant, M.; Trites, A.W.; Rosen, D.A.S.; Monestiez, P.; Guinet, C. Prey capture attempts can be detected in steller sea lions and other marine predators using accelerometers. Polar Biol. 2010, 33, 713–719. [Google Scholar] [CrossRef]
  21. Whitney, N.M.; Pratt, H.L.; Pratt, T.C.; Carrier, J.C. Identifying shark mating behaviour using three-dimensional acceleration loggers. Endanger. Species Res. 2010, 10, 71–82. [Google Scholar] [CrossRef]
  22. Brownscombe, J.W.; Gutowsky, L.F.G.; Danylchuk, A.J.; Cooke, S.J. Foraging behaviour and activity of a marine benthivorous fish estimated using tri-axial accelerometer biologgers. Mar. Ecol. Prog. Ser. 2014, 505, 241–251. [Google Scholar] [CrossRef] [Green Version]
  23. Brewster, L.R.; Dale, J.J.; Guttridge, T.L.; Gruber, S.H.; Hansell, A.C.; Elliott, M.; Cowx, I.G.; Whitney, N.M.; Gleiss, A.C. Development and application of a machine learning algorithm for classification of elasmobranch behaviour from accelerometry data. Mar. Biol. 2018, 165, 62. [Google Scholar] [CrossRef] [Green Version]
  24. Watanabe, S.; Izawa, M.; Kato, A.; Ropert-Coudert, Y.; Naito, Y. A new technique for monitoring the detailed behaviour of terrestrial animals: A case study with the domestic cat. Appl. Anim. Behav. Sci. 2005, 94, 117–131. [Google Scholar] [CrossRef]
  25. Soltis, J.; Wilson, R.P.; Douglas-Hamilton, I.; Vollrath, F.; King, L.E.; Savage, A. Accelerometers in collars identify behavioral states in captive African elephants Loxodonta africana. Endanger. Species Res. 2012, 18, 255–263. [Google Scholar] [CrossRef] [Green Version]
  26. McClune, D.W.; Marks, N.J.; Wilson, R.P.; Houghton, J.D.R.; Montgomery, I.W.; McGowan, N.E.; Gormley, E.; Scantlebury, M. Tri-axial accelerometers quantify behaviour in the Eurasian badger (Meles meles): Towards an automated interpretation of field data. Anim. Biotelem. 2014, 2, 5. [Google Scholar] [CrossRef] [Green Version]
  27. Graf, P.M.; Wilson, R.P.; Qasem, L.; Hackländer, K.; Rosell, F. The Use of Acceleration to Code for Animal Behaviours: A Case Study in Free-Ranging Eurasian Beavers Castor fiber. PLoS ONE 2015, 10, e0136751. [Google Scholar] [CrossRef] [Green Version]
  28. Kröschel, M.; Reineking, B.; Werwie, F.; Wildi, F.; Storch, I. Remote monitoring of vigilance behaviour in large herbivores using acceleration data. Anim. Biotelem. 2017, 5, 10. [Google Scholar] [CrossRef] [Green Version]
  29. Grünewälder, S.; Broekhuis, F.; Macdonald, D.W.; Wilson, A.M.; McNutt, J.W.; Shawe-Taylor, J.; Hailes, S. Movement Activity Based Classification of Animal Behaviour with an Application to Data from Cheetah (Acinonyx jubatus). PLoS ONE 2012, 7, e49120. [Google Scholar] [CrossRef] [Green Version]
  30. Campbell, H.A.; Gao, L.; Bidder, O.R.; Hunter, J.; Franklin, C.E. Creating a behavioural classification module for acceleration data: Using a captive surrogate for difficult to observe species. J. Exp. Biol. 2013, 216, 4501–4506. [Google Scholar] [CrossRef] [Green Version]
  31. Bidder, O.R.; Campbell, H.A.; Gómez-Laich, A.; Urgé, P.; Walker, J.; Cai, Y.; Gao, L.; Quintana, F.; Wilson, R.P. Love Thy Neighbour: Automatic Animal Behavioural Classification of Acceleration Data Using the K-Nearest Neighbour Algorithm. PLoS ONE 2014, 9, e88609. [Google Scholar] [CrossRef] [Green Version]
  32. Wilson, J.W.; Mills, M.G.L.; Wilson, R.P.; Peters, G.; Mills, M.E.J.; Speakman, J.R.; Durant, S.M.; Bennett, N.C.; Marks, N.J.; Scantlebury, M. Cheetahs, Acinonyx jubatus, balance turn capacity with pace when chasing prey. Biol. Lett. 2013, 9, 20130620. [Google Scholar] [CrossRef] [Green Version]
  33. Caro, T.M. Cheetahs of the Serengeti Plains; The University of Chicago Press: Chicago, IL, USA, 1994. [Google Scholar]
  34. Melzheimer, J.; Streif, S.; Wasiolka, B.; Fischer, M.; Thalwitzer, S.; Heinrich, S.K.; Weigold, A.; Hofer, H.; Wachter, B. Queuing, takeovers, and becoming a fat cat: Long-term data reveal two distinct male spatial tactics at different life-history stages in Namibian cheetahs. Ecosphere 2018, 9, e02308. [Google Scholar] [CrossRef]
  35. Melzheimer, J.; Heinrich, S.K.; Wasiolka, B.; Mueller, R.; Thalwitzer, S.; Palmegiani, I.; Weigold, A.; Portas, R.; Roeder, R.; Krofel, M.; et al. Communication hubs of an asocial cat are the source of a human–carnivore conflict and key to its solution. Proc. Natl. Acad. Sci. USA 2020, 117, 33325–33333. [Google Scholar] [CrossRef]
  36. Durant, S.; Mitchell, N.; Ipavec, A.; Groom, R. Acinonyx jubatus, Cheetah. IUCN Red List Threat. Species 2015, e.T219A50649567. [Google Scholar] [CrossRef]
  37. Durant, S.M.; Mitchell, N.; Groom, R.; Pettorelli, N.; Ipavec, A.; Jacobson, A.P.; Woodroffe, R.; Böhm, M.; Hunter, L.T.B.; Becker, M.S.; et al. The global decline of cheetah Acinonyx jubatus and what it means for conservation. Proc. Natl. Acad. Sci. USA 2017. [Google Scholar] [CrossRef] [Green Version]
  38. Weise, F.J.; Vijay, V.; Jacobson, A.P.; Schoonover, R.F.; Groom, R.J.; Horgan, J.; Keeping, D.; Klein, R.; Marnewick, K.; Maude, G.; et al. The distribution and numbers of cheetah (Acinonyx jubatus) in southern Africa. PeerJ 2017, 5, e4096. [Google Scholar] [CrossRef] [Green Version]
  39. Marker, L.L.; Mills, M.G.L.; Macdonald, D.W. Factors Influencing Perceptions of Conflict and Tolerance toward Cheetahs on Namibian Farmlands. Conserv. Biol. 2003, 17, 1290–1298. [Google Scholar] [CrossRef]
  40. Wachter, B.; Broekhuis, F.; Melzheimer, J.; Horgan, J.; Chelysheva, E.V.; Marker, L.; Mills, G.; Caro, T. Behaviour and communication of free-ranging cheetahs. In Cheetahs: Biology and Conservation; Marker, L., Boast, L., Schmidt-Küntzel, A., Eds.; Academic Press: San Diego, CA, USA, 2018; pp. 121–134. [Google Scholar]
  41. Menges, V.; Melzheimer, J. Making an uneasy co-existence easier—Leopards and farmers. Conserv. Environ. Namib. 2015, 23–25. [Google Scholar]
  42. Stein, A.B.; Bourquin, S.L.; McNutt, J.W. Avoiding intraguild competition: Leopard feeding ecology and prey caching in northern Botswana. Afr. J. Wildl. Res. 2015, 45, 1–11. [Google Scholar] [CrossRef]
  43. Gese, E.M.; Terletzky, P.A.; Cavalcanti, S.M.C. Identification of kill sites from GPS cluster for Jaguars (Panthera onca) in the southern Pantanal, Brazil. Wildl. Res. 2016, 43, 130–139. [Google Scholar] [CrossRef] [Green Version]
  44. Wang, Y.; Nickel, B.; Rutishauser, M.; Bryce, C.M.; Williams, T.M.; Elkaim, G.; Wilmers, C.C. Movement, resting, and attack behaviors of wild pumas are revealed by tri-axial accelerometer measurements. Mov. Ecol. 2015, 3, 2. [Google Scholar] [CrossRef] [Green Version]
  45. Fröhlich, M.; Berger, A.; Kramer-Schadt, S.; Heckmann, I.; Martins, Q. Complementing GPS Cluster Analysis with Activity Data for Studies of Leopard (Panthera pardus) Diet. South Afr. J. Wildl. Res. 2012, 42, 104–110. [Google Scholar] [CrossRef]
  46. Blecha, K.A.; Alldredge, M.W. Improvements on GPS Location Cluster Analysis for the Prediction of Large Carnivore Feeding Activities: Ground-Truth Detection Probability and Inclusion of Activity Sensor. PLoS ONE 2015, 10, e0138915. [Google Scholar] [CrossRef] [PubMed]
  47. e-obs GmbH System Manual. 2014: GPS-Acceleration-Tags. Grünwald, Germany. Available online: www.e-obs.de (accessed on 3 August 2021).
  48. Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental algorithms for scientific computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [Green Version]
  49. Van Rossum, G.; Drake, F.L. Python 3 Reference Manual; Createspace: Scotts Valley, CA, USA, 2009. [Google Scholar]
  50. Resheff, Y.S.; Rotics, S.; Harel, R.; Spiegel, O.; Nathan, R. AcceleRater: A web application for supervised learning of behavioral modes from acceleration measurements. Mov. Ecol. 2014, 2, 27. [Google Scholar] [CrossRef] [Green Version]
  51. Izenman, A.J. Linear Discriminant Analysis. In Modern Multivariate Statistical Techniques; Springer Texts in Statistics; Springer: New York, NY, USA, 2013. [Google Scholar] [CrossRef]
  52. Schliep, K.; Hechenbichler, K. kknn: Weighted k-Nearest Neighbors. R Package Version 1.3.1. 2016. Available online: https://CRAN.R-project.org/package=kknn (accessed on 3 August 2021).
  53. Therneau, T.; Atkinson, B. rpart: Recursive Partitioning and Regression Trees. R Package Version 4.1-15. 2019. Available online: https://CRAN.R-project.org/package=rpart (accessed on 3 August 2021).
  54. Martiskainen, P.; Järvinen, M.; Skön, J.P.; Tiirikainen, J.; Kolehmainen, M.; Mononen, J. Cow behaviour pattern recognition using three-dimensional accelerometer and support vector machines. Appl. Anim. Behav. Sci. 2009, 119, 32–38. [Google Scholar] [CrossRef]
  55. Meyer, D.; Dimitriadou, E.; Hornik, K.; Weingessel, A.; Leisch, F. e1071: Misc Functions of the Department of Statistics, Probability Theory Group (Formerly: E1071), TU Wien. R Package Version 1.7-7. 2021. Available online: https://CRAN.R-project.org/package=e1071 (accessed on 3 August 2021).
  56. Breimann, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  57. Liaw, A.; Wiener, M. Classification and Regression by randomForest. R News 2002, 2, 18–22. [Google Scholar]
  58. Mills, M.G.L.; Broomhall, L.S.; du Toit, J.T. Cheetah Acinonyx jubatus feeding ecology in the Kruger National Park and a comparison across African savanna habitats: Is the cheetah only a successful hunter on open grassland plains? Wildl. Biol. 2004, 10, 177–186. [Google Scholar] [CrossRef]
  59. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  60. Mills, M.G.L.; Mills, M.E.J. Kalahri Cheetahs, Adaptations to an Arid Region, 1st ed.; Oxford University Press: New York, NY, USA, 2017. [Google Scholar]
  61. Rast, W.; Kimmig, S.E.; Giese, L.; Berger, A. Machine learning goes wild: Using data from captive individuals to infer wildlife behaviours. PLoS ONE 2020, 15, e0227317. [Google Scholar] [CrossRef]
  62. Moreau, M.; Siebert, S.; Buerkert, A.; Schlecht, E. Use of a tri-axial accelerometer for automated recording and classification of goats’ grazing behaviour. Appl. Anim. Behav. Sci. 2009, 119, 158–170. [Google Scholar] [CrossRef]
  63. Painter, M.S.; Blanco, J.A.; Malkemper, E.P.; Anderson, C.; Sweeney, D.C.; Hewgley, C.W.; Červený, J.; Hart, V.; Topinka, V.; Belotti, E.; et al. Use of biologgers to characterize red fox behavior with implications for studies of magnetic alignment responses in free roaming animals. Anim. Biotelem. 2016, 4, 20. [Google Scholar] [CrossRef] [Green Version]
  64. Hayward, M.W.; Hofmeyr, M.; O’Brien, J.; Kerley, G.I.H. Prey preferences of the cheetah (Acinonyx jubatus) (Felidae: Carnivora): Morphological limitations or the need to capture rapidly consumable prey before kleptoparasites arrive? J. Zool. 2006, 270, 615–627. [Google Scholar] [CrossRef]
  65. Scantlebury, D.M.; Mills, M.G.L.; Wilson, R.P.; Wilson, J.W.; Mills, M.E.J.; Durant, S.M.; Bennett, N.C.; Bradford, P.; Marks, N.J.; Speakman, J.R. Flexible energetics of cheetah hunting strategies provide resistance against kleptoparasitism. Science 2014, 346, 79–81. [Google Scholar] [CrossRef] [Green Version]
  66. Hetem, R.S.; Mitchell, D.; de Witt, B.A.; Fick, L.G.; Maloney, S.K.; Meyer, L.C.R.; Fuller, A. Body temperature, activity patterns and hunting in free-living cheetah: Biologging reveals new insights. Integr. Zool. 2019, 14, 30–47. [Google Scholar] [CrossRef]
  67. Macdonald, D.W. Animal behaviour and its role in carnivore conservation: Examples of seven deadly threats. Anim. Behav. 2016. [Google Scholar] [CrossRef]
  68. Suraci, J.P.; Frank, L.G.; Oriol-Cotterill, A.; Ekwanga, S.; Williams, T.M.; Wilmers, C.C. Behavior-specific habitat selection by African lions may promote their persistence in a human-dominated landscape. Ecology 2019, 100, 1–11. [Google Scholar] [CrossRef]
  69. Smith, J.A.; Wang, Y.; Wilmers, C.C. Top carnivores increase their kill rates on prey as a response to human-induced fear. Proc. R. Soc. B Biol. Sci. 2015, 282, 20142711. [Google Scholar] [CrossRef]
  70. Weise, F.J.; Hauptmeier, H.; Stratford, K.J.; Hayward, M.W.; Aal, K.; Heuer, M.; Tomeletso, M.; Wulf, V.; Somers, M.J.; Stein, A.B. Lions at the Gates: Trans-disciplinary Design of an Early Warning System to Improve Human-Lion Coexistence. Front. Ecol. Evol. 2019, 6, 242. [Google Scholar] [CrossRef] [Green Version]
  71. Lamb, C.T.; Ford, A.T.; McLellan, B.N.; Proctor, M.F.; Mowat, G.; Ciarniello, L.; Nielsen, S.E.; Boutin, S. The ecology of human-carnivore coexistence. Proc. Natl. Acad. Sci. USA 2020, 117, 17876–17883. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Observer taking notes of a captive cheetah’s behaviour. Captive cheetahs allowed observers to approach to approximately 10–20 m. Collars were set to record acceleration data on the z-axis (yellow arrow in insert), which translates to up-down movements.
Figure 1. Observer taking notes of a captive cheetah’s behaviour. Captive cheetahs allowed observers to approach to approximately 10–20 m. Collars were set to record acceleration data on the z-axis (yellow arrow in insert), which translates to up-down movements.
Sensors 21 05426 g001
Figure 2. Movement paths from both males of coalition 2 for the same nine days of continuous GPS data with (a) collar taking a position every 3 min and (b) collar taking a position every 15 min. Blue stars mark the centre of a resting cluster. Green circles mark GPS coordinates representing walking events which could only be calculated for (a) since the GPS interval of 15 min in (b) was too long to detect these events. Red stars mark the centre of a feeding event that was also confirmed in the field. Orange stars mark the centre of a feeding event that was not confirmed in the field.
Figure 2. Movement paths from both males of coalition 2 for the same nine days of continuous GPS data with (a) collar taking a position every 3 min and (b) collar taking a position every 15 min. Blue stars mark the centre of a resting cluster. Green circles mark GPS coordinates representing walking events which could only be calculated for (a) since the GPS interval of 15 min in (b) was too long to detect these events. Red stars mark the centre of a feeding event that was also confirmed in the field. Orange stars mark the centre of a feeding event that was not confirmed in the field.
Sensors 21 05426 g002
Figure 3. Schematic overview of the methods and results of the study (EL = ensemble learning approach).
Figure 3. Schematic overview of the methods and results of the study (EL = ensemble learning approach).
Sensors 21 05426 g003
Table 1. Settings of the accelerometer recordings. The settings for the three captive brothers were programmed at 10.0 Hz and 4.0 s burst length during most of the 36 observation days and at 33.3 Hz and 3.3 s burst length during 7 to 15 days.
Table 1. Settings of the accelerometer recordings. The settings for the three captive brothers were programmed at 10.0 Hz and 4.0 s burst length during most of the 36 observation days and at 33.3 Hz and 3.3 s burst length during 7 to 15 days.
AnimalsTag FirmwareSampling Frequency [Hz]Burst Length [s]ACC Samples per Burst & Axis
Captive brothersA10.0 or 33.34.0 or 3.340 or 110
Captive maleA104.040
Captive femaleB10.543.840
Free-ranging malesB103.636
Table 2. Results of the leave-one-out cross-validations (LOOCV) with a probability threshold of 0.5 for behaviour classification for each of the six machine learning algorithms (MLA), i.e., linear discriminant analysis (LDA), quadratic discriminant analysis (QDA), the k-nearest neighbour (KNN) algorithm, the classification and regression tree (CART) algorithm, support vector machine (SVM) algorithm and the random forest (RF) algorithm, and the two ensemble-learning approaches, i.e., the majority and mean voting. Precisions (PR) across all acceleration data and for the behaviours drinking (D), feeding (F), grooming (G), resting (R), trotting/running (T) and walking (W), and recall (RE) for all behaviours are shown. PR and RE were calculated using the respective confusion matrices, with PR = TP/(TP+FP) and RE = TP/(TP+FN). TP = true positive, FP = false positive, FN = false negative.
Table 2. Results of the leave-one-out cross-validations (LOOCV) with a probability threshold of 0.5 for behaviour classification for each of the six machine learning algorithms (MLA), i.e., linear discriminant analysis (LDA), quadratic discriminant analysis (QDA), the k-nearest neighbour (KNN) algorithm, the classification and regression tree (CART) algorithm, support vector machine (SVM) algorithm and the random forest (RF) algorithm, and the two ensemble-learning approaches, i.e., the majority and mean voting. Precisions (PR) across all acceleration data and for the behaviours drinking (D), feeding (F), grooming (G), resting (R), trotting/running (T) and walking (W), and recall (RE) for all behaviours are shown. PR and RE were calculated using the respective confusion matrices, with PR = TP/(TP+FP) and RE = TP/(TP+FN). TP = true positive, FP = false positive, FN = false negative.
MLAOverallDFGRTW
PRREPRREPRREPRREPRREPRREPRRE
LDA71.555.5<0.10.081.654.871.76.894.087.310089.282.094.6
QDA75.576.241.856.278.271.350.044.699.090.797.899.286.194.9
KNN87.565.712.412.478.374.476.419.994.994.810099.295.093.7
CART57.559.7<0.10.074.475.40.00.092.193.298.097.780.191.9
RF87.474.482.437.880.477.970.844.497.696.898.898.594.191.3
SVM86.974.579.136.878.882.474.038.696.697.610099.292.892.2
Mean91.467.897.117.881.574.980.625.797.996.899.698.191.993.5
Majority87.970.984.426.579.479.477.529.096.597.599.699.290.094.0
Table 3. GPS-informed feeding events based on the 3-min schedule from the field study in 2015, results of the model application for each detected feeding event with the corresponding GPS cluster, information on whether bursts with high acceleration (ACC) variation were detected one hour before the feeding event and whether the two data sets matched in the identified feeding events. Matches are indicated by a check mark in the respective column, a cross indicates no match between the results. A minus indicates that the model did not detect a feeding event and an asterisk * indicates that the model detected a feeding event, but no prey remains were detected at the GPS cluster. The symbol # indicates instances for which 3-min GPS data were not available for male 3 and the 15-min GPS data of the coalition partner were used.
Table 3. GPS-informed feeding events based on the 3-min schedule from the field study in 2015, results of the model application for each detected feeding event with the corresponding GPS cluster, information on whether bursts with high acceleration (ACC) variation were detected one hour before the feeding event and whether the two data sets matched in the identified feeding events. Matches are indicated by a check mark in the respective column, a cross indicates no match between the results. A minus indicates that the model did not detect a feeding event and an asterisk * indicates that the model detected a feeding event, but no prey remains were detected at the GPS cluster. The symbol # indicates instances for which 3-min GPS data were not available for male 3 and the 15-min GPS data of the coalition partner were used.
GPS-Informed Feeding EventsModel ApplicationMatch
DateStartEndDuration [min]MaleDateStartEndDuration [min]GPS Schedule [min]Start GPS ClusterEnd GPS ClusterDuration GPS Cluster [min]High ACC before
COALITION 17 February19:0322:0618317 February19:2821:08100319:0322:06183no
27 February19:3020:44731519:1522:00165yes
11 February05:0007:30150111 February05:0207:34152305:0307:54171yes
211 February05:4607:16901505:0007:45165no
14 February19:3022:30180114 February19:4022:18158319:2722:33186yes
214 February19:5222:041321519:3022:30180no
17 February18:3920:0384117 February18:3020:0090318:3020:0090no
217 February18:4519:44591518:4520:0075no
21 March07:4609:2498121 March07:3808:5880307:3909:24105no
221 March00:1509:10551508:1509:1560yes
COALITION2 13 March17:0619:27141313 March17:0619:20134317:0619:30144no
413 March18:0418:56521517:1519:30135yes
16/17 March *22:48 *01:25 *157 *316/17 March22:5401:18144322:4801:21153nox
417 March00:5001:28381523:0001:15135nox
19 March20:1500:00225319 March20:1623:20184320:1500:06231yes
419 March20:2621:42761520:1500:00225yes
29/30 April23:3001:151053 #29/30 April23:2001:3413436023:3001:30120yes
4----15----x
6 May02:0005:001803 #6 May02:0005:0218236002:0005:00180no
46 May03:1804:24661502:0005:00180no
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Giese, L.; Melzheimer, J.; Bockmühl, D.; Wasiolka, B.; Rast, W.; Berger, A.; Wachter, B. Using Machine Learning for Remote Behaviour Classification—Verifying Acceleration Data to Infer Feeding Events in Free-Ranging Cheetahs. Sensors 2021, 21, 5426. https://doi.org/10.3390/s21165426

AMA Style

Giese L, Melzheimer J, Bockmühl D, Wasiolka B, Rast W, Berger A, Wachter B. Using Machine Learning for Remote Behaviour Classification—Verifying Acceleration Data to Infer Feeding Events in Free-Ranging Cheetahs. Sensors. 2021; 21(16):5426. https://doi.org/10.3390/s21165426

Chicago/Turabian Style

Giese, Lisa, Jörg Melzheimer, Dirk Bockmühl, Bernd Wasiolka, Wanja Rast, Anne Berger, and Bettina Wachter. 2021. "Using Machine Learning for Remote Behaviour Classification—Verifying Acceleration Data to Infer Feeding Events in Free-Ranging Cheetahs" Sensors 21, no. 16: 5426. https://doi.org/10.3390/s21165426

APA Style

Giese, L., Melzheimer, J., Bockmühl, D., Wasiolka, B., Rast, W., Berger, A., & Wachter, B. (2021). Using Machine Learning for Remote Behaviour Classification—Verifying Acceleration Data to Infer Feeding Events in Free-Ranging Cheetahs. Sensors, 21(16), 5426. https://doi.org/10.3390/s21165426

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