Moment estimation for statistics from marked point processes

In spatial statistics the data typically consist of measurements of some quantity at irregularly scattered locations; in other words, the data form a realization of a marked point process. In this paper, we formulate subsampling estimators of the moments of general statistics computed from marked point process data, and we establish their L ₂-consistency. The variance estimator in particular can be used for the construction of confidence intervals for estimated parameters. A practical data-based method for choosing a subsampling parameter is given and illustrated on a data set. Finite sample simulation examples are also presented.     

Local dimension reduction of summary statistics for likelihood-free inference

Statistics and Computing - Approximate Bayesian computation (ABC) and other likelihood-free inference methods have gained popularity in the last decade, as they allow rigorous statistical inference...     

Modified Wald statistics for generalized linear models

Summary: Wald statistics in generalized linear models are asymptotically ² distributed. The asymptotic chi–squared law of the corresponding quadratic form shows disadvantages with respect to the approximation of the finite–sample distribution. It is shown by means of a comprehensive simulation study that improvements can be achieved by applying simple finite–sample size approximations to the distribution of the quadratic form in generalized linear models. These approximations are based on a ² distribution with an estimated degree of freedom that generalizes an approach by Patnaik and Pearson. Simulation studies confirm that nominal level is maintained with higher accuracy compared to the Wald statistics.     

On special linear identities for order statistics

Linear identities for the distribution functions of order statistics from an iid sample are defined. It is shown that such identities are true for all distributions or to some discrete distributions taking a finite number of values.     

Considerate approaches to constructing summary statistics for ABC model selection

For nearly any challenging scientific problem evaluation of the likelihood is problematic if not impossible. Approximate Bayesian computation (ABC) allows us to employ the whole Bayesian formalism to problems where we can use simulations from a model, but cannot evaluate the likelihood directly. When summary statistics of real and simulated data are compared??rather than the data directly??information is lost, unless the summary statistics are sufficient. Sufficient statistics are, however, not common but without them statistical inference in ABC inferences are to be considered with caution. Previously other authors have attempted to combine different statistics in order to construct (approximately) sufficient statistics using search and information heuristics. Here we employ an information-theoretical framework that can be used to construct appropriate (approximately sufficient) statistics by combining different statistics until the loss of information is minimized. We start from a potentially large number of different statistics and choose the smallest set that captures (nearly) the same information as the complete set. We then demonstrate that such sets of statistics can be constructed for both parameter estimation and model selection problems, and we apply our approach to a range of illustrative and real-world model selection problems.     

Bayesian analysis of spatial point processes in the neighbourhood of Voronoi networks

A model for an inhomogeneous Poisson process with high intensity near the edges of a Voronoi tessellation in 2D or 3D is proposed. The model is analysed in a Bayesian setting with priors on nuclei of the Voronoi tessellation and other model parameters. An MCMC algorithm is constructed to sample from the posterior, which contains information about the unobserved Voronoi tessellation and the model parameters. A major element of the MCMC algorithm is the reconstruction of the Voronoi tessellation after a proposed local change of the tessellation. A simulation study and examples of applications from biology (animal territories) and material science (alumina grain structure) are presented.     

Optimal linear projections for enhancing desired data statistics

Problems involving high-dimensional data, such as pattern recognition, image analysis, and gene clustering, often require a preliminary step of dimension reduction before or during statistical analysis. If one restricts to a linear technique for dimension reduction, the remaining issue is the choice of the projection. This choice can be dictated by desire to maximize certain statistical criteria, including variance, kurtosis, sparseness, and entropy, of the projected data. Motivations for such criteria comes from past empirical studies of statistics of natural and urban images. We present a geometric framework for finding projections that are optimal for obtaining certain desired statistical properties. Our approach is to define an objective function on spaces of orthogonal linear projections—Stiefel and Grassmann manifolds, and to use gradient techniques to optimize that function. This construction uses the geometries of these manifolds to perform the optimization. Experimental results are presented to demonstrate these ideas for natural and facial images.     

A copula model for marked point processes

Many chronic diseases feature recurring clinically important events. In addition, however, there often exists a random variable which is realized upon the occurrence of each event reflecting the severity of the event, a cost associated with it, or possibly a short term response indicating the effect of a therapeutic intervention. We describe a novel model for a marked point process which incorporates a dependence between continuous marks and the event process through the use of a copula function. The copula formulation ensures that event times can be modeled by any intensity function for point processes, and any multivariate model can be specified for the continuous marks. The relative efficiency of joint versus separate analyses of the event times and the marks is examined through simulation under random censoring. An application to data from a recent trial in transfusion medicine is given for illustration.     

Two-step estimation for inhomogeneous spatial point processes

Summary.  The paper is concerned with parameter estimation for inhomogeneous spatial point processes with a regression model for the intensity function and tractable second-order properties ( K -function). Regression parameters are estimated by using a Poisson likelihood score estimating function and in the second step minimum contrast estimation is applied for the residual clustering parameters. Asymptotic normality of parameter estimates is established under certain mixing conditions and we exemplify how the results may be applied in ecological studies of rainforests.     

Bayesian inference for pairwise interacting point processes

