Practical Maximum Pseudolikelihood for Spatial Point Patterns(with Discussion)

This paper describes a technique for computing approximate maximum pseudolikelihood estimates of the parameters of a spatial point process. The method is an extension of Berman & Turner's (1992) device for maximizing the likelihoods of inhomogeneous spatial Poisson processes. For a very wide class of spatial point process models the likelihood is intractable, while the pseudolikelihood is known explicitly, except for the computation of an integral over the sampling region. Approximation of this integral by a finite sum in a special way yields an approximate pseudolikelihood which is formally equivalent to the (weighted) likelihood of a loglinear model with Poisson responses. This can be maximized using standard statistical software for generalized linear or additive models, provided the conditional intensity of the process takes an 'exponential family' form. Using this approach a wide variety of spatial point process models of Gibbs type can be fitted rapidly, incorporating spatial trends, interaction between points, dependence on spatial covariates, and mark information.     

A Note on the Innovation Distribution of a Gamma Distributed Autoregressive Process

A representation of the innovation random variable for a gamma distributed first-order autoregressive process was found by Lawrance (1982) in the form of a compound Poisson distribution, connected with a shot-noise process. In this note we simplify the representation of Lawrance by providing a direct representation in terms of density functions.     

Biomedical applications for a generalized linear functional Poisson model

Many medical and biological studies involve response in the form of Poisson counts which can bemodelled using explanatory variables which also arise from count data. If the explanatory variables are observable without error (also as Poisson counts) we have a generalized linear model with a logarithmic link function and Poisson error structure. If,however, some of the explanatory variables are not directly observable, but arise with superimposed errors (again of Poisson form), the model is of a new type:a generalised linear functional Poisson model. In this paper,maximum likelihood estimates of the parameters of this model are determined along with the information matrix which (on noting its particular patterned form) is amenable to inversion in explicit form. Methods are proposed of an iterative type for computing estimates of the parameters and of their variational properties (e.g. standard errors) for this model, which also has application in other fields such as road traffic studies.     

Optimal Maintenance Time for Repairable Systems Under Two Types of Failures

Determination of preventive maintenance is an important issue for systems under degradation. A typical maintenance policy calls for complete preventive repair actions at pre-scheduled times and minimal repair actions whenever a failure occurs. Under minimal repair, failures are modeled according to a non homogeneous Poisson process. A perfect preventive maintenance restores the system to the as good as new condition. The motivation for this article was a maintenance data set related to power switch disconnectors. Two different types of failures could be observed for these systems according to their causes. The major difference between these types of failures is their costs. Assuming that the system will be in operation for an infinite time, we find the expected cost per unit of time for each preventive maintenance policy and hence obtain the optimal strategy as a function of the processes intensities. Assuming a parametrical form for the intensity function, large sample estimates for the optimal maintenance check points are obtained and discussed.     

Local polynomial estimation of Poisson intensities in the presence of reporting delays

Summary.  The system for monitoring suicides in Hong Kong has considerable delays in reporting as the cause of death needs to be determined by a coroner's investigation. However, timely estimates of suicide rates are desirable to assist in the formulation of public health policies. This motivated us to develop a non-parametric procedure to estimate the intensity function of a Poisson process in the presence of reporting delays. We give closed form estimators of the Poisson intensity and the delay distribution, conduct simulation studies to evaluate the method proposed and derive their asymptotic properties. The method proposed is applied to estimate the intensity of suicide in Hong Kong.     

Approximating the distribution of a renewal process using generalized poisson distributions

The exact distribution of a renewal counting process is not easy to compute and is rarely of closed form. In this article, we approximate the distribution of a renewal process using families of generalized Poisson distributions. We first compute approximations to the first several moments of the renewal process. In some cases, a closed form approximation is obtained. It is found that each family considered has its own strengths and weaknesses. Some new families of generalized Poisson distributions are recommended. Theorems are obtained determining when these variance to mean ratios are less than (or exceed) one without having to find the mean and variance. Some numerical comparisons are also made.     

What is the effective sample size of a spatial point process?

