Barycentres and Hurricane Trajectories
The use of barycentres in data analysis is illustrated, using as example a dataset of hurricane trajectories.
Researcher profile
Wilfrid S. Kendall
Wilfrid S. Kendall contributes to research discovery and scholarly infrastructure.
ResearcherAffiliation not importedOpen to collaborate
Trust snapshot
Quick read
Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
6topics
4close collaborators
Actions
Decide how to stay connected
Identity and collaboration
How to connect with this researcher
Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.
Log in to claimDirect collaboration
Open a focused conversation when the fit is right
Claim this author entity first to unlock direct invitations.
Research graph
See the researcher in context
Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.
Building this graph slice
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
11 published item(s)
preprint2014arXiv
From Random Lines to Metric Spaces
Consider an improper Poisson line process, marked by positive speeds so as to satisfy a scale-invariance property (actually, scale-equivariance). The line process can be characterized by its intensity measure, which belongs to a one-parameter family if scale and Euclidean invariance are required. This paper investigates a proposal by Aldous, namely that the line process could be used to produce a scale-invariant random spatial network (SIRSN) by means of connecting up points using paths which follow segments from the line process at the stipulated speeds. It is shown that this does indeed produce a scale-invariant network, under suitable conditions on the parameter; indeed that this produces a parameter-dependent random geodesic metric for d-dimensional space ($d\geq2$), where geodesics are given by minimum-time paths. Moreover in the planar case it is shown that the resulting geodesic metric space has an almost-everywhere-unique-geodesic property, that geodesics are locally of finite mean length, and that if an independent Poisson point process is connected up by such geodesics then the resulting network places finite length in each compact region. It is an open question whether the result is a SIRSN (in Aldous' sense; so placing finite mean length in each compact region), but it may be called a pre-SIRSN.
0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv
Return to the Poissonian City
Consider the following random spatial network: in a large disk, construct a network using a stationary and isotropic Poisson line process of unit intensity. Connect pairs of points using the network, with initial / final segments of the connecting path formed by travelling off the network in the opposite direction to that of the destination / source. Suppose further that connections are established using "near-geodesics", constructed between pairs of points using the perimeter of the cell containing these two points and formed using only the Poisson lines not separating them. If each pair of points generates an infinitesimal amount of traffic divided equally between the two connecting near-geodesics, and if the Poisson line pattern is conditioned to contain a line through the centre, then what can be said about the total flow through the centre? In earlier work ("Geodesics and flows in a Poissonian city", Annals of Applied Probability, 21(3), 801--842, 2011) it was shown that a scaled version of this flow had asymptotic distribution given by the 4-volume of a region in 4-space, constructed using an improper anisotropic Poisson line process in an infinite planar strip. Here we construct a more amenable representation in terms of two "seminal curves" defined by the improper Poisson line process, and establish results which produce a framework for effective simulation from this distribution up to an L1 error which tends to zero with increasing computational effort.
0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv
Dynamic filtering of static dipoles in magnetoencephalography
We consider the problem of estimating neural activity from measurements of the magnetic fields recorded by magnetoencephalography. We exploit the temporal structure of the problem and model the neural current as a collection of evolving current dipoles, which appear and disappear, but whose locations are constant throughout their lifetime. This fully reflects the physiological interpretation of the model. In order to conduct inference under this proposed model, it was necessary to develop an algorithm based around state-of-the-art sequential Monte Carlo methods employing carefully designed importance distributions. Previous work employed a bootstrap filter and an artificial dynamic structure where dipoles performed a random walk in space, yielding nonphysical artefacts in the reconstructions; such artefacts are not observed when using the proposed model. The algorithm is validated with simulated data, in which it provided an average localisation error which is approximately half that of the bootstrap filter. An application to complex real data derived from a somatosensory experiment is presented. Assessment of model fit via marginal likelihood showed a clear preference for the proposed model and the associated reconstructions show better localisation.
0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv
Fibre-generated point processes and fields of orientations
This paper introduces a new approach to analyzing spatial point data clustered along or around a system of curves or "fibres." Such data arise in catalogues of galaxy locations, recorded locations of earthquakes, aerial images of minefields and pore patterns on fingerprints. Finding the underlying curvilinear structure of these point-pattern data sets may not only facilitate a better understanding of how they arise but also aid reconstruction of missing data. We base the space of fibres on the set of integral lines of an orientation field. Using an empirical Bayes approach, we estimate the field of orientations from anisotropic features of the data. We then sample from the posterior distribution of fibres, exploring models with different numbers of clusters, fitting fibres to the clusters as we proceed. The Bayesian approach permits inference on various properties of the clusters and associated fibres, and the results perform well on a number of very different curvilinear structures.
0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv
Geodesics and flows in a Poissonian city
The stationary isotropic Poisson line process was used to derive upper bounds on mean excess network geodesic length in Aldous and Kendall [Adv. in Appl. Probab. 40 (2008) 1-21]. The current paper presents a study of the geometry and fluctuations of near-geodesics in the network generated by the line process. The notion of a "Poissonian city" is introduced, in which connections between pairs of nodes are made using simple "no-overshoot" paths based on the Poisson line process. Asymptotics for geometric features and random variation in length are computed for such near-geodesic paths; it is shown that they traverse the network with an order of efficiency comparable to that of true network geodesics. Mean characteristics and limiting behavior at the center are computed for a natural network flow. Comparisons are drawn with similar network flows in a city based on a comparable rectilinear grid. A concluding section discusses several open problems.
0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv
Limit theorems for empirical Fréchet means of independent and non-identically distributed manifold-valued random variables
We prove weak laws of large numbers and central limit theorems of Lindeberg type for empirical centres of mass (empirical Fréchet means) of independent non-identically distributed random variables taking values in Riemannian manifolds. In order to prove these theorems we describe and prove a simple kind of Lindeberg-Feller central approximation theorem for vector-valued random variables, which may be of independent interest and is therefore the subject of a self-contained section. This vector-valued result allows us to clarify the number of conditions required for the central limit theorem for empirical Fréchet means, while extending its scope.
0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv
Coupling time distribution asymptotics for some couplings of the Levy stochastic area
We exhibit some explicit co-adapted couplings for n-dimensional Brownian motion and all its Levy stochastic areas. In the two-dimensional case we show how to derive exact asymptotics for the coupling time under various mixed coupling strategies, using Dufresne's formula for the distribution of exponential functionals of Brownian motion. This yields quantitative asymptotics for the distributions of random times required for certain simultaneous couplings of stochastic area and Brownian motion. The approach also applies to higher dimensions, but will then lead to upper and lower bounds rather than exact asymptotics.
0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv
Brownian couplings, convexity, and shy-ness
Benjamini, Burdzy and Chen (2007) introduced the notion of a shy coupling: a coupling of a Markov process such that, for suitable starting points, there is a positive chance of the two component processes of the coupling staying a positive distance away from each other for all time. Among other results, they showed no shy couplings could exist for reflected Brownian motions in C^2 bounded convex planar domains whose boundaries contain no line segments. Here we use potential-theoretic methods to extend this Benjamini et al. result (a) to all bounded convex domains (whether planar and smooth or not) whose boundaries contain no line segments, (b) to all bounded convex planar domains regardless of further conditions on the boundary.
0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv
Coupling all the Lévy stochastic areas of multidimensional Brownian motion
It is shown how to construct a successful co-adapted coupling of two copies of an $n$-dimensional Brownian motion $(B_1,...,B_n)$ while simultaneously coupling all corresponding copies of Lévy stochastic areas $\int B_i dB_j-\int B_j dB_i$. It is conjectured that successful co-adapted couplings still exist when the Lévy stochastic areas are replaced by a finite set of multiply iterated path- and time-integrals, subject to algebraic compatibility of the initial conditions.
0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv
Perfect simulation for a class of positive recurrent Markov chains
This paper generalizes the work of Kendall [Electron. Comm. Probab. 9 (2004) 140--151], which showed that perfect simulation, in the form of dominated coupling from the past, is always possible (although not necessarily practical) for geometrically ergodic Markov chains. Here, we consider the more general situation of positive recurrent chains and explore when it is possible to produce such a simulation algorithm for these chains. We introduce a class of chains which we name tame, for which we show that perfect simulation is possible.
0 citations0 reviewsNo code linkedOpen access