Bivariate exponentiated‐exponential geometric regression model

A bivariate exponentiated‐exponential geometric regression model that allows negative, zero, or positive correlation is defined and studied. The model can accommodate under‐ or over‐dispersed count data. The regression model is based on the univariate exponentiated‐exponential geometric distribution, and the marginal means of the bivariate model are functions of the explanatory variables. The parameters of the bivariate regression model are estimated by using the maximum likelihood method. Some test statistics including goodness of fit are discussed. A simulation study is conducted to compare the model with the bivariate generalized Poisson regression model. One numerical data set is used to illustrate the application of the regression model.     

A bivariate Poisson count data model using conditional probabilities

The applied econometrics of bivariate count data predominantly focus on a bivariate Poisson density with a correlation structure that is very restrictive. The main limitation is that this bivariate distribution excludes zero and negative correlation. This paper introduces a new model which allows for a more flexible correlation structure. To this end the joint density is decomposed by means of the multiplication rule in marginal and conditional densities. Simulation experiments and an application of the model to recreational data are presented.     

On three and five parameter bivariate beta distributions

Summary In this paper two bivariate beta distributions have been studied. The five parameter bivariate beta distribution is derived from the Morgenstern-system of curves while the three parameter distribution is the bivariate Dirichlet distribution. In both cases the distributions of the product and the quotient of random variables are derived and other properties are also studied.     

Limit theorems for multivariate discrete distributions

According to the usual law of small numbers a multivariate Poisson distribution is derived by defining an appropriate model for multivariate Binomial distributions and examining their behaviour for large numbers of trials and small probabilities of marginal and simultaneous successes. The weak limit law is a generalization of Poisson's distribution to larger finite dimensions with arbitrary dependence structure. Compounding this multivariate Poisson distribution by a Gamma distribution results in a multivariate Pascal distribution which is again asymptotically multivariate Poisson. These Pascal distributions contain a class of multivariate geometric distributions. Finally the bivariate Binomial distribution is shown to be the limit law of appropriate bivariate hypergeometric distributions. Proving the limit theorems mentioned here as well as understanding the corresponding limit distributions becomes feasible by using probability generating functions.     

A new bivariate distribution in natural exponential family

We propose a new bivariate distribution following a GLM form i.e., natural exponential family given the constantly correlated covariance matrix. The proposed distribution can represent an independent bivariate gamma distribution as a special case. In order to derive the distribution we utilize an integrating factor method to satisfy the integrability condition of the quasi-score function. The derived distribution becomes a mixture of discrete and absolute continuous distributions. The proposal of our new bivariate distribution will make it possible to develop some bivariate generalized linear models. Further the discrete correlated bivariate distribution will also arise from an independent bivariate Poisson mass function by compounding our proposed distribution (Iwasaki and Tsubaki, 2002).Received March 2003     

Bivariate generalized Poisson distribution with some applications

The univariate generalized Poisson probability model has many applications in various areas such as engineering, manufacturing, survival analysis, genetic, shunting accidents, queuing, and branching processes. A correlated bivariate version of the univariate generalized Poisson distribution is defined and studied. Estimation of its parameters and some of its properties are also discussed.     

Characterizations of bivariate distributions using concomitants of record values

In this paper, we consider a family of bivariate distributions which is a generalization of the Morgenstern family of bivariate distributions. We have derived some properties of concomitants of record values which characterize this generalized class of distributions. The role of concomitants of record values in the unique determination of the parent bivariate distribution has been established. We have also derived properties of concomitants of record values which characterize each of the following families viz Morgenstern family, bivariate Pareto family and a generalized Gumbel’s family of bivariate distributions. Some applications of the characterization results are discussed and important conclusions based on the characterization results are drawn.     

Bivariate semi-Pareto minification processes

Marshall-Olkin bivariate semi-Pareto distribution (MO-BSP) and Marshall-Olkin bivariate Pareto distribution (MO-BP) are introduced and studied. AR(1) and AR(k) time series models are developed with minification structure having MO-BSP stationary marginal distribution. Various characterizations are investigated.Acknowledgements. The authors thank the Editor and the referee for their valuable suggestions which led to an improved version of the original paper. The first author is grateful to the University Grants Commission of India for the support under Teacher Fellowship Scheme.     

