An empirical investigation of different operating characteristics of several estimators of the intraclass correlation in the analysis of binary data

This paper is concerned with properties (bias, standard deviation, mean square error and efficiency) of twenty six estimators of the intraclass correlation in the analysis of binary data. Our main interest is to study these properties when data are generated from different distributions. For data generation we considered three over-dispersed binomial distributions, namely, the beta-binomial distribution, the probit normal binomial distribution and a mixture of two binomial distributions. The findings regarding bias, standard deviation and mean squared error of all these estimators, are that (a) in general, the distributions of biases of most of the estimators are negatively skewed. The biases are smallest when data are generated from the beta-binomial distribution and largest when data are generated from the mixture distribution; (b) the standard deviations are smallest when data are generated from the beta-binomial distribution; and (c) the mean squared errors are smallest when data are generated from the beta-binomial distribution and largest when data are generated from the mixture distribution. Of the 26, nine estimators including the maximum likelihood estimator, an estimator based on the optimal quadratic estimating equations of Crowder (1987), and an analysis of variance type estimator is found to have least amount of bias, standard deviation and mean squared error. Also, the distributions of the bias, standard deviation and mean squared error for each of these estimators are, in general, more symmetric than those of the other estimators. Our findings regarding efficiency are that the estimator based on the optimal quadratic estimating equations has consistently high efficiency and least variability in the efficiency results. In the important range in which the intraclass correlation is small (≤0 5), on the average, this estimator shows best efficiency performance. The analysis of variance type estimator seems to do well for larger values of the intraclass correlation. In general, the estimator based on the optimal quadratic estimating equations seems to show best efficiency performance for data from the beta-binomial distribution and the probit normal binomial distribution, and the analysis of variance type estimator seems to do well for data from the mixture distribution.     

A statistical problem arising in ascertainment

In this paper, an alternative model is examined for the distribution arising out of ascertainment. The weighted beta-binomial is suggested for this purpose as it incorporates the variability in the parameter ø of the weighted binomial distribution. The latter distribution has been the model considered by Rao (1965, 1985) and Kocherlakota and Kocherlakota (1990). Techniques for separate families introduced in Kocherlakota and Kocherlakota (1986) are applied to demonstrate that the weighted beta-binomial model is more appropriate in this situation.     

Bootstrap goodness-of-fit test for the beta-binomial model

A common question in the analysis of binary data is how to deal with overdispersion. One widely advocated sampling distribution for overdispersed binary data is the beta-binomial model. For example, this distribution is often used to model litter effects in toxicological experiments. Testing the null hypothesis of a beta-binomial distribution against all other distributions is difficult, however, when the litter sizes vary greatly. Herein, we propose a test statistic based on combining Pearson statistics from individual litter sizes, and estimate the p-value using bootstrap techniques. A Monte Carlo study confirms the accuracy and power of the test against a beta-binomial distribution contaminated with a few outliers. The method is applied to data from environmental toxicity studies.     

Sensitivity of a bayesian inference to prior assumptions

The sensitivity of-a Bayesian inference to prior assumptions is examined by Monte Carlo simulation for the beta-binomial conjugate family of distributions. Results for the effect on a Bayesian probability interval of the binomial parameter indicate that the Bayesian inference is for the most part quite sensitive to misspecification of the prior distribution. The magnitude of the sensitivity depends primarily on the difference of assigned means and variances from the respective means and variances of the actually-sampled prior distributions. The effect of a disparity in form between the assigned prior and actually-sampled distributions was less important for the cases tested.     

A generalization of the logistic linear'model

Consider the logistic linear model, with some explanatory variables overlooked. Those explanatory variables may be quantitative or qualitative. In either case, the resulting true response variable is not a binomial or a beta-binomial but a sum of binomials. Hence, standard computer packages for logistic regression can be inappropriate even if an overdispersion factor is incorporated. Therefore, a discrete exponential family assumption is considered to broaden the class of sampling models. Likelihood and Bayesian analyses are discussed. Bayesian computation techniques such as Laplacian approximations and Markov chain simulations are used to compute posterior densities and moments. Approximate conditional distributions are derived and are shown to be accurate. The Markov chain simulations are performed effectively to calculate posterior moments by using the approximate conditional distributions. The methodology is applied to Keeler's hardness of winter wheat data for checking binomial assumptions and to Matsumura's Accounting exams data for detailed likelihood and Bayesian analyses.     

Discrete Linnik Weibull distribution

Abstract

