The non-central negative binomial distribution: Further properties and applications

Abstract

The non-central negative binomial distribution is both a mixed and compound Poisson distribution with applications in photon and neural counting, statistical optics, astronomy and a stochastic reversible counter system. In this paper various important probabilistic properties of the non-central negative binomial distribution in practical applications like log-concavity, discrete self-decomposability, unimodality, asymptotic behavior and tail length of the probability distribution have been derived. The construction as a mixed Poisson process by specifying a joint distribution for the inter-arrival times and its application is illustrated by a fit to real life data set.     

The finite-time ruin probability of a discrete-time risk model with GARCH discounted factors and dependent risks

The finite-time ruin probability of a discrete-time risk model with dependent stochastic discount factors and dependent insurance and financial risks is investigated in this paper. Assume that the stochastic discount factors follow a GARCH process and the one-period insurance and financial risks form a sequence of independent and identically distributed random pairs, which are the copies of a random pair with a bivariate Sarmanov dependent distribution. When the common distribution of claim-sizes is heavy-tailed, we establish an asymptotic estimate for the finite-time ruin probability. Applying the result to a special case, we also get conservative asymptotic bounds. A numerical simulation is given at the end of the paper.     

Stochastic differential equation systems for an SIS epidemic model with vaccination and immigration

Two Itô stochastic differential equation (SDE) systems are constructed for a Susceptible-Infected-Susceptible epidemic model with temporary vaccination. A constant number of new members enter the population and total size of the population is variable. Some conditions for disease extinction in the stochastic models are established and compared with conditions in deterministic one. It is shown that the two stochastic models are equivalent in the sense that their solutions come from same distribution. In addition, the SDE models are simulated and the equivalence of the two stochastic models is confirmed by numerical examples. The probability distribution for extinction is also obtained numerically, provided there exists a probability for disease persistence whereas the expected duration of epidemic is acquired when extinction occurs with probability 1.     

Analysis of a two incongruent arrivals and services M/G/1 queue with random feedback

A queuing system with two incongruent arrivals and services is considered. Two kinds of customers enter the system by Poisson process and the service times are assumed to have general distribution. After first kind service completion, it may feedback to repeat the first service, leave the system or go to give second service. The same policy is applied for the other kind of customer. All stochastic processes involved in this system are independent. We derive the probability generating function for each kind and for the system that yield the performance measures. Some numerical approaches examined the validity of the results.     

Wavelets for Nonparametric Stochastic Regression with Mixing Stochastic Process

We propose a wavelet based stochastic regression function estimator for the estimation of the regression function for a sequence of mixing stochastic process with a common one-dimensional probability density function. Some asymptotic properties of the proposed estimator are investigated. It is found that the estimators have similar properties to their counterparts studied earlier in literature.     

Bayesian Poisson models for the graphical combination of dependent expert information

A supra-Bayesian (SB) wants to combine the information from a group of k experts to produce her distribution of a probability θ. Each expert gives his counts of what he thinks are the numbers of successes and failures in a sequence of independent trials, each with probability θ of success. These counts, used as a surrogate for each expert's own individual probability assessment (together with his associated level of confidence in his estimate), allow the SB to build various plausible conjugate models. Such models reflect her beliefs about the reliability of different experts and take account of different possible patterns of overlap of information between them. Corresponding combination rules are then obtained and compared with other more established rules and their properties examined.     

On probabilistic parametric inference

This paper formulates a theory of probabilistic parametric inference and explores the limits of its applicability. Unlike Bayesian statistical models, the system does not comprise prior probability distributions. Objectivity is imposed on the theory: a particular direct probability density should always result in the same posterior probability distribution. For calibrated posterior probability distributions it is possible to construct credible regions with posterior-probability content equal to the coverage of the regions, but the calibration is not generally preserved under marginalization. As an application of the theory, the paper also constructs a filter for linear Gauss–Markov stochastic processes with unspecified initial conditions.     

Convolution of Binomial and Negative Binomial Variables

This paper considers a distribution formed by convolution of binomial and negative binomial variables. The distribution has the flexibility to adapt to the model under, equi, and over dispersion. Some properties of the proposed distribution are discussed, including characterization. Three stochastic processes leading to the distribution are also considered: (1) a three-dimensional random walk; (2) a birth, death, and immigration process; and (3) a thinned stochastic process.     

A stochastic graph process for epidemic modelling

A stochastic graph process with a Markov property is introduced to model the flow of an infectious disease over a known contact network. The model provides a probability distribution over unobserved infectious pathways. The basic reproductive number in compartmental models is generalized to a dynamic reproductive number based on the sequence of outdegrees in the graph process. The cumulative resistance and threat associated with each individual is also measured based on the cumulative indegree and outdegree of the graph process. The model is applied to the outbreak data from the 2001 foot‐and‐mouth (FMD) outbreak in the United Kingdom. The Canadian Journal of Statistics 40: 55–67; 2012 © 2012 Statistical Society of Canada     

