Parametric Bayesian inference for Y-linked two-sex branching models

A Y-linked two-sex branching process with blind choice is a suitable model for analyzing the evolution of the number of carriers of two alleles of a Y-linked gene in a two-sex monogamous population where each female chooses her partner from among the male population without caring about his type (i.e., the allele he carries). This work focuses on the development of Bayesian inference for this model, considering a parametric framework with the reproduction laws belonging to the power series family of distributions. A sample is considered given by the observation of the total number of females and males (regardless of their types) up to some generation as well as the number of each type of male in the last generation. Using a simulation method based on the Gibbs sampler, we approximate the posterior distributions of the main parameters of this model. The accuracy of the procedure based on this sample is illustrated by way of a simulated example.     

Estimation of the offspring mean in a branching process with non stationary immigration

ABSTRACT

In the paper, we consider a natural estimator of the offspring mean of a branching process with non stationary immigration based on observation of population sizes and number of immigrating individuals to each generation. We demonstrate that using a central limit theorem for multiple sums of dependent random variables it is possible to derive asymptotic distributions for the estimator without prior knowledge about the behavior (criticality) of the reproduction process. Before the three cases of criticality have been considered separately. Assuming that the immigration mean and variance vary regularly, conditions guaranteeing the strong consistency of the proposed estimator is also derived.     

LIMIT THEOREMS FOR SUBCRITICAL AGE-DEPENDENT BRANCHING PROCESSES WITH TWO TYPES OF IMMIGRATION

ABSTRACT

For the classical subcritical age-dependent branching process the effect of the following two-type immigration pattern is studied. At a sequence of renewal epochs a random number of immigrants enters the population. Each subpopulation stemming from one of these immigrants or one of the ancestors is revived by new immigrants and their offspring whenever it dies out, possibly after an additional delay period. All individuals have the same lifetime distribution and produce offspring according to the same reproduction law. This is the Bellman-Harris process with immigration at zero and immigration of renewal type (BHPIOR). We prove a strong law of large numbers and a central limit theorem for such processes. Similar conclusions are obtained for their discrete-time counterparts (lifetime per individual equals one), called Galton-Watson processes with immigration at zero and immigration of renewal type (GWPIOR). Our approach is based on the theory of regenerative processes, renewal theory and occupation measures and is quite different from those in earlier related work using analytic tools.     

The estimation of adopted mortality and morbidity rates using model and the phase type law: the Turkish case

Abstract

This paper aims to estimate mortality rate, morbidity-mortality rates of a chronic disease utilizing phase type law in the frame of two and three state processes. The application on commonly used mortality tables in Turkey are adopted to process to estimate the future mortalities with respect to phase type distribution for the purpose of justifying. Using one absorbing state, two and three state Models calculate the time until absorbing of the death and death by phase type distribution for each gender. Consequently, the 3-state probabilities in estimating the mortality-morbidity rates of IHD for Turkish population yield a significant information on the health management and pricing health insurance products.     

Model fitting and inference under latent equilibrium processes

This paper presents a methodology for model fitting and inference in the context of Bayesian models of the type f(Y | X,θ)f(X|θ)f(θ), where Y is the (set of) observed data, θ is a set of model parameters and X is an unobserved (latent) stationary stochastic process induced by the first order transition model f(X ^(t+1)|X ^(t),θ), where X ^(t) denotes the state of the process at time (or generation) t. The crucial feature of the above type of model is that, given θ, the transition model f(X ^(t+1)|X ^(t),θ) is known but the distribution of the stochastic process in equilibrium, that is f(X|θ), is, except in very special cases, intractable, hence unknown. A further point to note is that the data Y has been assumed to be observed when the underlying process is in equilibrium. In other words, the data is not collected dynamically over time. We refer to such specification as a latent equilibrium process (LEP) model. It is motivated by problems in population genetics (though other applications are discussed), where it is of interest to learn about parameters such as mutation and migration rates and population sizes, given a sample of allele frequencies at one or more loci. In such problems it is natural to assume that the distribution of the observed allele frequencies depends on the true (unobserved) population allele frequencies, whereas the distribution of the true allele frequencies is only indirectly specified through a transition model. As a hierarchical specification, it is natural to fit the LEP within a Bayesian framework. Fitting such models is usually done via Markov chain Monte Carlo (MCMC). However, we demonstrate that, in the case of LEP models, implementation of MCMC is far from straightforward. The main contribution of this paper is to provide a methodology to implement MCMC for LEP models. We demonstrate our approach in population genetics problems with both simulated and real data sets. The resultant model fitting is computationally intensive and thus, we also discuss parallel implementation of the procedure in special cases.     

The effect of autocorrelation on the ewma maxmin tolerance limits

