Modelling Animal Activity as Curves: An Approach Using Wavelet-Based Functional Data Analysis

Barbara Henning¹, Airton Kist², Alusio Pinheiro³, Rafael L. Camargo⁴, Thiago M. Batista⁴, Everardo M. Carneiro⁴, Sérgio F. dos Reis⁵

¹Programa de Pós-graduação em Ecologia, Universidade Estadual de Campinas, Campinas, São Paulo, Brazil.
²Departamento de Matemática e Estatstica, Universidade Estadual de Ponta Grossa, Ponta Grossa, Brazil.
³Departamento de Estatstica, Universidade Estadual de Campinas, Campinas, Brazil.
⁴Departamento de Biologia Estrutural e Funcional, Universidade Estadual de Campinas, Campinas, Brazil.
⁵Departamento de Biologia Animal, Universidade Estadual de Campinas, Campinas, Brazil.

DOI: 10.4236/ojs.2017.72016                

Abstract

Temporal activity patterns in animals emerge from complex interactions between choices made by organisms as responses to biotic interactions and challenges posed by external factors. Temporal activity pattern is an inherently continuous process, even being recorded as a time series. The discreteness of the data set is clearly due to data-acquisition limitations rather than a true underlying discrete nature of the phenomenon itself. Therefore, curves are a natural representation for high-frequency data. Here, we fully model temporal activity data as curves integrating wavelets and functional data analysis, allowing for testing hypotheses based on curves rather than on scalar and vector-valued data. Temporal activity data were obtained experimentally for males and females of a small-bodied marsupial and modelled as wavelets with independent and identically distributed errors and dependent errors. The null hypothesis of no difference in temporal activity pattern between male and female curves was tested with functional analysis of variance (FANOVA). The null hypothesis was rejected by FANOVA and we discussed the differences in temporal activity pattern curves between males and females in terms of ecological and life-history attributes of the reference species. We also performed numerical analysis that shed light on the regularity properties of the wavelet bases used and the thresholding parameters.

Keywords

Functional Analysis of Variance, High-Frequency Data, Temporal Activity Pattern, Shrinkage, Wavelet Thresholding

1. Introduction

Temporal activity of animals is the result of behavioral responses to external (environmental) factors, such as availability of food resources and predation risk, and the internal states of individuals, such as nutritional condition, aversion to risk and reproductive drive [1] [2] [3] [4] . The constraints imposed by external factors and internal states dictate that animal temporal activity pattern in ter- restrial, aquatic, and aerial environments is interrupted by pauses, angle turnings, and changes in speed [5] [6] [7] [8] [9] . Temporal activity pattern is marked with punctuations in the movement (pauses and changes in speed), temporally autocorrelated, and localized in nature [3] [9] [10] . Therefore, ecological and behavioral inference on patterns and processes of animal temporal activity are based on data with inherent autocorrelation and localization properties [3] [11] . Furthermore, technological advances are making high-dimensional data on ani- mal temporal activity increasingly available [10] [12] - [17] .

Animal temporal activity is recorded as a series of time points but one can successfully argue that the discreteness of the data set is due to a technological limitation (data-acquisition capabilities) rather than a true underlying discrete nature of the phenomenon itself. Thus, as opposed to the standard statistical cases in which the observations are often numbers (scalars) or vectors [18] , high- frequency data such as temporal activity records are continuous curves (functional data). Wavelets are functions that are able to represent a signal in a time series in both time and large and small scale domains [18] . Such decomposition into time-scale space allows the identification of the dominant modes of variability and how these modes vary with time [19] . Wavelets have become the method of choice to deal with such data because wavelet transforms are particularly suited to handle data that are high dimensional, autocorrelated, and localized [3] [18] [20] [21] . Whereas wavelet transforms have been recently used to characterize autocorrelative properties of animal temporal activity [3] [15] [22] , wavelets have not yet been explored to model animal temporal activity as curves. In fact, when- ever wavelets are used, their statistical analysis is based on statistical tests de- signed for scalar or vector random variables and not functional data [23] .

The versatility of wavelets is breathtaking. They have been used to model pro- blems as widely disconnected as resampling time series of surrogate data derived from random cascades on dyadic trees [24] , to establish a connection between discrete wavelet transforms, and in entanglement renormalization for quantum systems on the lattice [25] . Important applications of wavelets are also found in preserving motion discontinuities along the edges of weak textures and for dealing with rotations that exist in image sequences [26] , the combined analysis of ther- mal and visible light images of plants to detect early disease with high accuracy [27] , and early detection of melanomas from images of boundary irregularities of skin lesions [28] . Wavelets are also useful for high-speed detection of transient high impedance faults and power-quality disturbances [29] , feature extraction, discriminant analysis, and classification rules as crucial issues for face recogni- tion [30] , and in the development of a full-fledged theory of multiresolution sig- nal decomposition via wavelet representation [31] .

