Bayesian analysis for the two-parameter Pareto distribution based on record values and times

Doostparast and Balakrishnan (Pareto record-based analysis, Statistics, under review) recently developed optimal confidence intervals as well as uniformly most powerful tests for one- and two-sided hypotheses concerning shape and scale parameters, for the two-parameter Pareto distribution based on record data. In this paper, on the basis of record values and inter-record times from the two-parameter Pareto distribution, maximum-likelihood and Bayes estimators as well as credible regions are developed for the two parameters of the Pareto distribution. For illustrative purposes, a data set on annual wages of a sample of production-line workers in a large industrial firm is analysed using the proposed procedures.     

Statistical Inference for a New Class of Multivariate Pareto Distributions

Various solutions to the parameter estimation problem of a recently introduced multivariate Pareto distribution are developed and exemplified numerically. Namely, a density of the aforementioned multivariate Pareto distribution with respect to a dominating measure, rather than the corresponding Lebesgue measure, is specified and then employed to investigate the maximum likelihood estimation (MLE) approach. Also, in an attempt to fully enjoy the common shock origins of the multivariate model of interest, an adapted variant of the expectation-maximization (EM) algorithm is formulated and studied. The method of moments is discussed as a convenient way to obtain starting values for the numerical optimization procedures associated with the MLE and EM methods.     

Order Sampling Design with Prescribed Inclusion Probabilities

Order sampling with fixed distribution shape is a class of sampling schemes with inclusion probabilities approximately proportional to given size measures. In a recent article, methods were provided to compute the exact first and second order inclusion probabilities numerically when the distribution shape is of the Pareto type. In the same article, procedures were also provided for this case to adjust the parameters to get predetermined inclusion probabilities. In this paper we prove the existence and uniqueness of a solution for the latter problem, in general for any order sampling of fixed distribution shape.     

A Sieve model for extreme values

From the class of extreme value distributions, we focus on the set of heavy-tailed distributions which produce low-frequency, high-cost events. The regular Pareto distribution is the basic model of choice, being the simplest heavy-tailed distribution. Real data suggest that modifications of the Pareto distribution may be a better fit; an alternative model is the truncated Pareto distribution (TPD). For further study, this paper proposed a TPD Sieve class of distributions. The properties and estimation on the Sieve class are also discussed. We fit the models to the largest Black Sea bass caught in Buzzard's Bay, MA, USA and the costliest Atlantic hurricanes from 1900 to 2005. Using measures of model adequacy, the TPD Sieve model is generally found to be the best-fitting model.     

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.     

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.     

CONDITIONAL INFERENCE PROCEDURES FOR THE PARETO AND POWER FUNCTION DISTRIBUTIONS BASED ON TYPE-II RIGHT CENSORED SAMPLES

In this paper we consider conditional inference procedures for the Pareto and power function distributions. We develop procedures for obtaining confidence intervals for the location and scale parameters as well as upper and lower n probability tolerance intervals for a proportion g, given a Type-II right censored sample from the corresponding distribution. The intervals are exact, and are obtained by conditioning on the observed values of the ancillary statistics. Since, for each distribution, the procedures assume that a shape parameter x is known, a sensitivity analysis is also carried out to see how the procedures are affected by changes in x.     

Tests for the Validity of the Assumption that the Underlying Distribution of Life is Pareto

This article considers the problem of testing the validity of the assumption that the underlying distribution of life is Pareto. For complete and censored samples, the relationship between the Pareto and the exponential distributions could be of vital importance to test for the validity of this assumption. For grouped uncensored data the classical Pearson χ² test based on the multinomial model can be used. Attention is confined in this article to handle grouped data with withdrawals within intervals. Graphical as well as analytical procedures will be presented. Maximum likelihood estimators for the parameters of the Pareto distribution based on grouped data will be derived.     

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.     

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.     

Bayesian analysis for heteroscedastic normal-Pareto mixture model with application

ABSTRACT

In this article, we propose a new distribution by mixing normal and Pareto distributions, and the new distribution provides an unusual hazard function. We model the mean and the variance with covariates for heterogeneity. Estimation of the parameters is obtained by the Bayesian method using Markov Chain Monte Carlo (MCMC) algorithms. Proposal distribution in MCMC is proposed with a defined working variable related to the observations. Through the simulation, the method shows a dependable performance of the model. We demonstrate through establishing model under a real dataset that the proposed model and method can be more suitable than the previous report.     

Optimal progressive group censoring scheme under cost considerations for pareto distribution

