Simulations of Some Doubly Stochastic Poisson Point Processes

Computer simulations of point processes are important either to verify the results of certain theoretical calculations that can be very awkward at times or to obtain practical results when these calculations become almost impossible. One of the most common methods for the simulation of nonstationary Poisson processes is random thinning. Its extension when the intensity becomes random (doubly stochastic Poisson processes) depends on the structure of this intensity. If the random density takes only discrete values, which is a common situation in many physical problems where quantum mechanics introduces discrete states, it is shown that the thinning method can be applied without error. We study in particular the case of binary density and present the kind of theoretical calculations that then become possible. The results of various experiments realized with data obtained by simulation show a fairly good agreement with the theoretical calculations.     

Inversion of the Renewal Density with Dead-time

Stationary renewal point processes are defined by the probability distribution of the distances between successive points (lifetimes) that are independent and identically distributed random variables. For some applications it is also interesting to define the properties of a renewal process by using the renewal density. There are well-known expressions of this density in terms of the probability density of the lifetimes. It is more difficult to solve the inverse problem consisting in the determination of the density of the lifetimes in terms of the renewal density. Theoretical expressions between their Laplace transforms are available but the inversion of these transforms is often very difficult to obtain in closed form. We show that this is possible for renewal processes presenting a dead-time property characterized by the fact that the renewal density is zero in an interval including the origin. We present the principle of a recursive method allowing the solution of this problem and we apply this method to the case of some processes with input dead-time. Computer simulations on Poisson and Erlang (2) processes show quite good agreement between theoretical calculations and experimental measurements on simulated data.     

Point Processes with a Generalized Order Statistic Property

Mixed Poisson processes are characterized by a well-known order statistic property: their occurrence times are distributed like ordinary uniform order statistics, given the state of the process at a certain time. We study a generalization of this property using a generalized model of ordered random variables, including a sequence of real parameters.     

Output Dead-Time in Point Processes

When events of temporal point processes are too close to each other they can be erased by dead-time effects. Among various possible mechanisms of dead-time, the output dead-time is the most important. Dead-time effects modify the statistical properties of point processes and some of these modifications are analyzed in this article. To do so, we note that a point process is defined by the distance between its successive points called life-time which constitutes a discrete time positive signal. The dead-time mechanism is a system which transforms such a signal into another discrete time positive signal. Except in very specific cases this transformation cannot be expressed in closed form. We show, however, that it can be written in a recursive form analogous to the state representation of systems. By using this recursion, various statistical properties of point processes with dead-time are analyzed in computer experiments. In this study, we focus on the probability distribution of the intervals between points and the coincidence function which describes the second-order properties of the point process. For the rare processes where theoretical calculations are possible there is an excellent agreement between experiment and theory.     

On Random and Gibbsian Particle Motions for Point Processes Evolving in Space and Time

This article presents an analysis of space-time interdependencies of spatial point processes considering random and deterministic Gibbsian point motions caused by repulsion effects between particles. Two deterministic models of Gibbsian motions are considered by formulating a constant (i.e., Strauss-like) and a linear interaction motion functions. Given that theoretical development of continuous space-time stochastic processes are mathematically intractable, we have mainly based our analysis on numerical simulations. Our results suggest that to fully understand such complex dynamics, the analysis of purely spatial patterns should be combined with their interactions in the space-time domain. Otherwise, analysis of pure spacial patterns may not fully explain the real mechanism generating such dynamical configurations. We highlight that adding movement to sedentary points opens new areas of application and research to study biological phenomena, where particles not only evolve through time but also can change spatial positions in terms of their neighbor locations.     

Re-colouring the Intensity-Based Bootstrap for Point Processes

Abstract

A method for obtaining bootstrapping replicates for one-dimensional point processes is presented. The method involves estimating the conditional intensity of the process and computing residuals. The residuals are bootstrapped using a block bootstrap and used, together with the conditional intensity, to define the bootstrap realizations. The method is applied to the estimation of the cross-intensity function for data arising from a reaction time experiment.     

Statistical Inference for Transformation Inhomogeneous Point Processes

Abstract.  Statistical inference for exponential inhomogeneous Markov point processes by transformation is discussed. It is argued that the inhomogeneity parameter can be estimated, using a partial likelihood based on an inhomogeneous Poisson point process. The inhomogeneity parameter can thereby be estimated without taking the interaction into account, which simplifies the statistical analysis considerably. Data analysis and simulation experiments support the results.     

A Study for Missing Values in PINAR(1)T Processes

In this paper, we propose several approaches to estimate the parameters of the periodic first-order integer-valued autoregressive process with period T (PINAR(1)_T) in the presence of missing data. By using incomplete data, we propose two approaches that are based on the conditional expectation and conditional likelihood to estimate the parameters of interest. Then we study three kinds of imputation methods for the missing data. The performances of these approaches are compared via simulations.     

Intensity Estimation for Spatial Point Processes Observed with Noise

Abstract.  This article introduces a kernel estimator of the intensity function of spatial point processes taking into account location errors. The asymptotic properties of the estimator are derived and a bandwidth selection procedure is described. A simulation study compares our results with that of the classical kernel estimator and shows that the edge-corrected deconvoluting kernel estimator is more appropriate.     

A Note on Continuous Spatial-Temporal Dynamics of Stochastic Processes

We propose a general form to analyze the space-time interdependency of continuous space-time stochastic processes. We present a new space-time approach based on the intensity function of the underlying point process. These formulations can be, to some extent, analytically solved to obtain explicit formulae of interest. We define a general function that controls the space-time interaction and allows for closed forms depending on the particular choice of several mathematical tools playing a role in this interaction function. In particular, we make use of copulas and Laplace transforms to provide interesting examples of the dynamics of the random intensity function and, in turn, of the number of points contained in a given region.     

