On bivariate discrete Weibull distribution

Recently, Lee and Cha proposed two general classes of discrete bivariate distributions. They have discussed some general properties and some specific cases of their proposed distributions. In this paper we have considered one model, namely bivariate discrete Weibull distribution, which has not been considered in the literature yet. The proposed bivariate discrete Weibull distribution is a discrete analogue of the Marshall–Olkin bivariate Weibull distribution. We study various properties of the proposed distribution and discuss its interesting physical interpretations. The proposed model has four parameters, and because of that it is a very flexible distribution. The maximum likelihood estimators of the parameters cannot be obtained in closed forms, and we have proposed a very efficient nested EM algorithm which works quite well for discrete data. We have also proposed augmented Gibbs sampling procedure to compute Bayes estimates of the unknown parameters based on a very flexible set of priors. Two data sets have been analyzed to show how the proposed model and the method work in practice. We will see that the performances are quite satisfactory. Finally, we conclude the paper.     

A simple nonparametric test for bivariate symmetry about a line

Over the years many researchers have dealt with testing the hypotheses of symmetry in univariate and multivariate distributions in the parametric and nonparametric setup. In a multivariate setup, there are several formulations of symmetry, for example, symmetry about an axis, joint symmetry, marginal symmetry, radial symmetry, symmetry about a known point, spherical symmetry, and elliptical symmetry among others. In this paper, for the bivariate case, we formulate a concept of symmetry about a straight line passing through the origin in a plane and accordingly develop a simple nonparametric test for testing the hypothesis of symmetry about a straight line. The proposed test is based on a measure of deviance between observed counts of bivariate samples in suitably defined pairs of sets. The exact null distribution and non-null distribution, for specified classes of alternatives, of the test statistics are obtained. The null distribution is tabulated for sample size from n=5 up to n=30. The null mean, null variance and the asymptotic null distributions of the proposed test statistics are also obtained. The empirical power of the proposed test is evaluated by simulating samples from the suitable class of bivariate distributions. The empirical findings suggest that the test performs reasonably well against various classes of asymmetric bivariate distributions. Further, it is advocated that the basic idea developed in this work can be easily adopted to test the hypotheses of exchangeability of bivariate random variables and also bivariate symmetry about a given axis which have been considered by several authors in the past.     

Report of the ASA Technical Panel on the Census Undercount

A new model is proposed for the joint distribution of paired survival times generated from clinical trials and certain reliability settings. The new model can be considered an extension to the bivariate exponential models studied in the literature. Here, a more flexible bivariate Weibull model will be derived, and two exact parametric tests for testing the equality of marginal survival distributions are developed.     

Absolute continuous bivariate generalized exponential distribution

Generalized exponential distribution has been used quite effectively to model positively skewed lifetime data as an alternative to the well known Weibull or gamma distributions. In this paper we introduce an absolute continuous bivariate generalized exponential distribution by using a simple transformation from a well known bivariate exchangeable distribution. The marginal distributions of the proposed bivariate generalized exponential distributions are generalized exponential distributions. The joint probability density function and the joint cumulative distribution function can be expressed in closed forms. It is observed that the proposed bivariate distribution can be obtained using Clayton copula with generalized exponential distribution as marginals. We derive different properties of this new distribution. It is a five-parameter distribution, and the maximum likelihood estimators of the unknown parameters cannot be obtained in closed forms. We propose some alternative estimators, which can be obtained quite easily, and they can be used as initial guesses to compute the maximum likelihood estimates. One data set has been analyzed for illustrative purposes. Finally we propose some generalization of the proposed model.     

A generalization of the bivariate Beta-Binomial distribution

This paper presents a new bivariate discrete distribution that generalizes the bivariate Beta-Binomial distribution. It is generated by Appell hypergeometric function F₁ and can be obtained as a Binomial mixture with an Exton's Generalized Beta distribution. The model has different marginal distributions which are, together with the conditional distributions, more flexible than the Beta-Binomial distribution. It has non-linear regression curves and is useful for random variables with positive correlation. These features make the model very adequate to fit observed data as the two applications included show.     

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.     

Testing for departures from nominal dispersion in generalized nonlinear models with varying dispersion and/or additive random effects

