Multivariate Polya and Inverse Polya Distributions of Order k

Multivariate Polya and inverse Polya distributions of order k are derived by means of generalized urn models and by compounding the type II multinomial and multivariate negative binomial distributions of order k of PHILIPPOU , ANTZOULAKOS and TRIPSIANNIS (1990, 1988), respectively, with the Dirichlet distribution. It is noted that the above two distributions include as special cases a multivariate hypergeometric distribution of order k, a negative one, an inverse one, a negative inverse one and a discrete uniform of the same order. The probability generating functions, means, variances and covariances of the new distributions are obtained and five asymptotic results are established relating them to the above-mentioned multinomial and multivariate negative binomial distributions of order k, and to the type II negative binomial and the type I multivariate Poisson distributions of order k of PHILIPPOU (1983), and PHILIPPOU , ANTZOULAKOS and TRIPSIAN-NIS (1988), respectively. Potential applications are also indicated. The present paper extends to the multivariate case the work of PHILIPPOU , TRIPSIANNIS and ANTZOULAKOS (1989) on Polya and inverse Polya distributions of order k..     

On Bivariate Cumulative Damage Models With Application

Bivariate cumulative damage models are proposed where the responses given the damages are independent random variables. The bivariate damage process can be either bivariate Poisson or bivariate gamma. A bivariate continuous cumulative damage model is investigated in which the responses given the damages have gamma distributions. In this case evaluation of the joint density function and bivariate tail probability function is facilitated by expanding the gamma distributions of the conditional responses by Laguerre polynomials. This approach also leads to evaluation of associated survival models. Moments and estimating equations are discussed. In addition, a bivariate discrete cumulative damage model is investigated in which the responses given the damages have a distribution chosen from a class that includes the negative binomial, the Neyman Type‐A, the Polya‐Aeppli, and the Lagrangian Poisson. Probabilities are obtained from recursive formulas which do not involve cancellation error as all quantities are non‐negative. Moments and estimating equations are presented for these models also. The continuous and the discrete models are applied to describe the rise of systolic and diastolic blood pressure with age.     

Some Applications of a Bivariate BOREL-TANNER Distribution

This paper considers the probability distribution of the number of customers served in the busy period of a single server queue with POISSON input and constant service time. Applications to traffic flow, semi-infinite discrete dam and branching processes are described. It is shown that in many situations the bivariate BOREL-TANNER distribution provides a better fit than the bivariate negative binomial distribution.     

Quantitative trait linkage analysis using Gaussian copulas

Mapping and identifying variants that influence quantitative traits is an important problem for genetic studies. Traditional QTL mapping relies on a variance-components (VC) approach with the key assumption that the trait values in a family follow a multivariate normal distribution. Violation of this assumption can lead to inflated type I error, reduced power, and biased parameter estimates. To accommodate nonnormally distributed data, we developed and implemented a modified VC method, which we call the "copula VC method," that directly models the nonnormal distribution using Gaussian copulas. The copula VC method allows the analysis of continuous, discrete, and censored trait data, and the standard VC method is a special case when the data are distributed as multivariate normal. Through the use of link functions, the copula VC method can easily incorporate covariates. We use computer simulations to show that the proposed method yields unbiased parameter estimates, correct type I error rates, and improved power for testing linkage with a variety of nonnormal traits as compared with the standard VC and the regression-based methods.     

COMPARISON OF METHODS FOR ANALYZING REPLICATED PREFERENCE TESTS

Preference testing is commonly used in consumer sensory evaluation. Traditionally, it is done without replication, effectively leading to a single 0/1 (binary) measurement on each panelist. However, to understand the nature of the preference, replicated preference tests are a better approach, resulting in binomial counts of preferences on each panelist. Variability among panelists then leads to overdispersion of the counts when the binomial model is used and to an inflated Type I error rate for statistical tests of preference. Overdispersion can be adjusted by Pearson correction or by other models such as correlated binomial or beta‐binomial. Several methods are suggested or reviewed in this study for analyzing replicated preference tests and their Type I error rates and power are compared. Simulation studies show that all methods have reasonable Type I error rates and similar power. Among them, the binomial model with Pearson adjustment is probably the safest way to analyze replicated preference tests, while a normal model in which the binomial distribution is not assumed is the easiest.     

Random effects modeling of multiple binomial responses using the multivariate binomial logit-normal distribution