We introduce a new family of distributions using truncated discrete Linnik distribution. This family is a rich family of distributions which includes many important families of distributions such as Marshall–Olkin family of distributions, family of distributions generated through truncated negative binomial distribution, family of distributions generated through truncated discrete Mittag–Leffler distribution etc. Some properties of the new family of distributions are derived. A particular case of the family, a five parameter generalization of Weibull distribution, namely discrete Linnik Weibull distribution is given special attention. This distribution is a generalization of many distributions, such as extended exponentiated Weibull, exponentiated Weibull, Weibull truncated negative binomial, generalized exponential truncated negative binomial, Marshall-Olkin extended Weibull, Marshall–Olkin generalized exponential, exponential truncated negative binomial, Marshall–Olkin exponential and generalized exponential. The shape properties, moments, median, distribution of order statistics, stochastic ordering and stress–strength properties of the new generalized Weibull distribution are derived. The unknown parameters of the distribution are estimated using maximum likelihood method. The discrete Linnik Weibull distribution is fitted to a survival time data set and it is shown that the distribution is more appropriate than other competitive models.     

Exact Linear Restrictions on Parameters in the Classical Linear Regression Model

We consider a simple sampling scheme where after each drawing there is a “replacement” whose magnitude is a fixed real number. This gives a unified approach to a family of discrete distributions that includes the hypergeometric, binomial, and Polya distributions as well as their multivariate versions.     

PARAMETER ESTIMATION AND APPLICATIONS FOR A GENERALISATION OF THE BETA-BINOMIAL DISTRIBUTION

A three-parameter generalisation of the beta-binomial distribution (BBD) derived by Chandon (1976) is examined. We obtain the maximum likelihood estimates of the parameters and give the elements of the information matrix. To exhibit the applicability of the generalised distribution we show how it gives an improved fit over the BBD for magazine exposure and consumer purchasing data. Finally we derive an empirical Bayes estimate of a binomial proportion based on the generalised beta distribution used in this study.     

A generalized binomial distribution

A new generalization of the binomial distribution is introduced that allows dependence between trials, nonconstant probabilities of success from trial to trial, and which contains the usual binomial distribution as a special case. Along with the number of trials and an initial probability of ‘success’, an additional parameter that controls the degree of correlation between trials is introduced. The resulting class of distributions includes the binomial, unirnodal distributions, and bimodal distributions. Formulas for the moments, mean, and variance of this distribution are given along with a method for fitting the distribution to sample data.     

A family of bivariate binomial distributions generated by extreme bernoulli distributions

In this paper, bivariate binomial distributions generated by extreme bivariate Bernoulli distributions are obtained and studied. Representation of the bivariate binomial distribution generated by a convex combination of extreme bivariate Bernoulli distributions as a mixture of distributions in the class of bivariate binomial distribution generated by extreme bivariate Bernoulli distribution is obtained. A subfamily of bivariate binomial distributions exhibiting the property of positive and negative dependence is constructed. Some results on positive dependence notions as it relates to the bivariate binomial distribution generated by extreme bivariate Bernoulli distribution and a linear combination of such distributions are obtained.     

Comparison of dichotomized and distributional approaches in rare event clinical trial design: a fixed Bayesian design

This research was motivated by our goal to design an efficient clinical trial to compare two doses of docosahexaenoic acid supplementation for reducing the rate of earliest preterm births (ePTB) and/or preterm births (PTB). Dichotomizing continuous gestational age (GA) data using a classic binomial distribution will result in a loss of information and reduced power. A distributional approach is an improved strategy to retain statistical power from the continuous distribution. However, appropriate distributions that fit the data properly, particularly in the tails, must be chosen, especially when the data are skewed. A recent study proposed a skew-normal method. We propose a three-component normal mixture model and introduce separate treatment effects at different components of GA. We evaluate operating characteristics of mixture model, beta-binomial model, and skew-normal model through simulation. We also apply these three methods to data from two completed clinical trials from the USA and Australia. Finite mixture models are shown to have favorable properties in PTB analysis but minimal benefit for ePTB analysis. Normal models on log-transformed data have the largest bias. Therefore we recommend finite mixture model for PTB study. Either finite mixture model or beta-binomial model is acceptable for ePTB study.     

On a class of bivariate mixed Sarmanov distributions

Multivariate distributions are more and more used to model the dependence encountered in many fields. However, classical multivariate distributions can be restrictive by their nature, while Sarmanov's multivariate distribution, by joining different marginals in a flexible and tractable dependence structure, often provides a valuable alternative. In this paper, we introduce some bivariate mixed Sarmanov distributions with the purpose to extend the class of bivariate Sarmanov distributions and to obtain new dependency structures. Special attention is paid to the bivariate mixed Sarmanov distribution with Poisson marginals and, in particular, to the resulting bivariate Sarmanov distributions with negative binomial and with Poisson‐inverse Gaussian marginals; these particular types of mixed distributions have possible applications in, for example modelling bivariate count data. The extension to higher dimensions is also discussed. Moreover, concerning the dependency structure, we also present some correlation formulas.     

