Parameter estimation and model selection for Neyman-Scott point processes

This paper proposes an approximative method for maximum likelihood estimation of parameters of Neyman-Scott and similar point processes. It is based on the point pattern resulting from forming all difference points of pairs of points in the window of observation. The intensity function of this constructed point process can be expressed in terms of second-order characteristics of the original process. This opens the way to parameter estimation, if the difference pattern is treated as a non-homogeneous Poisson process. The computational feasibility and accuracy of this approach is examined by means of simulated data. Furthermore, the method is applied to two biological data sets. For these data, various cluster process models are considered and compared with respect to their goodness-of-fit.     

Nonparametric and parametric estimation for a linear germination-growth model

Seeds are planted on the interval [0, L] at various locations. Each seed has a location x and a potential germination time t epsilon [0, infinity), and it is assumed that the collection of such (x, t) pairs forms a Poisson process in [0, L] x [0, infinity) with intensity measure dxd lambda(t). From each seed that germinates, an inhibiting region grows bidirectionally at rate 2v. These regions inhibit germination of any seed in the region with a later potential germination time. Thus, seeds only germinate in the uninhibited part of [0, L]. We want to estimate lambda on the basis of one or more realizations of the process, the data being the locations and germination times of the germinated seeds. We derive the maximum likelihood estimator of v and a nonparametric estimator of lambda and describe methods of obtaining parametric estimates from it, illustrating these with reference to gamma densities. Simulation results are described and the methods applied to some neurobiological data. An Appendix outlines the S-PLUS code used.     

A continuous time version and a generalization of a Markov-recapture model for trapping experiments

Wileyto et al. [E.P. Wileyto, W.J. Ewens, M.A. Mullen, Markov-recapture population estimates: a tool for improving interpretation of trapping experiments, Ecology 75 (1994) 1109] propose a four-state discrete time Markov process, which describes the structure of a marking-capture experiment as a method of population estimation. They propose this method primarily for estimation of closed insect populations. Their method provides a mark-recapture estimate from a single trap observation by allowing subjects to mark themselves. The estimate of the unknown population size is based on the assumption of a closed population and a simple Markov model in which the rates of marking, capture, and recapture are assumed to be equal. Using the one step transition probability matrix of their model, we illustrate how to go from an embedded discrete time Markov process to a continuous time Markov process assuming exponentially distributed holding times. We also compute the transition probabilities after time t for the continuous time case and compare the limiting behavior of the continuous and discrete time processes. Finally, we generalize their model by relaxing the assumption of equal per capita rates for marking, capture, and recapture. Other questions about how their results change when using a continuous time Markov process are examined.