In this note, we use wavelets to model temporal activity data as curves in the context of functional data analysis (FDA [32] [33] [34] [35] ). One may argue the case for treating these temporal activity data as vector-valued data. Even though this may be possible, we can advance two main reasons why treating the data as curves is more appropriate in practice. First, the vector-based and the functional data paradigms interpret the augmentation of sample sizes in two opposing views. The former understands new observations as an elongation of the time interval in which data are being provided. For example, if the original time points are 0.1, 0.2, 0.3, 0.4, 0.5, a new data point will happen at some time larger than 0.5, say 0.6. On the other hand, functional data analysis is based on a enhanced capacity of finer sampling. On the same example, instead of augmenting the original sampling scope from [0.1, 0.5] to [0.1, 0.6], one captures more points on the same interval [0.1, 0.5], say 0.15. This increasing sampling capacity is sometimes referred to as micro dynamics in opposition to the macro dynamics of vector- based analysis. Because animals may live in a limited space in nature or they may be confined to a laboratory setting and have a finite life-time, increasing sample size as increasing space (or time) is less desirable than increasing sampling rate (timewise or spatially). This refinement can be easily related to infinitesimal calculus so functions are more appropriate than their finite dimensional vector counterparts.

Wavelets are therefore used to model temporal activity data as curves in the context of FDA, which deals with functional responses as when units are ob- served over time. FDA theory provides the means for testing hypotheses based on curves rather than on scalar and vector-valued data [18] [36] . In particular, we will use functional analysis of variance (FANOVA), which is the equivalent of ANOVA for functional data [37] . In FANOVA for two functional samples with a common covariance function one wishes to test the null hypothesis of equality of mean functions versus the alternative hypothesis that the mean func- tions are different [37] . Our approach is exemplified by activity data obtained experimentally for a small neotropical marsupial, Gracilinanus microtarsus. We chose this species as a reference system because there is currently relevant eco- logical data available, including demography and life history [38] - [43] . Survival rate for females remains fairly constant in the pre- and post-mating periods, whereas male survival rate decreases significantly in the post-mating period due to stress from aggressive interactions between males [41] [44] . Adult males (30 - 45 g) are much heavier than female gracile opossum (20 - 30 g [41] [42] ). Body size has a positive, linear association with the area traversed by an individual in its normal activities of food gathering, mating and caring for the young [45] [46] [47] . As expected, the area traversed in nature is larger and more variable in male than in female gracile opossum [48] . Therefore, we expect that activity patterns will differ between male and female gracile opossum. Specifically, we will test the null hypothesis of homogeneity of mean curves of activity patterns from male and female gracile opossum. To the best of our knowledge this is the first attempt to integrate wavelets to model curves derived from activity data with functional data analysis in the framework of functional analysis of variance.

2. Theory

Ramsay et al. (2005) [36] provide the reader with a generalized framework for the analysis of functional data, which basically depends on regularity conditions of the underlying curve. Wavelets are especially built to provide regular esti- mates through multiscale shrinkage [18] . We refer to Kist and Pinheiro (2015) [20] for a detailed development of the wavelet functional data analysis for de- pendent errors.

Wavelets are basically elements of some specially built basis of the space of square-integrable functions. This means the following. If $f$ is a function so that ${\displaystyle \int }f{\left(x\right)}^2\text{d}x$ is finite and $\left\{{\varphi }_{0k};{\psi }_{jk}\left(x\right),,j\ge 0,k=0,1,2,\cdots \right\}$ is a wavelet basis for the space of square-integrable functions, there are constants $\left\{c_{0k}\right\}$ and $\left\{d_{jk}\right\}$ such that $f$ can be written as:

$f\left(x\right)={\displaystyle \underset{k}{\sum }}{c}_{0k}{\varphi }_{0k}\left(x\right)+{\displaystyle \underset{j,k}{\sum }}{d}_{jk}{\psi }_{jk}\left(x\right).$

There are many different wavelet bases. Their main characteristic is that all elements of a wavelet basis have the same form, being different only in location and/or dilation. Therefore, wavelet analysis is immediately equipped with a fast transform algorithm.

In what follows, we give a general presentation of this procedure as needed for the animal temporal activity data set.

Each observation is composed of $K$ time-point evaluations of a function of interest $f_g$ , $g=1,2$ , such that

$\text{d}{X}_{gl}\left(t\right)=f_g\left(t\right)\text{d}t+ϵ\text{d}{W}_{gl}\left(t\right),t\in \left[0,1\right],$ (1)

where $ϵ$ is the diffusion parameter, $f_g\left(t\right)$ is a deterministic unknown func- tion and $\left\{W_{gl}\left(t\right):t\in ℝ\right\}$ are either independent standard Brownian motions, - case independent and identically distributed, or independent Continuous Time Autoregressive Moving Average (CTARMA) processes [20] case dependent, for $g=1,\cdots ,G$ , $l=1,\cdots ,n_g$ . We call $X_{11},\cdots ,X_{1n_1},\cdots ,X_{G1},\cdots ,X_{G{n}_G}$ the $n_g$ $K$ -di- mensional observations in group $g=1,2$ , respectively.

Suppose $f$ belongs to a convenient Besov space [20] . We write the non- linear wavelet estimators as

${\hat{f}}_g(\cdot )={\displaystyle \underset{k}{\sum }}{\hat{c}}_{g,0k}{\varphi }_{0k}\left(x\right)+{\displaystyle \underset{j,k}{\sum }}{\hat{d}}_{g,jk}^{thr}{\psi }_{jk}\left(x\right).$ (2)