The multivariate binomial logit-normal distribution is a mixture distribution for which, (i) conditional on a set of success probabilities and sample size indices, a vector of counts is independent binomial variates, and (ii) the vector of logits of the parameters has a multivariate normal distribution. We use this distribution to model multivariate binomial-type responses using a vector of random effects. The vector of logits of parameters has a mean that is a linear function of explanatory variables and has an unspecified or partly specified covariance matrix. The model generalizes and provides greater flexibility than the univariate model that uses a normal random effect to account for positive correlations in clustered data. The multivariate model is useful when different elements of the response vector refer to different characteristics, each of which may naturally have its own random effect. It is also useful for repeated binary measurement of a single response when there is a nonexchangeable association structure, such as one often expects with longitudinal data or when negative association exists for at least one pair of responses. We apply the model to an influenza study with repeated responses in which some pairs are negatively associated and to a developmental toxicity study with continuation-ratio logits applied to an ordinal response with clustered observations.     

Zero‐Inflated Negative Binomial Mixed Regression Modeling of Over‐Dispersed Count Data with Extra Zeros

In many biometrical applications, the count data encountered often contain extra zeros relative to the Poisson distribution. Zero‐inflated Poisson regression models are useful for analyzing such data, but parameter estimates may be seriously biased if the nonzero observations are over‐dispersed and simultaneously correlated due to the sampling design or the data collection procedure. In this paper, a zero‐inflated negative binomial mixed regression model is presented to analyze a set of pancreas disorder length of stay (LOS) data that comprised mainly same‐day separations. Random effects are introduced to account for inter‐hospital variations and the dependency of clustered LOS observations. Parameter estimation is achieved by maximizing an appropriate log‐likelihood function using an EM algorithm. Alternative modeling strategies, namely the finite mixture of Poisson distributions and the non‐parametric maximum likelihood approach, are also considered. The determination of pertinent covariates would assist hospital administrators and clinicians to manage LOS and expenditures efficiently.     

On Countable Mixtures of Bivariate Binomial Distributions

A unified treatment is given for mixtures of bivariate binomial distributions with respect to their index parameter(s). The use of probability generating functions is employed and a number of interesting properties including probabilities, factorial moments, factorial cumulants and conditional distributions are derived. Five classes of such mixtures are examined and several well known bivariate discrete distributions are used as illustrative examples. Biological applications are indicated including the fit of three bivariate distributions to an actual set of human family data.     

Multinomial N-mixture models for removal sampling

Multinomial $N$-mixture models are commonly used to fit data from a removal sampling protocol. If the mixing distribution is negative binomial, the distribution of the counts does not appear to have been identified, and practitioners approximate the requisite likelihood by placing an upper bound on the embedded infinite sum. In this paper, the distribution which underpins the multinomial $N$-mixture model with a negative binomial mixing distribution is shown to belong to the broad class of multivariate negative binomial distributions. Specifically, the likelihood can be expressed in closed form as the product of conditional and marginal likelihoods and the information matrix shown to be block diagonal. As a consequence, the nature of the maximum likelihood estimates of the unknown parameters and their attendant standard errors can be examined and tests of the hypothesis of the Poisson against the negative binomial mixing distribution formulated. In addition, appropriate multinomial $N$-mixture models for data sets which include zero site totals can also be constructed. Two illustrative examples are provided.     

One‐inflation and unobserved heterogeneity in population size estimation

We present the one‐inflated zero‐truncated negative binomial (OIZTNB) model, and propose its use as the truncated count distribution in Horvitz–Thompson estimation of an unknown population size. In the presence of unobserved heterogeneity, the zero‐truncated negative binomial (ZTNB) model is a natural choice over the positive Poisson (PP) model; however, when one‐inflation is present the ZTNB model either suffers from a boundary problem, or provides extremely biased population size estimates. Monte Carlo evidence suggests that in the presence of one‐inflation, the Horvitz–Thompson estimator under the ZTNB model can converge in probability to infinity. The OIZTNB model gives markedly different population size estimates compared to some existing truncated count distributions, when applied to several capture–recapture data that exhibit both one‐inflation and unobserved heterogeneity.     

A score test for testing a zero-inflated Poisson regression model against zero-inflated negative binomial alternatives

Count data often show a higher incidence of zero counts than would be expected if the data were Poisson distributed. Zero-inflated Poisson regression models are a useful class of models for such data, but parameter estimates may be seriously biased if the nonzero counts are overdispersed in relation to the Poisson distribution. We therefore provide a score test for testing zero-inflated Poisson regression models against zero-inflated negative binomial alternatives.     

Genetic analysis of life-history constraint and evolution in a wild ungulate population

