Source author record

Jesper Møller

Jesper Møller appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

18 published item(s)

preprint2024arXiv

Coupling results and Markovian structures for number representations of continuous random variables

A general setting for nested subdivisions of a bounded real set into intervals defining the digits $X_1,X_2,...$ of a random variable $X$ with a probability density function $f$ is considered. Under the weak condition that $f$ is almost everywhere lower semi-continuous, a coupling between $X$ and a non-negative integer-valued random variable $N$ is established so that $X_1,...,X_N$ have an interpretation as the ``sufficient digits'', since the distribution of $R=(X_{N+1},X_{N+2},...)$ conditioned on $S=(X_1,...,X_N)$ does not depend on $f$. Adding a condition about a Markovian structure of the lengths of the intervals in the nested subdivisions, $R\,|\,S$ becomes a Markov chain of a certain order $s\ge0$. If $s=0$ then $X_{N+1},X_{N+2},...$ are IID with a known distribution. When $s>0$ and the Markov chain is uniformly geometric ergodic, a coupling is established between $(X,N)$ and a random time $M$ so that the chain after time $\max\{N,s\}+M-s$ is stationary and $M$ follows a simple known distribution. The results are related to several examples of number representations generated by a dynamical system, including base-$q$ expansions, generalized Lüroth series, $β$-expansions, and continued fraction representations. The importance of the results and some suggestions and open problems for future research are discussed.

preprint2022arXiv

Characterization of random variables with stationary digits

Let $q\ge2$ be an integer, $\{X_n\}_{n\geq 1}$ a stochastic process with state space $\{0,\ldots,q-1\}$, and $F$ the cumulative distribution function (CDF) of $\sum_{n=1}^\infty X_n q^{-n}$. We show that stationarity of $\{X_n\}_{n\geq 1}$ is equivalent to a functional equation obeyed by $F$ and use this to characterize the characteristic function of $X$ and the structure of $F$ in terms of its Lebesgue decomposition. More precisely, while the absolutely continuous component of $F$ can only be the uniform distribution on the unit interval, its discrete component can only be a countable convex combination of certain explicitly computable CDFs for probability distributions with finite support. We also show that $\mathrm{d} F$ is a Rajchman measure if and only if $F $ is the uniform CDF on $[0,1]$.

preprint2022arXiv

Determinantal shot noise Cox processes

We present a new class of cluster point process models, which we call determinantal shot noise Cox processes (DSNCP), with repulsion between cluster centres. They are the special case of generalized shot noise Cox processes where the cluster centres are determinantal point processes. We establish various moment results and describe how these can be used to easily estimate unknown parameters in two particularly tractable cases, namely when the offspring density is isotropic Gaussian and the kernel of the determinantal point process of cluster centres is Gaussian or like in a scaled Ginibre point process. Through a simulation study and the analysis of a real point pattern data set we see that when modelling clustered point patterns, a much lower intensity of cluster centres may be needed in DSNCP models as compared to shot noise Cox processes.

preprint2022arXiv

Fitting three-dimensional Laguerre tessellations by hierarchical marked point process models

We present a general statistical methodology for analysing a Laguerre tessellation data set viewed as a realization of a marked point process model. In the first step, for the points we use a nested sequence of multiscale processes which constitute a flexible parametric class of pairwise interaction point process models. In the second step, for the marks/radii conditioned on the points we consider various exponential family models where the canonical sufficient statistic is based on tessellation characteristics. For each step parameter estimation based on maximum pseudolikelihood methods is tractable. Model checking is performed using global envelopes and corresponding tests in the first step and by comparing observed and simulated tessellation characteristics in the second step. We apply our methodology for a 3D Laguerre tessellation data set representing the microstructure of a polycrystalline metallic material, where simulations under a fitted model may substitute expensive laboratory experiments.

preprint2022arXiv

Should we condition on the number of points when modelling spatial point patterns?

We discuss the practice of directly or indirectly assuming a model for the number of points when modelling spatial point patterns even though it is rarely possible to validate such a model in practice since most point pattern data consist of only one pattern. We therefore explore the possibility to condition on the number of points instead when fitting and validating spatial point process models. In a simulation study with different popular spatial point process models, we consider model validation using global envelope tests based on functional summary statistics. We find that conditioning on the number of points will for some functional summary statistics lead to more narrow envelopes and that it can also be useful for correcting for some conservativeness in the tests when testing composite hypothesis. However, for other functional summary statistics, it makes little or no difference to condition on the number of points. When estimating parameters in popular spatial point process models, we conclude that for mathematical and computational reasons it is convenient to assume a distribution for the number of points.

preprint2020arXiv

Approximate Bayesian inference for a spatial point process model exhibiting regularity and random aggregation

In this paper, we propose a doubly stochastic spatial point process model with both aggregation and repulsion. This model combines the ideas behind Strauss processes and log Gaussian Cox processes. The likelihood for this model is not expressible in closed form but it is easy to simulate realisations under the model. We therefore explain how to use approximate Bayesian computation (ABC) to carry out statistical inference for this model. We suggest a method for model validation based on posterior predictions and global envelopes. We illustrate the ABC procedure and model validation approach using both simulated point patterns and a real data example.

preprint2016arXiv

A tutorial on Palm distributions for spatial point processes

This tutorial provides an introduction to Palm distributions for spatial point processes. Initially, in the context of finite point processes , we give an explicit definition of Palm distributions in terms of their density functions. Then we review Palm distributions in the general case. Finally we discuss some examples of Palm distributions for specific models and some applications.

