Calculation Methods for Wallenius' Noncentral Hypergeometric Distribution

Two different probability distributions are both known in the literature as “the” noncentral hypergeometric distribution. Wallenius' noncentral hypergeometric distribution can be described by an urn model without replacement with bias. Fisher's noncentral hypergeometric distribution is the conditional distribution of independent binomial variates given their sum. No reliable calculation method for Wallenius' noncentral hypergeometric distribution has hitherto been described in the literature. Several new methods for calculating probabilities from Wallenius' noncentral hypergeometric distribution are derived. Range of applicability, numerical problems, and efficiency are discussed for each method. Approximations to the mean and variance are also discussed. This distribution has important applications in models of biased sampling and in models of evolutionary systems.     

An alternative representation of noncentral beta and F distributions

In this paper the doubly noncentral beta and F distributions are represented alternatively by using the results on the product of two hypergeometric functions. Their moments and the cumulative distribution functions are also given in terms of hypergeometric functions, which can be easily calculated by the Mathematica package.     

Characterization of hypergeometric type distributions by regression

The conditional distribution of Y given X is proved to be of the positive or the negative hypergeometric law when X follows a positive or a negative hypergeometric distribution, respectively and the regression function of X on Y is of a known linear form     

Large sample results for tests of association ii. multinomial and stratified sampling

This paper deals with the asymptotics of a class of tests for association in 2-way contingency tables based on square forms in cell frequencies, given the total number of observations (multinomial sampling) or one set of marginal totals (stratified sampling). The case when both row and column marginal totals are fixed (hypergeometric sampling) was studied in Kulinskaya (1994), The class of tests under consideration includes a number of classical measures for association, Its two subclasses are the tests based on statistics using centralized cell frequencies (asymptotically distributed as weighted sums of central chi-squares) and those using the non-centralized cell frequencies (asymptotically normal). The parameters of asymptotic distributions depend on the sampling model and on true marginal probabilities. Maximum efficiency for asymptotically normal statistics is achieved under hypergeometric sampling, If the cell frequencies or the statistic as a whole are centralized using marginal proportions as estimates for marginal probabilities, the asymptotic distribution does not differ much between models and it is equivalent to that under hypergeometric sampling. These findings give an extra justification for the use of permutation tests for association (which are based on hypergeometric sampling). As an application, several well known measures of association are analysed.     

On umvu estimators for the multivariate lognormal distribution and their variances

Uniformly minimum variance unbiased estimators of several parameters of the multivariate lognormal distribution are expressed by using the hypergeometric functions of matrix argument. And the variances are given in special cases.     

A Family Which Integrates the Generalized Charlier Series and Extended Noncentral Negative Binomial Distributions

The generalized Charlier series distribution includes the binomial distribution, and the noncentral negative binomial distribution extends the negative binomial distribution. The present article proposes a family of counting distributions, which contains both the generalized Charlier series and extended noncentral negative binomial distributions. Compound and mixture formulations of the proposed distribution are given. The probability mass function is expressible in terms of the confluent hypergeometric function as well as the Gauss hypergeometric function. Recursive formulae for probability mass function have been studied by Panjer, Sundt and Jewell, Schröter, Sundt, and Kitano et al. in the context of insurance risk. This article explores horizontal, vertical, triangular, and diagonal recursions. Recursive formulae as well as exact expressions for descending factorial moments are studied. The proposed distribution allows overdispersion or underdispersion relative to a Poisson distribution. An illustrative example of data fitting is given.     

On the Computation of Gauss Hypergeometric Functions

The pioneering study undertaken by Liang et al. in 2008 (Journal of the American Statistical Association, 103, 410–423) and the hundreds of papers citing that work make use of certain hypergeometric functions. Liang et al. and many others claim that the computation of the hypergeometric functions is difficult. Here, we show that the hypergeometric functions can in fact be reduced to simpler functions that can often be computed using a pocket calculator.     

A Phase I nonparametric Shewhart-type control chart based on the median