A unified approach to bivariate discrete distributions

Here we introduce a bivariate generalized hypergeometric factorial moment distribution (BGHFMD) through its probability generating function (p.g.f.) whose marginal distributions are the generalized hypergeometric factorial moment distributions introduced by Kemp and Kemp (Bull Int Stat Inst 43:336–338,1969). Well-known bivariate versions of distributions such as binomial, negative binomial and Poisson are special cases of this distribution. A genesis of the distribution and explicit closed form expressions for the probability mass function of the BGHFMD, its factorial moments and the p.g.f.’s of its conditional distributions are derived here. Certain recurrence relations for probabilities, moments and factorial moments of the bivariate distribution are also established.     

Bivariate non-central negative binomial distribution: Another generalisation

Summary Another bivariate generalisation (Type V) of the non-central negative binomial distribution is considered. This generalisation is constructed (i) as a latent structure model; (ii) as an extension of an accident proneness model investigated by Edwards/Gurland (1961); and (iii) as a reversible stochastic counter model. The third construction gives, as a result, an apparently new formulation of the Edwards/Gurland model. The probabilities, moments, recurrence formulas and some properties are given. An application to the data used by Holgate (1966) is considered.     

Dose–response modelling for bivariate covariates with and without a spike at zero: theory and application to binary outcomes

In epidemiology and clinical research, there is often a proportion of unexposed individuals resulting in zero values of exposure, meaning that some individuals are not exposed and those exposed have some continuous distribution. Examples are smoking or alcohol consumption. We will call these variables with a spike at zero (SAZ). In this paper, we performed a systematic investigation on how to model covariates with a SAZ and derived theoretical odds ratio functions for selected bivariate distributions. We consider the bivariate normal and bivariate log normal distribution with a SAZ. Both confounding and effect modification can be elegantly described by formalizing the covariance matrix given the binary outcome variable Y. To model the effect of these variables, we use a procedure based on fractional polynomials first introduced by Royston and Altman (1994, Applied Statistics 43: 429–467) and modified for the SAZ situation (Royston and Sauerbrei, 2008, Multivariable model‐building: a pragmatic approach to regression analysis based on fractional polynomials for modelling continuous variables, Wiley; Becher et al., 2012, Biometrical Journal 54: 686–700). We aim to contribute to theory, practical procedures and application in epidemiology and clinical research to derive multivariable models for variables with a SAZ. As an example, we use data from a case–control study on lung cancer.     

A new family of bivariate max-infinitely divisible distributions

In this article we discuss the asymptotic behaviour of the componentwise maxima for a specific bivariate triangular array. Its components are given in terms of linear transformations of bivariate generalised symmetrised Dirichlet random vectors introduced in Fang and Fang (Statistical inference in elliptically contoured and related distributions. Allerton Press, New York, 1990). We show that the componentwise maxima of such triangular arrays is attracted by a bivariate max-infinitely divisible distribution function, provided that the associated random radius is in the Weibull max-domain of attraction.     

The Asymptotic Efficacies and Relative Efficiencies of Various Linear Rank Tests for Independence

In this paper we consider eight general models of independence, the Hájek–Šidák model, the Janssen–Mason model, Konijn’s model, Steffensen’s model, the Farlie model, the bivariate Gamma distribution, the Mardia model and the Frechet model. The asymptotic efficacies and relative efficiencies of various linear rank tests are computed. It turns out that the asymptotic power depends heavily on the underlying model. However, for the vast majority of considered models, the Spearman test is, asymptotically, a good choice.     

A generalisation of Gumbel's bivariate logistic distribution

Summary The bivariate distributionF(x, y)=1/[1+exp(–x)+exp(–y)] was examined byGumbel. We have generalised this expression by raising it to an arbitarary power. Such a distribution may occur as a mixture of bivariate extreme-value distribution. As well as giving its basic properties, we have paid special attention to measures of correlation alternative to the product-moment, namely, Kendall's and Spearman's rank correlations and the product-moment correlation calculated after transforming the marginal distributions into Normal ones. An application to the multifactorial model of disease transmission is outlined.     

Compound gamma,beta and F distributions