Trade-offs among life-history traits are central to evolutionary theory. In quantitative genetic terms, trade-offs may be manifested as negative genetic covariances relative to the direction of selection on phenotypic traits. Although the expression and selection of ecologically important phenotypic variation are fundamentally multivariate phenomena, the in situ quantification of genetic covariances is challenging. Even for life-history traits, where well-developed theory exists with which to relate phenotypic variation to fitness variation, little evidence exists from in situ studies that negative genetic covariances are an important aspect of the genetic architecture of life-history traits. In fact, the majority of reported estimates of genetic covariances among life-history traits are positive. Here we apply theory of the genetics and selection of life histories in organisms with complex life cycles to provide a framework for quantifying the contribution of multivariate genetically based relationships among traits to evolutionary constraint. We use a Bayesian framework to link pedigree-based inference of the genetic basis of variation in life-history traits to evolutionary demography theory regarding how life histories are selected. Our results suggest that genetic covariances may be acting to constrain the evolution of female life-history traits in a wild population of red deer Cervus elaphus: genetic covariances are estimated to reduce the rate of adaptation by about 40%, relative to predicted evolutionary change in the absence of genetic covariances. Furthermore, multivariate phenotypic (rather than genetic) relationships among female life-history traits do not reveal this constraint.     

Analyzing Short-Term Noise Dependencies of Spike-Counts in Macaque Prefrontal Cortex Using Copulas and the Flashlight Transformation

Simultaneous spike-counts of neural populations are typically modeled by a Gaussian distribution. On short time scales, however, this distribution is too restrictive to describe and analyze multivariate distributions of discrete spike-counts. We present an alternative that is based on copulas and can account for arbitrary marginal distributions, including Poisson and negative binomial distributions as well as second and higher-order interactions. We describe maximum likelihood-based procedures for fitting copula-based models to spike-count data, and we derive a so-called flashlight transformation which makes it possible to move the tail dependence of an arbitrary copula into an arbitrary orthant of the multivariate probability distribution. Mixtures of copulas that combine different dependence structures and thereby model different driving processes simultaneously are also introduced. First, we apply copula-based models to populations of integrate-and-fire neurons receiving partially correlated input and show that the best fitting copulas provide information about the functional connectivity of coupled neurons which can be extracted using the flashlight transformation. We then apply the new method to data which were recorded from macaque prefrontal cortex using a multi-tetrode array. We find that copula-based distributions with negative binomial marginals provide an appropriate stochastic model for the multivariate spike-count distributions rather than the multivariate Poisson latent variables distribution and the often used multivariate normal distribution. The dependence structure of these distributions provides evidence for common inhibitory input to all recorded stimulus encoding neurons. Finally, we show that copula-based models can be successfully used to evaluate neural codes, e.g., to characterize stimulus-dependent spike-count distributions with information measures. This demonstrates that copula-based models are not only a versatile class of models for multivariate distributions of spike-counts, but that those models can be exploited to understand functional dependencies.     

A quantile count model of water depth constraints on Cape Sable seaside sparrows

1. A quantile regression model for counts of breeding Cape Sable seaside sparrows Ammodramus maritimus mirabilis (L.) as a function of water depth and previous year abundance was developed based on extensive surveys, 1992-2005, in the Florida Everglades. The quantile count model extends linear quantile regression methods to discrete response variables, providing a flexible alternative to discrete parametric distributional models, e.g. Poisson, negative binomial and their zero-inflated counterparts. 2. Estimates from our multiplicative model demonstrated that negative effects of increasing water depth in breeding habitat on sparrow numbers were dependent on recent occupation history. Upper 10th percentiles of counts (one to three sparrows) decreased with increasing water depth from 0 to 30 cm when sites were not occupied in previous years. However, upper 40th percentiles of counts (one to six sparrows) decreased with increasing water depth for sites occupied in previous years. 3. Greatest decreases (-50% to -83%) in upper quantiles of sparrow counts occurred as water depths increased from 0 to 15 cm when previous year counts were 1, but a small proportion of sites (5-10%) held at least one sparrow even as water depths increased to 20 or 30 cm. 4. A zero-inflated Poisson regression model provided estimates of conditional means that also decreased with increasing water depth but rates of change were lower and decreased with increasing previous year counts compared to the quantile count model. Quantiles computed for the zero-inflated Poisson model enhanced interpretation of this model but had greater lack-of-fit for water depths > 0 cm and previous year counts 1, conditions where the negative effect of water depths were readily apparent and fitted better with the quantile count model.     

