Multivariate t and beta distributions associated with the multivariate F distribution

Relationships between F, skew t and beta distributions in the univariate case are in this paper extended in a natural way to the multivariate case. The result is two new distributions: a multivariate t/skew t distribution (on ℜ^m) and a multivariate beta distribution (on (0,1)^m). A special case of the former distribution is a new multivariate symmetric t distribution. The new distributions have a natural relationship to the standard multivariate F distribution (on (ℜ⁺)^m) and many of their properties run in parallel. We look at: joint distributions, mathematically and graphically; marginal and conditional distributions; moments; correlations; local dependence; and some limiting cases. Received: March 2001     

A Skew‐normal copula‐driven GLMM

This paper presents a method for fitting a copula‐driven generalized linear mixed models. For added flexibility, the skew‐normal copula is adopted for fitting. The correlation matrix of the skew‐normal copula is used to capture the dependence structure within units, while the fixed and random effects coefficients are estimated through the mean of the copula. For estimation, a Monte Carlo expectation–maximization algorithm is developed. Simulations are shown alongside a real data example from the Framingham Heart Study.     

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.     

Characterizations and Modelling of Multivariate Lack of Memory Property

In this article we establish characterizations of multivariate lack of memory property in terms of the hazard gradient (whenever exists), the survival function and the cumulative hazard function. Based on one of these characterizations we establish a method of generating bivariate lifetime distributions possessing bivariate lack of memory property (BLMP) with specified marginals. It is observed that the marginal distributions have to satisfy certain conditions to be stated. The method generates absolutely continuous bivariate distributions as well as those containing a singular component. Bivariate exponential distributions due to Proschan and Sullo (Reliability and biometry, pp 423–440, 1974), Freund (in J Am Stat Assoc 56:971–977, 1961), Block and Basu (J Am Stat Assoc 89:1091–1097, 1974) and Marshall and Olkin (J Am Math Assoc 62:30–44, 1967) are generated as particular cases among others using the proposed method. Some other distributions generated using the method may be of practical importance. Shock models leading to bivariate distributions possessing BLMP are given. Some closure properties of a class of univariate failure rate functions that can generate distributions possessing BLMP and of the class of bivariate survival functions having BLMP are studied.     

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.     

Description of purchase incidence by multivariate heterogeneous Poisson processes

In this paper the application of bivariate Poisson heterogeneous models to budget data is studied. This study was motivated from inconsistencies that we encountered when univariate Poisson based models were applied to cumulative data sets. Application of a multivariate Poisson based model is a possible solution to this problem. In this paper we will study the feasibility of estimators based on these models.     

Modelling security market events in continuous time: Intensity based,multivariate point process models

A continuous time econometric modelling framework for multivariate financial market event (or ‘transactions’) data is developed in which the model is specified via the vector conditional intensity. Generalised Hawkes models are introduced that incorporate inhibitory events and dependence between trading days. Novel omnibus specification tests based on a multivariate random time change theorem are proposed. A bivariate point process model of the timing of trades and mid-quote changes is then presented for a New York Stock Exchange stock and related to the market microstructure literature. The two-way interaction of trades and quote changes in continuous time is found to be important empirically.     

Probabilistic forecasting in day-ahead electricity markets: Simulating peak and off-peak prices

In this article we include dependency structures for electricity price forecasting and forecasting evaluation. We work with off-peak and peak time series from the German-Austrian day-ahead price; hence, we analyze bivariate data. We first estimate the mean of the two time series, and then in a second step we estimate the residuals. The mean equation is estimated by ordinary least squares and the elastic net, and the residuals are estimated by maximum likelihood. Our contribution is to include a bivariate jump component in a mean reverting jump diffusion model in the residuals. The models’ forecasts are evaluated with use of four different criteria, including the energy score to measure whether the correlation structure between the time series is properly included. It is observed that the models with bivariate jumps provide better results with the energy score, which means that it is important to consider this structure to properly forecast correlated time series.     

Corrections for coarsely categorized measures: LISREL's polyserial and polychoric correlations

Versions 5 and 6 of LISREL (Joreskog and Sorbom, 1981) contain procedures that estimate the underlying correlation between continuous variables on the basis of crude rank category measures. The procedures assume that the distribution of the measured variables would have been bivariate normal if they had not been categorized. Using survey data and simulations, the accuracy of these polyserial/polychoric (P/P) based estimates of the underlying correlations are compared with those based on simple equal distance scoring of the categories. The results indicate that under some conditions, e.g., nearly normally distributed variables and moderate to high correlations, the polyserial/polychoric based estimates are better. Under other conditions, e.g., a moderate to high degree of skew and kurtosis and low correlations, the equal distance score based estimates are better. Under all conditions, the amount of error decreases fairly rapidly as the number of categories is increased from two to five.     

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 and multivariate scaled association models. An application to homogamy of social origin and education in Hungary between 1930 and 1979

