Characteristic properties of multivariate survival functions in terms of residual life distributions

This paper discusses the relationships among some characteristic properties of the multivariate survival function based on the residual life distribution, and provides the conditions for their equivalence. In the meanwhile, the corrected version of Ma (1996, Theorem 1 (ii) and (iii)) is given.     

Multivariate gamma distributions with one-factorial accompanying correlation matrices and applications to the distribution of the multivariate range

Summary A new representation for the characteristic function of the joint distribution of the Mahalanobis distances betweenk independentN(μ, Σ)-distributed points is given. Especially fork=3 the corresponding distribution function is obtained as a special case of multivariate gamma distributions whose accompanying normal distribution has a positive semidefinite correlation matrix with correlationsϱ _ij=−a _i a _j. These gamma distribution functions are given here by one-dimensional parameter integrals. With some further trivariate gamma distributions third order Bonferroni inequalities are derived for the upper tails of the distribution function of the multivariate range ofk independentN(μ, I)-distributed points. From these inequalities very accurate (conservative) approximations to upperα-level bounds can also be computed for studentized multivariate ranges.     

An Appraisal and Bibliography of Tests for Multivariate Normality

This paper is a review of many of the dozens of procedures currently available for testing a data set for goodness-of-fit to the multivariate normal distribution. A majority of the procedures can be placed into one of four basic categories. Most procedures are multivariate extensions or adaptations of procedures used for testing univariate normality. Results of several power studies are summarized, and an extensive bibliography of literature pertaining to testing for multivariate normality is provided.     

Percentage Points of the Multivariate t Distribution

The known methods for computing percentage points of multivariate t distributions are reviewed. We believe that this review will serve as an important reference and encourage further research activities in the area.     

The t Copula and Related Copulas

The t copula and its properties are described with a focus on issues related to the dependence of extreme values. The Gaussian mixture representation of a multivariate t distribution is used as a starting point to construct two new copulas, the skewed t copula and the grouped t copula, which allow more heterogeneity in the modelling of dependent observations. Extreme value considerations are used to derive two further new copulas: the t extreme value copula is the limiting copula of componentwise maxima of t distributed random vectors; the t lower tail copula is the limiting copula of bivariate observations from a t distribution that are conditioned to lie below some joint threshold that is progressively lowered. Both these copulas may be approximated for practical purposes by simpler, better-known copulas, these being the Gumbel and Clayton copulas respectively.     

Increasing failure rate and decreasing reversed hazard rate properties of the minimum and maximum of multivariate distributions with log-concave densities

For a multivariate random vector X = (X ₁,...,X _n ) with a log-concave density function, it is shown that the minimum min{X ₁,...,X _n } has an increasing failure rate, and the maximum max{X ₁,...,X _n } has a decreasing reversed hazard rate. As an immediate consequence, the result of Gupta and Gupta (in Metrika 53:39–49, 2001) on the multivariate normal distribution is obtained. One error in Gupta and Gupta method is also pointed out.      

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     

The association between two random elements: A complete characterization and odds ratio models

For random elements X and Y (e.g. vectors) a complete characterization of their association is given in terms of an odds ratio function. The main result establishes for any odds ratio function and any pre-specified marginals the unique existence of a corresponding joint distribution (the joint density is obtained as a limit of an iterative procedure of marginal fittings). Restricting only the odds ratio function but not the marginals leads to semi-parmetric association models for which statistical inference is available for samples drawn conditionally on either X or Y. Log-bilinear association models for random vectors X and Y are introduced which generalize standard (regression) models by removing restrictions on the marginals. In particular, the logistic regression model is recognized as a log-bilinear association model. And the joint distribution of X and Y is shown to be multivariate normal if and only if both marginals are normal and the association is log-bilinear.Acknowledgements The author thanks both referees for their helpful comments which improved the first draft of the paper.     

The emperor's new clothes: a critique of the multivariate t regression model