Pairwise interacting point processes are commonly used to model spatial point patterns. To perform inference, the established frequentist methods can produce good point estimates when the interaction in the data is moderate, but some methods may produce severely biased estimates when the interaction in strong. Furthermore, because the sampling distributions of the estimates are unclear, interval estimates are typically obtained by parametric bootstrap methods. In the current setting however, the behavior of such estimates is not well understood. In this article we propose Bayesian methods for obtaining inferences in pairwise interacting point processes. The requisite application of Markov chain Monte Carlo (MCMC) techniques is complicated by an intractable function of the parameters in the likelihood. The acceptance probability in a Metropolis-Hastings algorithm involves the ratio of two likelihoods evaluated at differing parameter values. The intractable functions do not cancel, and hence an intractable ratio r must be estimated within each iteration of a Metropolis-Hastings sampler. We propose the use of importance sampling techniques within MCMC to address this problem. While r may be estimated by other methods, these, in general, are not readily applied in a Bayesian setting. We demonstrate the validity of our importance sampling approach with a small simulation study. Finally, we analyze the Swedish pine sapling dataset (Strand 1972) and contrast the results with those in the literature.     

Moment expressions and summary statistics for the complete and truncated weibull distribution

Traditionally, the moments of the Weibull distribution have been calculated using the standard Weibull (Johnson and Kotz, 1970) . This article will expand on that idea and cover the truncated cases for the standard Weibull distributions. Also, the same techniques used for the standard form will be used to derive the moment expressions for the three-parameter complete and truncated Weibull distributions. The summary statistics are then calculated from the moment expressions. Weibull moments involve the gamma and incomplete gamma functions.     

Control charts based on linear combinations of order statistics

The last 20 years have seen an increasing emphasis on statistical process control as a practical approach to reducing variability in industrial applications. Control charts are used to detect problems such as outliers or excess variability in subgroup means that may have a special cause. We describe an approach to the computation of control limits for exponentially weighted moving average control charts where the usual statistics in classical charts are replaced by linear combinations of order statistics; in particular, the trimmed mean and Gini's mean difference instead of the mean and range, respectively. Control limits are derived, and simulated average run length experiments show the trimmed control charts to be less influenced by extreme observations than their classical counterparts, and lead to tighter control limits. An example is given that illustrates the benefits of the proposed charts. parameters; see, for example, Hunter (1986) and Montgomery (1996). On the other hand, EWMA charts have been shown to be more efficient than Shewharttype charts in detecting small shifts in the process mean; see, for example, Ng & Case (1989), Crowder (1989), Lucas & Saccucci (1990), Amin & Searcy (1991) and Wetherill & Brown (1991). In fact, the EWMA control chart has become popular for monitoring a process mean; see Hunter (1986) for a good discussion. More recently, EWMA charts have been developed for monitoring process variability;     

Control charts based on linear combinations of order statistics

The last 20 years have seen an increasing emphasis on statistical process control as a practical approach to reducing variability in industrial applications. Control charts are used to detect problems such as outliers or excess variability in subgroup means that may have a special cause. We describe an approach to the computation of control limits for exponentially weighted moving average control charts where the usual statistics in classical charts are replaced by linear combinations of order statistics; in particular, the trimmed mean and Gini's mean difference instead of the mean and range, respectively. Control limits are derived, and simulated average run length experiments show the trimmed control charts to be less influenced by extreme observations than their classical counterparts, and lead to tighter control limits. An example is given that illustrates the benefits of the proposed charts. parameters; see, for example, Hunter (1986) and Montgomery (1996). On the other hand, EWMA charts have been shown to be more efficient than Shewharttype charts in detecting small shifts in the process mean; see, for example, Ng & Case (1989), Crowder (1989), Lucas & Saccucci (1990), Amin & Searcy (1991) and Wetherill & Brown (1991). In fact, the EWMA control chart has become popular for monitoring a process mean; see Hunter (1986) for a good discussion. More recently, EWMA charts have been developed for monitoring process variability;     

Edgeworth expansion for signed linear rank statistics with regression constants

An Edgeworth expansion with remainder o(N^?1) is obtained for signed linear rank statistics under suitable assumptions. The theorem is proved for a wide class of score generating functions including the Chi-quantile function by adapting van Zwet's methodand Does's conditioning arguments.     

Change-of-variance problem for linear processes with long memory

In this paper we investigate the asymptotic properties of the test statistics for detecting change-points in the variance of infinite moving average sequences with long memory. This research is partly supported by NSFC Grants and SRF for ROCS, SEM.     

Location estimators based on linear combinations of modified order statistics

Let Xj,…, X be i.i.d random variables with common distribution function F(x-0), and let a(u) be a function defined on [0,1], For each t$$$R define the t-order statistics as: X. (t) = X. if there in k exist exactly (i- 1) X. !s such that |X.-t|greater|X.-t| and define the variable T (t) = n 2. , a(i/n) X, (t) . We consider estimates of 8 defined as solutions of the eauaticn T (8) =6 , and     

Testing for a change point in linear regression models

In this paper, the Schwarz Information Criterion (SIC) is proposed to locate a change point in the simple linear regression model, as well as in the multiple linear regression model. The method is then applied to a financial data set, and a change point is successfully detected.