This paper discusses the tests for departures from nominal dispersion in the framework of generalized nonlinear models with varying dispersion and/or additive random effects. We consider two classes of exponential family distributions. The first is discrete exponential family distributions, such as Poisson, binomial, and negative binomial distributions. The second is continuous exponential family distributions, such as normal, gamma, and inverse Gaussian distributions. Correspondingly, we develop a unifying approach and propose several tests for testing for departures from nominal dispersion in two classes of generalized nonlinear models. The score test statistics are constructed and expressed in simple, easy to use, matrix formulas, so that the tests can easily be implemented using existing statistical software. The properties of test statistics are investigated through Monte Carlo simulations.     

Approximation of the Bivariate Renewal Function

Recently, Kambo and his co-researchers (2012) proposed a method of approximation for evaluating the one-dimensional renewal function based on the first three moments. Their method is simple and elegant, which gives exact values for well-known distributions. In this article, we propose an analogous method for the evaluation of bivariate renewal function based on the first two moments of the variables and their joint moment. The proposed method yields exact results for certain widely used bivariate distributions like bivariate exponential distribution, bivariate Weibull distributions, and bivariate Pareto distributions. An illustrative example in the form of a two-dimensional warranty problem is considered and comparisons of our method are made with the results of other models.     

Discrimination among bivariate beta-generated distributions

Some properties of the general families of bivariate distributions generated by beta dependent random variables are derived and discussed here. Some classic measures of dependence and information are derived, and their behaviours and properties are discussed as well. Finally, a discrimination procedure within this general family of bivariate distributions is proposed based on Shannon entropy. A real-life example is presented to illustrate the model as well as the inferential results developed here.     

Discrete and continuous bivariate lifetime models in presence of cure rate: a comparative study under Bayesian approach

The modeling and analysis of lifetime data in which the main endpoints are the times when an event of interest occurs is of great interest in medical studies. In these studies, it is common that two or more lifetimes associated with the same unit such as the times to deterioration levels or the times to reaction to a treatment in pairs of organs like lungs, kidneys, eyes or ears. In medical applications, it is also possible that a cure rate is present and needed to be modeled with lifetime data with long-term survivors. This paper presented a comparative study under a Bayesian approach among some existing continuous and discrete bivariate distributions such as the bivariate exponential distributions and the bivariate geometric distributions in presence of cure rate, censored data and covariates. In presence of lifetimes related to cured patients, it is assumed standard mixture cure rate models in the data analysis. The posterior summaries of interest are obtained using Markov Chain Monte Carlo methods. To illustrate the proposed methodology two real medical data sets are considered.     

On the Moments of Multivariate Discrete Distributions Using Finite Difference Operators

This article discusses a general approach to finding the moments of two classes of multivariate discrete distributions, which include those widely used in applied and theoretical statistics. The two classes of multivariate discrete distributions are the multivariate generalized power series distributions (GPSD) and the unified multivariate hypergeometric (UMH) Distributions. The results of Link (1981) follow as special cases.     

Weighted bivariate logarithmic series distributions

In this paper, we are interested in the weighted distributions of a bivariate three parameter logarithmic series distribution studied by Kocherlakota and Kocherlakota (1990). The weighted versions of the model are derived with weight W(x,y) = x^[r] y^[s]. Explicit expressions for the probability mass function and probability generating functions are derived in the case r = s = l. The marginal and conditional distributions are derived in the general case. The maximum likelihood estimation of the parameters, in both two parameter and three parameter cases, is studied. A procedure for computer generation of bivariate data from a discrete distribution is described. This enables us to present two examples, in order to illustrate the methods developed, for finding the maximum likelihood estimates.     

Modelling Poisson marked point processes using bivariate mixture transition distributions

A class of bivariate continuous-discrete distributions is proposed to fit Poisson dynamic models in a single unified framework via bivariate mixture transition distributions (BMTDs). Potential advantages of this class over the current models include its ability to capture stretches, bursts and nonlinear patterns characterized by Internet network traffic, high-frequency financial data and many others. It models the inter-arrival times and the number of arrivals (marks) in a single unified model which benefits from the dependence structure of the data. The continuous marginal distributions of this class include as special cases the exponential, gamma, Weibull and Rayleigh distributions (for the inter-arrival times), whereas the discrete marginal distributions are geometric and negative binomial. The conditional distributions are Poisson and Erlang. Maximum-likelihood estimation is discussed and parameter estimates are obtained using an expectation–maximization algorithm, while the standard errors are estimated using the missing information principle. It is shown via real data examples that the proposed BMTD models appear to capture data features better than other competing models.     