A Note on Empirical Process Methods in the Theory of Poisson Point Processes

For a sample taken from an i.i.d. sequence of Poisson point processes with not necessarily finite unknown intensity measure the arithmetic mean is shown to be an estimator which is consistent uniformly on certain classes of functions. The method is a reduction to the case of finite intensity measure, which in turn can be dealt with using empirical process methods. A functional central limit theorem is also established in this context.     

Gamma Kernel Intensity Estimation in Temporal Point Processes

In this article, we propose a nonparametric approach for estimating the intensity function of temporal point processes based on kernel estimators. In particular, we use asymmetric kernel estimators characterized by the gamma distribution, in order to describe features of observed point patterns adequately. Some characteristics of these estimators are analyzed and discussed both through simulated results and applications to real data from different seismic catalogs.     

On Some Applications of Poisson Point Processes

Poisson point processes play important role in various domains of Probability Theory and Mathematical Statistics. In this article, we investigate only two applications of Poisson point processes: a generated white noise problem and parameters estimation problem. This work continues the investigations started in paper Egorov and Kondybaev (2009 Egorov , V. A. , Kondybaev , N. S. ( 2009 ). On the estimation of a signal covered by background om Poisson noise . Methods and programs of data processing 4 : 75 – 81 . [Google Scholar]).     

Optimal Alarm Systems for Count Processes

In many phenomena described by stochastic processes, the implementation of an alarm system becomes fundamental to predict the occurrence of future events. In this work we develop an alarm system to predict whether a count process will upcross a certain level and give an alarm whenever the upcrossing level is predicted. We consider count models with parameters being functions of covariates of interest and varying on time. This article presents classical and Bayesian methodology for producing optimal alarm systems. Both methodologies are illustrated and their performance compared through a simulation study. The work finishes with an empirical application to a set of data concerning the number of sunspot on the surface of the sun.     

Takacs–Fiksel Method for Stationary Marked Gibbs Point Processes

Abstract. This article studies a method to estimate the parameters governing the distribution of a stationary marked Gibbs point process. This procedure, known as the Takacs–Fiksel method, is based on the estimation of the left and right hand sides of the Georgii–Nguyen–Zessin formula and leads to a family of estimators due to the possible choices of test functions. We propose several examples illustrating the interest and flexibility of this procedure. We also provide sufficient conditions based on the model and the test functions to derive asymptotic properties (consistency and asymptotic normality) of the resulting estimator. The different assumptions are discussed for exponential family models and for a large class of test functions. A short simulation study is proposed to assess the correctness of the methodology and the asymptotic results.     

Smoothed Extreme Value Estimators of Non-Uniform Point Processes Boundaries with Application to Star-Shaped Supports Estimation

We address the problem of estimating the edge of a bounded set in ? ^d given a random set of points drawn from the interior. Our method is based on a transformation of estimators dedicated to uniform point processes and obtained by smoothing some of its bias corrected extreme points. An application to the estimation of star-shaped supports is presented.     

A Pointwise Estimator for the k-Fold Convolution of a Distribution Function

ABSTRACT

This article is concerned with some parametric and nonparametric estimators for the k-fold convolution of a distribution function. An alternative estimator is proposed and its unbiasedness, asymptotic unbiasedness, and consistency properties are investigated. The asymptotic normality of this estimator is established. Some applications of the estimator are given in renewal processes. Finally, the computational procedures are described and the relative performance of these estimators for small sample sizes is investigated by a simulation study.     

Estimation for Continuous Branching Processes

The maximum-likelihood estimator for the curved exponential family given by continuous branching processes with immigration is investigated. These processes originated from population biology but also model the dynamics of interest rates and development of the state of technology in economics. It is proved that in contrast to branching processes with discrete space and/or time the MLE gives a unified approach to the inference. In order to include singular subdomains of the parameter space we modify the MLE slightly. Consistency and asymptotic normality for the MLE are considered. Concerning the asymptotic theory of the experiments, all three properties LAQ, LAN, and LAMN occur for different submodels     

A semiparametric additive rates model for recurrent event data

Recurrent event data often arise in biomedical studies, with examples including hospitalizations, infections, and treatment failures. In observational studies, it is often of interest to estimate the effects of covariates on the marginal recurrent event rate. The majority of existing rate regression methods assume multiplicative covariate effects. We propose a semiparametric model for the marginal recurrent event rate, wherein the covariates are assumed to add to the unspecified baseline rate. Covariate effects are summarized by rate differences, meaning that the absolute effect on the rate function can be determined from the regression coefficient alone. We describe modifications of the proposed method to accommodate a terminating event (e.g., death). Proposed estimators of the regression parameters and baseline rate are shown to be consistent and asymptotically Gaussian. Simulation studies demonstrate that the asymptotic approximations are accurate in finite samples. The proposed methods are applied to a state-wide kidney transplant data set.     

The Input Evaluation of Generalized Bernoulli Processes for Salary Lines Construction by Means of Continuous Time Generalized Non-Homogeneous Semi-Markov Processes

Salary line forecasting assumes a relevant role in manpower and in pension funds Previously, the authors presented a generalized Bernoulli process, useful for forecasting the evolution of salary lines, taking into account the salary costs, the number of workers at each rank and the probability transitions between the ranks. The problem with applying this model is constructing the probability of transition between the grades. In this article, we will present a model that allows obtaining these probabilities by means of the solution of the evolution equation of a generalization of continuous time non-homogeneous semi-Markov processes.