Zellner (1976) proposed a regression model in which the data vector (or the error vector) is represented as a realization from the multivariate Student t distribution. This model has attracted considerable attention because it seems to broaden the usual Gaussian assumption to allow for heavier-tailed error distributions. A number of results in the literature indicate that the standard inference procedures for the Gaussian model remain appropriate under the broader distributional assumption, leading to claims of robustness of the standard methods. We show that, although mathematically the two models are different, for purposes of statistical inference they are indistinguishable. The empirical implications of the multivariate t model are precisely the same as those of the Gaussian model. Hence the suggestion of a broader distributional representation of the data is spurious, and the claims of robustness are misleading. These conclusions are reached from both frequentist and Bayesian perspectives.     

A new continuous distribution and two new families of distributions based on the exponential

Recent work on social status led to derivation of a new continuous distribution based on the exponential. The new variate, termed the ring(2)-exponential, in turn leads to derivation of two closely related new families of continuous distributions, the mirror-exponential and the ring-exponential. Both the standard exponential and the ring(2)-exponential are special cases of both the new families. In this paper, we first focus on the ring(2)-exponential, describing its derivation and examining its properties, and next introduce the two new families, describing their derivation and initiating exploration of their properties. The mirror-exponential arises naturally in the study of status; the ring-exponential arises from the mathematical structure of the ring(2)-exponential. Both have the potential for broad application in diverse contexts across science and engineering. Within sociobehavioral contexts, the new mirror-exponential may have application to the problem of approximating the form and inequality of the wage distribution.     

Closed-form likelihood expansions for multivariate time-inhomogeneous diffusions

The aim of this paper is to find approximate log-transition density functions for multivariate time-inhomogeneous diffusions in closed form. There are many empirical evidences supporting that the data generating process governing dynamics of many economics variables might vary over time because of economic climate changes or time effects. One possible way to explain the time-dependent dynamics of state variables is to model the drift or volatility terms as functions of time t$t$ as well as state variables. A way to find closed-form likelihood expansion for a multivariate time-homogeneous diffusion has been developed by Aït-Sahalia (2008). This research is built on his work and extends his results to time-inhomogeneous cases. We conduct Monte Carlo simulation studies to examine performance of the approximate transition density function when it is used to obtain ML estimates. The results reveal that our method yields a very accurate approximate likelihood function, which can be a good candidate when the true likelihood function is unavailable as is often the case.     

Probability Integrals of the Multivariate t Distribution

Results on probability integrals of multivariate t distributions are reviewed. The results discussed include: Dunnett and Sobel's probability integrals, Gupta and Sobel's probability integrals, John's probability integrals, Amos and Bulgren's probability integrals, Steffens' non‐central probabilities, Dutt's probability integrals, Amos' probability integral, Fujikoshi's probability integrals, probabilities of cone, probabilities of convex polyhedra, probabilities of linear inequalities, maximum probability content, and Monte Carlo evaluation.     

Forecasting correlated time series with exponential smoothing models

This paper presents the Bayesian analysis of a general multivariate exponential smoothing model that allows us to forecast time series jointly, subject to correlated random disturbances. The general multivariate model, which can be formulated as a seemingly unrelated regression model, includes the previously studied homogeneous multivariate Holt-Winters’ model as a special case when all of the univariate series share a common structure. MCMC simulation techniques are required in order to approach the non-analytically tractable posterior distribution of the model parameters. The predictive distribution is then estimated using Monte Carlo integration. A Bayesian model selection criterion is introduced into the forecasting scheme for selecting the most adequate multivariate model for describing the behaviour of the time series under study. The forecasting performance of this procedure is tested using some real examples.     

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.     

Selection of the best population: an information theoretic approach

This paper is devoted to the statistical problem of ranking and selection populations by using the subset selection formulation. The interest is focused (i) on the selection of the best population among k independent populations and (ii) on the selection of the best population, which is closest to an additional standard or control population. With respect to the first problem the populations are ranked in terms of entropies of their distributions and the population whose distribution has maximum entropy is selected. For the second problem the populations are ranked in terms of divergences between their distributions and the distribution of the standard or control population and the population with the minimum divergence is selected. In each case the populations are assumed to have general parametric densities satisfying the classical regularity conditions of asymptotic statistic. Large sample properties of the estimators of entropies and divergences of the populations will be studied and used in order to determine the probabilities of correct selection of the proposed asymptotic selection rules. Illustrative examples, including a numerical example using real medical data appeared in the literature, will be given for multivariate homoscedastic normal populations and populations described by the regular exponential family of distributions. Received December 2001     

