The Distribution of Sums of Certain I.I.D. Pareto Variates

ABSTRACT

Though the Pareto distribution is important to actuaries and economists, an exact expression for the distribution of the sum of n i.i.d. Pareto variates has been difficult to obtain in general. This article considers Pareto random variables with common probability density function (pdf) f(x) = (α/β) (1 + x/β)^α+1 for x > 0, where α = 1,2,… and β > 0 is a scale parameter. To date, explicit expressions are known only for a few special cases: (i) α = 1 and n = 1,2,3; (ii) 0 < α < 1 and n = 1,2,…; and (iii) 1 < α < 2 and n = 1,2,…. New expressions are provided for the more general case where β > 0, and α and n are positive integers. Laplace transforms and generalized exponential integrals are used to derive these expressions, which involve integrals of real valued functions on the positive real line. An important attribute of these expressions is that the integrands involved are non oscillating.     

Limit Distributions of Random Sums of Z+-Valued Random Variables

In this article, we formulate a transfer theorem in terms of probability generating functions and discuss two approaches to limit distributions of random sums of Z ₊-valued random variables. We then develop Z ₊-valued N-ID and ?-ID laws.     

Shrinkage Testimators for the Shape Parameter of Pareto Distribution Using LINEX Loss Function

In this article, shrinkage testimators for the shape parameter of a Pareto distribution are considered, when its prior guess value is available. The choices of shrinkage factor are also suggested. The proposed testimators are compared with the minimum risk estimator among the class of unbiased estimators with the LINEX loss function.     

A Simple Estimator for the Shape Parameter of the Pareto Distribution with Economics and Medical Applications

In the present paper, an estimator of the shape parameter of the Pareto failure model is presented using grouped data. This estimator is based on obtaining the parameter in terms of the hazard rate, then replacing the unknown hazard rate by a grouped data estimator available in the literature. Death records are given as a numerical illustration in the medical context. The relation between the hazard rate and the income elasticity is derived. This relation allows the presentation of the same estimator in terms of the income elasticity so that it could be used in an economic context. Two illustrations are presented using income data. Simulated data are generated to compare the estimator with the maximum likelihood estimator.     

A Note on the Complementary Mixture Pareto II Distribution

We introduce a new survival distribution, of Pareto type, that arises from a cure-mixture frailty model. We describe its properties and demonstrate connections with familiar distributions including the Pareto and exponential. We derive its characteristic function and moments.     

Estimation of a Parameter of Bivariate Pareto Distribution by Ranked Set Sampling

Ranked set sampling is applicable whenever ranking of a set of sampling units can be done easily by a judgement method or based on the measurement of an auxiliary variable on the units selected. In this work, we derive different estimators of a parameter associated with the distribution of the study variate Y, based on a ranked-set sample obtained by using an auxiliary variable X correlated with Y for ranking the sample units, when (X, Y) follows a bivariate Pareto distribution. Efficiency comparisons among these estimators are also made. Real-life data have been used to illustrate the application of the results obtained.     

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.     

An Improvement of Poisson Approximation for Sums of Dependent Bernoulli Random Variables

This article uses the Stein-Chen method to obtain new non uniform bounds on the error of the cumulative distribution function of sums of dependent Bernoulli random variables and the Poisson cumulative distribution function. The bounds obtained in the present study are sharper than those reported in Teerapabolarn and Neammanee (2006 Teerapabolarn, K., Neammanee, K. (2006). Poisson approximation for sums of dependent Bernoulli random variables. Acta Mathematica Academiae Paedagogicae Nyiregyhaziensis. 22:87–99. [Google Scholar]). Examples are provided to illustrate applications of the obtained results.     

Parameter Estimation for the Tail Distribution of a Random Sequence

The tail Y_t = X_t – u of a random sequence {X_t ∈ , t ∈ } with identically distributed X_t is approximated by the generalized Pareto distribution according to the extreme value theory, wherein Y_t occurs in clusters because of the dependence in the random sequence. Nevertheless, the parameters of the generalized Pareto distribution are estimated by the same methods as in the case of independent and identically distributed Y_t, provided that there is independence between the clusters of Y_t. The estimation variances and confidence intervals can be estimated by the jackknife method. The approaches are theoretically discussed and verified by extensive numerical researches.     

Density of the Ratio of Two Normal Random Variables and Applications

In reply to a question raised in the literature, and to settle an argument debated in the last decades, we give the exact closed form expression of the density of X/Y, where X and Y are normal random variables, in terms of Hermite and confluent hypergeometric functions. All cases will be considered: standardized and nonstandardized variables, independent or correlated variables. Examples in applied disciplines are presented, and generalizations to ratios of variables from scale mixtures of bivariate normal distributions show the potential of further new applications in applied statistics and operations research.     