Amin et al. (1999) developed an exponentially weighted moving average (EWMA) control chart, based on the smallest and largest observations in each sample. The resulting plot of the extremes suggests that the MaxMin EWMA may also be viewed as smoothed tolerance limits. Tolerance limits are limits that include a specific proportion of the population at a given confidence level. In the context of process control, they are used to make sure that production will not be outside specifications. Amin and Li (2000) provided the coverages of the MaxMin EWMA tolerance limits for independent data. In this article, it is shown how autocorrelation affects the confidence level of MaxMin tolerance limits, for a specified level of coverage of the population, and modified smoothed tolerance limits are suggested for autocorrelated processes.     

Bayesian analysis of Cannabis offences using generalized Poisson geometric process model with flexible dispersion

ABSTRACT

This paper derives models to analyse Cannabis offences count series from New South Wales, Australia. The data display substantial overdispersion as well as underdispersion for a subset, trend movement and population heterogeneity. To describe the trend dynamic in the data, the Poisson geometric process model is first adopted and is extended to the generalized Poisson geometric process model to capture both over- and underdispersion. By further incorporating mixture effect, the model accommodates population heterogeneity and enables classification of homogeneous units. The model is implemented using Markov chain Monte Carlo algorithms via the user-friendly WinBUGS software and its performance is evaluated through a simulation study.     

Bayesian Estimation of the Size of a Closed Population Using Photo-ID Data with Part of the Population Uncatchable

Abstract

We develop a Bayesian statistical model for estimating bowhead whale population size from photo-identification data when most of the population is uncatchable. The proposed conditional likelihood function is a product of Darroch's model, formulated as a function of the number of good photos, and a binomial distribution of captured whales given the total number of good photos at each occasion. The full Bayesian model is implemented via adaptive rejection sampling for log concave densities. We apply the model to data from 1985 and 1986 bowhead whale photographic studies and the results compare favorably with the ones obtained in the literature. Also, a comparison with the maximum likelihood procedure with bootstrap simulation is considered using different vague priors for the capture probabilities.     

POISSON REGRESSION WITH A PERIODIC FUNCTION

ABSTRACT

Let {y_t } be a Poisson-like process with the mean μ _t which is a periodic function of time t. We discuss how to fit this type of data set using quasi-likelihood method. Our method provides a new avenue to fit a time series data when the usual assumption of stationarity and homogeneous residual variances are invalid. We show that the estimators obtained are strongly consistent and also asymptotically normal.     

Generalized run tests for statistical process control

ABSTRACT

In a sequence of elements, a run is defined as a maximal subsequence of like elements. The number of runs or the length of the longest run has been widely used to test the randomness of an ordered sequence. Based on two different sampling methods and two types of test statistics used, run tests can be classified into one of four cases. Numerous researchers have derived the probability distributions in many different ways, treating each case separately. In the paper, we propose a unified approach which is based on recurrence arguments of two mutually exclusive sub-sequences. We also consider the sequence of nominal data that has more than two classes. Thus, the traditional run tests for a binary sequence are special cases of our generalized run tests. We finally show that the generalized run tests can be applied to many quality management areas, such as testing changes in process variation, developing non-parametric multivariate control charts, and comparing the shapes and locations of more than two process distributions.     

A Resampling Approach for Under-estimating a Finite Population Total from a Censored Sample

Abstract

Suppose a finite population of N objects each of which has an unknown value μ _i  ≥ 0, i = 1, … , N of a nonnegative characteristic of interest. A random sample has been drawn, but only for a selected subset of the sample the μ-values have been observed. The subset selection procedure has been somewhat obscure, and thus the subsample is censorized rather than random. Despite that, a reliable lower bound for the population total (the sum of all μ _i ) is required which uses the statistical information contained in the data. We propose a resampling procedure to construct an under-estimate of the population total. We also consider the case when the objects of the population have unequal sampling probabilities, in particular when the population is divided into a few number of strata with constant probabilities within each stratum. A real data example illustrates the method.     

Model Order Estimation of a Multivariable Stochastic Process

ABSTRACT

This article presents a procedure allowing us to estimate the minimal order of a state-space representation, for a multivariable stochastic process, from a sequence of observations. The method proposes a statistical rule for testing the rank of a block Hankel matrix of data, since this rank is related to the order of the process. A new information criterion is then developed and used to decide upon the order of the model. In this article we generalize the Aoki C-test. Using two representative data sets as the basis for a Monte Carlo experiment and real data based on Danish economy, we estimate the order of multivariable stochastic processes.     

Ultimate Extinction of the Promiscuous Bisexual Galton—Watson Metapopulation

