Modelling and inference for the movement of interacting animals

1. Statistical modelling of animal movement data is a rapidly growing area of research. Typically though, these models have been developed for analysing the tracks of individual animals and we lose sight of the impact animals have on each other with regards to their movement behaviours. We aim to develop a model with a flexible social framework that allows us to capture that information. 2. Our approach is based on the concept of social hierarchies, and this is embedded in a multivariate diffusion process which models the movement of a group of animals. The possibility of switching between behavioural states facilitates dynamic social behaviours and we augment the observed data with sampled state switching times in order to model the animals' behaviour naturally in continuous time. In addition, this enables us to carry out exact inference in a Bayesian setting with the benefits of being able to handle regular, irregular and missing data. All movement and behaviour parameters are estimated with Markov chain Monte Carlo methods. 3. We examine the capability of our model with simulated data before fitting it to GPS locations of five wild olive baboons Papio anubis. The results enable us to identify which animals are influencing the movement of others and when, which provides both a dynamic and long-term static insight into the group's social behaviours. 4. Our model offers a flexible method in continuous time with which to model the network of social interactions within animal movement. Doing so avoids the limitations caused by a discrete-time approach and it allows us to capture rich information with regards to a group's social structure, leading to constructive applications in conservation and management decisions. However, currently it is a computationally expensive task to fit the model to data, which in turns limits extending the model to more fruitful but complex cases such as heterogeneity in space or individual characteristics. Furthermore, our social hierarchy approach assumes all relevant animals are tracked and that any interactions have some ordering, both of which narrow the scope within which this approach is appropriate

Higher physical fitness levels are associated with less language decline in healthy ageing

publishedVersio

Modelling group movement with behaviour switching in continuous time

This article presents a new method for modelling collective movement in continuous time with behavioural switching, motivated by simultaneous tracking of wild or semi‐domesticated animals. Each individual in the group is at times attracted to a unobserved leading point. However, the behavioural state of each individual can switch between ‘following’ and ‘independent’. The ‘following’ movement is modelled through a linear stochastic differential equation, while the ‘independent’ movement is modelled as Brownian motion. The movement of the leading point is modelled either as an Ornstein‐Uhlenbeck (OU) process or as Brownian motion (BM), which makes the whole system a higher‐dimensional Ornstein‐Uhlenbeck process, possibly an intrinsic non‐stationary version. An inhomogeneous Kalman filter Markov chain Monte Carlo algorithm is developed to estimate the diffusion and switching parameters and the behaviour states of each individual at a given time point. The method successfully recovers the true behavioural states in simulated data sets , and is also applied to model a group of simultaneously tracked reindeer (Rangifer tarandus)

Momentum state engineering and control in Bose-Einstein condensates

We demonstrate theoretically the use of genetic learning algorithms to coherently control the dynamics of a Bose-Einstein condensate. We consider specifically the situation of a condensate in an optical lattice formed by two counterpropagating laser beams. The frequency detuning between the lasers acts as a control parameter that can be used to precisely manipulate the condensate even in the presence of a significant mean-field energy. We illustrate this procedure in the coherent acceleration of a condensate and in the preparation of a superposition of prescribed relative phase.Comment: 9 pages incl. 6 PostScript figures (.eps), LaTeX using RevTeX, submitted to Phys. Rev. A, incl. small modifications, some references adde

A-dependence of nuclear transparency in quasielastic A(e,e'p) at high Q^2

The A-dependence of the quasielastic A(e,e'p) reaction has been studied at SLAC with H-2, C, Fe, and Au nuclei at momentum transfers Q^2 = 1, 3, 5, and 6.8 (GeV/c)^2. We extract the nuclear transparency T(A,Q^2), a measure of the average probability that the struck proton escapes from the nucleus A without interaction. Several calculations predict a significant increase in T with momentum transfer, a phenomenon known as Color Transparency. No significant rise within errors is seen for any of the nuclei studied.Comment: 5 pages incl. 2 figures, Caltech preprint OAP-73

Determination of the neutron electric form factor in quasielastic scattering of polarized electrons from polarized 3He

We report a measurement of the asymmetry in spin-dependent quasielastic scattering of longitudinally polarized electrons from a polarized 3He gas target. The asymmetry is measured at kinematics sensitive to the transverse-longitudinal response function RTL(Q2,ω). The value of the neutron electric form factor GEn(Q2=0.16 (GeV/c2))=+0.070±0.100±0.035 is extracted from the asymmetry using a Faddeev calculation of the 3He wave function

Observation of a Coherence Length Effect in Exclusive Rho^0 Electroproduction

Exclusive incoherent electroproduction of the rho^0(770) meson from 1H, 2H, 3He, and 14N targets has been studied by the HERMES experiment at squared four-momentum transfer Q**2>0.4 GeV**2 and positron energy loss nu from 9 to 20 GeV. The ratio of the 14N to 1H cross sections per nucleon, known as the nuclear transparency, was found to decrease with increasing coherence length of quark-antiquark fluctuations of the virtual photon. The data provide clear evidence of the interaction of the quark- antiquark fluctuations with the nuclear medium.Comment: RevTeX, 5 pages, 3 figure

Determination of the Deep Inelastic Contribution to the Generalised Gerasimov-Drell-Hearn Integral for the Proton and Neutron

The virtual photon absorption cross section differences [sigma_1/2-sigma_3/2] for the proton and neutron have been determined from measurements of polarised cross section asymmetries in deep inelastic scattering of 27.5 GeV longitudinally polarised positrons from polarised 1H and 3He internal gas targets. The data were collected in the region above the nucleon resonances in the kinematic range nu < 23.5 GeV and 0.8 GeV**2 < Q**2 < 12 GeV**2. For the proton the contribution to the generalised Gerasimov-Drell-Hearn integral was found to be substantial and must be included for an accurate determination of the full integral. Furthermore the data are consistent with a QCD next-to-leading order fit based on previous deep inelastic scattering data. Therefore higher twist effects do not appear significant.Comment: 6 pages, 3 figures, 1 table, revte

Measurement of the Neutron Spin Structure Function g1ng_1^n with a Polarized ^3He Target

Results are reported from the HERMES experiment at HERA on a measurement of the neutron spin structure function $g_1^n(x,Q^2)$ in deep inelastic scattering using 27.5 GeV longitudinally polarized positrons incident on a polarized $^3$He internal gas target. The data cover the kinematic range $0.023<x<0.6$ and $1 (GeV/c)^2 < Q^2 <15 (GeV/c)^2$. The integral $\int_{0.023}^{0.6} g_1^n(x) dx$ evaluated at a fixed $Q^2$ of $2.5 (GeV/c)^2$ is $-0.034\pm 0.013(stat.)\pm 0.005(syst.)$. Assuming Regge behavior at low $x$, the first moment $\Gamma_1^n=\int_0^1 g_1^n(x) dx$ is $-0.037\pm 0.013(stat.)\pm 0.005(syst.)\pm 0.006(extrapol.)$.Comment: 4 pages TEX, text available at http://www.krl.caltech.edu/preprints/OAP.htm