One can write the $L-2$ norm as ${\Vert {\hat{f}}_g-{\hat{f}}_{g^{\prime }}\Vert }^2={\displaystyle {\sum }_k}{\left({\hat{c}}_{g,0k}-{\hat{c}}_{g^{\prime },0k}\right)}^2+{\displaystyle {\sum }_{j,k}}{\left({\hat{d}}_{g,jk}^{thr}-{\hat{d}}_{g^{\prime },jk}^{thr}\right)}^2$ , where ${\hat{c}}_{g,jk}$ is the esti- mated approximation coefficient for the $g$ -th group, 0-th scale and $k$ -th position, while ${\hat{d}}_{g^{\prime },jk}^{thr}$ is the thresholded wavelet coefficient for the $g$ -th group, $j$ -th scale and $k$ -th position. The hypotheses of interest are

$H_0:f_1(\cdot )=\cdots =f_G(\cdot )\text{vs}{H}_1:\exists 1\le g\ne g^{\prime }\le G:f_g(\cdot )\ne f_{g^{\prime }}(\cdot ).$

The dual relationship between Besov spaces and wavelet bases leads to a na- tural change in the hypotheses being tested. Instead of testing $H_0$ vs $H_1$ , a slight formal change is made, as proposed in [49] for independent errors. A two- step procedure is developed by Kist and Pinheiro (2015) [20] for dependent errors. The test statistic is compared to the cut-off point and the decision (reject or not $H_0$ ) can then be made.

Although these tests and hypotheses are mathematically different from the aforementioned hypotheses, for all applied purposes, they all yield the same inter- pretation as follows. Whenever the empirical evidences lead to the rejection of $H_0$ one can conclude that the data provide statistical evidence that at least two of the functions $f_1,\cdots ,f_G$ are different. For instance, for our case, $G=2$ , $f_1$ and $f_2$ are the underlying functional behavior of female and male specimens, respectively. Thus, rejecting $H_0$ means that the average time-curves are taken with respect to observations of female and male behaviors are statistically diffe- rent. On the other hand, whenever $H_0$ is not rejected, one understands that there is not enough empirical evidence for each group to have different underlying functions. Again for our case, this can be interpreted as the data not providing statistically significant evidence that male and female specimens differ in their temporal activity behavior.

3. Materials and Methods

The specimens of the gracile mouse opossum (G. microtarsus) were live trapped in a savannah-like habitat in the city of Mogi-Guaçu in the state of São Paulo using Sherman live-traps (dimensions 7.5 ´ 9.0 ´ 23.5 cm) baited with banana and peanut butter. Individuals captured were marked with a numbered ear tag and their sex and age were recorded. In the laboratory, gracile mouse opossums were housed individually in acrylic boxes (44 cm width ´ 33 cm length ´ 20 cm height) with ad libitum access to food (commercial cat and dog chow) and water, and kept under artificial conditions of light (12-h light/dark cycle, light on at 7:00 A.M.) and temperature $\left(22{}^{\circ }\text{C}\right)$ .

To assess spontaneous locomotor activity of the gracile mouse opossum we used an automated motor activity monitor (Acti-Track v2.7.10, PanLab, S.L. Instrument, Barcelona, Spain [50] ). The apparatus consists of a transparent Perspex box (45 ´ 45 cm base, 35 cm height) connected to a photoelectric cell and locomotor activity is detected by light beam breaks. Thirty-two infrared beam breaks, 16 each on perpendicular walls, were mounted 3 cm above the box frame floor and connected to an interface (LE 8811, LSI Letica Scientific In- struments, Barcelona, Spain), and data were sent to a computer. Thence, loco- motor activity was assessed as a rate of light beam breaks during the period of the experiment [51] . This means that whenever the mouse opossum was moving vertically or horizontally the rate of light breaks would be recorded as the activity variable. Hence, the higher the rate of light beam breaks the higher the activity of the individual.

At the beginning of the experiment, each gracile mouse opossum was placed in the Perspex box and allowed to freely explore for 24 hours to habituate the individuals before conducting the experiments. Testing in the actimeter was done in an isolated room between 18:00 and 06:00, which corresponds to the activity of gracile mouse opossums in the wild. After each experiment, the Perspex box was carefully cleaned with a 5% ethanol cleaning solution.

Permission for animal collection was provided by SISBIO (Sistema de Auto- rização de Informação em Biodiversidade). Animal housing and experimental procedures were approved by Comissão de Ética no Uso de Animais, Univer- sidade Estadual de Campinas.

Experimental settings such as those used in our study provide relevant infor- mation not only for the study of animal temporal activity pattern, but also for several areas of ecology and behavior including for example the association be- tween social and sexual preferences and genetic variation at microsattelite loci [52] , modulation of vocalization by hormones [53] , and the link between heri- table neuroendocrine variation and male sexual behavior [54] .