Multivariate semi-nonparametric distributions with dynamic conditional correlations

This paper generalizes the Dynamic Conditional Correlation (DCC) model of Engle (2002), incorporating a flexible non-Gaussian distribution based on Gram-Charlier expansions. The resulting semi-nonparametric-DCC (SNP-DCC) model allows estimation in two stages and deals with the negativity problem which is inherent in truncated SNP densities. We test the performance of a SNP-DCC model with respect to the (Gaussian)-DCC through an empirical application of density forecasting for portfolio returns. Our results show that the proposed multivariate model provides a better in-sample fit and forecast of the portfolio returns distribution, and thus is useful for financial risk forecasting and evaluation.     

Exact Skewness–Kurtosis Tests for Multivariate Normality and Goodness‐of‐Fit in Multivariate Regressions with Application to Asset Pricing Models*

We study the problem of testing the error distribution in a multivariate linear regression (MLR) model. The tests are functions of appropriately standardized multivariate least squares residuals whose distribution is invariant to the unknown cross‐equation error covariance matrix. Empirical multivariate skewness and kurtosis criteria are then compared with a simulation‐based estimate of their expected value under the hypothesized distribution. Special cases considered include testing multivariate normal and stable error distributions. In the Gaussian case, finite‐sample versions of the standard multivariate skewness and kurtosis tests are derived. To do this, we exploit simple, double and multi‐stage Monte Carlo test methods. For non‐Gaussian distribution families involving nuisance parameters, confidence sets are derived for the nuisance parameters and the error distribution. The tests are applied to an asset pricing model with observable risk‐free rates, using monthly returns on New York Stock Exchange (NYSE) portfolios over 5‐year subperiods from 1926 to 1995.     

Global and regional integration of the Middle East and North African (MENA) stock markets

We investigate financial integration of MENA region to facilitate a more in-depth exploration of the structure of interdependence and transmission mechanism of stock returns and volatility between MENA and world stock markets. The EGARCH-M models with a generalized error distribution are employed to consider both leverage effect of negative shocks and leptokurtosis prevalent in the MENA stock markets. The estimation results of multivariate AR-GARCH models indicate that there are large and predominantly positive volatility spillovers and volatility persistence in conditional volatility between MENA and world stock markets. Own-volatility spillovers are generally higher than cross-volatility spillovers for all the markets.     

Parametric bootstrap tests for continuous and discrete distributions

Statistical procedures based on the estimated empirical process are well known for testing goodness of fit to parametric distribution families. These methods usually are not distribution free, so that the asymptotic critical values of test statistics depend on unknown parameters. This difficulty may be overcome by the utilization of parametric bootstrap procedures. The aim of this paper is to prove a weak approximation theorem for the bootstrapped estimated empirical process under very general conditions, which allow both the most important continuous and discrete distribution families, along with most parameter estimation methods. The emphasis is on families of discrete distributions, and simulation results for families of negative binomial distributions are also presented.     

Approximation for the cumulative and inverse gamma distribution

The gamma distribution function can be expressed in terms of the Normal distribution and density functions with sufficient accuracy for most practical purposes.
The distribution function for the density x^Λ-1e^-x/μ/μ^ΛΓ(A) on 0 -R(Λ){(1 + 1/1 2Λ) φ(z) + 11 -z/4Λ^1/2+2(z²+ 2)/45Λ] φ(z) /3 Λ^1/2} where φ(z)≅1/[1 +_e^{-2z(√2/π+z2 /28)}] and φ(z) = e^{-z2 /2}/√2π are the Normal distribution and density functions, y is the appropriate root of y-y²/6+y³/36-y⁴/270= In (x/Λμ), z= Λ^1/2 y, and R( Λ) is the remainder term in Stirling's approximation for In Γ(Λ).