A variant of a sexual Gallon–Watson process is considered. At each generation the population is partitioned among n‘hosts’ (population patches) and individual members mate at random only with others within the same host. This is appropriate for many macroparasite systems, and at low parasite loads it gives rise to a depressed rate of reproduction relative to an asexual system, due to the possibility that females are unmated. It is shown that stochasticity mitigates against this effect, so that for small initial populations the probability of ultimate extinction (the complement of an ‘epidemic’) displays a tradeoff as a function of n between the strength of fluctuations which overcome this ‘mating’ probability, and the probability of the subpopulation in one host being ‘rescued’ by that in another. Complementary approximations are developed for the extinction probability: an asymptotically exact approximation at large n, and for small n a short‐time probability that is exact in the limit where the mean number of offspring per parent is large.     

Nonparametric Growth Curve Model with Local Linear Approximation

The growth curve model has been developed for longitudinal data, and its time trend is usually described by polynomials. However, it is difficult to interpret each coefficient of the polynomials with higher degrees, even when the number of repetitions is sufficiently large. We propose herein an alternative growth curve model having time-varying coefficients.     

The pitfall of instrumental variables in big data: What the rule of thumb can't give you

ABSTRACT

Background: Instrumental variables (IVs) have become much easier to find in the “Big data era” which has increased the number of applications of the Two-Stage Least Squares model (TSLS). With the increased availability of IVs, the possibility that these IVs are weak has increased. Prior work has suggested a ‘rule of thumb’ that IVs with a first stage F statistic at least ten will avoid a relative bias in point estimates greater than 10%. We investigated whether or not this threshold was also an efficient guarantee of low false rejection rates of the null hypothesis test in TSLS applications with many IVs.

Objective: To test how the ‘rule of thumb’ for weak instruments performs in predicting low false rejection rates in the TSLS model when the number of IVs is large.

Method: We used a Monte Carlo approach to create 28 original data sets for different models with the number of IVs varying from 3 to 30. For each model, we generated 2000 observations for each iteration and conducted 50,000 iterations to reach convergence in rejection rates. The point estimate was set to 0, and probabilities of rejecting this hypothesis were recorded for each model as a measurement of false rejection rate. The relationship between the endogenous variable and IVs was carefully adjusted to let the F statistics for the first stage model equal ten, thus simulating the ‘rule of thumb.’

Results: We found that the false rejection rates (type I errors) increased when the number of IVs in the TSLS model increased while holding the F statistics for the first stage model equal to 10. The false rejection rate exceeds 10% when TLSL has 24 IVs and exceed 15% when TLSL has 30 IVs.

Conclusion: When more instrumental variables were applied in the model, the ‘rule of thumb’ was no longer an efficient guarantee for good performance in hypothesis testing. A more restricted margin for F statistics is recommended to replace the ‘rule of thumb,’ especially when the number of instrumental variables is large.     

The discretely observed immigration-death process: Likelihood inference and spatiotemporal applications

ABSTRACT

We consider a stochastic process, the homogeneous spatial immigration-death (HSID) process, which is a spatial birth-death process with as building blocks (i) an immigration-death (ID) process (a continuous-time Markov chain) and (ii) a probability distribution assigning iid spatial locations to all events. For the ID process, we derive the likelihood function, reduce the likelihood estimation problem to one dimension, and prove consistency and asymptotic normality for the maximum likelihood estimators (MLEs) under a discrete sampling scheme. We additionally prove consistency for the MLEs of HSID processes. In connection to the growth-interaction process, which has a HSID process as basis, we also fit HSID processes to Scots pine data.     

A semi-Markov decision model for the optimal control of a simple immigration-birth-death process through the introduction of a predator

ABSTRACT

This article is concerned with the problem of controlling a simple immigration-birth-death process, which represents a pest population, by the introduction of a predator in the habitat of the pests. The optimization criterion is the minimization of the expected long-run average cost per unit time. It is possible to construct an appropriate semi-Markov decision model with a finite set of states if and only if the difference between the per capita birth rate and the per capita death rate of the pests is smaller than half of the rate at which the predator is introduced in the habitat.     

THE BEHAVIOR OF HEGY TESTS FOR QUARTERLY TIME SERIES WITH SEASONAL MEAN SHIFTS

ABSTRACT

This paper studies the behavior of the HEGY statistics for quarterly data, for seasonal autoregressive unit roots, when the analyzed time series is deterministic seasonal stationary but exhibits a change in the seasonal pattern. We analyze also the HEGY test for the nonseasonal unit root. the data generation process being trend stationary too. Our results show that when the break magnitudes are finite, the HEGY test statistics are not asymptotically biased toward the nonrejection of the seasonal and nonseasonal unit root hypotheses. However, the finite sample power properties may be substantially affected, the behavior of the tests depending on the type of the break.