In this article we compare bivariate and multivariate models for homogamy of social origin and education to test whether bivariate models of homogamy lead to biased results. We use data on Hungarian couples married between 1930 and 1979 and loglinear models of scaled association. The results indicate some differences between bivariate and multivariate analyses. At each point of time bivariate models overestimate homogamy, both with respect to education and social origin. However, results on trends in time do not differ much between the two analyses. The exception is the period 1940–1959, in which bivariate analysis showed decreasing educational homogamy, and multivariate analysis showed an increasing trend. The latter finding can be explained by declining homogamy of social origin, as well as the weaker reproduction and cross-effects in this period.     

Replication of the results in ‘learning about heterogeneity in returns to schooling’

A recent article (Koop and Tobias, 2004 ) proposes a direct way to characterize the extent of heterogeneity in returns to education. They investigate the adequacy of several competing models and conclude that returns to schooling are heterogeneous and are best modelled as a bivariate normal distribution. The results of this replication paper basically agree with the authors. Copyright © 2005 John Wiley & Sons, Ltd.     

THE ROLE OF CONDITIONAL HETEROSKEDASTICITY IN IDENTIFYING AND ESTIMATING LINEAR TRIANGULAR SYSTEMS,WITH APPLICATIONS TO ASSET PRICING MODELS THAT INCLUDE A MISMEASURED FACTOR

A new estimator is proposed for linear triangular systems, where identification results from the model errors following a bivariate and diagonal GARCH(1,1) process with potentially time‐varying error covariances. This estimator applies when traditional instruments are unavailable. I demonstrate its usefulness on asset pricing models like the capital asset pricing model and Fama–French three‐factor model. In the context of a standard two‐pass cross‐sectional regression approach, this estimator improves the pricing performance of both models. Set identification bounds and an associated estimator are also provided for cases where the conditions supporting point identification fail. Copyright © 2013 John Wiley & Sons, Ltd.     

A bivariate count data model for household tourism demand

Households' choice of the number of leisure trips and the total number of overnight stays is empirically studied using Swedish tourism data. A bivariate hurdle approach separating the participation (to travel and stay the night or not) from the quantity (the number of trips and nights) decision is employed. The quantity decision is modelled with a bivariate mixed Poisson lognormal model allowing for both positive as well as negative correlation between count variables. The observed endogenous variables are drawn from a truncated density and estimation is pursued by simulated maximum likelihood. The estimation results indicate a negative correlation between the number of trips and nights. In most cases own price effects are as expected negative, while estimates of cross‐price effects vary between samples. Copyright © 2005 John Wiley & Sons, Ltd.     

A negative binomial thinning-based bivariate INAR(1) process

This article considers a bivariate INAR(1) process based on an extension of the negative binomial thinning operator by prespecifying the distribution of the innovations. The dependence is introduced through the innovation components. The existence, uniqueness, strict stationarity, ergodicity, and some probabilistic properties of the process are derived. The estimation methods of conditional least squares and conditional maximum likelihood are considered. Some numerical results of the estimates are presented by simulation study. An application to crime data set is provided.     

Distribution of correlation coefficient under mixture normal model

The present paper obtains the nonnull distribution of the product moment correlation coefficient r when sample is drawn from a mixture of two bivariate Gaussian distributions. The moments of 1−r ² have been used to derive the nonnull density of r. Received September 2000     

Semiparametric estimation of a bivariate Tobit model

The existing semiparametric estimation literature has mainly focused on univariate Tobit models and no semiparametric estimation has been considered for bivariate Tobit models. In this paper, we consider semiparametric estimation of the bivariate Tobit model proposed by Amemiya (1974), under the independence condition without imposing any parametric restriction on the error distribution. Our estimator is shown to be consistent and asymptotically normal, and simulation results show that our estimator performs well in finite samples. It is also worth noting that while Amemiya’s (1974) instrumental variables estimator (IV) requires the normality assumption, our semiparametric estimator actually outperforms his IV estimator even when normality holds. Our approach can be extended to higher dimensional multivariate Tobit models.     

A new bivariate generalized Poisson distribution

In this paper, a new bivariate generalized Poisson distribution (GPD) that allows any type of correlation is defined and studied. The marginal distributions of the bivariate model are the univariate GPDs. The parameters of the bivariate distribution are estimated by using the moment and maximum likelihood methods. Some test statistics are discussed and one numerical data set is used to illustrate the applications of the bivariate model.