Summary In this paper a compound gamma distribution has been derived by compounding a gamma distribution with another gamma distribution. The resulting compound gamma distribution has been reduced to the Beta distributions of the first kind and the second kind and to theF distribution by suitable transformations. This includes theLomax distribution as a special case which enjoys a useful property. Moment estimators for two of its parameters are explicitly obtained, which tend to a bivariate normal distribution. The paper contains expressions for a bivariate probability density function, its conditional expectation, conditional variance and the product moment correlation coefficient. Finally, all the parameters of the compound gamma distribution are explicitly expressed in terms of the functions of the moments of the functions of random variables in two different ways. This note is based on a technical report prepared by the author while he was with the Procter and Gamble Company.     

Estimation and model adequacy checking for multivariate seasonal autoregressive time series models with periodically varying parameters

We introduce a class of multivariate seasonal time series models with periodically varying parameters, abbreviated by the acronym SPVAR. The model is suitable for multivariate data, and combines a periodic autoregressive structure and a multiplicative seasonal time series model. The stationarity conditions (in the periodic sense) and the theoretical autocovariance functions of SPVAR stochastic processes are derived. Estimation and checking stages are considered. The asymptotic normal distribution of the least squares estimators of the model parameters is established, and the asymptotic distributions of the residual autocovariance and autocorrelation matrices in the class of SPVAR time series models are obtained. In order to check model adequacy, portmanteau test statistics are considered and their asymptotic distributions are studied. A simulation study is briefly discussed to investigate the finite-sample properties of the proposed test statistics. The methodology is illustrated with a bivariate quarterly data set on travelers entering in to Canada.     

Non-causality in bivariate binary time series

In this paper we develop a dynamic discrete-time bivariate probit model, in which the conditions for Granger non-causality can be represented and tested. The conditions for simultaneous independence are also worked out. The model is extended in order to allow for covariates, representing individual as well as time heterogeneity. The proposed model can be estimated by Maximum Likelihood. Granger non-causality and simultaneous independence can be tested by Likelihood Ratio or Wald tests. A specialized version of the model, aimed at testing Granger non-causality with bivariate discrete-time survival data is also discussed. The proposed tests are illustrated in two empirical applications.     

Combining competing forecasts of inflation using a bivariate arch model

This paper sets out the basic structure of the bivariate generalization of Engle's ARCH model. Conditions which guarantee that the conditional covariance matrix is well defined are summarized, as are estimation and hypothesis testing.The process is used to combine forecasts where the weights are allowed to vary over time. Forecast errors from competing models are treated as a bivariate ARCH process so that the conditional covariance matrix adapts over time. At each point in time these conditional estimates of the variances and covariances are used to construct the optimal weights for combining the forecasts. Consequently, when one model is fitting well, its variance will be reduced and its weight will be increased.Two models of US inflation are constructed; one is a stylized monetarist model while the other is a mark-up model. The forecast errors are modeled as a simple bivariate ARCH process. Diagnostic tests reveal that this has overly restricted the parameterization of the covariance matrix. An approximation to the theoretically anticipated factor structure model is then estimated. The results in both cases show the weights varying over the sample period in moderately interpretable fashion.     

THE BIVARIATE EXPONENTIAL DISTRIBUTION; A REVIEW AND SOME NEW RESULTS

Various models have been proposed as bivariate forms of the exponential distribution. A brief but comprehensive review is presented which classifies, interrelates and contrasts the different models and outlines what is known about distributional properties, applicability and estimation and testing of parameters (particularly the association parameter). Some new results are presented for one particular model. Maximum likelihood, and moment–type, estimators of the association parameter are examined. Asymptotic variances are derived and attention is given to the relative efficiency of the estimators and to problems of their evaluation.     

The M/M/1 queue with gated random order of service

We analyse an M / M /1 queueing model with gated random order of service discipline. In this service discipline there is a waiting room, in which arriving customers are collected, and a service queue. Each time the service queue becomes empty, all customers in the waiting room are instantaneously put in random order in the service queue. We find the joint stationary distribution of the number of customers in the waiting room and the service queue. Furthermore, we obtain the bivariate Laplace–Stieltjes transform of the joint distribution of the sojourn times of a customer in the waiting room and the service queue.