BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Wilfrid S. Kendall

Wilfrid S. Kendall appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.PRApplicationsMethodologyComputationmath.STPerformanceStatistics Theory

Catalog footprint

What is connected

15works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2016arXiv

Coupling the Kolmogorov Diffusion: maximality and efficiency considerations

This is a case study concerning the rate at which probabilistic coupling occurs for nilpotent diffusions. We focus on the simplest case of Kolmogorov diffusion (Brownian motion together with its time integral, or, slightly more generally, together with a finite number of iterated time integrals). In this case there can be no Markovian maximal coupling. Indeed, Markovian couplings cannot even be efficient (extending the terminology of Burdzy and Kendall, Efficient Markovian couplings: examples and counterexamples; Annals of Applied Probability, 2000). Finally, at least in the classical case of a single time integral, it is not possible to choose a Markovian coupling that is optimal in the sense of simultaneously minimizing the probability of failing to couple by time t for all positive t. In recompense for all these negative results, we exhibit a simple efficient non-Markovian coupling.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rigidity for Markovian Maximal Couplings of Elliptic Diffusions

Maximal couplings are (probabilistic) couplings of Markov processes such that the tail probabilities of the coupling time attain the total variation lower bound (Aldous bound) uniformly for all time. Markovian (or immersion) couplings are couplings defined by strategies where neither process is allowed to look into the future of the other before making the next transition. Markovian couplings are typically easier to construct and analyze than general couplings, and play an important role in many branches of probability and analysis. Hsu and Sturm (2013) proved that the reflection-coupling of Brownian motion is the unique Markovian maximal coupling (MMC) of Brownian motions starting from two different points. Later, Kuwada (2009) proved that the existence of a MMC for Brownian motions on a Riemannian manifold enforces existence of a reflection structure on the manifold. In this work, we investigate suitably regular elliptic diffusions on manifolds, and show how consideration of the diffusion geometry (including dimension of the isometry group and flows of isometries) is fundamental in classification of the space and the generator of the diffusion for which an MMC exists, especially when the MMC also holds under local perturbations of the starting points for the coupled diffusions. We also describe such diffusions in terms of Killing vectorfields (generators of isometry groups) and dilation vectorfields (generators of scaling symmetry groups).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coupling, local times, immersions

This paper answers a question of Émery [In Séminaire de Probabilités XLII (2009) 383-396 Springer] by constructing an explicit coupling of two copies of the Beneš et al. [In Applied Stochastic Analysis (1991) 121-156 Gordon & Breach] diffusion (BKR diffusion), neither of which starts at the origin, and whose natural filtrations agree. The paper commences by surveying probabilistic coupling, introducing the formal definition of an immersed coupling (the natural filtration of each component is immersed in a common underlying filtration; such couplings have been described as co-adapted or Markovian in older terminologies) and of an equi-filtration coupling (the natural filtration of each component is immersed in the filtration of the other; consequently the underlying filtration is simultaneously the natural filtration for each of the two coupled processes). This survey is followed by a detailed case-study of the simpler but potentially thematic problem of coupling Brownian motion together with its local time at $0$. This problem possesses its own intrinsic interest as well as being closely related to the BKR coupling construction. Attention focusses on a simple immersed (co-adapted) coupling, namely the reflection/synchronized coupling. It is shown that this coupling is optimal amongst all immersed couplings of Brownian motion together with its local time at $0$, in the sense of maximizing the coupling probability at all possible times, at least when not started at pairs of initial points lying in a certain singular set. However numerical evidence indicates that the coupling is not a maximal coupling, and is a simple but non-trivial instance for which this distinction occurs. It is shown how the reflection/synchronized coupling can be converted into a successful equi-filtration coupling, by modifying the coupling using a deterministic time-delay and then by concatenating an infinite sequence of such modified couplings. The construction of an explicit equi-filtration coupling of two copies of the BKR diffusion follows by a direct generalization, although the proof of success for the BKR coupling requires somewhat more analysis than in the local time case.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Perfect Simulation of $M/G/c$ Queues

In this paper we describe a perfect simulation algorithm for the stable $M/G/c$ queue. Sigman (2011: Exact Simulation of the Stationary Distribution of the FIFO M/G/c Queue. Journal of Applied Probability, 48A, 209--213) showed how to build a dominated CFTP algorithm for perfect simulation of the super-stable $M/G/c$ queue operating under First Come First Served discipline, with dominating process provided by the corresponding $M/G/1$ queue (using Wolff's sample path monotonicity, which applies when service durations are coupled in order of initiation of service), and exploiting the fact that the workload process for the $M/G/1$ queue remains the same under different queueing disciplines, in particular under the Processor Sharing discipline, for which a dynamic reversibility property holds. We generalize Sigman's construction to the stable case by comparing the $M/G/c$ queue to a copy run under Random Assignment. This allows us to produce a naive perfect simulation algorithm based on running the dominating process back to the time it first empties. We also construct a more efficient algorithm that uses sandwiching by lower and upper processes constructed as coupled $M/G/c$ queues started respectively from the empty state and the state of the $M/G/c$ queue under Random Assignment. A careful analysis shows that appropriate ordering relationships can still be maintained, so long as service durations continue to be coupled in order of initiation of service. We summarize statistical checks of simulation output, and demonstrate that the mean run-time is finite so long as the second moment of the service duration distribution is finite.

0 citations0 reviewsNo code linkedOpen access
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