On maximum likelihood estimation of the binomial parameter n

This paper gives a necessary and sufficient condition for the existence of a finite conditional maximum-likelihood estimate for the binomial parameter n. Upper bounds for the conditional and beta-binomial maximum-likelihood estimators are derived. An example is given to show that the conditional likelihood function and the beta-binomial likelihood function may not be unimodal.     

A new class of modified binomial distributions with applications to certain toxicological experiments ∗

A large class of distributions is proposed to fit the binary data obtained from certain toxicological experiments in which, for example, the outcome of interest is the occurrence of dead or malformed fetuses in a litter. This class of distribution includes the additive model proposed by Altham (1978) as a special case. The fits to three real-life data sets using this new distribution are shown to be much better than those provided by beta-binomial distribution used by Williams (1975) and by the correlated-binomial distribution proposed by Kupper and Haseman (1978).     

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.     

Properties and applications of the Poisson-reciprocal inverse Gaussian distribution

The barely known continuous reciprocal inverse Gaussian distribution is used in this paper to introduce the Poisson-reciprocal inverse Gaussian discrete distribution. Several of its most relevant statistical properties are examined, some of them directly inherited from the reciprocal of the inverse Gaussian distribution. Furthermore, a mixed Poisson regression model that uses the reciprocal inverse Gaussian as mixing distribution is presented. Parameters estimation in this regression model is performed via an EM type algorithm. In light of the numerical results displayed in the paper, the distributions introduced in this work are competitive with the classical negative binomial and Poisson-inverse Gaussian distributions.     

A full likelihood procedure of exchangeable negative binomials for modelling correlated and overdispersed count data

This paper introduces an exchangeable negative binomial distribution resulting from relaxing the independence of the Bernoulli sequence associated with a negative binomial distribution to exchangeability. It is demonstrated that the introduced distribution is a mixture of negative binomial distributions and can be characterized by infinitely many parameters that form a completely monotone sequence. The moments of the distribution are derived and a small simulation is conducted to illustrate the distribution. For data analytic purposes, two methods, truncation and completely-monotone links, are given for converting the saturated distribution of infinitely many parameters to parsimonious distributions of finitely many parameters. A full likelihood procedure is described which can be used to investigate correlated and overdispersed count data common in biomedical sciences and teratology. In the end, the introduced distribution is applied to analyze a real clinical data of burn wounds on patients.     

Modeling asset returns with alternative stable distributions

In the 1960's Benoit Mandelbrot and Eugene Fama argued strongly in favor of the stable Paretian distribution as a model for the unconditional distribution of asset returns. Although a substantial body of subsequent empirical studies supported this position, the stable Paretian model plays a minor role in current empirical work.

While in the economics and finance literature stable distributions are virtually exclusively associated with stable Paretian distributions, in this paper we adopt a more fundamental view and extend the concept of stability to a variety of probabilistic schemes. These schemes give rise to alternative stable distributions, which we compare empirically using S&P 500 stock return data. In this comparison the Weibull distribution, associated with both the nonrandom-minimum and geometric-random summation schemes dominates the other stable distributions considered-including the stable Paretian model.     

Modeling asset returns with alternative stable distributions *

In the 1960's Benoit Mandelbrot and Eugene Fama argued strongly in favor of the stable Paretian distribution as a model for the unconditional distribution of asset returns. Although a substantial body of subsequent empirical studies supported this position, the stable Paretian model plays a minor role in current empirical work.

While in the economics and finance literature stable distributions are virtually exclusively associated with stable Paretian distributions, in this paper we adopt a more fundamental view and extend the concept of stability to a variety of probabilistic schemes. These schemes give rise to alternative stable distributions, which we compare empirically using S&P 500 stock return data. In this comparison the Weibull distribution, associated with both the nonrandom-minimum and geometric-random summation schemes dominates the other stable distributions considered-including the stable Paretian model.     

Local efficiency of moment estimators in beta-binomial model

In this paper, the beta-binomial model is introduced as a Markov chain. It is shown that the correlated binomial model of Kupper and Haseman (1978) is identical to the additive binomial model of AItham(1978) and both are a first order approximation of the beta-binomial model. For small γ, the local efficiency of the moment estimators for the mean ρ and the extra-binomial variation γ is examined analytically. It is shown that, locally, the moment estimator for p is efficient up to the second order of y. Exact formulae for the relative efficiency are obtained for both the cases with γ known and unknown. Generalization to the unequal sample size case is also carried out. In particular, the gain in efficiency by using the quasi-likelihood estimator instead of the ratio estimator for p is studied when γ is known. These results are in agreement with the Monte Carlo results of Kleinman(1973) and Crowder(1985).