A nonparametric Shewhart-type control chart is proposed for monitoring the location of a continuous variable in a Phase I process control setting. The chart is based on the pooled median of the available Phase I samples and the charting statistics are the counts (number of observations) in each sample that are less than the pooled median. An exact expression for the false alarm probability (FAP) is given in terms of the multivariate hypergeometric distribution and this is used to provide tables for the control limits for a specified nominal FAP value (of 0.01, 0.05 and 0.10, respectively) and for some values of the sample size (n) and the number of Phase I samples (m). Some approximations are discussed in terms of the univariate hypergeometric and the normal distributions. A simulation study shows that the proposed chart performs as well as, and in some cases better than, an existing Shewhart-type chart based on the normal distribution. Numerical examples are given to demonstrate the implementation of the new chart.     

A generalization for two-sided power distributions and adjusted method of moments

A new rich class of generalized two-sided power (TSP) distributions, where their density functions are expressed in terms of the Gauss hypergeometric functions, is introduced and studied. In this class, the symmetric distributions are supported by finite intervals and have normal shape densities. Our study on TSP distributions also leads us to a new class of discrete distributions on {0, 1, …, k}. In addition, a new numerical method for parameter estimation using moments is given.     

A weighted likelihood ratio of two related negative hypergeometric distributions

In this paper we consider some related negative hypergeometric distributions arising from the problem of sampling without replacement from an urn containing balls of different colours and in different proportions but stopping only after some specific number of balls of different colours have been obtained. With the aid of some simple recurrence relations and identities we obtain in the case of two colours the moments for the maximum negative hypergeometric distribution, the minimum negative hypergeometric distribution, the likelihood ratio negative hypergeometric distribution and consequently the likelihood proportional negative hypergeometric distribution. To the extent that the sampling scheme is applicable to modelling data as illustrated with a biological example and, in fact, many situations of estimating Bernoulli parameters for binary traits within a finite population, these are important first-step results.     

Further results on the trace of a noncentral wishart matrix