Point process models are a natural approach for modelling data that arise as point events. In the case of Poisson counts, these may be fitted easily as a weighted Poisson regression. Point processes lack the notion of sample size. This is problematic for model selection, because various classical criteria such as the Bayesian information criterion (BIC) are a function of the sample size, n, and are derived in an asymptotic framework where n tends to infinity. In this paper, we develop an asymptotic result for Poisson point process models in which the observed number of point events, m, plays the role that sample size does in the classical regression context. Following from this result, we derive a version of BIC for point process models, and when fitted via penalised likelihood, conditions for the LASSO penalty that ensure consistency in estimation and the oracle property. We discuss challenges extending these results to the wider class of Gibbs models, of which the Poisson point process model is a special case.     

Hierarchical Bayesian analysis of a discrete time series of Poisson counts

A problem involving non-stationary, discrete-time series of counts from a Poisson process with a varying but smooth intensity function is studied. A smoothness prior for the underlying intensity process is modelled using the hierarchical Bayesian approach, which is shown to provide an AR(1) representation for the intensity process. Since conjugate priors are not assumed, analytic derivation of estimates and predictions of the Poisson series are not available. Some reasonably good approximations are given and illustrated using data on British road casualties before and after the introduction of the seatbelt law.     

Bridging the Gap: A Generalized Stochastic Process for Count Data

The Bernoulli and Poisson processes are two popular discrete count processes; however, both rely on strict assumptions. We instead propose a generalized homogenous count process (which we name the Conway–Maxwell–Poisson or COM-Poisson process) that not only includes the Bernoulli and Poisson processes as special cases, but also serves as a flexible mechanism to describe count processes that approximate data with over- or under-dispersion. We introduce the process and an associated generalized waiting time distribution with several real-data applications to illustrate its flexibility for a variety of data structures. We consider model estimation under different scenarios of data availability, and assess performance through simulated and real datasets. This new generalized process will enable analysts to better model count processes where data dispersion exists in a more accommodating and flexible manner.     

Likelihood Inference for Unions of Interacting Discs

Abstract. This is probably the first paper which discusses likelihood inference for a random set using a germ‐grain model, where the individual grains are unobservable, edge effects occur and other complications appear. We consider the case where the grains form a disc process modelled by a marked point process, where the germs are the centres and the marks are the associated radii of the discs. We propose to use a recent parametric class of interacting disc process models, where the minimal sufficient statistic depends on various geometric properties of the random set, and the density is specified with respect to a given marked Poisson model (i.e. a Boolean model). We show how edge effects and other complications can be handled by considering a certain conditional likelihood. Our methodology is illustrated by analysing Peter Diggle's heather data set, where we discuss the results of simulation‐based maximum likelihood inference and the effect of specifying different reference Poisson models.     

Semiparametric estimation for count data through weighted distributions

This paper is concerned with semiparametric discrete kernel estimators when the unknown count distribution can be considered to have a general weighted Poisson form. The estimator is constructed by multiplying the Poisson estimate with a nonparametric discrete kernel-type estimate of the Poisson weight function. Comparisons are then carried out with the ordinary discrete kernel probability mass function estimators. The Poisson weight function is thus a local multiplicative correction factor, and is considered as the uniform measure to detect departures from the equidispersed Poisson distribution. In this way, the effects of dispersion and zero-proportion with respect to the standard Poisson distribution are also minimized. This method of estimation is also applied to the weighted binomial form for the count distribution having a finite support. The proposed estimators, in addition to being simple, easy-to-implement and effective, also outperform the competing nonparametric and parametric estimators in finite-sample situations. Two examples illustrate this new semiparametric estimation.     

Exact Bayesian modeling for bivariate Poisson data and extensions

Bivariate count data arise in several different disciplines (epidemiology, marketing, sports statistics just to name a few) and the bivariate Poisson distribution being a generalization of the Poisson distribution plays an important role in modelling such data. In the present paper we present a Bayesian estimation approach for the parameters of the bivariate Poisson model and provide the posterior distributions in closed forms. It is shown that the joint posterior distributions are finite mixtures of conditionally independent gamma distributions for which their full form can be easily deduced by a recursively updating scheme. Thus, the need of applying computationally demanding MCMC schemes for Bayesian inference in such models will be removed, since direct sampling from the posterior will become available, even in cases where the posterior distribution of functions of the parameters is not available in closed form. In addition, we define a class of prior distributions that possess an interesting conjugacy property which extends the typical notion of conjugacy, in the sense that both prior and posteriors belong to the same family of finite mixture models but with different number of components. Extension to certain other models including multivariate models or models with other marginal distributions are discussed.     