Simple bounds for terminating Poisson and renewal shock processes

A system subject to a point process of shocks is considered. The shocks occur in accordance with a renewal process or a nonhomogeneous Poisson process. Each shock independently of the previous history leads to a system failure with probability θ and is survived with a complimentary probability $\text{θ}\text{̄}$. A number of problems in reliability and safety analysis can be interpreted by means of this model. The exact solution for the probability of survival $\text{W}\text{̄}\text{(t,θ)}$ can be obtained only in the form of infinite series (renewal process of shocks). Approximate solutions and new simple bounds for the probability of survival are obtained. The introduced method is based on the notion of a stochastic hazard rate process. A supplementary characteristic in this analysis is the mean of the hazard rate process. This method makes it possible to consider a generalization important in practical applications when the probability of a system failure under the effect of a current shock depends on the time since the previous one.     

Fast Inference for Network Models of Infectious Disease Spread

Models of infectious disease over contact networks offer a versatile means of capturing heterogeneity in populations during an epidemic. Highly connected individuals tend to be infected at a higher rate early during an outbreak than those with fewer connections. A powerful approach based on the probability generating function of the individual degree distribution exists for modelling the mean field dynamics of outbreaks in such a population. We develop the same idea in a stochastic context, by proposing a comprehensive model for 1‐week‐ahead incidence counts. Our focus is inferring contact network (and other epidemic) parameters for some common degree distributions, in the case when the network is non‐homogeneous ‘at random’. Our model is initially set within a susceptible–infectious–removed framework, then extended to the susceptible–infectious–removed–susceptible scenario, and we apply this methodology to influenza A data.     

Model-based clustering of multivariate ordinal data relying on a stochastic binary search algorithm

We design a probability distribution for ordinal data by modeling the process generating data, which is assumed to rely only on order comparisons between categories. Contrariwise, most competitors often either forget the order information or add a non-existent distance information. The data generating process is assumed, from optimality arguments, to be a stochastic binary search algorithm in a sorted table. The resulting distribution is natively governed by two meaningful parameters (position and precision) and has very appealing properties: decrease around the mode, shape tuning from uniformity to a Dirac, identifiability. Moreover, it is easily estimated by an EM algorithm since the path in the stochastic binary search algorithm can be considered as missing values. Using then the classical latent class assumption, the previous univariate ordinal model is straightforwardly extended to model-based clustering for multivariate ordinal data. Parameters of this mixture model are estimated by an AECM algorithm. Both simulated and real data sets illustrate the great potential of this model by its ability to parsimoniously identify particularly relevant clusters which were unsuspected by some traditional competitors.     

Process monitoring of exponentially distributed characteristics through an optimal normalizing transformation

Many process characteristics follow an exponential distribution, and control charts based on such a distribution have attracted a lot of attention. However, traditional control limits may be not appropriate because of the lack of symmetry. In this paper, process monitoring through a normalizing power transformation is studied. The traditional individual measurement control charts can be used based on the transformed data. The properties of this control chart are investigated. A comparison with the chart when using probability limits is also carried out for cases of known and estimated parameters. Without losing much accuracy, even compared with the exact probability limits, the power transformation approach can easily be used to produce charts that can be interpreted when the normality assumption is valid.     

A stochastic process approach to multilayer neutron detectors

The sparsity of the isotope Helium‐3, ongoing since 2009, has initiated a new generation of neutron detectors. One particularly promising development line for detectors is the multilayer gaseous detector. In this paper, a stochastic process approach is used to determine the neutron energy from the additional data afforded by the multilayer nature of these novel detectors. The data from a multilayer detector consist of counts of the number of absorbed neutrons along the sequence of the detector's layers, in which the neutron absorption probability is unknown. We study the maximum likelihood estimator for the intensity and absorption probability and show its consistency and asymptotic normality, as the number of incoming neutrons goes to infinity. We combine these results with known results on the relation between the absorption probability and the wavelength to derive an estimator of the wavelength and to show its consistency and asymptotic normality.     

First hitting probabilities for semi-Markov chains and estimation

We first consider a stochastic system described by an absorbing semi-Markov chain (SMC) with finite state space, and we introduce the absorption probability to a class of recurrent states. Afterwards, we study the first hitting probability to a subset of states for an irreducible SMC. In the latter case, a non-parametric estimator for the first hitting probability is proposed and the asymptotic properties of strong consistency and asymptotic normality are proven. Finally, a numerical application on a five-state system is presented to illustrate the performance of this estimator.     

An Interest-rate Model Analysis Based on Data Augmentation Bayesian Forecasting

