-
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.
-
The Spatial Structure of Young Stellar Clusters. I. Subclusters
Authors:
Michael A. Kuhn,
Eric D. Feigelson,
Konstantin V. Getman,
Adrian J. Baddeley,
Patrick S. Broos,
Alison Sills,
Matthew R. Bate,
Matthew S. Povich,
Kevin L. Luhman,
Heather A. Busk,
Tim Naylor,
Robert R. King
Abstract:
The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of…
▽ More
The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of young stars using finite mixture models -- collections of isothermal ellipsoids that model individual subclusters. Maximum likelihood estimation is used to estimate the model parameters, and the Akaike Information Criterion is used to determine the number of subclusters. This procedure often successfully finds famous subclusters, such as the BN/KL complex behind the Orion Nebula Cluster and the KW-object complex in M17. A catalog of 142 subclusters is presented, with 1 to 20 subclusters per region. The subcluster core radius distribution for this sample is peaked at 0.17 pc with a standard deviation of 0.43 dex, and subcluster core radius is negatively correlated with gas/dust absorption of the stars -- a possible age effect. Based on the morphological arrangements of subclusters, we identify four classes of spatial structure: long chains of subclusters, clumpy structures, isolated clusters with a core-halo structure, and isolated clusters well fit by a single isothermal ellipsoid.
△ Less
Submitted 17 March, 2014;
originally announced March 2014.
-
Overview of the Massive Young Star-Forming Complex Study in Infrared and X-ray (MYStIX) Project
Authors:
Eric D. Feigelson,
Leisa K. Townsley,
Patrick S. Broos,
Heather A. Busk,
Konstantin V. Getman,
Robert King,
Michael A. Kuhn,
Tim Naylor,
Matthew Povich,
Adrian Baddeley,
Matthew Bate,
Remy Indebetouw,
Kevin Luhman,
Mark McCaughrean,
Julian Pittard,
Ralph Pudritz,
Alison Sills,
Yong Song,
James Wadsley
Abstract:
MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined…
▽ More
MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.
△ Less
Submitted 17 September, 2013;
originally announced September 2013.
-
Variational estimators for the parameters of Gibbs point process models
Authors:
Adrian Baddeley,
David Dereudre
Abstract:
This paper proposes a new estimation technique for fitting parametric Gibbs point process models to a spatial point pattern dataset. The technique is a counterpart, for spatial point processes, of the variational estimators for Markov random fields developed by Almeida and Gidas. The estimator does not require the point process density to be hereditary, so it is applicable to models which do not h…
▽ More
This paper proposes a new estimation technique for fitting parametric Gibbs point process models to a spatial point pattern dataset. The technique is a counterpart, for spatial point processes, of the variational estimators for Markov random fields developed by Almeida and Gidas. The estimator does not require the point process density to be hereditary, so it is applicable to models which do not have a conditional intensity, including models which exhibit geometric regularity or rigidity. The disadvantage is that the intensity parameter cannot be estimated: inference is effectively conditional on the observed number of points. The new procedure is faster and more stable than existing techniques, since it does not require simulation, numerical integration or optimization with respect to the parameters.
△ Less
Submitted 23 July, 2013;
originally announced July 2013.
-
MYStIX First Results: Spatial Structures of Massive Young Stellar Clusters
Authors:
Michael A. Kuhn,
Adrian Baddeley,
Eric D. Feigelson,
Konstantin V. Getman,
Patrick S. Broos,
Leisa K. Townsley,
Matthew S. Povich,
Tim Naylor,
Robert R. King,
Heather A. Busk,
Kevin L. Luhman
Abstract:
Observations of the spatial distributions of young stars in star-forming regions can be linked to the theory of clustered star formation using spatial statistical methods. The MYStIX project provides rich samples of young stars from the nearest high-mass star-forming regions. Maps of stellar surface density reveal diverse structure and subclustering. Young stellar clusters and subclusters are fit…
▽ More
Observations of the spatial distributions of young stars in star-forming regions can be linked to the theory of clustered star formation using spatial statistical methods. The MYStIX project provides rich samples of young stars from the nearest high-mass star-forming regions. Maps of stellar surface density reveal diverse structure and subclustering. Young stellar clusters and subclusters are fit with isothermal spheres and ellipsoids using the Bayesian Information Criterion to estimate the number of subclusters. Clustering is also investigated using Cartwright and Whitworth's Q statistic and the inhomogeneous two-point correlation function. Mass segregation is detected in several cases, in both centrally concentrated and fractally structured star clusters, but a few clusters are not mass segregated.
△ Less
Submitted 16 August, 2012;
originally announced August 2012.
-
Score, Pseudo-Score and Residual Diagnostics for Spatial Point Process Models
Authors:
Adrian Baddeley,
Ege Rubak,
Jesper Møller
Abstract:
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…
▽ More
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.
△ Less
Submitted 17 May, 2012;
originally announced May 2012.
-
Estimating the J function without edge correction
Authors:
Adrian Baddeley,
Martin Kerscher,
Katja Schladitz,
Bryan T. Scott
Abstract:
The interaction between points in a spatial point process can be measured by its empty space function F, its nearest-neighbour distance distribution function G, and by combinations such as the J-function $J = (1-G)/(1-F)$. The estimation of these functions is hampered by edge effects: the uncorrected, empirical distributions of distances observed in a bounded sampling window W give severely bias…
▽ More
The interaction between points in a spatial point process can be measured by its empty space function F, its nearest-neighbour distance distribution function G, and by combinations such as the J-function $J = (1-G)/(1-F)$. The estimation of these functions is hampered by edge effects: the uncorrected, empirical distributions of distances observed in a bounded sampling window W give severely biased estimates of F and G. However, in this paper we show that the corresponding {\em uncorrected} estimator of the function $J=(1-G)/(1-F)$ is approximately unbiased for the Poisson case, and is useful as a summary statistic. Specifically, consider the estimate $\hat{J}_W$ of J computed from uncorrected estimates of F and G. The function $J_W(r)$, estimated by $\hat{J}_W$, possesses similar properties to the J function, for example $J_W(r)$ is identically 1 for Poisson processes. This enables direct interpretation of uncorrected estimates of J, something not possible with uncorrected estimates of either F, G or K. We propose a Monte Carlo test for complete spatial randomness based on testing whether $J_W(r)(r)\equiv 1$. Computer simulations suggest this test is at least as powerful as tests based on edge corrected estimators of J.
△ Less
Submitted 4 October, 1999;
originally announced October 1999.