Bacterial Density in Water Determined by Poisson or Negative Binomial Distributions

The question of how to characterize the bacterial density in a body of water when data are available as counts from a number of small-volume samples was examined for cases where either the Poisson or negative binomial probability distributions could be used to describe the bacteriological data. The suitability of the Poisson distribution when replicate analyses were performed under carefully controlled conditions and of the negative binomial distribution for samples collected from different locations and over time were illustrated by two examples. In cases where the negative binomial distribution was appropriate, a procedure was given for characterizing the variability by dividing the bacterial counts into homogeneous groups. The usefulness of this procedure was illustrated for the second example based on survey data for Lake Erie. A further illustration of the difference between results based on the Poisson and negative binomial distributions was given by calculating the probability of obtaining all samples sterile, assuming various bacterial densities and sample sizes.     

Smoothing of the bivariate LOD score for non-normal quantitative traits

Variance component analysis provides an efficient method for performing linkage analysis for quantitative traits. However, type I error of variance components-based likelihood ratio testing may be affected when phenotypic data are non-normally distributed (especially with high values of kurtosis). This results in inflated LOD scores when the normality assumption does not hold. Even though different solutions have been proposed to deal with this problem with univariate phenotypes, little work has been done in the multivariate case. We present an empirical approach to adjust the inflated LOD scores obtained from a bivariate phenotype that violates the assumption of normality. Using the Collaborative Study on the Genetics of Alcoholism data available for the Genetic Analysis Workshop 14, we show how bivariate linkage analysis with leptokurtotic traits gives an inflated type I error. We perform a novel correction that achieves acceptable levels of type I error.     

An estimator of the mutant frequency in assays using transgenic animals

The Poisson distribution is a fundamental probability model for count data, and is a natural model for the observed plaque counts in mutation assays using animals with lambda or PhiX174 transgenes. The Poisson likelihood for observed counts is a function of the mutant fraction, and it is straightforward to derive the associated maximum likelihood estimate of the mutant fraction and its variance. The estimate is easy to calculate, and if not the same, very similar to ad hoc estimates in current use. The model indicates the proper way to combine data from a number of plates, possibly prepared with different sample dilutions. The estimator of the mutant fraction is biased as a consequence of dividing by a random variable, the plaque count used to calculate the total recovered plaque-forming units. Fortunately, the bias becomes negligible as this count becomes large. On the other hand, increasing this count can increase the variance by decreasing the amount of sample assayed for mutant phages. Concurrent heed to the bias and the variance provides some guidance as to the optimum allocation of a sample into portions assayed for mutant phages and total recovered phages. The distribution of the estimate of the mutant fraction is related to the binomial distribution. This relationship implies a binomial distribution for the mutant count conditional on an overall count (either the sum of mutant and counted total plaques or the sum of counted mutant and non-mutant plaques). A special but important case occurs when each plate can be evaluated for mutant plaques and non-mutant plaques. Then, the observed proportion of mutants estimates the mutant fraction. More generally, the relationship to a binomial distribution provides a procedure for calculating a confidence interval.     

Overestimation of Yield Loss of Tobacco Caused by the Aggregated Spatial Pattern of Meloidogyne incognita

Overestimation of yield loss caused by Meloidogyne incognita on tobacco was calculated as a function of the statistical frequency distribution of sample counts. Sampling frequency distributions were described by a negative binomial model, with parameter k, and the resulting probability generating function was used to calculate discrete damage probabilities. Negative binomial damage predictions were compared to mean-density estimates of damage. Predictions based on mean density alone overestimate yield loss by values ranging from 300% at a k of 0.1 to less than 10% at a k of 1.0. Damage overestimation was described as an exponential function of k and mean density. Preplant sampling data for M. incognita were used to derive a linear model for the estimation of k from mean density, allowing the calculation of yield-loss overestimation based on one parameter, the field mean density. Overestimation of damage ranged from 288% at a density of 50 juveniles/500 cm³ soil, to 5% at a density of 1,000 juvelfiles/500 cm³ soil.     

One‐inflation and unobserved heterogeneity in population size estimation by Ryan T. Godwin

In this study, we would like to show that the one‐inflated zero‐truncated negative binomial (OIZTNB) regression model can be easily implemented in R via built‐in functions when we use mean‐parameterization feature of negative binomial distribution to build OIZTNB regression model. From the practitioners' point of view, we believe that this approach presents a computationally convenient way for implementation of the OIZTNB regression model.     

Calculation of narrower confidence intervals for tree mortality rates when we know nothing but the location of the death/survival events

相似文献