preprint2016arXiv

Determinantal point process models on the sphere

We consider determinantal point processes on the $d$-dimensional unit sphere $\mathbb S^d$. These are finite point processes exhibiting repulsiveness and with moment properties determined by a certain determinant whose entries are specified by a so-called kernel which we assume is a complex covariance function defined on $\mathbb S^d\times\mathbb S^d$. We review the appealing properties of such processes, including their specific moment properties, density expressions and simulation procedures. Particularly, we characterize and construct isotropic DPPs models on $\mathbb{S}^d$, where it becomes essential to specify the eigenvalues and eigenfunctions in a spectral representation for the kernel, and we figure out how repulsive isotropic DPPs can be. Moreover, we discuss the shortcomings of adapting existing models for isotropic covariance functions and consider strategies for developing new models, including a useful spectral approach.

preprint2016arXiv

Some recent developments in statistics for spatial point patterns

This paper reviews developments in statistics for spatial point processes obtained within roughly the last decade. These developments include new classes of spatial point process models such as determinantal point processes, models incorporating both regularity and aggregation, and models where points are randomly distributed around latent geometric structures. Regarding parametric inference the main focus is on various types of estimating functions derived from so-called innovation measures. Optimality of such estimating functions is discussed as well as computational issues. Maximum likelihood inference for determinantal point processes and Bayesian inference are briefly considered too. Concerning non-parametric inference, we consider extensions of functional summary statistics to the case of inhomogeneous point processes as well as new approaches to simulation based inference.

preprint2016arXiv

The cylindrical K-function and Poisson line cluster point processes

Analyzing point patterns with linear structures has recently been of interest in e.g. neuroscience and geography. To detect anisotropy in such cases, we introduce a functional summary statistic, called the cylindrical $K$-function, since it is a directional $K$-function whose structuring element is a cylinder. Further we introduce a class of anisotropic Cox point processes, called Poisson line cluster point processes. The points of such a process are random displacements of Poisson point processes defined on the lines of a Poisson line process. Parameter estimation based on moment methods or Bayesian inference for this model is discussed when the underlying Poisson line process and the cluster memberships are treated as hidden processes. To illustrate the methodologies, we analyze a two and a three-dimensional point pattern data set. The 3D data set is of particular interest as it relates to the minicolumn hypothesis in neuroscience, claiming that pyramidal and other brain cells have a columnar arrangement perpendicular to the pial surface of the brain.

preprint2015arXiv

Modelling aggregation on the large scale and regularity on the small scale in spatial point pattern datasets

We consider a dependent thinning of a regular point process with the aim of obtaining aggregation on the large scale and regularity on the small scale in the resulting target point process of retained points. Various parametric models for the underlying processes are suggested and the properties of the target point process are studied. Simulation and inference procedures are discussed when a realization of the target point process is observed, depending on whether the thinned points are observed or not. The paper extends previous work by Dietrich Stoyan on interrupted point processes.

preprint2014arXiv

Determinantal point process models and statistical inference : Extended version

Statistical models and methods for determinantal point processes (DPPs) seem largely unexplored. We demonstrate that DPPs provide useful models for the description of spatial point pattern datasets where nearby points repel each other. Such data are usually modelled by Gibbs point processes, where the likelihood and moment expressions are intractable and simulations are time consuming. We exploit the appealing probabilistic properties of DPPs to develop parametric models, where the likelihood and moment expressions can be easily evaluated and realizations can be quickly simulated. We discuss how statistical inference is conducted using the likelihood or moment properties of DPP models, and we provide freely available software for simulation and statistical inference.

preprint2014arXiv

Fingerprint Analysis with Marked Point Processes

We present a framework for fingerprint matching based on marked point process models. An efficient Monte Carlo algorithm is developed to calculate the marginal likelihood ratio for the hypothesis that two observed prints originate from the same finger against the hypothesis that they originate from different fingers. Our model achieves good performance on an NIST-FBI fingerprint database of 258 matched fingerprint pairs.

preprint2014arXiv

Variational approach for spatial point process intensity estimation

We introduce a new variational estimator for the intensity function of an inhomogeneous spatial point process with points in the $d$-dimensional Euclidean space and observed within a bounded region. The variational estimator applies in a simple and general setting when the intensity function is assumed to be of log-linear form $β+{θ}^{\top}z(u)$ where $z$ is a spatial covariate function and the focus is on estimating ${θ}$. The variational estimator is very simple to implement and quicker than alternative estimation procedures. We establish its strong consistency and asymptotic normality. We also discuss its finite-sample properties in comparison with the maximum first order composite likelihood estimator when considering various inhomogeneous spatial point process models and dimensions as well as settings were $z$ is completely or only partially known.

preprint2012arXiv

Score, Pseudo-Score and Residual Diagnostics for Spatial Point Process Models

We develop new tools for formal inference and informal model validation in the analysis of spatial point pattern data. The score test is generalized to a "pseudo-score" test derived from Besag's pseudo-likelihood, and to a class of diagnostics based on point process residuals. The results lend theoretical support to the established practice of using functional summary statistics, such as Ripley's $K$-function, when testing for complete spatial randomness; and they provide new tools such as the compensator of the $K$-function for testing other fitted models. The results also support localization methods such as the scan statistic and smoothed residual plots. Software for computing the diagnostics is provided.