In this article, optimal design under the restriction of pre-determined budget of experiment is developed for the Pareto distribution when the life test is progressively group censored. We use the maximum-likelihood method to obtain the point estimator of the Pareto parameter. We propose two approaches to decide the number of test units, the number of inspections, and the length of inspection interval under limited budget such that the asymptotic variance of estimator of Pareto parameter is minimum. A numerical example is given to illustrate the proposed method. Some sensitivity analysis is also studied.     

Estimation of the Lognormal-Pareto Distribution Using Probability Weighted Moments and Maximum Likelihood

This article deals with the estimation of the lognormal-Pareto and the lognormal-generalized Pareto distributions, for which a general result concerning asymptotic optimality of maximum likelihood estimation cannot be proved. We develop a method based on probability weighted moments, showing that it can be applied straightforwardly to the first distribution only. In the lognormal-generalized Pareto case, we propose a mixed approach combining maximum likelihood and probability weighted moments. Extensive simulations analyze the relative efficiencies of the methods in various setups. Finally, the techniques are applied to two real datasets in the actuarial and operational risk management fields.     

On the identification of extreme outliers and dragon-kings mechanisms in the upper tail of income distribution

The presence of extreme outliers in the upper tail data of income distribution affects the Pareto tail modeling. A simulation study is carried out to compare the performance of three types of boxplot in the detection of extreme outliers for Pareto data, including standard boxplot, adjusted boxplot and generalized boxplot. It is found that the generalized boxplot is the best method for determining extreme outliers for Pareto distributed data. For the application, the generalized boxplot is utilized for determining the exreme outliers in the upper tail of Malaysian income distribution. In addition, for this data set, the confidence interval method is applied for examining the presence of dragon-kings, extreme outliers which are beyond the Pareto or power-laws distribution.     

Projection Mean–Variance Bounds on Expectations of kth Record Values from Restricted Families

We present sharp mean–variance bounds for expectations of kth record values based on distributions coming from restricted families of distributions. These families are defined in terms of convex or star ordering with respect to generalized Pareto distribution. The bounds for expectations of kth record values from DD, DFR, DDA, and DFRA families are special cases of our results. The bounds are derived by application of the projection method.     

Bayesian prediction of order statistics with fixed and random sample sizes based on k-record values from Pareto distribution

In this paper, the two-parameter Pareto distribution is considered and the problem of prediction of order statistics from a future sample and that of its geometric mean are discussed. The Bayesian approach is applied to construct predictors based on observed k-record values for the cases when the future sample size is fixed and when it is random. Several Bayesian prediction intervals are derived. Finally, the results of a simulation study and a numerical example are presented for illustrating all the inferential procedures developed here.     

Pareto Regression: A Bayesian Analysis

Abstract

In this article, a new model is presented that is based on the Pareto distribution of the second kind, when the location parameter depends on covariates as well as unobserved heterogeneity. Bayesian analysis of the model can be performed using Markov Chain Monte Carlo techniques. The new procedures are illustrated in the context of artificial data as well as international output data.     

General location transform of the order statistics from the exponential,pareto and weibull,with application to maximum likelihood estimation

Based one some common distribution properties of the order statistics and the transformation theory by Efron(1982), we determine unified explicit general location transformations, which map the distributions of the order statistics from the Exponential, Pareto and Weibull to a standard normal distribution. This result is used to derive analytical formulas for the maximum likelihood estimators of the shape parameter of these distributions of order statistics. The presented exact method is applied to catastrophe earthquake life reinsurance.     

Bayesian estimation and prediction based on generalized Type-II hybrid censored sample

ABSTRACT

In this paper, we consider a general form for the underlying distribution and a general conjugate prior, and develop a general procedure for deriving the maximum likelihood and Bayesian estimators based on an observed generalized Type-II hybrid censored sample. The problems of predicting the future order statistics from the same sample and that from a future sample are also discussed from a Bayesian viewpoint. For the illustration of the developed results, the exponential and Pareto distributions are used as examples. Finally, two numerical examples are presented for illustrating all the inferential procedures developed here.     

Inference and prediction for pareto progressively censored data

Progressively censored data from a classical Pareto distribution are to be used to make inferences about its shape and precision parameters and the reliability function. An approximation form due to Tierney and Kadane (1986) is used for obtaining the Bayes estimates. Bayesian prediction of further observations from this distribution is also considered. When the Bayesian approach is concerned, conjugate priors for either the one or the two parameters cases are considered. To illustrate the given procedures, a numerical example and a simulation study are given.