-
Composite likelihood inference for space-time point processes
Authors:
Abdollah Jalilian,
Francisco Cuevas-Pacheco,
Ganggang Xu,
Rasmus Waagepetersen
Abstract:
The dynamics of a rain forest is extremely complex involving births, deaths and growth of trees with complex interactions between trees, animals, climate, and environment. We consider the patterns of recruits (new trees) and dead trees between rain forest censuses. For a current census we specify regression models for the conditional intensity of recruits and the conditional probabilities of death…
▽ More
The dynamics of a rain forest is extremely complex involving births, deaths and growth of trees with complex interactions between trees, animals, climate, and environment. We consider the patterns of recruits (new trees) and dead trees between rain forest censuses. For a current census we specify regression models for the conditional intensity of recruits and the conditional probabilities of death given the current trees and spatial covariates. We estimate regression parameters using conditional composite likelihood functions that only involve the conditional first order properties of the data. When constructing assumption lean estimators of covariance matrices of parameter estimates we only need mild assumptions of decaying conditional correlations in space while assumptions regarding correlations over time are avoided by exploiting conditional centering of composite likelihood score functions. Time series of point patterns from rain forest censuses are quite short while each point pattern covers a fairly big spatial region. To obtain asymptotic results we therefore use a central limit theorem for the fixed timespan - increasing spatial domain asymptotic setting. This also allows us to handle the challenge of using stochastic covariates constructed from past point patterns. Conveniently, it suffices to impose weak dependence assumptions on the innovations of the space-time process. We investigate the proposed methodology by simulation studies and applications to rain forest data.
△ Less
Submitted 19 February, 2024;
originally announced February 2024.
-
Pairwise interaction function estimation of Gibbs point processes using basis expansion
Authors:
Ismaïla Ba,
Jean-François Coeurjolly,
Francisco Cuevas-Pacheco
Abstract:
The class of Gibbs point processes (GPP) is a large class of spatial point processes able to model both clustered and repulsive point patterns. They are specified by their conditional intensity, which for a point pattern $\mathbf{x}$ and a location $u$, is roughly speaking the probability that an event occurs in an infinitesimal ball around $u$ given the rest of the configuration is $\mathbf{x}$.…
▽ More
The class of Gibbs point processes (GPP) is a large class of spatial point processes able to model both clustered and repulsive point patterns. They are specified by their conditional intensity, which for a point pattern $\mathbf{x}$ and a location $u$, is roughly speaking the probability that an event occurs in an infinitesimal ball around $u$ given the rest of the configuration is $\mathbf{x}$. The most simple and natural class of models is the class of pairwise interaction point processes where the conditional intensity depends on the number of points and pairwise distances between them. This paper is concerned with the problem of estimating the pairwise interaction function non parametrically. We propose to estimate it using an orthogonal series expansion of its logarithm. Such an approach has numerous advantages compared to existing ones. The estimation procedure is simple, fast and completely data-driven. We provide asymptotic properties such as consistency and asymptotic normality and show the efficiency of the procedure through simulation experiments and illustrate it with several datasets.
△ Less
Submitted 6 January, 2023; v1 submitted 11 October, 2021;
originally announced October 2021.
-
The $\mathcal{F}$-family of covariance functions: A Matérn analogue for modeling random fields on spheres
Authors:
Alfredo Alegría,
Francisco Cuevas-Pacheco,
Peter Diggle,
Emilio Porcu
Abstract:
The Mat{é}rn family of isotropic covariance functions has been central to the theoretical development and application of statistical models for geospatial data. For global data defined over the whole sphere representing planet Earth, the natural distance between any two locations is the great circle distance. In this setting, the Mat{é}rn family of covariance functions has a restriction on the smo…
▽ More
The Mat{é}rn family of isotropic covariance functions has been central to the theoretical development and application of statistical models for geospatial data. For global data defined over the whole sphere representing planet Earth, the natural distance between any two locations is the great circle distance. In this setting, the Mat{é}rn family of covariance functions has a restriction on the smoothness parameter, making it an unappealing choice to model smooth data. Finding a suitable analogue for modelling data on the sphere is still an open problem. This paper proposes a new family of isotropic covariance functions for random fields defined over the sphere. The proposed family has a parameter that indexes the mean square differentiability of the corresponding Gaussian field, and allows for any admissible range of fractal dimension. Our simulation study mimics the fixed domain asymptotic setting, which is the most natural regime for sampling on a closed and bounded set. As expected, our results support the analogous results (under the same asymptotic scheme) for planar processes that not all parameters can be estimated consistently. We apply the proposed model to a dataset of precipitable water content over a large portion of the Earth, and show that the model gives more precise predictions of the underlying process at unsampled locations than does the Mat{é}rn model using chordal distances.
△ Less
Submitted 13 January, 2021;
originally announced January 2021.
-
A convolution type model for the intensity of spatial point processes applied to eye-movement data
Authors:
Francisco Cuevas-Pacheco,
Jean-François Coeurjolly,
Marie-Hélène Descary
Abstract:
Estimating the first-order intensity function in point pattern analysis is an important problem, and it has been approached so far from different perspectives: parametrically, semiparametrically or nonparametrically. Our approach is close to a semiparametric one. Motivated by eye-movement data, we introduce a convolution type model where the log-intensity is modelled as the convolution of a functi…
▽ More
Estimating the first-order intensity function in point pattern analysis is an important problem, and it has been approached so far from different perspectives: parametrically, semiparametrically or nonparametrically. Our approach is close to a semiparametric one. Motivated by eye-movement data, we introduce a convolution type model where the log-intensity is modelled as the convolution of a function $β(\cdot)$, to be estimated, and a single spatial covariate (the image an individual is looking at for eye-movement data). Based on a Fourier series expansion, we show that the proposed model \rev{can be viewed as a} log-linear model with an infinite number of coefficients, which correspond to the spectral decomposition of $β(\cdot)$. After truncation, we estimate these coefficients through a penalized Poisson likelihood. We illustrate the efficiency of the proposed methodology on simulated data and on eye-movement data.
△ Less
Submitted 5 October, 2021; v1 submitted 17 December, 2020;
originally announced December 2020.
-
Karhunen-Loève Expansions for Axially Symmetric Gaussian Processes: Modeling Strategies and $L^2$ Approximations
Authors:
Alfredo Alegría,
Francisco Cuevas-Pacheco
Abstract:
Axially symmetric processes on spheres, for which the second-order dependency structure may substantially vary with shifts in latitude, are a prominent alternative to model the spatial uncertainty of natural variables located over large portions of the Earth. In this paper, we focus on Karhunen-Loève expansions of axially symmetric Gaussian processes. First, we investigate a parametric family of K…
▽ More
Axially symmetric processes on spheres, for which the second-order dependency structure may substantially vary with shifts in latitude, are a prominent alternative to model the spatial uncertainty of natural variables located over large portions of the Earth. In this paper, we focus on Karhunen-Loève expansions of axially symmetric Gaussian processes. First, we investigate a parametric family of Karhunen-Loève coefficients that allows for versatile spatial covariance functions. The isotropy as well as the longitudinal independence can be obtained as limit cases of our proposal. Second, we introduce a strategy to render any longitudinally reversible process irreversible, which means that its covariance function could admit certain types of asymmetries along longitudes. Then, finitely truncated Karhunen-Loève expansions are used to approximate axially symmetric processes. For such approximations, bounds for the $L^2$-error are provided. Numerical experiments are conducted to illustrate our findings.
△ Less
Submitted 3 July, 2020;
originally announced July 2020.
-
Regularized estimation for highly multivariate log Gaussian Cox processes
Authors:
Achmad Choiruddin,
Francisco Cuevas-Pacheco,
Jean-François Coeurjolly,
Rasmus Waagepetersen
Abstract:
Statistical inference for highly multivariate point pattern data is challenging due to complex models with large numbers of parameters. In this paper, we develop numerically stable and efficient parameter estimation and model selection algorithms for a class of multivariate log Gaussian Cox processes. The methodology is applied to a highly multivariate point pattern data set from tropical rain for…
▽ More
Statistical inference for highly multivariate point pattern data is challenging due to complex models with large numbers of parameters. In this paper, we develop numerically stable and efficient parameter estimation and model selection algorithms for a class of multivariate log Gaussian Cox processes. The methodology is applied to a highly multivariate point pattern data set from tropical rain forest ecology.
△ Less
Submitted 4 May, 2019;
originally announced May 2019.
-
Second-order variational equations for spatial point processes with a view to pair correlation function estimation
Authors:
Jean-François Coeurjolly,
Francisco Cuevas-Pacheco,
Rasmus Waagepetersen
Abstract:
Second-order variational type equations for spatial point processes are established. In case of log linear parametric models for pair correlation functions, it is demonstrated that the variational equations can be applied to construct estimating equations with closed form solutions for the parameter estimates. This result is used to fit orthogonal series expansions of log pair correlation function…
▽ More
Second-order variational type equations for spatial point processes are established. In case of log linear parametric models for pair correlation functions, it is demonstrated that the variational equations can be applied to construct estimating equations with closed form solutions for the parameter estimates. This result is used to fit orthogonal series expansions of log pair correlation functions of general form.
△ Less
Submitted 15 January, 2019;
originally announced January 2019.
-
Structured space-sphere point processes and $K$-functions
Authors:
Jesper Møller,
Heidi S. Christensen,
Francisco Cuevas-Pacheco,
Andreas D. Christoffersen
Abstract:
This paper concerns space-sphere point processes, that is, point processes on the product space of $\mathbb R^d$ (the $d$-dimensional Euclidean space) and $\mathbb S^k$ (the $k$-dimen\-sional sphere). We consider specific classes of models for space-sphere point processes, which are adaptations of existing models for either spherical or spatial point processes. For model checking or fitting, we pr…
▽ More
This paper concerns space-sphere point processes, that is, point processes on the product space of $\mathbb R^d$ (the $d$-dimensional Euclidean space) and $\mathbb S^k$ (the $k$-dimen\-sional sphere). We consider specific classes of models for space-sphere point processes, which are adaptations of existing models for either spherical or spatial point processes. For model checking or fitting, we present the space-sphere $K$-function which is a natural extension of the inhomogeneous $K$-function for point processes on $\mathbb R^d$ to the case of space-sphere point processes. Under the assumption that the intensity and pair correlation function both have a certain separable structure, the space-sphere $K$-function is shown to be proportional to the product of the inhomogeneous spatial and spherical $K$-functions. For the presented space-sphere point process models, we discuss cases where such a separable structure can be obtained. The usefulness of the space-sphere $K$-function is illustrated for real and simulated datasets with varying dimensions $d$ and $k$.
△ Less
Submitted 29 October, 2019; v1 submitted 21 December, 2018;
originally announced December 2018.
-
Log Gaussian Cox processes on the sphere
Authors:
Jesper Møller,
Francisco Cuevas-Pacheco
Abstract:
A log Gaussian Cox process (LGCP) is a doubly stochastic construction consisting of a Poisson point process with a random log-intensity given by a Gaussian random field. Statistical methodology have mainly been developed for LGCPs defined in the $d$-dimensional Euclidean space. This paper concerns the case of LGCPs on the $d$-dimensional sphere, with $d=2$ of primary interest. We discuss the exist…
▽ More
A log Gaussian Cox process (LGCP) is a doubly stochastic construction consisting of a Poisson point process with a random log-intensity given by a Gaussian random field. Statistical methodology have mainly been developed for LGCPs defined in the $d$-dimensional Euclidean space. This paper concerns the case of LGCPs on the $d$-dimensional sphere, with $d=2$ of primary interest. We discuss the existence problem of such LGCPs, provide sufficient existence conditions, and establish further useful theoretical properties. The results are applied for the description of sky positions of galaxies, in comparison with previous analysis based on a Thomas process, using simple estimation procedures and making a careful model checking. We account for inhomogeneity in our models, and as the model checking is based on a thinning procedure which produces homogeneous/isotropic LGCPs, we discuss its sensitivity.
△ Less
Submitted 5 May, 2018; v1 submitted 8 March, 2018;
originally announced March 2018.