Expectation formulas for integer valued multivariate random variables

Nadarajah and Mitov [Communications in Statistics—Theory and Methods, 32, 2003, 47–60] derived an expectation formula for continuous multivariate random variables involving the joint survival function. Their result is extended here for discrete multivariate random variables. Examples proposing new discrete bivariate distributions are given.     

Global prior distributions for the analysis of discrete graphical models

Summary We propose a new class of prior distributions for the analysis of discrete graphical models. Such a class, obtained following a conditional approach, generalizes the hyper Dirichlet distributions of Dawid and Lauritzen (1993), since it can be extended to non decomposable graphical models. The two classes are compared in terms of model selection, with an application to a medical data-set illustrating the performance of the two resulting procedures. The proposed class turns out to select simpler, more par-simonious structures.     

Modelling conditional distributions in bivariate survival

Conditional distributions for bivariate survival can be obtained via a model for the joint distribution, or, as has sometimes been suggested, by modelling the conditioned variable directly, with the conditioning variable included as a covariate. A quantitative comparison of estimated covariate effects and predictive distributions under the two approaches is given. The results are illustrated in a novel frailty application.     

Parameter estimation in nonlinear regression models 1

Fisher information contained in record values, inter-record times and their concomitants from a sample of fixed size is derived in general and explicit expressions are deduced for some specific known bivariate classes of distributions. A comparison between fixed sampling and inverse sampling schemes with equal number of records and concomitants is also carried out. We also consider parameter estimation based on bivariate records and a small simulation study is done.     

On the monotonic properties of discrete failure rates

As is well known, the monotonicity of failure rate of a life distribution plays an important role in modeling failure time data. In this paper, we develop techniques for the determination of increasing failure rate (IFR) and decreasing failure rate (DFR) property for a wide class of discrete distributions. Instead of using the failure rate, we make use of the ratio of two consecutive probabilities. The method developed is applied to various well known families of discrete distributions which include the binomial, negative binomial and Poisson distributions as special cases. Finally, a formula is presented to determine explicitly the failure rate of the families considered. This formula is used to determine the failure rate of various classes of discrete distributions. These formulas are explicit but complicated and cannot normally be used to determine the monotonicity of the failure rates.     

A discrete analogue of odd Weibull-G family of distributions: properties,classical and Bayesian estimation with applications to count data

In the statistical literature, several discrete distributions have been developed so far. However, in this progressive technological era, the data generated from different fields is getting complicated day by day, making it difficult to analyze this real data through the various discrete distributions available in the existing literature. In this context, we have proposed a new flexible family of discrete models named discrete odd Weibull-G (DOW-G) family. Its several impressive distributional characteristics are derived. A key feature of the proposed family is its failure rate function that can take a variety of shapes for distinct values of the unknown parameters, like decreasing, increasing, constant, J-, and bathtub-shaped. Furthermore, the presented family not only adequately captures the skewed and symmetric data sets, but it can also provide a better fit to equi-, over-, under-dispersed data. After producing the general class, two particular distributions of the DOW-G family are extensively studied. The parameters estimation of the proposed family, are explored by the method of maximum likelihood and Bayesian approach. A compact Monte Carlo simulation study is performed to assess the behavior of the estimation methods. Finally, we have explained the usefulness of the proposed family by using two different real data sets.     

Bayesian latent variable models for clustered mixed outcomes

A general framework is proposed for modelling clustered mixed outcomes. A mixture of generalized linear models is used to describe the joint distribution of a set of underlying variables, and an arbitrary function relates the underlying variables to be observed outcomes. The model accommodates multilevel data structures, general covariate effects and distinct link functions and error distributions for each underlying variable. Within the framework proposed, novel models are developed for clustered multiple binary, unordered categorical and joint discrete and continuous outcomes. A Markov chain Monte Carlo sampling algorithm is described for estimating the posterior distributions of the parameters and latent variables. Because of the flexibility of the modelling framework and estimation procedure, extensions to ordered categorical outcomes and more complex data structures are straightforward. The methods are illustrated by using data from a reproductive toxicity study.