Parameter estimation of the generalized Pareto distribution—Part II

This is the second part of a paper which focuses on reviewing methods for estimating the parameters of the generalized Pareto distribution (GPD). The GPD is a very important distribution in the extreme value context. It is commonly used for modeling the observations that exceed very high thresholds. The ultimate success of the GPD in applications evidently depends on the parameter estimation process. Quite a few methods exist in the literature for estimating the GPD parameters. Estimation procedures, such as the maximum likelihood (ML), the method of moments (MOM) and the probability weighted moments (PWM) method were described in Part I of the paper. We shall continue to review methods for estimating the GPD parameters, in particular methods that are robust and procedures that use the Bayesian methodology. As in Part I, we shall focus on those that are relatively simple and straightforward to be applied to real world data.     

Approximate studentization for Pareto distribution with application to censored data

The use of the Pareto distribution as a model for various socio-economic phenomena dates back to the late nineteenth century. Recently, it has also been recognized as a useful model for the analysis of lifetime data. In this paper, we apply the approximate studentization method to obtain inference for the scale parameter of the Pareto distribution, and also for the strong Pareto law. Moreover, we extend the method to construct prediction limits for thejth smallest future observation based on the firstk observed data.     

Estimation for the Parameters of the Pareto Distribution Under Progressive Censoring

In this article, progressive Type-II right censored sample from Pareto distribution is considered. Exact confidence region is derived for the parameters of the corresponding distribution under progressive censoring. Simulation study is performed to investigate the coverage probabilities of the proposed confidence region. Illustrative example is also given.     

Improved Estimation of the Pareto Location Parameter Under Squared Error Loss

Suppose two independent observations are drawn from Pareto distributions with known shape parameters and an order restriction on the unknown location parameters. An isotonic regression estimator of the smaller location parameter dominates a preferred marginal estimator under squared error loss, but fails to dominate under stochastic domination. The results expressed herein advance the theory of order restricted inference.     

Shared frailty models with baseline generalized Pareto distribution

Abstract

In this article, we have considered three different shared frailty models under the assumption of generalized Pareto Distribution as baseline distribution. Frailty models have been used in the survival analysis to account for the unobserved heterogeneity in an individual risks to disease and death. These three frailty models are with gamma frailty, inverse Gaussian frailty and positive stable frailty. Then we introduce the Bayesian estimation procedure using Markov chain Monte Carlo (MCMC) technique to estimate the parameters. We applied these three models to a kidney infection data and find the best fitted model for kidney infection data. We present a simulation study to compare true value of the parameters with the estimated values. Model comparison is made using Bayesian model selection criterion and a well-fitted model is suggested for the kidney infection data.     

Sum,product and ratio of Pareto and gamma variables

The exact distributions of X+Y, X Y and X/(X+Y) are studied when X and Y are independent Pareto and gamma random variables. Applications are discussed, to real problems in clinical trials, computer networks and economics.     

Estimation of the exponential Pareto II distribution parameters

In this article, we estimate the parameters of exponential Pareto II distribution by two new methods. The first one is based on the principle of maximum entropy (POME) and the second is by Kullback–Leibler divergence of survival function (KLS). Monte Carlo simulated data are used to evaluate these methods and compare them with the maximum likelihood method. Finally, we fit this distribution to a set of real data by estimation procedures.     

Simultaneous Confidence Intervals for the Parameters of Pareto Distribution under Progressive Censoring

In this article, we consider the progressive Type II right censored sample from Pareto distribution. We introduce a new approach for constructing the simultaneous confidence interval of the unknown parameters of this distribution under progressive censoring. A Monte Carlo study is also presented for illustration. It is shown that this confidence region has a smaller area than that introduced by Ku? and Kaya (2007 Ku? , C. , Kaya , M. F. ( 2007 ). Estimation for the parameters of the Pareto distribution under progressive censoring . Commun. Statist. Theor. Meth. 36 : 1359 – 1365 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]).     

The Joint Distribution of the Sum and the Maximum of IID Exponential Random Variables

The authors establish the joint distribution of the sum X and the maximum Y of IID exponential random variables. They derive exact formuli describing the random vector (X, Y), including its joint PDF, CDF, and other characteristics; marginal and conditional distributions; moments and related parameters; and stochastic representations leading to further properties of infinite divisibility and self-decomposability. The authors also discuss parameter estimation and include an example from climatology that illustrates the modeling potential of this new bivariate model.