-
Function-valued marked spatial point processes on linear networks: application to urban cycling profiles
Authors:
Matthias Eckardt,
Jorge Mateu,
Mehdi Moradi
Abstract:
In the literature on spatial point processes, there is an emerging challenge in studying marked point processes with points being labelled by functions. In this paper, we focus on point processes living on linear networks and, from distinct points of view, propose several marked summary characteristics that are of great use in studying the average association and dispersion of the function-valued…
▽ More
In the literature on spatial point processes, there is an emerging challenge in studying marked point processes with points being labelled by functions. In this paper, we focus on point processes living on linear networks and, from distinct points of view, propose several marked summary characteristics that are of great use in studying the average association and dispersion of the function-valued marks. Through a simulation study, we evaluate the performance of our proposed marked summary characteristics, both when marks are independent and when some sort of spatial dependence is evident among them. Finally, we employ our proposed mark summary characteristics to study the spatial structure of urban cycling profiles in Vancouver, Canada.
△ Less
Submitted 10 July, 2024;
originally announced July 2024.
-
On the semi-additivity of the $1/2$-symmetric caloric capacity in the plane
Authors:
Joan Hernández,
Joan Mateu,
Laura Prat
Abstract:
In this paper we study properties of a variant of the $1/2$-caloric capacity, called $1/2$-symmetric caloric capacity. The latter is associated simultaneously with the $1/2$-fractional heat equation and its conjugate. We establish its semi-additivity in $\mathbb{R}^2$ and, moreover, we compute explicitly the $1/2$-symmetric caloric capacity of rectangles, which illustrates its anisotropic behaviou…
▽ More
In this paper we study properties of a variant of the $1/2$-caloric capacity, called $1/2$-symmetric caloric capacity. The latter is associated simultaneously with the $1/2$-fractional heat equation and its conjugate. We establish its semi-additivity in $\mathbb{R}^2$ and, moreover, we compute explicitly the $1/2$-symmetric caloric capacity of rectangles, which illustrates its anisotropic behaviour.
△ Less
Submitted 2 May, 2024;
originally announced May 2024.
-
Semi-parametric profile pseudolikelihood via local summary statistics for spatial point pattern intensity estimation
Authors:
Nicoletta D'Angelo,
Giada Adelfio,
Jorge Mateu,
Ottmar Cronie
Abstract:
Second-order statistics play a crucial role in analysing point processes. Previous research has specifically explored locally weighted second-order statistics for point processes, offering diagnostic tests in various spatial domains. However, there remains a need to improve inference for complex intensity functions, especially when the point process likelihood is intractable and in the presence of…
▽ More
Second-order statistics play a crucial role in analysing point processes. Previous research has specifically explored locally weighted second-order statistics for point processes, offering diagnostic tests in various spatial domains. However, there remains a need to improve inference for complex intensity functions, especially when the point process likelihood is intractable and in the presence of interactions among points. This paper addresses this gap by proposing a method that exploits local second-order characteristics to account for local dependencies in the fitting procedure. Our approach utilises the Papangelou conditional intensity function for general Gibbs processes, avoiding explicit assumptions about the degree of interaction and homogeneity. We provide simulation results and an application to real data to assess the proposed method's goodness-of-fit. Overall, this work contributes to advancing statistical techniques for point process analysis in the presence of spatial interactions.
△ Less
Submitted 16 April, 2024;
originally announced April 2024.
-
Analiticity of the Flow for the Aggregation Equation
Authors:
J. M. Burgués,
J. Mateu
Abstract:
Let $Ω$ be a bounded domain in $\R^n$ whose boundary is $\ka{1,\,γ}$ for $γ\in(0,\,1)$. Consider the aggregation equation in the case of the initial condition being a positive multiple of the characteristic function of $Ω$. In this paper we prove global in time analyticity of the flow generated by the velocity field which propagates the density solution of this equation.
Let $Ω$ be a bounded domain in $\R^n$ whose boundary is $\ka{1,\,γ}$ for $γ\in(0,\,1)$. Consider the aggregation equation in the case of the initial condition being a positive multiple of the characteristic function of $Ω$. In this paper we prove global in time analyticity of the flow generated by the velocity field which propagates the density solution of this equation.
△ Less
Submitted 7 February, 2024; v1 submitted 9 November, 2023;
originally announced November 2023.
-
Summary characteristics for multivariate function-valued spatial point process attributes
Authors:
Matthias Eckardt,
Carles Comas,
Jorge Mateu
Abstract:
Prompted by modern technologies in data acquisition, the statistical analysis of spatially distributed function-valued quantities has attracted a lot of attention in recent years. In particular, combinations of functional variables and spatial point processes yield a highly challenging instance of such modern spatial data applications. Indeed, the analysis of spatial random point configurations, w…
▽ More
Prompted by modern technologies in data acquisition, the statistical analysis of spatially distributed function-valued quantities has attracted a lot of attention in recent years. In particular, combinations of functional variables and spatial point processes yield a highly challenging instance of such modern spatial data applications. Indeed, the analysis of spatial random point configurations, where the point attributes themselves are functions rather than scalar-valued quantities, is just in its infancy, and extensions to function-valued quantities still remain limited. In this view, we extend current existing first- and second-order summary characteristics for real-valued point attributes to the case where in addition to every spatial point location a set of distinct function-valued quantities are available. Providing a flexible treatment of more complex point process scenarios, we build a framework to consider points with multivariate function-valued marks, and develop sets of different cross-function (cross-type and also multi-function cross-type) versions of summary characteristics that allow for the analysis of highly demanding modern spatial point process scenarios. We consider estimators of the theoretical tools and analyse their behaviour through a simulation study and two real data applications.
△ Less
Submitted 11 July, 2023;
originally announced July 2023.
-
Time periodic solutions close to localized radial monotone profiles for the 2D Euler equations
Authors:
Claudia García,
Taoufik Hmidi,
Joan Mateu
Abstract:
In this paper, we address for the 2D Euler equations the existence of rigid time periodic solutions close to stationary radial vortices of type $f_0(|x|){\bf 1}_{\mathbb{D}}(x)$, with $\mathbb{D}$ the unit disc and $f_0$ being a strictly monotonic profile with constant sign. We distinguish two scenarios according to the sign of the profile: defocusing and focusing. In the first regime, we have sca…
▽ More
In this paper, we address for the 2D Euler equations the existence of rigid time periodic solutions close to stationary radial vortices of type $f_0(|x|){\bf 1}_{\mathbb{D}}(x)$, with $\mathbb{D}$ the unit disc and $f_0$ being a strictly monotonic profile with constant sign. We distinguish two scenarios according to the sign of the profile: defocusing and focusing. In the first regime, we have scarcity of the bifurcating curves associated with lower symmetry. However in the focusing case we get a countable family of bifurcating solutions associated with large symmetry. The approach developed in this work is new and flexible, and the explicit expression of the radial profile is no longer required as in [41] with the quadratic shape. The alternative for that is a refined study of the associated spectral problem based on Sturm-Liouville differential equation with a variable potential that changes the sign depending on the shape of the profile and the location of the time period. Deep hidden structure on positive definiteness of some intermediate integral operators are also discovered and used in a crucial way. Notice that a special study will be performed for the linear problem associated with the first mode founded on Prüfer transformation and Kneser's Theorem on the non-oscillation phenomenon.
△ Less
Submitted 2 February, 2023;
originally announced February 2023.
-
Explicit minimisers for anisotropic Coulomb energies in 3D
Authors:
Joan Mateu,
Maria Giovanna Mora,
Luca Rondi,
Lucia Scardia,
Joan Verdera
Abstract:
In this paper we consider a general class of anisotropic energies in three dimensions and give a complete characterisation of their minimisers. We show that, depending on the Fourier transform of the interaction potential, the minimiser is either the normalised characteristic function of an ellipsoid or a measure supported on a two-dimensional ellipse. In particular, it is always an ellipsoid if t…
▽ More
In this paper we consider a general class of anisotropic energies in three dimensions and give a complete characterisation of their minimisers. We show that, depending on the Fourier transform of the interaction potential, the minimiser is either the normalised characteristic function of an ellipsoid or a measure supported on a two-dimensional ellipse. In particular, it is always an ellipsoid if the transform is strictly positive, while when the Fourier transform is degenerate both cases can occur. Finally, we show an explicit example where loss of dimensionality of the minimiser does occur.
△ Less
Submitted 13 October, 2022;
originally announced October 2022.
-
Locally weighted minimum contrast estimation for spatio-temporal log-Gaussian Cox processes
Authors:
Nicoletta D'Angelo,
Giada Adelfio,
Jorge Mateu
Abstract:
We propose a local version of spatio-temporal log-Gaussian Cox processes using Local Indicators of Spatio-Temporal Association (LISTA) functions into the minimum contrast procedure to obtain space as well as time-varying parameters.
We resort to the joint minimum contrast fitting method to estimate the set of second-order parameters. This approach has the advantage of being suitable in both sepa…
▽ More
We propose a local version of spatio-temporal log-Gaussian Cox processes using Local Indicators of Spatio-Temporal Association (LISTA) functions into the minimum contrast procedure to obtain space as well as time-varying parameters.
We resort to the joint minimum contrast fitting method to estimate the set of second-order parameters. This approach has the advantage of being suitable in both separable and non-separable parametric specifications of the correlation function of the underlying Gaussian Random Field.
We present simulation studies to assess the performance of the proposed fitting procedure, and show an application to seismic spatio-temporal point pattern data.
△ Less
Submitted 15 September, 2022;
originally announced September 2022.
-
Local inhomogeneous weighted summary statistics for marked point processes
Authors:
Nicoletta D'Angelo,
Giada Adelfio,
Jorge Mateu,
Ottmar Cronie
Abstract:
We introduce a family of local inhomogeneous mark-weighted summary statistics, of order two and higher, for general marked point processes. Depending on how the involved weight function is specified, these summary statistics capture different kinds of local dependence structures. We first derive some basic properties and show how these new statistical tools can be used to construct most existing s…
▽ More
We introduce a family of local inhomogeneous mark-weighted summary statistics, of order two and higher, for general marked point processes. Depending on how the involved weight function is specified, these summary statistics capture different kinds of local dependence structures. We first derive some basic properties and show how these new statistical tools can be used to construct most existing summary statistics for (marked) point processes. We then propose a local test of random labelling. This procedure allows us to identify points, and consequently regions, where the random labelling assumption does not hold, e.g.~when the (functional) marks are spatially dependent. Through a simulation study we show that the test is able to detect local deviations from random labelling. We also provide an application to an earthquake point pattern with functional marks given by seismic waveforms.
△ Less
Submitted 19 August, 2022;
originally announced August 2022.
-
Clustering constrained on linear networks
Authors:
Asael Fabian Martínez,
Somnath Chaudhuri,
Carlos Díaz-Avalos,
Pablo Juan,
Jorge Mateu,
Ramsés H. Mena
Abstract:
An unsupervised classification method for point events occurring on a network of lines is proposed. The idea relies on the distributional flexibility and practicality of random partition models to discover the clustering structure featuring observations from a particular phenomenon taking place on a given set of edges. By incorporating the spatial effect in the random partition distribution, induc…
▽ More
An unsupervised classification method for point events occurring on a network of lines is proposed. The idea relies on the distributional flexibility and practicality of random partition models to discover the clustering structure featuring observations from a particular phenomenon taking place on a given set of edges. By incorporating the spatial effect in the random partition distribution, induced by a Dirichlet process, one is able to control the distance between edges and events, thus leading to an appealing clustering method. A Gibbs sampler algorithm is proposed and evaluated with a sensitivity analysis. The proposal is motivated and illustrated by the analysis of crime and violence patterns in Mexico City.
△ Less
Submitted 19 July, 2022;
originally announced July 2022.
-
A Spatio-Temporal Dirichlet Process Mixture Model for Coronavirus Disease-19
Authors:
Jaewoo Park,
Seorim Yi,
Won Chang,
Jorge Mateu
Abstract:
Understanding the spatio-temporal patterns of the coronavirus disease 2019 (COVID-19) is essential to construct public health interventions. Spatially referenced data can provide richer opportunities to understand the mechanism of the disease spread compared to the more often encountered aggregated count data. We propose a spatio-temporal Dirichlet process mixture model to analyze confirmed cases…
▽ More
Understanding the spatio-temporal patterns of the coronavirus disease 2019 (COVID-19) is essential to construct public health interventions. Spatially referenced data can provide richer opportunities to understand the mechanism of the disease spread compared to the more often encountered aggregated count data. We propose a spatio-temporal Dirichlet process mixture model to analyze confirmed cases of COVID-19 in an urban environment. Our method can detect unobserved cluster centers of the epidemics, and estimate the space-time range of the clusters that are useful to construct a warning system. Furthermore, our model can measure the impact of different types of landmarks in the city, which provides an intuitive explanation of disease spreading sources from different time points. To efficiently capture the temporal dynamics of the disease patterns, we employ a sequential approach that uses the posterior distribution of the parameters for the previous time step as the prior information for the current time step. This approach enables us to incorporate time dependence into our model in a computationally efficient manner without complicating the model structure. We also develop a model assessment by comparing the data with theoretical densities, and outline the goodness-of-fit of our fitted model.
△ Less
Submitted 13 July, 2022;
originally announced July 2022.
-
Time periodic doubly connected solutions for the 3D quasi-geostrophic model
Authors:
C. García,
T. Hmidi,
J. Mateu
Abstract:
In this paper, we construct time periodic doubly connected solutions for the 3D quasi-geostrophic model in the patch setting. More specifically, we prove the existence of nontrivial $m$-fold doubly connected rotating patches bifurcating from a generic doubly connected revolution shape domain with higher symmetry $m\geq m_0$ and $m_0$ is large enough. The linearized matrix operator at the equilibri…
▽ More
In this paper, we construct time periodic doubly connected solutions for the 3D quasi-geostrophic model in the patch setting. More specifically, we prove the existence of nontrivial $m$-fold doubly connected rotating patches bifurcating from a generic doubly connected revolution shape domain with higher symmetry $m\geq m_0$ and $m_0$ is large enough. The linearized matrix operator at the equilibrium state is with variable and singular coefficients and its spectral analysis is performed via the approach devised in [27] where a suitable symmetrization has been introduced. New difficulties emerge due to the interaction between the surfaces making the spectral problem richer and involved.
△ Less
Submitted 21 June, 2022;
originally announced June 2022.
-
Removable singularities for solutions of the fractional Heat Equation in time varying domains
Authors:
Joan Mateu,
Laura Prat
Abstract:
In this paper we study removable singularities for solutions of the fractional heat equation in time varying domains. We introduce associated capacities and we study some of its metric and geometric properties.
In this paper we study removable singularities for solutions of the fractional heat equation in time varying domains. We introduce associated capacities and we study some of its metric and geometric properties.
△ Less
Submitted 5 May, 2022;
originally announced May 2022.
-
ANOVA for Data in Metric Spaces, with Applications to Spatial Point Patterns
Authors:
Raoul Müller,
Dominic Schuhmacher,
Jorge Mateu
Abstract:
We give a review of recent ANOVA-like procedures for testing group differences based on data in a metric space and present a new such procedure. Our statistic is based on the classic Levene's test for detecting differences in dispersion. It uses only pairwise distances of data points and and can be computed quickly and precisely in situations where the computation of barycenters ("generalized mean…
▽ More
We give a review of recent ANOVA-like procedures for testing group differences based on data in a metric space and present a new such procedure. Our statistic is based on the classic Levene's test for detecting differences in dispersion. It uses only pairwise distances of data points and and can be computed quickly and precisely in situations where the computation of barycenters ("generalized means") in the data space is slow, only by approximation or even infeasible. We show the asymptotic normality of our test statistic and present simulation studies for spatial point pattern data, in which we compare the various procedures in a 1-way ANOVA setting. As an application, we perform a 2-way ANOVA on a data set of bubbles in a mineral flotation process.
△ Less
Submitted 18 February, 2022; v1 submitted 21 January, 2022;
originally announced January 2022.
-
Generalised functional additive mixed models with compositional covariates for areal Covid-19 incidence curves
Authors:
Matthias Eckardt,
Jorge Mateu,
Sonja Greven
Abstract:
We extend the generalised functional additive mixed model to include (functional) compositional covariates carrying relative information of a whole. Relying on the isometric isomorphism of the Bayes Hilbert space of probability densities with a subspace of the $L^2$, we include functional compositions as transformed functional covariates with constrained effect function. The extended model allows…
▽ More
We extend the generalised functional additive mixed model to include (functional) compositional covariates carrying relative information of a whole. Relying on the isometric isomorphism of the Bayes Hilbert space of probability densities with a subspace of the $L^2$, we include functional compositions as transformed functional covariates with constrained effect function. The extended model allows for the estimation of linear, nonlinear and time-varying effects of scalar and functional covariates, as well as (correlated) functional random effects, in addition to the compositional effects. We use the model to estimate the effect of the age, sex and smoking (functional) composition of the population on regional Covid-19 incidence data for Spain, while accounting for climatological and socio-demographic covariate effects and spatial correlation.
△ Less
Submitted 20 January, 2022;
originally announced January 2022.
-
Measurement Error Models for Spatial Network Lattice Data: Analysis of Car Crashes in Leeds
Authors:
Andrea Gilardi,
Riccardo Borgoni,
Luca Presicce,
Jorge Mateu
Abstract:
Road casualties represent an alarming concern for modern societies. During the last years, several authors proposed sophisticated approaches to help authorities implement new policies. These models were usually developed considering a set of socioeconomic variables and ignoring the measurement error, which can bias the statistical inference. This paper presents a Bayesian model to analyse car cras…
▽ More
Road casualties represent an alarming concern for modern societies. During the last years, several authors proposed sophisticated approaches to help authorities implement new policies. These models were usually developed considering a set of socioeconomic variables and ignoring the measurement error, which can bias the statistical inference. This paper presents a Bayesian model to analyse car crashes occurrences at the network-lattice level, taking into account measurement error in the spatial covariate. The suggested methodology is exemplified by considering the collisions in the road network of Leeds (UK) during 2011-2019. Traffic volumes are approximated using an extensive set of counts obtained from mobile devices and the estimates are adjusted using a spatial measurement error correction.
△ Less
Submitted 10 March, 2023; v1 submitted 7 January, 2022;
originally announced January 2022.
-
Stability of ellipsoids as the energy minimisers of perturbed Coulomb energies
Authors:
Joan Mateu,
Maria Giovanna Mora,
Luca Rondi,
Lucia Scardia,
Joan Verdera
Abstract:
In this paper we characterise the minimiser for a class of nonlocal perturbations of the Coulomb energy. We show that the minimiser is the normalised characteristic function of an ellipsoid, under the assumption that the perturbation kernel has the same homogeneity as the Coulomb potential, is even, smooth off the origin and sufficiently small. This result can be seen as the stability of ellipsoid…
▽ More
In this paper we characterise the minimiser for a class of nonlocal perturbations of the Coulomb energy. We show that the minimiser is the normalised characteristic function of an ellipsoid, under the assumption that the perturbation kernel has the same homogeneity as the Coulomb potential, is even, smooth off the origin and sufficiently small. This result can be seen as the stability of ellipsoids as energy minimisers, since the minimiser of the Coulomb energy is the normalised characteristic function of a ball.
△ Less
Submitted 28 December, 2021;
originally announced December 2021.
-
Mapping the intensity function of a non-stationary point process in unobserved areas
Authors:
Edith Gabriel,
Francisco J. Rodríguez-Cortes,
Jérôme Coville,
Jorge Mateu,
Joël Chadoeuf
Abstract:
Seismic networks provide data that are used as basis both for public safety decisions and for scientific research. Their configuration affects the data completeness, which in turn, critically affects several seismological scientific targets (e.g., earthquake prediction, seismic hazard...). In this context, a key aspect is how to map earthquakes density in seismogenic areas from censored data or ev…
▽ More
Seismic networks provide data that are used as basis both for public safety decisions and for scientific research. Their configuration affects the data completeness, which in turn, critically affects several seismological scientific targets (e.g., earthquake prediction, seismic hazard...). In this context, a key aspect is how to map earthquakes density in seismogenic areas from censored data or even in areas that are not covered by the network. We propose to predict the spatial distribution of earthquakes from the knowledge of presence locations and geological relationships, taking into account any interaction between records. Namely, in a more general setting, we aim to estimate the intensity function of a point process, conditional to its censored realization, as in geostatistics for continuous processes. We define a predictor as the best linear unbiased combination of the observed point pattern. We show that the weight function associated to the predictor is the solution of a Fredholm equation of second kind. Both the kernel and the source term of the Fredholm equation are related to the first-and second-order characteristics of the point process through the intensity and the pair correlation function. Results are presented and illustrated on simulated non-stationary point processes and real data for mapping Greek Hellenic seismicity in a region with unreliable and incomplete records.
△ Less
Submitted 3 June, 2022; v1 submitted 29 November, 2021;
originally announced November 2021.
-
A new class of spatial covariance functions generated by higher-order kernels
Authors:
Mohammad Ghorbani,
Jorge Mateu
Abstract:
Covariance functions and variograms play a fundamental role in exploratory analysis and statistical modelling of spatial and spatio-temporal datasets. In this paper, we construct a new class of spatial covariance functions using the Fourier transform of some higher-order kernels. Further, we extend this class of the spatial covariance functions to the spatio-temporal setting by using the idea used…
▽ More
Covariance functions and variograms play a fundamental role in exploratory analysis and statistical modelling of spatial and spatio-temporal datasets. In this paper, we construct a new class of spatial covariance functions using the Fourier transform of some higher-order kernels. Further, we extend this class of the spatial covariance functions to the spatio-temporal setting by using the idea used in Ma (2003).
△ Less
Submitted 12 November, 2021;
originally announced November 2021.
-
Non-stationary spatio-temporal point process modeling for high-resolution COVID-19 data
Authors:
Zheng Dong,
Shixiang Zhu,
Yao Xie,
Jorge Mateu,
Francisco J. Rodríguez-Cortés
Abstract:
Most COVID-19 studies commonly report figures of the overall infection at a state- or county-level. This aggregation tends to miss out on fine details of virus propagation. In this paper, we analyze a high-resolution COVID-19 dataset in Cali, Colombia, that records the precise time and location of every confirmed case. We develop a non-stationary spatio-temporal point process equipped with a neura…
▽ More
Most COVID-19 studies commonly report figures of the overall infection at a state- or county-level. This aggregation tends to miss out on fine details of virus propagation. In this paper, we analyze a high-resolution COVID-19 dataset in Cali, Colombia, that records the precise time and location of every confirmed case. We develop a non-stationary spatio-temporal point process equipped with a neural network-based kernel to capture the heterogeneous correlations among COVID-19 cases. The kernel is carefully crafted to enhance expressiveness while maintaining model interpretability. We also incorporate some exogenous influences imposed by city landmarks. Our approach outperforms the state-of-the-art in forecasting new COVID-19 cases with the capability to offer vital insights into the spatio-temporal interaction between individuals concerning the disease spread in a metropolis.
△ Less
Submitted 9 March, 2023; v1 submitted 18 September, 2021;
originally announced September 2021.
-
Phenotypic-dependent variability and the emergence of tolerance in bacterial populations
Authors:
José Camacho Mateu,
Matteo Sireci,
Miguel A. Muñoz
Abstract:
Ecological and evolutionary dynamics have been historically regarded as unfolding at broadly separated timescales. However, these two types of processes are nowadays well documented to much more tightly than traditionally assumed, especially in communities of microorganisms. With this motivation in mind, here we scrutinize recent experimental results showing evidence of rapid evolution of toleranc…
▽ More
Ecological and evolutionary dynamics have been historically regarded as unfolding at broadly separated timescales. However, these two types of processes are nowadays well documented to much more tightly than traditionally assumed, especially in communities of microorganisms. With this motivation in mind, here we scrutinize recent experimental results showing evidence of rapid evolution of tolerance by lag in bacterial populations that are periodically exposed to antibiotic stress in laboratory conditions. In particular, the distribution of single-cell lag times evolves its average value to approximately fit the antibiotic-exposure time. Moreover, the distribution develops right skewed heavy tails, revealing the presence of individuals with anomalously large lag times. Here, we develop a parsimonious individual-based model mimicking the actual demographic processes of the experimental setup. Individuals are characterized by a single phenotypic trait: their intrinsic lag time, which is transmitted with variation to the progeny. The model (in a version in which the amplitude of phenotypic variations grows with the parentś lag time) is able to reproduce quite well the key empirical observations. Furthermore, we develop a general mathematical framework allowing us to describe with good accuracy the properties of the stochastic model by means of a macroscopic equation, which generalizes the Crow Kimura equation in population genetics. From a broader perspective, this work represents a benchmark for the mathematical framework designed to tackle much more general eco-evolutionary problems, thus paving the road to further research avenues.
△ Less
Submitted 6 October, 2021; v1 submitted 17 September, 2021;
originally announced September 2021.
-
A non-separable first-order spatio-temporal intensity for events on linear networks: an application to ambulance interventions
Authors:
Andrea Gilardi,
Riccardo Borgoni,
Jorge Mateu
Abstract:
The algorithms used for the optimal management of an ambulance fleet require an accurate description of the spatio-temporal evolution of the emergency events. In the last years, several authors have proposed sophisticated statistical approaches to forecast ambulance dispatches, typically modelling the data as a point pattern occurring on a planar region. Nevertheless, ambulance interventions can b…
▽ More
The algorithms used for the optimal management of an ambulance fleet require an accurate description of the spatio-temporal evolution of the emergency events. In the last years, several authors have proposed sophisticated statistical approaches to forecast ambulance dispatches, typically modelling the data as a point pattern occurring on a planar region. Nevertheless, ambulance interventions can be more appropriately modelled as a realisation of a point process occurring on a linear network. The constrained spatial domain raises specific challenges and unique methodological problems that cannot be ignored when developing a proper statistical approach. Hence, this paper proposes a spatio-temporal model to analyse ambulance dispatches focusing on the interventions that occurred in the road network of Milan (Italy) from 2015 to 2017. We adopt a non-separable first-order intensity function with spatial and temporal terms. The temporal dimension is estimated semi-parametrically using a Poisson regression model, while the spatial dimension is estimated non-parametrically using a network kernel function. A set of weights is included in the spatial term to capture space-time interactions, inducing non-separability in the intensity function. A series of tests show that our approach successfully models the ambulance interventions and captures the space-time patterns more accurately than planar or separable point process models.
△ Less
Submitted 28 April, 2023; v1 submitted 1 June, 2021;
originally announced June 2021.
-
Bayesian Effect Selection for Additive Quantile Regression with an Analysis to Air Pollution Thresholds
Authors:
Nadja Klein,
Jorge Mateu
Abstract:
Statistical techniques used in air pollution modelling usually lack the possibility to understand which predictors affect air pollution in which functional form; and are not able to regress on exceedances over certain thresholds imposed by authorities directly. The latter naturally induce conditional quantiles and reflect the seriousness of particular events. In the present paper we focus on this…
▽ More
Statistical techniques used in air pollution modelling usually lack the possibility to understand which predictors affect air pollution in which functional form; and are not able to regress on exceedances over certain thresholds imposed by authorities directly. The latter naturally induce conditional quantiles and reflect the seriousness of particular events. In the present paper we focus on this important aspect by developing quantile regression models further. We propose a general Bayesian effect selection approach for additive quantile regression within a highly interpretable framework. We place separate normal beta prime spike and slab priors on the scalar importance parameters of effect parts and implement a fast Gibbs sampling scheme. Specifically, it enables to study quantile-specific covariate effects, allows these covariates to be of general functional form using additive predictors, and facilitates the analysts' decision whether an effect should be included linearly, non-linearly or not at all in the quantiles of interest. In a detailed analysis on air pollution data in Madrid (Spain) we find the added value of modelling extreme nitrogen dioxide (NO2) concentrations and how thresholds are driven differently by several climatological variables and traffic as a spatial proxy. Our results underpin the need of enhanced statistical models to support short-term decisions and enable local authorities to mitigate or even prevent exceedances of NO2 concentration limits.
△ Less
Submitted 23 May, 2021;
originally announced May 2021.
-
The regularity of the boundary of vortex patches for some non-linear transport equations
Authors:
J. C. Cantero,
J. Mateu,
J. Orobitg,
J. Verdera
Abstract:
We prove the persistence of boundary smoothness of vortex patches for a non-linear transport equation in $\mathbb{R}^n$ with velocity field given by convolution of the density with an odd kernel, homogeneous of degree $-(n-1)$ and of class $C^2(\mathbb{R}^n\setminus\{0\}, \mathbb{R}^n).$ This allows the velocity field to have non-trivial divergence. The quasi-geostrophic equation in…
▽ More
We prove the persistence of boundary smoothness of vortex patches for a non-linear transport equation in $\mathbb{R}^n$ with velocity field given by convolution of the density with an odd kernel, homogeneous of degree $-(n-1)$ and of class $C^2(\mathbb{R}^n\setminus\{0\}, \mathbb{R}^n).$ This allows the velocity field to have non-trivial divergence. The quasi-geostrophic equation in $\mathbb{R}^3$ and the Cauchy transport equation in the plane are examples.
△ Less
Submitted 26 September, 2023; v1 submitted 9 March, 2021;
originally announced March 2021.
-
Multivariate hierarchical analysis of car crashes data considering a spatial network lattice
Authors:
Andrea Gilardi,
Jorge Mateu,
Riccardo Borgoni,
Robin Lovelace
Abstract:
Road traffic casualties represent a hidden global epidemic, demanding evidence-based interventions. This paper demonstrates a network lattice approach for identifying road segments of particular concern, based on a case study of a major city (Leeds, UK), in which 5,862 crashes of different severities were recorded over an eight-year period (2011-2018). We consider a family of Bayesian hierarchical…
▽ More
Road traffic casualties represent a hidden global epidemic, demanding evidence-based interventions. This paper demonstrates a network lattice approach for identifying road segments of particular concern, based on a case study of a major city (Leeds, UK), in which 5,862 crashes of different severities were recorded over an eight-year period (2011-2018). We consider a family of Bayesian hierarchical models that include spatially structured and unstructured random effects, to capture the dependencies between the severity levels. Results highlight roads that are more prone to collisions, relative to estimated traffic volumes, in the northwest and south of city-centre. We analyse the Modifiable Areal Unit Problem (MAUP), proposing a novel procedure to investigate the presence of MAUP on a network lattice. We conclude that our methods enable a reliable estimation of road safety levels to help identify "hotspots" on the road network and to inform effective local interventions.
△ Less
Submitted 25 June, 2021; v1 submitted 25 November, 2020;
originally announced November 2020.
-
Spatio-temporal small area surveillance of the Covid-19 pandemics
Authors:
Miguel A. Martinez-Beneito,
Jorge Mateu,
Paloma Botella-Rocamora
Abstract:
The emergence of Covid-19 requires new effective tools for epidemiological surveillance. Spatio-temporal disease mapping models, which allow dealing with highly disaggregated spatial and temporal units of analysis, are a priority in this sense. Spatio-temporal models provide a geographically detailed and temporally updated overview of the current state of the pandemics, making public health interv…
▽ More
The emergence of Covid-19 requires new effective tools for epidemiological surveillance. Spatio-temporal disease mapping models, which allow dealing with highly disaggregated spatial and temporal units of analysis, are a priority in this sense. Spatio-temporal models provide a geographically detailed and temporally updated overview of the current state of the pandemics, making public health interventions to be more effective. Moreover, the use of spatio-temporal disease mapping models in the new Covid-19 epidemic context, facilitates estimating newly demanded epidemiological indicators, such as the instantaneous reproduction number (R_t), even for small areas. This, in turn, allows to adapt traditional disease mapping models to these new circumstancies and make their results more useful in this particular context.
In this paper we propose a new spatio-temporal disease mapping model, particularly suited to Covid-19 surveillance. As an additional result, we derive instantaneous reproduction number estimates for small areas, enabling monitoring this parameter with a high spatial disaggregation. We illustrate the use of our proposal with the separate study of the disease pandemics in two Spanish regions. As a result, we illustrate how touristic flows could haved shaped the spatial distribution of the disease. In these real studies, we also propose new surveillance tools that can be used by regional public health services to make a more efficient use of their resources.
△ Less
Submitted 8 November, 2020;
originally announced November 2020.
-
A spatial functional count model for heterogeneity analysis in time
Authors:
Torres-Signes,
M. P. Frías,
J. Mateu,
M. D. Ruiz-Medina
Abstract:
A spatial curve dynamical model framework is adopted for functional prediction of counts in a spatiotemporal log-Gaussian Cox process model. Our spatial functional estimation approach handles both wavelet-based heterogeneity analysis in time, and spectral analysis in space. Specifically, model fitting is achieved by minimising the information divergence or relative entropy between the multiscale m…
▽ More
A spatial curve dynamical model framework is adopted for functional prediction of counts in a spatiotemporal log-Gaussian Cox process model. Our spatial functional estimation approach handles both wavelet-based heterogeneity analysis in time, and spectral analysis in space. Specifically, model fitting is achieved by minimising the information divergence or relative entropy between the multiscale model underlying the data and the corresponding candidates in the spatial spectral domain. A simulation study is carried out within the family of log-Gaussian Spatial Autoregressive l2-valued processes (SARl2 processes) to illustrate the asymptotic properties of the proposed spatial functional estimators. We apply our modelling strategy to spatiotemporal prediction of respiratory disease mortality.
△ Less
Submitted 8 October, 2020;
originally announced October 2020.
-
A hierarchical spatio-temporal model to analyze relative risk variations of COVID-19: a focus on Spain, Italy and Germany
Authors:
Abdollah Jalilian,
Jorge Mateu
Abstract:
The novel coronavirus disease (COVID-19) has spread rapidly across the world in a short period of time and with a heterogeneous pattern. Understanding the underlying temporal and spatial dynamics in the spread of COVID-19 can result in informed and timely public health policies. In this paper, we use a spatio-temporal stochastic model to explain the temporal and spatial variations in the daily num…
▽ More
The novel coronavirus disease (COVID-19) has spread rapidly across the world in a short period of time and with a heterogeneous pattern. Understanding the underlying temporal and spatial dynamics in the spread of COVID-19 can result in informed and timely public health policies. In this paper, we use a spatio-temporal stochastic model to explain the temporal and spatial variations in the daily number of new confirmed cases in Spain, Italy and Germany from late February to mid September 2020. Using a hierarchical Bayesian framework, we found that the temporal trend of the epidemic in the three countries rapidly reached their peaks and slowly started to decline at the beginning of April and then increased and reached their second maximum in August. However decline and increase of the temporal trend seems to be sharper in Spain and smoother in Germany. The spatial heterogeneity of the relative risk of COVID-19 in Spain is also more pronounced than Italy and Germany.
△ Less
Submitted 28 September, 2020;
originally announced September 2020.
-
Removable singularities for Lipschitz caloric functions in time varying domains
Authors:
Joan Mateu,
Laura Prat,
Xavier Tolsa
Abstract:
In this paper we study removable singularities for regular $(1,1/2)$-Lipschitz solutions of the heat equation in time varying domains. We introduce an associated Lipschitz caloric capacity and we study its metric and geometric properties and the connection with the $L^2$ boundedness of the singular integral whose kernel is given by the gradient of the fundamental solution of the heat equation.
In this paper we study removable singularities for regular $(1,1/2)$-Lipschitz solutions of the heat equation in time varying domains. We introduce an associated Lipschitz caloric capacity and we study its metric and geometric properties and the connection with the $L^2$ boundedness of the singular integral whose kernel is given by the gradient of the fundamental solution of the heat equation.
△ Less
Submitted 24 December, 2020; v1 submitted 7 May, 2020;
originally announced May 2020.
-
Nonparametric testing of the dependence structure among points-marks-covariates in spatial point patterns
Authors:
Jiří Dvořák,
Tomáš Mrkvička,
Jorge Mateu,
Jonatan González
Abstract:
We investigate testing of the hypothesis of independence between a covariate and the marks in a marked point process. It would be rather straightforward if the (unmarked) point process were independent of the covariate and the marks. In practice, however, such an assumption is questionable and possible dependence between the point process and the covariate or the marks may lead to incorrect conclu…
▽ More
We investigate testing of the hypothesis of independence between a covariate and the marks in a marked point process. It would be rather straightforward if the (unmarked) point process were independent of the covariate and the marks. In practice, however, such an assumption is questionable and possible dependence between the point process and the covariate or the marks may lead to incorrect conclusions. Therefore, we propose to investigate the complete dependence structure in the triangle points--marks--covariates together. We take advantage of the recent development of the nonparametric random shift methods, namely the new variance correction approach, and propose tests of the null hypothesis of independence between the marks and the covariate and between the points and the covariate. We present a detailed simulation study showing the performance of the methods and provide two theorems establishing the appropriate form of the correction factors for the variance correction. Finally, we illustrate the use of the proposed methods in two real applications.
△ Less
Submitted 13 May, 2022; v1 submitted 3 May, 2020;
originally announced May 2020.
-
Time periodic solutions for 3D quasi-geostrophic model
Authors:
Claudia García,
Taoufik Hmidi,
Joan Mateu
Abstract:
This paper aims to study time periodic solutions for 3D inviscid quasi-geostrophic model. We show the existence of non trivial rotating patches by suitable perturbation of stationary solutions given by generic revolution shapes around the vertical axis. The construction of those special solutions are done through bifurcation theory. In general, the spectral problem is very delicate and strongly de…
▽ More
This paper aims to study time periodic solutions for 3D inviscid quasi-geostrophic model. We show the existence of non trivial rotating patches by suitable perturbation of stationary solutions given by generic revolution shapes around the vertical axis. The construction of those special solutions are done through bifurcation theory. In general, the spectral problem is very delicate and strongly depends on the shape of the initial stationary solutions. More specifically, the spectral study can be related to an eigenvalue problem of a self-adjoint compact operator. We are able to implement the bifurcation only from the largest eigenvalues of the operator, which are simple. Additional difficulties generated by the singularities of the poles are solved through the use of suitable function spaces with Dirichlet boundary condition type and refined potential theory with anisotropic kernels.
△ Less
Submitted 13 December, 2021; v1 submitted 3 April, 2020;
originally announced April 2020.
-
Explicit minimisers of some nonlocal anisotropic energies: a short proof
Authors:
J. Mateu,
M. G. Mora,
l. Rondi,
L. Scardia,
J. Verdera
Abstract:
In this paper we consider nonlocal energies defined on probability measures in the plane, given by a convolution interaction term plus a quadratic confinement. The interaction kernel is $-\log|z|+α\, x^2/|z|^2, \; z=x+iy,$ with $-1 < α< 1.$ This kernel is anisotropic except for the Coulombic case $α=0.$ We present a short compact proof of the known surprising fact that the unique minimiser of the…
▽ More
In this paper we consider nonlocal energies defined on probability measures in the plane, given by a convolution interaction term plus a quadratic confinement. The interaction kernel is $-\log|z|+α\, x^2/|z|^2, \; z=x+iy,$ with $-1 < α< 1.$ This kernel is anisotropic except for the Coulombic case $α=0.$ We present a short compact proof of the known surprising fact that the unique minimiser of the energy is the normalised characteristic function of the domain enclosed by an ellipse with horizontal semi-axis $\sqrt{1-α}$ and vertical semi-axis $\sqrt{1+α}.$ Letting $α\to 1^-$ we find that the semicircle law on the vertical axis is the unique minimiser of the corresponding energy, a result related to interacting dislocations, and previously obtained by some of the authors. We devote the first sections of this paper to presenting some well-known background material in the simplest way possible, so that readers unfamiliar with the subject find the proofs accessible
△ Less
Submitted 9 October, 2020; v1 submitted 30 March, 2020;
originally announced March 2020.
-
Graphical modelling and partial characteristics for multitype and multivariate-marked spatio-temporal point processes
Authors:
Matthias Eckardt,
Jonatan A. González,
Jorge Mateu
Abstract:
This paper contributes to the multivariate analysis of marked spatio-temporal point process data by introducing different partial point characteristics and extending the spatial dependence graph model formalism. Our approach yields a unified framework for different types of spatio-temporal data including both, purely qualitatively (multivariate) cases and multivariate cases with additional quantit…
▽ More
This paper contributes to the multivariate analysis of marked spatio-temporal point process data by introducing different partial point characteristics and extending the spatial dependence graph model formalism. Our approach yields a unified framework for different types of spatio-temporal data including both, purely qualitatively (multivariate) cases and multivariate cases with additional quantitative marks. The proposed graphical model is defined through partial spectral density characteristics, it is highly computationally efficient and reflects the conditional similarity among sets of spatio-temporal sub-processes of either points or marked points with identical discrete marks. The paper considers three applications, two on crime data and a third one on forestry.
△ Less
Submitted 5 March, 2020;
originally announced March 2020.
-
$L^2$-bounded singular integrals on a purely unrectifiable set in $\mathbb{R}^d$
Authors:
Joan Mateu,
Laura Prat
Abstract:
We construct an example of a purely unrectifiable measure $μ$ in $\mathbb{R}^d$ for which the singular integrals associated to the kernels $\displaystyle{K(x)=\frac{P_{2k+1}(x)}{|x|^{2k+d}}}$, with $k\geq 1$ and $P_{2k+1}$ a homogeneous harmonic polynomial of degree $2k+1$, are bounded in $L^2(μ)$. This contrasts starkly with the results concerning the Riesz kernel…
▽ More
We construct an example of a purely unrectifiable measure $μ$ in $\mathbb{R}^d$ for which the singular integrals associated to the kernels $\displaystyle{K(x)=\frac{P_{2k+1}(x)}{|x|^{2k+d}}}$, with $k\geq 1$ and $P_{2k+1}$ a homogeneous harmonic polynomial of degree $2k+1$, are bounded in $L^2(μ)$. This contrasts starkly with the results concerning the Riesz kernel $\displaystyle{\frac{x}{|x|^{d}}}$ in $\mathbb{R}^d$.
△ Less
Submitted 24 December, 2019;
originally announced December 2019.
-
Functional marked point processes -- A natural structure to unify spatio-temporal frameworks and to analyse dependent functional data
Authors:
Ottmar Cronie,
Mohammad Ghorbani,
Jorge Mateu,
Jun Yu
Abstract:
This paper treats functional marked point processes (FMPPs), which are defined as marked point processes where the marks are random elements in some (Polish) function space. Such marks may represent e.g. spatial paths or functions of time. To be able to consider e.g. multivariate FMPPs, we also attach an additional, Euclidean, mark to each point. We indicate how FMPPs quite naturally connect the p…
▽ More
This paper treats functional marked point processes (FMPPs), which are defined as marked point processes where the marks are random elements in some (Polish) function space. Such marks may represent e.g. spatial paths or functions of time. To be able to consider e.g. multivariate FMPPs, we also attach an additional, Euclidean, mark to each point. We indicate how FMPPs quite naturally connect the point process framework with both the functional data analysis framework and the geostatistical framework. We further show that various existing models fit well into the FMPP framework. In addition, we introduce a new family of summary statistics, weighted marked reduced moment measures, together with their non-parametric estimators, in order to study features of the functional marks. We further show how they generalise other summary statistics and we finally apply these tools to analyse population structures, such as demographic evolution and sex ratio over time, in Spanish provinces.
△ Less
Submitted 29 November, 2019;
originally announced November 2019.
-
Revisiting the random shift approach for testing in spatial statistics
Authors:
Tomas Mrkvicka,
Jiri Dvorak,
Jonatan A. Gonzalez,
Jorge Mateu
Abstract:
We consider the problem of non-parametric testing of independence of two components of a stationary bivariate spatial process. In particular, we revisit the random shift approach that has become a standard method for testing the independent superposition hypothesis in spatial statistics, and it is widely used in a plethora of practical applications. However, this method has a problem of liberality…
▽ More
We consider the problem of non-parametric testing of independence of two components of a stationary bivariate spatial process. In particular, we revisit the random shift approach that has become a standard method for testing the independent superposition hypothesis in spatial statistics, and it is widely used in a plethora of practical applications. However, this method has a problem of liberality caused by breaking the marginal spatial correlation structure due to the toroidal correction. This indeed causes that the assumption of exchangability, which is essential for the Monte Carlo test to be exact, is not fulfilled.
We present a number of permutation strategies and show that the random shift with the variance correction brings a suitable improvement compared to the torus correction in the random field case. It reduces the liberality and achieves the largest power from all investigated variants. To obtain the variance for the variance correction method, several approaches were studied. The best results were achieved, for the sample covariance as the test statistics, with the correction factor $1/n$. This corresponds to the asymptotic order of the variance of the test statistics.
In the point process case, the problem of deviations from exchangeability is far more complex and we propose an alternative strategy based on the mean cross nearest-neighbor distance and torus correction. It reduces the liberality but achieves slightly lower power than the usual cross $K$-function. Therefore we recommend it, when the point patterns are clustered, where the cross $K$-function achieves liberality.
△ Less
Submitted 1 November, 2019;
originally announced November 2019.
-
Inhomogeneous higher-order summary statistics for linear network point processes
Authors:
Ottmar Cronie,
Mehdi Moradi,
Jorge Mateu
Abstract:
We introduce the notion of intensity reweighted moment pseudostationary point processes on linear networks. Based on arbitrary general regular linear network distances, we propose geometrically corrected versions of different higher-order summary statistics, including the inhomogeneous empty space function, the inhomogeneous nearest neighbour distance distribution function and the inhomogeneous…
▽ More
We introduce the notion of intensity reweighted moment pseudostationary point processes on linear networks. Based on arbitrary general regular linear network distances, we propose geometrically corrected versions of different higher-order summary statistics, including the inhomogeneous empty space function, the inhomogeneous nearest neighbour distance distribution function and the inhomogeneous $J$-function. We also discuss their non-parametric estimators. Through a simulation study, considering models with different types of spatial interaction, we study the performance of our proposed summary statistics. Finally, we make use of our methodology to analyse two datasets: motor vehicle traffic accidents and spider data.
△ Less
Submitted 8 October, 2019;
originally announced October 2019.
-
Metrics and barycenters for point pattern data
Authors:
Raoul Müller,
Dominic Schuhmacher,
Jorge Mateu
Abstract:
We introduce the transport-transform (TT) and the relative transport-transform (RTT) metrics between finite point patterns on a general space, which provide a unified framework for earlier point pattern metrics, in particular the generalized spike time and the normalized and unnormalized OSPA metrics. Our main focus is on barycenters, i.e. minimizers of a $q$-th order Fréchet functional with respe…
▽ More
We introduce the transport-transform (TT) and the relative transport-transform (RTT) metrics between finite point patterns on a general space, which provide a unified framework for earlier point pattern metrics, in particular the generalized spike time and the normalized and unnormalized OSPA metrics. Our main focus is on barycenters, i.e. minimizers of a $q$-th order Fréchet functional with respect to these metrics.
We present a heuristic algorithm that terminates in a local minimum and is shown to be fast and reliable in a simulation study. The algorithm serves as an umbrella method that can be applied on any state space where an appropriate algorithm for solving the location problem for individual points is available. We present applications to geocoded data of crimes in Euclidean space and on a street network, illustrating that barycenters serve as informative summary statistics. Our work is a first step towards statistical inference in covariate-based models of repeated point pattern observations.
△ Less
Submitted 16 September, 2019;
originally announced September 2019.
-
On the analyticity of the trajectories of the particles in the patch problem for $2$D Euler and aggregation equations
Authors:
J. M. Burgués,
J. Mateu
Abstract:
We give a proof of the analiticity in time for the particle trajectories associated with the solutions of some transport equations when the initial datum is a patch. These results are obtained from a precise study of the Beurling transform, which provides estimates for the solutions of some new equations satisfied by the lagrangian flow.
We give a proof of the analiticity in time for the particle trajectories associated with the solutions of some transport equations when the initial datum is a patch. These results are obtained from a precise study of the Beurling transform, which provides estimates for the solutions of some new equations satisfied by the lagrangian flow.
△ Less
Submitted 14 September, 2020; v1 submitted 31 July, 2019;
originally announced July 2019.
-
The equilibrium measure for an anisotropic nonlocal energy
Authors:
J. A. Carrillo,
J. Mateu,
M. G. Mora,
L. Rondi,
L. Scardia,
J. Verdera
Abstract:
In this paper we characterise the minimisers of a one-parameter family of nonlocal and anisotropic energies $I_α$ defined on probability measures in $\R^n$, with $n\geq 3$. The energy $I_α$ consists of a purely nonlocal term of convolution type, whose interaction kernel reduces to the Coulomb potential for $α=0$ and is anisotropic otherwise, and a quadratic confinement. The two-dimensional case ar…
▽ More
In this paper we characterise the minimisers of a one-parameter family of nonlocal and anisotropic energies $I_α$ defined on probability measures in $\R^n$, with $n\geq 3$. The energy $I_α$ consists of a purely nonlocal term of convolution type, whose interaction kernel reduces to the Coulomb potential for $α=0$ and is anisotropic otherwise, and a quadratic confinement. The two-dimensional case arises in the study of defects in metals and has been solved by the authors by means of complex-analysis techniques. We prove that for $α\in (-1, n-2]$, the minimiser of $I_α$ is unique and is the (normalised) characteristic function of a spheroid. This result is a paradigmatic example of the role of the anisotropy of the kernel on the shape of minimisers. In particular, the phenomenon of loss of dimensionality, observed in dimension $n=2$, does not occur in higher dimension at the value $α=n-2$ corresponding to the sign change of the Fourier transform of the interaction potential.
△ Less
Submitted 30 June, 2019;
originally announced July 2019.
-
A spatial dependence graph model for multivariate spatial hybrid processes
Authors:
Matthias Eckardt,
Jorge Mateu
Abstract:
This paper is concerned with the joint analysis of multivariate mixed-type spatial data, where some components are point processes and some are of lattice-type by nature. After a survey of statistical methods for marked spatial point and lattice processes, the class of multivariate spatial hybrid processes is defined and embedded within the framework of spatial dependence graph models. In this mod…
▽ More
This paper is concerned with the joint analysis of multivariate mixed-type spatial data, where some components are point processes and some are of lattice-type by nature. After a survey of statistical methods for marked spatial point and lattice processes, the class of multivariate spatial hybrid processes is defined and embedded within the framework of spatial dependence graph models. In this model, the point and lattice sub-processes are identified with nodes of a graph whereas missing edges represent conditional independence among the components. This finally leads to a general framework for any type of spatial data in a multivariate setting. We demonstrate the application of our method in the analysis of a multivariate point-lattice pattern on crime and ambulance service call-out incidents recorded in London, where the points are the locations of different pre-classified crime events and the lattice components report different aggregated incident rates at ward level.
△ Less
Submitted 18 June, 2019;
originally announced June 2019.
-
Partial and semi-partial measures of spatial associations for multivariate lattice data
Authors:
Matthias Eckardt,
Jorge Mateu
Abstract:
This paper concerns the development of partial and semi-partial measures of spatial associations in the context of multivariate spatial lattice data which describe global or local associations among spatially aggregated measurements for pairs of different components conditional on all remaining components. The new measures are illustrated using aggregated data on crime counts at ward level.
This paper concerns the development of partial and semi-partial measures of spatial associations in the context of multivariate spatial lattice data which describe global or local associations among spatially aggregated measurements for pairs of different components conditional on all remaining components. The new measures are illustrated using aggregated data on crime counts at ward level.
△ Less
Submitted 4 June, 2019;
originally announced June 2019.
-
Spatial Cox processes in an infinite-dimensional framework
Authors:
M. P. Frías,
A. Torres-Signes,
M. D. Ruiz-Medina,
J. Mateu
Abstract:
We introduce a new class of spatial Cox processes driven by a Hilbert--valued random log--intensity. We adopt a parametric framework in the spectral domain, to estimate its spatial functional correlation structure. Specifically, we consider a spectral functional, based on the periodogram operator, inspired on Whittle estimation methodology. Strong-consistency of the parametric estimator is proved…
▽ More
We introduce a new class of spatial Cox processes driven by a Hilbert--valued random log--intensity. We adopt a parametric framework in the spectral domain, to estimate its spatial functional correlation structure. Specifically, we consider a spectral functional, based on the periodogram operator, inspired on Whittle estimation methodology. Strong-consistency of the parametric estimator is proved in the linear case. We illustrate this property in a simulation study under a Gaussian first order Spatial Autoregressive Hilbertian scenario for the log--intensity model. Our method is applied to the spatial functional prediction of respiratory disease mortality in the Spanish Iberian Peninsula, in the period 1980--2015.
△ Less
Submitted 10 March, 2021; v1 submitted 27 November, 2018;
originally announced November 2018.
-
Resample-smoothing of Voronoi intensity estimators
Authors:
M. Mehdi Moradi,
Ottmar Cronie,
Ege Rubak,
Raphael Lachieze-Rey,
Jorge Mateu,
Adrian Baddeley
Abstract:
Voronoi intensity estimators, which are non-parametric estimators for intensity functions of point processes, are both parameter-free and adaptive; the intensity estimate at a given location is given by the reciprocal size of the Voronoi/Dirichlet cell containing that location. Their major drawback, however, is that they tend to under-smooth the data in regions where the point density of the obser…
▽ More
Voronoi intensity estimators, which are non-parametric estimators for intensity functions of point processes, are both parameter-free and adaptive; the intensity estimate at a given location is given by the reciprocal size of the Voronoi/Dirichlet cell containing that location. Their major drawback, however, is that they tend to under-smooth the data in regions where the point density of the observed point pattern is high and over-smooth in regions where the point density is low. To remedy this problem, i.e. to find some middle-ground between over- and under-smoothing, we propose an additional smoothing technique for Voronoi intensity estimators for point processes in arbitrary metric spaces, which is based on repeated independent thinnings of the point process/pattern. Through a simulation study we show that our resample-smoothing technique improves the estimation significantly. In addition, we study statistical properties such as unbiasedness and variance, and propose a rule-of-thumb and a data-driven cross-validation approach to choose the amount of thinning/smoothing to apply. We finally apply our proposed intensity estimation scheme to two datasets: locations of pine saplings (planar point pattern) and motor vehicle traffic accidents (linear network point pattern).
△ Less
Submitted 6 July, 2018;
originally announced July 2018.
-
A family of singular integral operators which control the Cauchy transform
Authors:
Petr Chunaev,
Joan Mateu,
Xavier Tolsa
Abstract:
We study the behaviour of singular integral operators $T_{k_t}$ of convolution type on $\mathbb{C}$ associated with the parametric kernels $$ k_t(z):=\frac{(\Re z)^{3}}{|z|^{4}}+t\cdot \frac{\Re z}{|z|^{2}}, \quad t\in \mathbb{R},\qquad k_\infty(z):=\frac{\Re z}{|z|^{2}}\equiv \Re \frac{1}{z},\quad z\in \mathbb{C}\setminus\{0\}. $$ It is shown that for any positive locally finite Borel measure wit…
▽ More
We study the behaviour of singular integral operators $T_{k_t}$ of convolution type on $\mathbb{C}$ associated with the parametric kernels $$ k_t(z):=\frac{(\Re z)^{3}}{|z|^{4}}+t\cdot \frac{\Re z}{|z|^{2}}, \quad t\in \mathbb{R},\qquad k_\infty(z):=\frac{\Re z}{|z|^{2}}\equiv \Re \frac{1}{z},\quad z\in \mathbb{C}\setminus\{0\}. $$ It is shown that for any positive locally finite Borel measure with linear growth the corresponding $L^2$-norm of $T_{k_0}$ controls the $L^2$-norm of $T_{k_\infty}$ and thus of the Cauchy transform. As a corollary, we prove that the $L^2(\mathcal{H}^1\lfloor E)$-boundedness of $T_{k_t}$ with a fixed $t\in (-t_0,0)$, where $t_0>0$ is an absolute constant, implies that $E$ is rectifiable. This is so in spite of the fact that the usual curvature method fails to be applicable in this case. Moreover, as a corollary of our techniques, we provide an alternative and simpler proof of the bi-Lipschitz invariance of the $L^2$-boundedness of the Cauchy transform, which is the key ingredient for the bi-Lipschitz invariance of analytic capacity.
△ Less
Submitted 15 September, 2018; v1 submitted 7 March, 2018;
originally announced March 2018.
-
Prediction of spatial functional random processes: Comparing functional and spatio-temporal kriging approaches
Authors:
Johan Strandberg,
Sara Sjöstedt de Luna,
Jorge Mateu
Abstract:
In this paper, we present and compare functional and spatio-temporal (Sp.T.) kriging approaches to predict spatial functional random processes (which can also be viewed as Sp.T. random processes). Comparisons with respect to computational time and prediction performance via functional cross-validation is evaluated, mainly through a simulation study but also on two real data sets. We restrict compa…
▽ More
In this paper, we present and compare functional and spatio-temporal (Sp.T.) kriging approaches to predict spatial functional random processes (which can also be viewed as Sp.T. random processes). Comparisons with respect to computational time and prediction performance via functional cross-validation is evaluated, mainly through a simulation study but also on two real data sets. We restrict comparisons to Sp.T. kriging versus ordinary kriging for functional data (OKFD), since the more flexible functional kriging approaches, pointwise functional kriging (PWFK) and functional kriging total model, coincide with OKFD in several situations. We contribute with new knowledge by proving that OKFD and PWFK coincide under certain conditions. From the simulation study, it is concluded that the prediction performance for the two kriging approaches in general is rather equal for stationary Sp.T. processes, with a tendency for functional kriging to work better for small sample sizes and Sp.T. kriging to work better for large sample sizes. For non-stationary Sp.T. processes, with a common deterministic time trend and/or time varying variances and dependence structure, OKFD performs better than Sp.T. kriging irrespective of sample size. For all simulated cases, the computational time for OKFD was considerably lower compared to those for the Sp.T. kriging methods.
△ Less
Submitted 15 February, 2018;
originally announced February 2018.
-
Second-order and local characteristics of network intensity functions
Authors:
Matthias Eckardt,
Jorge Mateu
Abstract:
The last decade has witnessed an increase of interest in the spatial analysis of structured point patterns over networks whose analysis is challenging because of geometrical complexities and unique methodological problems. In this context, it is essential to incorporate the network specificity into the analysis as the locations of events are restricted to areas covered by line segments. Relying on…
▽ More
The last decade has witnessed an increase of interest in the spatial analysis of structured point patterns over networks whose analysis is challenging because of geometrical complexities and unique methodological problems. In this context, it is essential to incorporate the network specificity into the analysis as the locations of events are restricted to areas covered by line segments. Relying on concepts originating from graph theory, we extend the notions of first-order network intensity functions to second-order and local network intensity functions. We consider two types of local indicators of network association functions which can be understood as adaptations of the primary ideas of local analysis on the plane. We develop the node-wise and cross-hierarchical type of local functions. A real dataset on urban disturbances is also presented.
△ Less
Submitted 5 December, 2017;
originally announced December 2017.
-
Substationarity in Spatial Point Processes
Authors:
Tonglin Zhang,
Jorge Mateu
Abstract:
The goal of the article is to develop the approach of substationarity to spatial point processes (SPPs). Substationarity is a new concept, which has never been studied in the literature. It means that the distribution of SPPs can only be invariant under location shifts within a linear subspace of the domain. Theoretically, substationarity is a concept between stationariy and nonstationarity, but i…
▽ More
The goal of the article is to develop the approach of substationarity to spatial point processes (SPPs). Substationarity is a new concept, which has never been studied in the literature. It means that the distribution of SPPs can only be invariant under location shifts within a linear subspace of the domain. Theoretically, substationarity is a concept between stationariy and nonstationarity, but it belongs to nonstationarity. To formally propose the approach, the article provides the definition of substationarity and an estimation method for the first-order intensity function. As the linear subspace may be unknown, it recommends using a parametric way to estimate the linear subspace and a nonparametric way to estimate the first-order intensity function, indicating that it is a semiparametric approach. The simulation studies show that both the estimators of the linear subspace and the first-order intensity function are reliable. In an application to a forest wildfire data set, the article concludes that substationarity of wildfire occurrences may be assumed along the longitude, indicating that latitude is a more important factor than longitude in forest wildfire studies.
△ Less
Submitted 8 October, 2017;
originally announced October 2017.
-
The ellipse law: Kirchhoff meets dislocations
Authors:
J. A. Carrillo,
J. Mateu,
M. G. Mora,
L. Rondi,
L. Scardia,
J. Verdera
Abstract:
In this paper we consider a nonlocal energy $I_α$ whose kernel is obtained by adding to the Coulomb potential an anisotropic term weighted by a parameter $α\in \R$. The case $α=0$ corresponds to purely logarithmic interactions, minimised by the celebrated circle law for a quadratic confinement; $α=1$ corresponds to the energy of interacting dislocations, minimised by the semi-circle law. We show t…
▽ More
In this paper we consider a nonlocal energy $I_α$ whose kernel is obtained by adding to the Coulomb potential an anisotropic term weighted by a parameter $α\in \R$. The case $α=0$ corresponds to purely logarithmic interactions, minimised by the celebrated circle law for a quadratic confinement; $α=1$ corresponds to the energy of interacting dislocations, minimised by the semi-circle law. We show that for $α\in (0,1)$ the minimiser can be computed explicitly and is the normalised characteristic function of the domain enclosed by an \emph{ellipse}. To prove our result we borrow techniques from fluid dynamics, in particular those related to Kirchhoff's celebrated result that domains enclosed by ellipses are rotating vortex patches, called \emph{Kirchhoff ellipses}. Therefore we show a surprising connection between vortices and dislocations.
△ Less
Submitted 20 March, 2017;
originally announced March 2017.
-
Covariance Functions for Multivariate Gaussian Fields Evolving Temporally over Planet Earth
Authors:
Alfredo Alegría,
Emilio Porcu,
Reinhard Furrer,
Jorge Mateu
Abstract:
The construction of valid and flexible cross-covariance functions is a fundamental task for modeling multivariate space-time data arising from climatological and oceanographical phenomena. Indeed, a suitable specification of the covariance structure allows to capture both the space-time dependencies between the observations and the development of accurate predictions. For data observed over large…
▽ More
The construction of valid and flexible cross-covariance functions is a fundamental task for modeling multivariate space-time data arising from climatological and oceanographical phenomena. Indeed, a suitable specification of the covariance structure allows to capture both the space-time dependencies between the observations and the development of accurate predictions. For data observed over large portions of planet Earth it is necessary to take into account the curvature of the planet. Hence the need for random field models defined over spheres across time. In particular, the associated covariance function should depend on the geodesic distance, which is the most natural metric over the spherical surface. In this work, we propose a flexible parametric family of matrix-valued covariance functions, with both marginal and cross structure being of the Gneiting type. We additionally introduce a different multivariate Gneiting model based on the adaptation of the latent dimension approach to the spherical context. Finally, we assess the performance of our models through the study of a bivariate space-time data set of surface air temperatures and precipitations.
△ Less
Submitted 21 November, 2017; v1 submitted 21 January, 2017;
originally announced January 2017.