The data analysis was performed as follows. Twelve hours were selected from the continuous observed curves. Three families of wavelet bases were employed on the data: Symmlets, Coiflets and Daubechies. Preliminary analyses led to one smoothness parameter for each family: Symmlets 8, Coiflets 3 and Daubechies 6. These are different wavelet bases. This means the aforementioned functions ${\varphi }_{0k}$ and ${\psi }_{jk}$ will be different if we use each of these bases. We could write $\left\{{\varphi }_{0k}^{S8},{\psi }_{jk}^{S8}\right\}$ , $\left\{{\varphi }_{0k}^{C3},{\psi }_{jk}^{C3}\right\}$ and $\left\{{\varphi }_{0k}^{D6},{\psi }_{jk}^{D6}\right\}$ for each case, respectively.

The aforementioned data set was composed of temporal activity curves from 12 hours of data acquisition, in which data were taken every second for 6 males and 7 females of G. microtarsus. This data set was analyzed using the proposed wavelet model for which there were $G=2$ groups: females and males. We then estimated $f_1$ and $f_2$ .

4. Results

Figure 1 shows three estimators for each gender, based on three different wave- let bases: Symmlets 8, Coiflets 3, and Daubechies 6. The data is shown in gray. The respective estimators for the FANOVA model based on independent errors are shown in blue, and the estimators for the FANOVA model based on depen- dent errors are shown in red. With no loss of generality, time is transformed to $\left[0,1\right]$ . One can notice the differences between ${\hat{f}}_1$ and ${\hat{f}}_2$ (see Equations (1) and (2) above). Moreover, the proposed dependent estimators are much more regularized then the previously available wavelet estimators.

Numerical results did not differ much among the bases chosen (shown in Figure 1). However, the visual results for Coifflets were in general coarser than the other two bases. Daubechies’ bases were interesting because of their theoretical and numerical properties, whereas Symmlets were the bases which more closely relate to Daubechies, and look the most symmetrical. We should point out that there are no wavelets bases that are both compactly supported and symmetrical

[21] .

Figure 1. Temporal activity curve curves analyzed using the proposed wavelet model. Curves obtained for three estimators for each gender, each onE based on a different wavelet base: Symmlets 8, Coiflets 3, and Daubechies 6. Data are shown in gray; the respective estimators for the FANOVA model based on independent errors are shown in blue; and the estimators for the FANOVA model based on dependent errors are shown in red. (a) Regularized Wavelet Mean 12-hour temporal activity curve of female Gracilianus microtarsus. The wavelet filter is the Symmlet 8, and the regularizing parameters are $j_s=5$ and $J_{\eta }=9$ . (b) Regularized Wavelet Mean 12-hour temporal activity curve of female G. microtarsus. The wavelet filter is the Coiflets 3, and the regularizing parameters are $j_s=5$ and $J_{\eta }=9$ . (c) Regularized Wavelet Mean 12-hour temporal activity curve of female G. microtarsus. The wavelet filter is the Daubechies 6, and the regularizing parameters are $j_s=5$ and $J_{\eta }=9$ . (d) Regularized Wavelet Mean 12-hour temporal activity curve of male G. microtarsus. The wavelet filter is the Symmlet 8, and the regularizing parameters are $j_s=5$ and $J_{\eta }=9$ . (e) Regularized Wavelet Mean 12-hour temporal activity curve of male G. microtarsus. The wavelet filter is the Coiflets 3, and the regularizing parameters are $j_s=5$ and $J_{\eta }=9$ . (f) Regularized Wavelet Mean 12-hour temporal activity curve of male G. microtarsus. The wavelet filter is the Daubechies 6, and the regularizing parameters are $j_s=5$ and $J_{\eta }=9$ .

The robust mean absolute deviation (MAD) or the standard deviation may be employed for the estimates of the measure of the noise variability. Standard deviation is in general superior to MAD, since the latter yields less regular estimated curves [20] . The choice of the wavelet basis is usually quite unim- portant. The use of any such basis leads to the same inferential results. Some local characteristics of the estimated curves are highlighted or shadowed by each basis, but the results are the same.

The data were comprised of a total of 43,200 observations for each specimen, which means one observation at every second for 12 hours of experiment. The curves were estimated for each sex, and several choices of wavelet basis and/or thresholding were employed. The individual autocorrelation estimates varied from 0.38 to 0.66, whilst the overall autocorrelation estimates for females and males were 0.46 and 0.55, respectively. The maximum difference between any estimates for the same data set given by the choice of the wavelet basis was not greater than .02.

One should note that, as expected, the estimated curve with independent errors was much less regular than its dependent counterpart. This happens considering each pair of curves for fixed wavelet basis and thresholding pro- cedure. Levels 6 - 9 were thresholded to estimate $\rho$ , while levels 8 - 12 were thresholded in the final curve estimate for either independent or dependent models. Figure 1 displays the preliminary results for three families of bases. From this, Coifflets do present visually less appealing results.

The FANOVA test results using Daubechies 6 and Symmlets 8, and the aforementioned thresholding configuration were all favourable to the inequality of $f_1$ and $f_2$ . For db6, the test statistics were $2.6\times {10}^5$ and $1.2\times {10}^6$ against respective critical values of 10.54 and 5.8. For sym8, the test statistics were $5.2\times {10}^5$ and $4.2\times {10}^6$ against respective critical values of 10.67 and 5.78.