In this paper, the author presents an efficient method of analyzing an interest-rate model using a new approach called 'data augmentation Bayesian forecasting.' First, a dynamic linear model estimation was constructed with a hierarchically-incorporated model. Next, an observational replication was generated based on the one-step forecast distribution derived from the model. A Markov-chain Monte Carlo sampling method was conducted on it as a new observation and unknown parameters were estimated. At that time, the EM algorithm was applied to establish initial values of unknown parameters while the 'quasi Bayes factor' was used to appreciate parameter candidates. 'Data augmentation Bayesian forecasting' is a method of evaluating the transition and history of 'future,' 'present' and 'past' of an arbitrary stochastic process by which an appropriate evaluation is conducted based on the probability measure that has been sequentially modified with additional information. It would be possible to use future prediction results for modifying the model to grasp the present state or re-evaluate the past state. It would be also possible to raise the degree of precision in predicting the future through the modification of the present and the past. Thus, 'data augmentation Bayesian forecasting' is applicable not only in the field of financial data analysis but also in forecasting and controlling the stochastic process.     

Double-blind deconvolution: the analysis of post-synaptic currents in nerve cells

This paper is concerned with the analysis of observations made on a system that is being stimulated at fixed time intervals but where the precise nature and effect of any individual stimulus is unknown. The realized values are modelled as a stochastic process consisting of a random signal embedded in noise. The aim of the analysis is to use the data to unravel the unknown structure of the system and to ascertain the probabilistic behaviour of the stimuli. A method of parameter estimation based on quasi-profile likelihood is presented and the statistical properties of the estimates are established while recognizing that there will be a discrepancy between the model and the true data-generating mechanism. A method of model validation and determination is also advanced and kernel smoothing techniques are proposed as a basis for identifying the amplitude distribution of the stimuli. The data processing techniques described have a direct application to the investigation of excitatory post-synaptic currents recorded from nerve cells in the central nervous system and their use in quantal analysis of such data is illustrated.     

The Distribution of Family Sizes Under a Time-Homogeneous Birth and Death Process

The number of extant individuals within a lineage, as exemplified by counts of species numbers across genera in a higher taxonomic category, is known to be a highly skewed distribution. Because the sublineages (such as genera in a clade) themselves follow a random birth process, deriving the distribution of lineage sizes involves averaging the solutions to a birth and death process over the distribution of time intervals separating the origin of the lineages. In this article, we show that the resulting distributions can be represented by hypergeometric functions of the second kind. We also provide approximations of these distributions up to the second order, and compare these results to the asymptotic distributions and numerical approximations used in previous studies. For two limiting cases, one with a relatively high rate of lineage origin, one with a low rate, the cumulative probability densities and percentiles are compared to show that the approximations are robust over a wide range of parameters. It is proposed that the probability distributions of lineage size may have a number of relevant applications to biological problems such as the coalescence of genetic lineages and in predicting the number of species in living and extinct higher taxa, as these systems are special instances of the underlying process analyzed in this article.     

Bayesian semiparametric modeling for stochastic precedence,with applications in epidemiology and survival analysis

We propose a prior probability model for two distributions that are ordered according to a stochastic precedence constraint, a weaker restriction than the more commonly utilized stochastic order constraint. The modeling approach is based on structured Dirichlet process mixtures of normal distributions. Full inference for functionals of the stochastic precedence constrained mixture distributions is obtained through a Markov chain Monte Carlo posterior simulation method. A motivating application involves study of the discriminatory ability of continuous diagnostic tests in epidemiologic research. Here, stochastic precedence provides a natural restriction for the distributions of test scores corresponding to the non-infected and infected groups. Inference under the model is illustrated with data from a diagnostic test for Johne’s disease in dairy cattle. We also apply the methodology to the comparison of survival distributions associated with two distinct conditions, and illustrate with analysis of data on survival time after bone marrow transplantation for treatment of leukemia.     

Matrix‐Exponential Distributions: Calculus and Interpretations via Flows

By considering randomly stopped deterministic flow models, we develop an intuitively appealing way to generate probability distributions with rational Laplace–Stieltjes transforms on [0,∞). That approach includes and generalizes the formalism of PH-distributions. That generalization results in the class of matrix-exponential probability distributions. To illustrate the novel way of thinking that is required to use these in stochastic models, we retrace the derivations of some results from matrix-exponential renewal theory and prove a new extension of a result from risk theory. Essentially the flow models allows for keeping track of the dynamics of a mechanism that generates matrix-exponential distributions in a similar way to the probabilistic arguments used for phase-type distributions involving transition rates. We also sketch a generalization of the Markovian arrival process (MAP) to the setting of matrix-exponential distribution. That process is known as the Rational arrival process (RAP).