Multivariate Pareto Distributions: Inference and Financial Applications

Univariate Pareto distributions are extensively studied. In this article, we propose a Bayesian inference methodology in the context of multivariate Pareto distributions of the second kind (Mardia's type). Computational techniques organized around Gibbs sampling with data augmentation are proposed to implement Bayesian inference in practice. The new methods are shown to work well in artificial examples involving a trivariate distribution, and to an empirical application involving daily exchange rate data for four major currencies.     

Properties and Inference for a New Class of Skew-t Distributions

In this article we introduce a new generalization of skew-t distributions, which contains the standard skew-t distribution, as a special case. This new class of distributions is an adequate model for modeling some dataset rather than the standard skew-t distributions. This kind of distributions can be represented as a scale-shape mixture of the extended skew-normal distributions. The main properties of this family of distributions are studied and a recurrence relation for the cumulative distribution functions (cdf) of them is presented. We derive the distribution of the order statistics from the trivariate exchangeable t-distribution in terms of our distribution and then an exact expression for the cdf of order statistics is derived. Likelihood inference for this distribution is also examined. The method is illustrated with a numerical example via a simulation study.     

Parameter estimation of the generalized Pareto distribution—Part I

The generalized Pareto distribution (GPD) has been widely used in the extreme value framework. The success of the GPD when applied to real data sets depends substantially on the parameter estimation process. Several methods exist in the literature for estimating the GPD parameters. Mostly, the estimation is performed by maximum likelihood (ML). Alternatively, the probability weighted moments (PWM) and the method of moments (MOM) are often used, especially when the sample sizes are small. Although these three approaches are the most common and quite useful in many situations, their extensive use is also due to the lack of knowledge about other estimation methods. Actually, many other methods, besides the ones mentioned above, exist in the extreme value and hydrological literatures and as such are not widely known to practitioners in other areas. This paper is the first one of two papers that aim to fill in this gap. We shall extensively review some of the methods used for estimating the GPD parameters, focusing on those that can be applied in practical situations in a quite simple and straightforward manner.     

On the Parameter Estimation of the Asymmetric Multivariate Laplace Distribution

This article examines a family of three-parameter multivariate Laplace distributions ML _p (a, μ, Σ) which is closed under constant shifts. Parameter vectors a and μ are called shift and shape parameter, respectively, positive definite p × p-matrix Σ is a scale parameter. The first three moments are derived and used for estimating the parameters. The behavior of the obtained estimates is explored in a simulation experiment.     

On the mixtures of Weibull and Pareto (IV) distribution: An alternative to Pareto distribution

Finite mixture models have provided a reasonable tool to model various types of observed phenomena, specially those which are random in nature. In this article, a finite mixture of Weibull and Pareto (IV) distribution is considered and studied. Some structural properties of the resulting model are discussed including estimation of the model parameters via expectation maximization (EM) algorithm. A real-life data application exhibits the fact that in certain situations, this mixture model might be a better alternative than the rival popular models.     

Inference for a Hidden Truncated (Both-Sided) Bivariate Pareto (II) Distribution

The Pareto distribution is a simple model for non negative data with a power law probability tail. Income and wealth data are typically modeled using some variant of the classical Pareto distribution. In practice, it is frequently likely that the observed data have been truncated with respect to some unobserved covariable. In this paper, a hidden truncation formulation of this scenario is proposed and analyzed. A bivariate Pareto (II) distribution is assumed for the variable of interest and the unobserved covariable. Distributional properties of the resulting model are investigated. A variety of parameter estimation strategies (under the classical set up) are investigated.     

Entropy-based goodness-of-fit tests for the Pareto I distribution

In this article, having observed the generalized order statistics in a sample, we construct a test for the hypothesis that the underlying distribution is the Pareto I distribution. The Shannon entropy of generalized order statistics is used to test the null hypothesis.     

On the Estimation of Parameters of Kumaraswamy-G Distributions

Cordeiro and de Castro proposed a new family of generalized distributions based on the Kumaraswamy distribution (denoted as Kw-G). Nadarajah et al. showed that the density function of the new family of distributions can be expressed as a linear combination of the density of exponentiated family of distributions. They derived some properties of Kw-G distributions and discussed estimation of parameters using the maximum likelihood (ML) method. Cheng and Amin and Ranneby introduced a new method of estimating parameters based on Kullback–Leibler divergence (the maximum spacing (MSP) method). In this article, the estimates of parameters of Kw-G distributions are obtained using the MSP method. For some special Kw-G distributions, the new estimators are compared with ML estimators. It is shown by simulations and a real data application that MSP estimators have better properties than ML estimators.     

Parametric Estimation Procedures in Multivariate Generalized Pareto Models

Abstract.  Modelling the tails of a multivariate distribution can be reasonably done by multivariate generalized Pareto distributions (GPDs). We present several methods of parametric estimation in these models, which use decompositions of the corresponding random vectors with the help of different versions of Pickands coordinates. The estimators are compared to each other with simulated data sets. To show the practical value of the methods, they are applied to a real hydrological data set.     

Statistical Inference for Damaged Poisson Distribution

Abstract

For non-negative integer-valued random variables, the concept of “damaged” observations was introduced, for the first time, by Rao and Rubin [Rao, C. R., Rubin, H. (1964). On a characterization of the Poisson distribution. Sankhya 26:295–298] in 1964 on a paper concerning the characterization of Poisson distribution. In 1965, Rao [Rao, C. R. (1965). On discrete distribution arising out of methods of ascertainment. Sankhya Ser. A. 27:311–324] discusses some results related with inferences for parameters of a Poisson Model when it has occurred partial destruction of observations. A random variable is said to be damaged if it is unobservable, due to a damage mechanism which randomly reduces its magnitude. In subsequent years, considerable attention has been given to characterizations of distributions of such random variables that satisfy the “Rao–Rubin” condition. This article presents some inference aspects of a damaged Poisson distribution, under reasonable assumption that, when an observation on the random variable is made, it is also possible to determine whether or not some damage has occurred. In other words, we do not know how many items are damaged, but we can identify the existence of damage. Particularly it is illustrated the situation in which it is possible to identify the occurrence of some damage although it is not possible to determine the amount of items damaged. Maximum likelihood estimators of the underlying parameters and their asymptotic covariance matrix are obtained. Convergence of the estimates of parameters to the asymptotic values are studied through Monte Carlo simulations.     

Estimation of the Reliability of a Stress-Strength System from Power Lindley Distributions

In this paper, we are interested in the estimation of the reliability parameter R = P(X > Y) where X, a component strength, and Y, a component stress, are independent power Lindley random variables. The point and interval estimation of R, based on maximum likelihood, nonparametric and parametric bootstrap methods, are developed. The performance of the point estimate and confidence interval of R under the considered estimation methods is studied through extensive simulation. A numerical example, based on a real data, is presented to illustrate the proposed procedure.     

Adjusted Likelihood Inference in an Elliptical Multivariate Errors-in-Variables Model

In this paper, we obtain an adjusted version of the likelihood ratio (LR) test for errors-in-variables multivariate linear regression models. The error terms are allowed to follow a multivariate distribution in the class of the elliptical distributions, which has the multivariate normal distribution as a special case. We derive a modified LR statistic that follows a chi-squared distribution with a high degree of accuracy. Our results generalize those in Melo and Ferrari (Advances in Statistical Analysis, 2010, 94, pp. 75–87) by allowing the parameter of interest to be vector-valued in the multivariate errors-in-variables model. We report a simulation study which shows that the proposed test displays superior finite sample behavior relative to the standard LR test.     

Bayesian Statistical Inference for Weighted Exponential Distribution

Exponential distribution has an extensive application in reliability. Introducing shape parameter to this distribution have produced various distribution functions. In their study in 2009, Gupta and Kundu brought another distribution function using Azzalini's method, which is applicable in reliability and named as weighted exponential (WE) distribution. The parameters of this distribution function have been recently estimated by the above two authors in classical statistics. In this paper, Bayesian estimates of the parameters are derived. To achieve this purpose we use Lindley's approximation method for the integrals that cannot be solved in closed form. Furthermore, a Gibbs sampling procedure is used to draw Markov chain Monte Carlo samples from the posterior distribution indirectly and then the Bayes estimates of parameters are derived. The estimation of reliability and hazard functions are also discussed. At the end of the paper, some comparisons between classical and Bayesian estimation methods are studied by using Monte Carlo simulation study. The simulation study incorporates complete and Type-II censored samples.     

Statistical Inference for a Relative Risk Measure

ABSTRACT

For monitoring systemic risk from regulators’ point of view, this article proposes a relative risk measure, which is sensitive to the market comovement. The asymptotic normality of a nonparametric estimator and its smoothed version is established when the observations are independent. To effectively construct an interval without complicated asymptotic variance estimation, a jackknife empirical likelihood inference procedure based on the smoothed nonparametric estimation is provided with a Wilks type of result in case of independent observations. When data follow from AR-GARCH models, the relative risk measure with respect to the errors becomes useful and so we propose a corresponding nonparametric estimator. A simulation study and real-life data analysis show that the proposed relative risk measure is useful in monitoring systemic risk.     

Inference based on progressively censored sample from Pareto population

Abstract

In this paper, we assume that the lifetimes have a two-parameter Pareto distribution and discuss some results of progressive Type-II censored sample. We obtain maximum likelihood estimators and Bayes estimators of the unknown parameters under squared error loss and a precautionary loss functions in progressively Type-II censored sample. Robust Bayes estimation of unknown parameters over three different classes of priors under progressively Type-II censored sample, squared error loss, and precautionary loss functions are obtained. We discuss estimation of unknown parameters on competing risks progressive Type-II censoring. Finally, we consider the problem of estimating the common scale parameter of two Pareto distributions when samples are progressively Type-II censored.     

Multivariate extension of modified Sarhan-Balakrishnan bivariate distribution

Recently Kundu and Gupta [2010, Modified Sarhan-Balakrishnan singular bivariate distribution, Journal of Statistical Planning and Inference, 140, 526-538] introduced the modified Sarhan-Balakrishnan bivariate distribution and established its several properties. In this paper we provide a multivariate extension of the modified Sarhan-Balakrishnan bivariate distribution. It is a distribution with a singular part. Different ageing and dependence properties of the proposed multivariate distribution have been established. The moment generating function, the product moments can be obtained in terms of infinite series. The multivariate hazard rate has been obtained. We provide the EM algorithm to compute the maximum likelihood estimators and an illustrative example is performed to see the effectiveness of the proposed method.     

Trace Class Markov Chains for Bayesian Inference with Generalized Double Pareto Shrinkage Priors

Bayesian shrinkage methods have generated a lot of interest in recent years, especially in the context of high‐dimensional linear regression. In recent work, a Bayesian shrinkage approach using generalized double Pareto priors has been proposed. Several useful properties of this approach, including the derivation of a tractable three‐block Gibbs sampler to sample from the resulting posterior density, have been established. We show that the Markov operator corresponding to this three‐block Gibbs sampler is not Hilbert–Schmidt. We propose a simpler two‐block Gibbs sampler and show that the corresponding Markov operator is trace class (and hence Hilbert–Schmidt). Establishing the trace class property for the proposed two‐block Gibbs sampler has several useful consequences. Firstly, it implies that the corresponding Markov chain is geometrically ergodic, thereby implying the existence of a Markov chain central limit theorem, which in turn enables computation of asymptotic standard errors for Markov chain‐based estimates of posterior quantities. Secondly, because the proposed Gibbs sampler uses two blocks, standard recipes in the literature can be used to construct a sandwich Markov chain (by inserting an appropriate extra step) to gain further efficiency and to achieve faster convergence. The trace class property for the two‐block sampler implies that the corresponding sandwich Markov chain is also trace class and thereby geometrically ergodic. Finally, it also guarantees that all eigenvalues of the sandwich chain are dominated by the corresponding eigenvalues of the Gibbs sampling chain (with at least one strict domination). Our results demonstrate that a minor change in the structure of a Markov chain can lead to fundamental changes in its theoretical properties. We illustrate the improvement in efficiency resulting from our proposed Markov chains using simulated and real examples.     

Some simple estimators for the two-parameter gamma distribution

ABSTRACT

In this paper, we propose two new simple estimation methods for the two-parameter gamma distribution. The first one is a modified version of the method of moments, whereas the second one makes use of some key properties of the distribution. We then derive the asymptotic distributions of these estimators. Also, bias-reduction methods are suggested to reduce the bias of these estimators. The performance of the estimators are evaluated through a Monte Carlo simulation study. The probability coverages of confidence intervals are also discussed. Finally, two examples are used to illustrate the proposed methods.     

Simulation-based Inference in a Zero-inflated Bernoulli Regression Model

The logistic regression model has become a standard tool to investigate the relationship between a binary outcome and a set of potential predictors. When analyzing binary data, it often arises that the observed proportion of zeros is greater than expected under the postulated logistic model. Zero-inflated binomial (ZIB) models have been developed to fit binary data that contain too many zeros. Maximum likelihood estimators in these models have been proposed and their asymptotic properties established. Several aspects of ZIB models still deserve attention however, such as the estimation of odds-ratios and event probabilities. In this article, we propose estimators of these quantities and we investigate their properties both theoretically and via simulations. Based on these results, we provide recommendations about the range of conditions (minimum sample size, maximum proportion of zeros in excess) under which a reliable statistical inference on the odds-ratios and event probabilities can be obtained in a ZIB regression model. A real-data example illustrates the proposed estimators.