An extensive numerical study was performed varying the regularity of the Daubechies and Symmlets bases and the thresholding parameters. The general results were that:

1) Test results were robust with respect to the choice of the wavelet basis and thresholding parameters.

2) The chosen basis was the least relevant aspect on the procedure, both in its inferential (and numerical) aspects, as well as the visual characteristics of the estimates.

3) The choice of the thresholded levels bears importance for the test results - one should understand that the differences between males and females were statistically significant for any estimation procedure, but the value of the test statistics and the cut-off points changed considerably with the thresholded levels.

4) The choice of the thresholded levels influenced the visual characteristics of the estimated curves.

5) Standard deviation was superior to MAD.

5. Discussion

In this study, we used wavelets to model activity data as curves and functional data analysis, in particular functional analysis of variance, to test the hypothesis that the mean activity curves differed between male and female of the small marsupial G. microtarsus. We rejected the null hypothesis by FANOVA showing that despite the common laboratory environment, male and female G. micro- tarsus did differ in their temporal activity pattern. The statistical differences in the temporal activity curves may be attributable to endogenous factors associated with the sexes. In fact, male and female G. microtarsus differ in many aspects of their life history and ecology [55] . Male and female G. microtarsus differ most remarkably in their demographics. Survival rate for females is constant during the pre- and post-mating periods. On the other hand, survival rate in males de- creases sharply in the post-mating period and is significantly lower than that of females [41] . The decrease in survival rate in males is probably explained by post-mating mortality associated with stress that results from aggressive beha- vior and fighting between males during the mating period [41] [44] . This striking difference in life-history is possibly associated with differences in activity be- tween male and female. Other aspects of the ecology and life-history of G. micro- tarsus may also play a role in the observed statistical difference between male and female curves as revealed by FANOVA. Sexual dimorphism is very pro- nounced with adult males (30 - 45 g) being heavier than adult females (20 - 30 g [41] [42] ). The dimorphism in body size arises from males growing much faster than females, as inferred from a Gompertz growth model [48] . Body size has a positive, linear association with home range size, which is the area traversed by an individual in its normal activities of food gathering, mating and caring for the young [45] [46] [47] . In fact, observations in nature showed that home-range size in G. microtarsus was larger and more variable in males (0.14 ± 0.18 ha) than that in females (0.12 ± 0.09 ha [48] ).

Whereas wavelets have long been applied to problems in behavior and ecology (e.g. [9] [56] ), the use of fuctional data analysis is very promising in this context. For example, functional principal components analysis (FPCA) was recently used to investigate the relationship between a prey species, the sandeel, and its pre- dator, the black-legged kittiwake, in a dynamic marine ecosystem [57] . They de- monstrated that FPCA was an useful tool to assess spatio-temporal patterns in natural ecosystems and their study revealed the fine scale details of the inte- raction between environmental factors and prey behavior and predator foraging behavior. Just recently, [58] developed a Bayesian model to study animal move- ment patterns at different temporal scales within the context of functional data analysis. They applied this model to estimate movement paths and associated movement descriptors of the canadian lynx reintroduced to Colorado. In this application, B-splines were used but the model was general enough to incor- porate other basis functions such as Fourier series and wavelets. The approach by [58] seems extremely promising to reveal details of animal movements with important implications for population dynamics.

We believe that functional data analysis when applied in conjunction with wavelets to model curves derived from temporal activity data has the potential to provide an important methodological breakthrough in the study of ecology and animal behavior. One final observation regarding the applicability of wavelet tech- niques such as those proposed here for animal temporal activity data concerns regular sampling. Some studies of important phenomena do not allow for equally spaced data registers. Several options are available for adapting the proposed procedure for these cases. Among them padding [21] , lifted wavelets, and other second-generation wavelets can be employed [59] .

Acknowledgements

We are very grateful to João Del Giudice Neto and Marcos Mecca Pinto of Mogi- Guaçu Biological Reserve for logistical support. We are indebted to Eduardo Guimarães Martins for critical comments that improved the quality of the manuscript.

Funding