A New Goodness-of-Fit Test Based on the Laplace Statistic for a Large Class of NHPP Models

ABSTRACT

Nonhomogeneous Poisson processes (NHPP) provide many models for hardware and software reliability analysis. In order to get an appropriate NHPP model, goodness-of-Fit (GOF for short) tests have to be carried out. For the power-law processes, lots of GOF tests have been developed. For other NHPP models, only the Conditional Probability Integral Transformation (CPIT) test has been proposed. However, the CPIT test is less powerful and cannot be applied to some NHPP models. This article proposes a general GOF test based on the Laplace statistic for a large class of NHPP models with intensity functions of the form αλ(t, β). The simulation results show that this test is more powerful than CPIT test.     

Joint modelling of event counts and survival times

Summary.  In studies of recurrent events, such as epileptic seizures, there can be a large amount of information about a cohort over a period of time, but current methods for these data are often unable to utilize all of the available information. The paper considers data which include post-treatment survival times for individuals experiencing recurring events, as well as a measure of the base-line event rate, in the form of a pre-randomization event count. Standard survival analysis may treat this pre-randomization count as a covariate, but the paper proposes a parametric joint model based on an underlying Poisson process, which will give a more precise estimate of the treatment effect.     

A class of correlated weighted Poisson processes

In this paper, we propose new classes of correlated Poisson processes and correlated weighted Poisson processes on the interval [0,1], which generalize the class of weighted Poisson processes defined by Balakrishnan and Kozubowski (2008), by incorporating a dependence structure between the standard uniform variables used in the construction. In this manner, we obtain another process that we refer to as correlated weighted Poisson process. Various properties of this process such as marginal and joint distributions, stationarity of the increments, moments, and the covariance function, are studied. The results are then illustrated through some examples, which include processes with length-biased Poisson, exponentially weighted Poisson, negative binomial, and COM-Poisson distributions.     

Bias reduction for the maximum likelihood estimator of the doubly-truncated Poisson distribution

ABSTRACT

We derive an analytic expression for the bias of the maximum likelihood estimator of the parameter in a doubly-truncated Poisson distribution, which proves highly effective as a means of bias correction. For smaller sample sizes, our method outperforms the alternative of bias correction via the parametric bootstrap. Bias is of little concern in the positive Poisson distribution, the most common form of truncation in the applied literature. Bias appears to be the most severe in the doubly-truncated Poisson distribution, when the mean of the distribution is close to the right (upper) truncation.     

A density function connected with a non-negative self-decomposable random variable

The innovation random variable for a non-negative self-decomposable random variable can have a compound Poisson distribution. In this case, we provide the density function for the compounded variable. When it does not have a compound Poisson representation, there is a straightforward and easily available compound Poisson approximation for which the density function of the compounded variable is also available. These results can be used in the simulation of Ornstein–Uhlenbeck type processes with given marginal distributions. Previously, simulation of such processes used the inverse of the corresponding tail Lévy measure. We show this approach corresponds to the use of an inverse cdf method of a certain distribution. With knowledge of this distribution and hence density function, the sampling procedure is open to direct sampling methods.     

Property of bivariate poisson distribution and its application to stochastic processes

Assuming that the random vectors X ₁ and X ₂ have independent bivariate Poisson distributions, the conditional distribution of X ₁ given X ₁?+?X ₂?=?n is obtained. The conditional distribution turns out to be a finite mixture of distributions involving univariate binomial distributions and the mixing proportions are based on a bivariate Poisson (BVP) distribution. The result is used to establish two properties of a bivariate Poisson stochastic process which are the bivariate extensions of the properties for a Poisson process given by Karlin, S. and Taylor, H. M. (1975). A First Course in Stochastic Processes, Academic Press, New York.     

A Third Order Point Process Characteristic

Second order characteristics, in particular Ripley's K -function, are widely used for the statistical analysis of point patterns. We examine a third order analogue, namely the mean number T ( r ) of r -close triples of points per unit area. Equivalently this is the expected number of r -close point pairs in an r -neighbourhood of the typical point. Various estimators for this function are proposed and compared, and we give an explicit formula for the isotropic edge correction. The theoretical value of T seems to be unobtainable for most point process models apart from the homogeneous Poisson process. However, simulation studies show that the function T discriminates well between different types of point processes. In particular it detects a clear difference between the Poisson process and the Baddeley–Silverman cell process whereas the K -functions for these processes coincide.