In an earlier article Mathai (1980) has given compact representations for the moments and cumulants of the trace of a noncentral Wishart matrix. He has also shown that (trA-ntr;∑)/(2ntri∑²)172. is asymptotically standard normal where A is a noncentral wishart matrix with n degrees of freedom and covariance matrix [0, In the present article explicit expressions for the exact density of the trace are given in terms of confluent hypergeometric functions and in terms of zonal polynomials for the general case and as finite sums when the sample size is odd. As a consequence of some of these representations some summation formulae for zonal polynomials are also given     

A general non-central hypergeometric distribution

We construct a general non-central hypergeometric distribution, which models biased sampling without replacement. Our distribution is constructed from the combined order statistics of two samples: one of independent and identically distributed random variables with absolutely continuous distribution F and the other of independent and identically distributed random variables with absolutely continuous distribution G. The distribution depends on F and G only through F○G^{( ? 1)} (F composed with the quantile function of G), and the standard hypergeometric distribution and Wallenius’ non-central hypergeometric distribution arise as special cases. We show in efficient economic markets the quantity traded has a general non-central hypergeometric distribution.     

A Lagrangian Non Central Negative Binomial Distribution of the First Kind

A Lagrangian probability distribution of the first kind is proposed. Its probability mass function is expressed in terms of generalized Laguerre polynomials or, equivalently, a generalized hypergeometric function. The distribution may also be formulated as a Charlier series distribution generalized by the generalizing Consul distribution and a non central negative binomial distribution generalized by the generalizing Geeta distribution. This article studies formulation and properties of the distribution such as mixture, dispersion, recursive formulas, conditional distribution and the relationship with queuing theory. Two illustrative examples of application to fitting are given.     

New polya and inverse polya distributions of order k

New Polya and inverse Polya distributions of order k are derived by means of generalized urn models and by compounding the binomial and negative binomial distributions of order k of Philippou (1986, 1983) with the beta distribution. It i s noted that the present Polpa distribution of order k includes as special cases a new hypergeometric distribution of order k, a negative one,an inverse one, and a discrete uniform of the same order. The probability generating functions, means and variances of the new distributions are obtained, and five asymptotic results are established relating them to the abovedmentioned binomial and negative binomial distributions of order k, and to the Poisson distribution of the same order of Philippou (1983).Moment estimates are also given and applications are indicated.     

An introduction to hypergeometric functions for economists

Hypergeometric functions are a generalization of exponential functions. They are explicit, computable functions that can also be manipulated analytically. The functions and series we use in quantitative economics are all special cases of them. In this paper, a unified approach to hypergeometric functions is given. As a result, some potentially useful general applications emerge in a number of areas such as in econometrics and economic theory. The greatest benefit from using these functions stems from the fact that they provide parsimonious explicit (and interpretable) solutions to a wide range of general problems.     

EXACT AND APPROXIMATED RELATIONS BETWEEN NEGATIVE HYPERGEOMETRIC AND NEGATIVE BINOMIAL PROBABILITIES

We derive orthogonal expansions in terms of the Meixner polynomials of the first kind for hypergeometric probabilities. We show how these expansions can be used to obtain negative binomial approximations to negative hypergeometric probabilities. Some limit properties of these approximations are studied and also the extension of these results to cumulative probabilities.     

A Conway Maxwell Poisson type generalization of the negative hypergeometric distribution

Abstract

Negative hypergeometric distribution arises as a waiting time distribution when we sample without replacement from a finite population. It has applications in many areas such as inspection sampling and estimation of wildlife populations. However, as is well known, the negative hypergeometric distribution is over-dispersed in the sense that its variance is greater than the mean. To make it more flexible and versatile, we propose a modified version of negative hypergeometric distribution called COM-Negative Hypergeometric distribution (COM-NH) by introducing a shape parameter as in the COM-Poisson and COMP-Binomial distributions. It is shown that under some limiting conditions, COM-NH approaches to a distribution that we call the COM-Negative binomial (COMP-NB), which in turn, approaches to the COM Poisson distribution. For the proposed model, we investigate the dispersion characteristics and shape of the probability mass function for different combinations of parameters. We also develop statistical inference for this model including parameter estimation and hypothesis tests. In particular, we investigate some properties such as bias, MSE, and coverage probabilities of the maximum likelihood estimators for its parameters by Monte Carlo simulation and likelihood ratio test to assess shape parameter of the underlying model. We present illustrative data to provide discussion.     

Optimal domain estimation under summation restriction

The domain estimators that do not sum up to the population total (estimated or known) are considered. In order to achieve their additivity, the theory of the general restriction (GR)-estimator [Knottnerus P., 2003. Sample Survey Theory: Some Pythagorean Perspectives. Springer, New York] is used. The elaborated domain GR-estimators are optimal, they have the minimum variance in a class of estimators that satisfy summation restriction. Furthermore, their variances are smaller than the variances of the corresponding initial domain estimators. The variance/covariance formulae of the domain GR-estimators are explicitly given.The ratio estimators as representatives of the non-additive domain estimators are considered. Their design-based covariance matrix, being crucial for the GR-estimator, is presented. Its structure simplifies under certain assumptions on sampling design (and population model). The corresponding simpler forms of the domain GR-estimators are elaborated as well. The hypergeometric [Traat I., Ilves M., 2007. The hypergeometric sampling design, theory and practice. Acta Appl. Math. 97, 311–321] and the simple random sampling designs are considered in more detail. The results are illustrated in a simulation study where the optimal domain estimator displays its superiority among other meaningful domain estimators. It is noteworthy that due to the imposed restrictions also these other estimators, though not optimal, can be much more precise than the initial estimators.     

On unbiased estimation of gini index and yntema-pietra index of lognormal distribution and their variances

Here we derive the uniformly minimum variance unbiased (best) estimator and a strongly consistent, asymptotically normal, unbiased estimator of each of Gini index and Yntema-Pietra index of lognormal distribution . These estimators are in terms of generalized hypergeometric functions ₁F₂. Further, the variances of these estimators and the best estimators of variances of best estimators are found out. They are in terms of Kempé de Fériet's hypergeometric functions.     

Construction on non-adaptive hypergeometric group testing desings for identifying at most two defectives

In the literature a systematic method of obtaining a group testing design is not available at present. Weideman and Raghavarao (1987a, b) gave methods for the construction of non - adaptive hypergeometric group testing designs for identifying at most two defectives by using a dual method. In the present investigation we have developed a method of construction of group testing designs from (i) Hypercubic Designs for t ≡ 3 (mod 6) and (ii) Balanced Incomplete Block Designs for t ≡ 1 (mod 6) and t ≡ 3 (mod 6). These constructions are accomplished by the use of dual designs. The designs so constructed satisfy specified properties and attained an optimal bound as discussed by Weidman and Raghavarao (1987a,b). Here it is also shown that the condition for pairwise disjoint sets of BIBD for t ≡ 1 (mod 6) given by Weideman and Raghavarao (1987b) is not true for all such designs.