Research supported by Fundação de Amparo à Pesquisa do Estado de São Paulo, Brazil (FAPESP) [Alusio Pinheiro-13/00506-1 and Sérgio F. dos Reis-2005/513- 53-4]; Conselho Nacional de Desenvolvimento Cientfico e Tecnológico, Brazil (CNPq) [Alusio Pinheiro-304512/2011-7, Sérgio F. dos Reis-303544/2011-2 and Barbara Henning-140773/2013-4]; and CAPES [Babara Henning].

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Getz, W.M. and Saltz, D. (2008) A Framework for Generating and Analyzing Movement Paths on Ecological Landscapes. Proceedings of the National Academy of Sciences of the United States of America, 105, 19066-19071. https://doi.org/10.1073/pnas.0801732105
[2]	Nathan, R., Getz, W.M., Revilla, E., Holyoak, M., Kadmon, R., Saltz, D. and Smouse, P.E. (2008) A Movement Ecology Paradigm for Unifying Organismal Movement Research. Proceedings of the National Academy of Sciences of the United States of America, 105, 19052-19059. https://doi.org/10.1073/pnas.0800375105
[3]	Wittemyer, G., Polansky, L., Douglas-Hamilton, I. and Getz, W.M. (2008) Disentangling the Effects of Forage, Social Rank, and Risk on Movement Autocorrelation of Elephants Using Fourier and Wavelet Analyses. Proceedings of the National Academy of Sciences of the United States of America, 105, 19108-19113. https://doi.org/10.1073/pnas.0801744105
[4]	Gurarie, E., Bracis, C., Delgado, M., Meckley, T.D., Kojola, I. and Wagner, C.M. (2015) What Is the Animal Doing? Tools for Exploring Behavioral Structure in Animal Movements. Journal of Animal Ecology, 85, 69-84. https://doi.org/10.1111/1365-2656.12379
[5]	Kramer, D.L. and McLaughlin, R.L. (2001) The Behavioral Ecology of Intermittent Locomotion. Integrative and Comparative Biology, 41, 137-153. https://doi.org/10.1093/icb/41.2.137
[6]	McAdam, A.G. and Kramer, D.L. (1998) Vigilance as a Benefit of Intermittent Locomotion in Small Mammals. Animal Behaviour, 55, 109-117.
[7]	Mashanova, A., Oliver, T.H. and Jansen, V.A. (2010) Evidence for Intermittency and a Truncated Power Law from Highly Resolved Aphid Movement Data. Journal of the Royal Society Interface, 7, 199-208. https://doi.org/10.1098/rsif.2009.0121
[8]	O’Brien, W.J., Browman, H.I. and Evans, B.I. (1990) Search Strategies of Foraging Animals. American Scientist, 78, 152-160.
[9]	Bazazi, S., Bartumeus, F., Hale, J.J. and Couzin, I.D. (2012) Intermittent Motion in Desert Locusts: Behavioural Complexity in Simple Environments. PLoS Computational Biology, 8, e1002498. https://doi.org/10.1371/journal.pcbi.1002498
[10]	Hart, T., Coulson, T. and Trathan, P.N. (2010) Time Series Analysis of Biologging Data: Autocorrelation Reveals Periodicity of Diving Behaviour in Macaroni Penguins. Animal Behaviour, 79, 845-855.
[11]	Bartumeus, F. and Levin, S.A. (2008) Fractal Reorientation Clocks: Linking Animal Behavior to statistical Patterns of Search. Proceedings of the National Academy of Sciences of the United States of America, 105, 19072-19077. https://doi.org/10.1073/pnas.0801926105
[12]	Bagniewska, J.M., Hart, T., Harrington, L.A. and MacDonald, D.W. (2013) Hidden Markov Analysis Describes Dive Patterns in Semiaquatic Animals. Behavioral Ecology, 24, 659-667.
[13]	Cagnacci, F., Boitani, L., Powell, R.A. and Boyce, M.S. (2010) Animal Ecology Meets GPS-Based Radiotelemetry: A Perfect Storm of Opportunities and Challenges. Philosophical Transactions of the Royal Society B: Biological Sciences, 365, 2157-2162. https://doi.org/10.1098/rstb.2010.0107
[14]	Rutz, C. and Hays, G.C. (2009) New Frontiers in Biologging Science. Biology Letters, 5, 289-292. https://doi.org/10.1098/rsbl.2009.0089
[15]	Polansky, L., Wittemyer, G., Cross, P.C., Tambling, C.J. and Getz, W.M. (2010) From Moonlight to Movement and Synchronized Randomness: Fourier and Wavelet Analyses of Animal Location Time Series Data. Ecology, 91, 1506-1518. https://doi.org/10.1890/08-2159.1
[16]	Benhamou, S. (2014) Of Scales and Stationarity in Animal Movements. Ecology Letters, 17, 261-272. https://doi.org/10.1111/ele.12225
[17]	Pettorelli, N., Laurance, W.F., O’Brien, T.G., Wegmann, M., Nagendra, H. and Turner, W. (2014) Satellite Remote Sensing for Applied Ecologists: Opportunities and Challenges. Journal of Applied Ecology, 51, 839-848. https://doi.org/10.1111/1365-2664.12261
[18]	Pinheiro, A. and Vidakovic, B. (2009) Wavelets in Functional Data Analysis. ICMC-USP, Sao Carlos.
[19]	Torrence, C. and Compo, G.P. (1998) A Practical Guide to Wavelet Analysis. Bulletin of the American Meteorological Society, 79, 61-78. https://doi.org/10.1175/1520-0477(1998)079<0061:APGTWA>2.0.CO;2
[20]	Kist, A. and Pinheiro, A. (2015) Wavelet Functional Data Analysis for FANOVA Models under Dependent Errors. https://arxiv.org/pdf/1511.00119
[21]	Vidakovic, B. (2009) Statistical Modeling by Wavelets. John Wiley & Sons, New York.
[22]	Riotte-Lambert, L., Benhamou, S. and Chamaillé-Jammes, S. (2013) Periodicity Analysis of Movement Recursions. Journal of Theoretical Biology, 317, 238-243.
[23]	Horváth, L. and Rice, G. (2015) Testing Equality of Means When the Observations Are from Functional Time Series. Journal of Time Series Analysis, 36, 84-108. https://doi.org/10.1111/jtsa.12095
[24]	Palus, M. (2008) Bootstrapping Multifractals: Surrogate Data from Random Cascades on Wavelet Dyadic Trees. Physical Review Letters, 101, Article ID: 134101. https://doi.org/10.1103/PhysRevLett.101.134101
[25]	Evenbly, G. and White, S.R. (2016) Entanglement Renormalization and Wavelets. Physical Review Letters, 116, Article ID: 140403. https://doi.org/10.1103/PhysRevLett.116.140403
[26]	Ali, S., Daul, C., Galbrun, E., Guillemin, F. and Blondel, W. (2016) Anisotropic Motion Estimation on Edge Preserving Riesz Wavelets for Robust Video Mosaicing. Pattern Recognition, 51, 425-442.
[27]	Raza, S.E.A., Silva, V.S., Prince, G., Clarkson, J.P. and Rajpoot, N.M. (2015) Registration of Thermal and Visible Light Images of Diseased Plants Using Silhouette Extraction in the Wavelet Domain. Pattern Recognition, 48, 2119-2128.
[28]	Ma, L. and Staunton, R.C. (2013) Analysis of the Contour Structural Irregularity of Skin Lesions Using Wavelet Decomposition. Pattern Recognition, 46, 98-106.
[29]	Costa, F.B. (2014) Boundary Wavelet Coefficients for Real-Time Detection of Transients Induced by Faults and Power-Quality Disturbances. IEEE Transactions on Power Delivery, 29, 2674-2687. https://doi.org/10.1109/TPWRD.2014.2321178
[30]	Chien, J.T. and Wu, C.C. (2002) Discriminant Waveletfaces and Nearest Feature Classifiers for Face Recognition. IEEE Transactions on Pattern Analysis and Machine Intelligence, 24, 1644-1649. https://doi.org/10.1109/TPAMI.2002.1114855
[31]	Mallat, S.G. (1989) A Theory for Multiresolution Signal Decomposition: The Wavelet Representation. IEEE Transactions on Pattern Analysis and Machine Intelligence, 11, 674-693. https://doi.org/10.1109/34.192463
[32]	Graves, S., Hooker, G. and Ramsay, J. (2009) Functional Data Analysis with R and MATLAB. Springer, New York.
[33]	Ainsworth, L.M., Routledge, R. and Cao, J. (2011) Functional Data Analysis in Ecosystem Research: The Decline of Oweekeno Lake Sockeye Salmon and Wannock River Flow. Journal of Agricultural, Biological, and Environmental Statistics, 16, 282-300. https://doi.org/10.1007/s13253-010-0049-z
[34]	Chang, C., Chen, Y. and Ogden, R.T. (2014) Functional Data Classification: A Wavelet Approach. Computational Statistics, 29, 1497-1513. https://doi.org/10.1007/s00180-014-0503-4
[35]	Horváth, L. and Rice, G. (2015) An Introduction to Functional Data Analysis and a Principal Component Approach for Testing the Equality of Mean Curves. Revista Matemática Complutense, 28, 505-548. https://doi.org/10.1007/s13163-015-0169-7
[36]	Ramsay, J.O. and Silverman, B.W. (2005) Functional Data Analysis. Springer, New York. https://doi.org/10.1002/0470013192.bsa239
[37]	Zhang, J.T. (2013) Analysis of Variance for Functional Data. CRC Press, Boca Raton.
[38]	da Silva, C.Q., Gomes, A., Martins, E., Bonato, V. and dos Reis, S.F. (2008) A Bayesian Change-Point Model for Describing Partial Semelparity of a Neotropical Didelphid Marsupial. Review of Brazilian Biometry, 26, 31-44.
[39]	da Silva, C.Q., Martins, E.G., Bonato, V. and dos Reis, S.F. (2008) Bayesian Capture-Recapture Analysis: An Application in Modeling Semelparity of a Neotropical Didelphid Marsupial. Communications in Statistics-Simulation and Computation, 37, 816-828. https://doi.org/10.1080/03610910701884310
[40]	Martins, E.G., Bonato, V., Pinheiro, H.P. and dos Reis, S.F. (2006) Diet of the Gracile Mouse Opossum (Gracilinanus microtarsus) (Didelphimorphia: Didelphidae) in a Brazilian Cerrado: Patterns of Food Consumption and Intrapopulation Variation. Journal of Zoology, 269, 21-28. https://doi.org/10.1111/j.1469-7998.2006.00052.x
[41]	Martins, E.G., Bonato, V., da Silva, C.Q. and dos Reis, S.F. (2006) Partial Semelparity in the Neotropical Didelphid Marsupial Gracilinanus microtarsus. Journal of Mammalogy, 87, 915-920. https://doi.org/10.1644/05-MAMM-A-403R1.1
[42]	Martins, E.G., Bonato, V., da Silva, C.Q. and dos Reis, S.F. (2006) Seasonality in Reproduction, Age Structure and Density of the Gracile Mouse Opossum Gracilinanus microtarsus (Marsupialia: Didelphidae) in a Brazilian Cerrado. Journal of Tropical Ecology, 22, 461-468. https://doi.org/10.1017/S0266467406003269
[43]	Martins, E.G., Araújo, M.S., Bonato, V. and dos Reis, S.F. (2008) Sex and Season Affect Individual-Level Diet Variation in the Neotropical Marsupial Gracilinanus microtarsus (Didelphidae). Biotropica, 40, 132-135.
[44]	Boonstra, R. (2005) Equipped for Life: The Adaptive Role of the Stress Axis in Male Mammals. Journal of Mammalogy, 86, 236-224. https://doi.org/10.1644/BHE-001.1
[45]	Lindstedt, S.L., Miller, B.J. and Buskirk, S.W. (1986) Home Range, Time, and Body Size in Mammals. Ecology, 67, 413-418. https://doi.org/10.2307/1938584
[46]	Burt, W.H. (1943) Territoriality and Home Range Concepts as Applied to Mammals. Journal of Mammalogy, 24, 346-352. https://doi.org/10.2307/1374834
[47]	Kelt, D.A. and Van Vuren, D. (1999) Energetic Constraints and the Relationship between Body Size and Home Range Area in Mammals. Ecology, 80, 337-340. https://doi.org/10.1890/0012-9658(1999)080[0337:ECATRB]2.0.CO;2
[48]	Fernandes, F.R., Cruz, L.D., Martins, E.G. and dos Reis, S.F. (2010) Growth and Home Range Size of the Gracile Mouse Opossum Gracilinanus microtarsus (Marsupialia: Didelphidae) in Brazilian Cerrado. Journal of Tropical Ecology, 26, 185-192. https://doi.org/10.1017/S0266467409990526
[49]	Abramovich, F., Antoniadis, A., Sapatinas, T. and Vidakovic, B. (2004) Optimal Testing in a Fixed-Effects Functional Analysis of Variance Model. International Journal of Wavelets, Multiresolution and Information Processing, 2, 323-349. https://doi.org/10.1142/S0219691304000639
[50]	Abedi, P.M., Delaville, C., De Deurwaerdere, P., Benjelloun, W. and Benazzouz, A. (2013) Intrapallidal Administration of 6-Hydroxydopamine Mimics in Large Part the Electrophysiological and Behavioral Consequences of Major Dopamine Depletion in the Rat. Neuroscience, 236, 289-297. https://doi.org/10.1016/j.neuroscience.2013.01.043
[51]	Dulawa, S.C., Grandy, D.K., Low, M.J., Paulus, M.P. and Geyer, M.A. (1999) Dopamine D4 Receptor-Knock-Out Mice Exhibit Reduced Exploration of Novel Stimuli. The Journal of Neuroscience, 19, 9550-9556.
[52]	Castelli, F.R., Kelley, R.A., Keane, B. and Solomon, N.G. (2011) Female Prairie Voles Show Social and Sexual Preferences for Males with Longer avpr1a Microsatellite Alleles. Animal Behaviour, 82, 1117-1126.
[53]	Pasch, B., George, A.S., Campbell, P. and Phelps, S.M. (2011) Androgen-Dependent Male Vocal Performance Influences Female Preference in Neotropical Singing Mice. Animal Behaviour, 82, 177-183.
[54]	Sharp, K., Bucci, D., Zelensky, P.K., Chesney, A., Tidhar, W., Broussard, D.R. and Heideman, P.D. (2015) Genetic Variation in Male Sexual Behaviour in a Population of White-Footed Mice in Relation to Photoperiod. Animal Behaviour, 104, 203-212.
[55]	Araújo, M.S., Martins, E.G., Cruz, L.D., Fernandes, F.R., Linhares, A.X., dos Reis, S.F. and Guimaraes, P.R. (2010) Nested Diets: A Novel Pattern of Individual-Level Resource Use. Oikos, 119, 81-88. https://doi.org/10.1111/j.1600-0706.2009.17624.x
[56]	Detto, M. and Muller-Landau, H.C. (2013) Fitting Ecological Process Models to Spatial Patterns Using Scalewise Variances and Moment Equations. The American Naturalist, 181, E68-E82. https://doi.org/10.1086/669678
[57]	Embling, C.B., Illian, J., Armstrong, E., van der Kooij, J., Sharples, J., Camphuysen, K.C. and Scott, B.E. (2012) Investigating Fine-Scale Spatio-Temporal Predator-Prey Patterns in Dynamic Marine Ecosystems: A Functional Data Analysis Approach. Journal of Applied Ecology, 49, 481-492. https://doi.org/10.1111/j.1365-2664.2012.02114.x
[58]	Buderman, F.E., Hooten, M.B., Ivan, J.S. and Shenk, T.M. (2015) A Functional Model for Characterizing Long-Distance Movement Behaviour. Methods in Ecology and Evolution, 7, 264-273. https://doi.org/10.1111/2041-210X.12465
[59]	Sweldens, W. (1998) The Lifting Scheme: A Construction of Second Generation Wavelets. SIAM Journal on Mathematical Analysis, 29, 511-546. https://doi.org/10.1137/S0036141095289051