The complementary beta distribution

The complementary beta distribution is proposed as a new distribution on the unit interval. It results from reversing the roles of the distribution and quantile functions of the beta distribution. It has some attractive properties that are complementary to those of the beta distribution. In particular, the complementary beta distribution is much more amenable than the beta distribution to exact computations involving expectations of order statistics, including L-moments. At least for a wide range of parameter values, complementary beta and beta distributions with parameters that are reciprocals of the other's parameters are good approximations to one another. We also note the position of the complementary beta distribution in a wider family of distributions defined through the same simple form for their quantile density functions.     

The model identification of beta distribution based on quantiles

This paper addresses the modelling of a beta distribution based on quantile pairs. The elicited beta distribution is proved to exist uniquely for any given two-quantile pairs. By combining the results of this study with the research documented in the literature, we concluded that for any given two-quantile pairs, there is one and only one beta distribution fits the two specified quantile pairs. The numerical results that verify this finding were obtained by comprehensively examining 1300 two-quantile pairs. The first four moments of the fitted beta distribution were also derived numerically using these two-quantile pairs. The variations in the first four moments, which reflect the perception of uncertainty, were plotted and were explored.     

The McDonald arcsine distribution: a new model to proportional data

We propose a new three-parameter continuous model called the McDonald arcsine distribution, which is a very competitive model to the beta, beta type I and Kumaraswamy distributions for modelling rates and proportions. We provide a mathematical treatment of the new distribution including explicit expressions for the density function, moments, generating and quantile functions, mean deviations, two probability measures based on the Bonferroni and Lorenz curves, Shannon entropy, Rényi entropy and cumulative residual entropy. Maximum likelihood is used to estimate the model parameters and the expected information matrix is determined. An application of the proposed model to real data shows that it can give consistently a better fit than other important statistical models.     

A simple formula based on quantiles for the moments of beta generalized distributions

In this article, we derive explicit expansions for the moments of beta generalized distributions from power series expansions for the quantile functions of the baseline distributions. We apply our formula to the beta normal, beta Student t, beta gamma and beta beta generalized distributions. We propose a simple way to express the quantile function of any beta generalized distribution as a power series expansion with known coefficients.     

On quantile residuals in beta regression

Beta regression is often used to model the relationship between a dependent variable that assumes values on the open interval (0, 1) and a set of predictor variables. An important challenge in beta regression is to find residuals whose distribution is well approximated by the standard normal distribution. Two previous works compared residuals in beta regression, but the authors did not include the quantile residual. Using Monte Carlo simulation techniques, this article studies the behavior of certain residuals in beta regression in several scenarios. Overall, the results suggest that the distribution of the quantile residual is better approximated by the standard normal distribution than that of the other residuals in most scenarios. Three applications illustrate the effectiveness of the quantile residual.     

Saddlepoint Approximations for the Difference of Order Statistics and Studentized Sample Quantiles

For a sample from a given distribution the difference of two order statistics and the Studentized quantile are statistics whose distribution is needed to obtain tests and confidence intervals for quantiles and quantile differences. This paper gives saddlepoint approximations for densities and saddlepoint approximations of the Lugannani–Rice form for tail probabilities of these statistics. The relative errors of the approximations are n ⁻¹ uniformly in a neighbourhood of the parameters and this uniformity is global if the densities are log-concave.     

A matching prior for extreme quantile estimation of the generalized Pareto distribution

Extreme quantile estimation plays an important role in risk management and environmental statistics among other applications. A popular method is the peaks-over-threshold (POT) model that approximate the distribution of excesses over a high threshold through generalized Pareto distribution (GPD). Motivated by a practical financial risk management problem, we look for an appropriate prior choice for Bayesian estimation of the GPD parameters that results in better quantile estimation. Specifically, we propose a noninformative matching prior for the parameters of a GPD so that a specific quantile of the Bayesian predictive distribution matches the true quantile in the sense of Datta et al. (2000).     

The McDonald extended distribution: properties and applications

We study a five-parameter lifetime distribution called the McDonald extended exponential model to generalize the exponential, generalized exponential, Kumaraswamy exponential and beta exponential distributions, among others. We obtain explicit expressions for the moments and incomplete moments, quantile and generating functions, mean deviations, Bonferroni and Lorenz curves and Gini concentration index. The method of maximum likelihood and a Bayesian procedure are adopted for estimating the model parameters. The applicability of the new model is illustrated by means of a real data set.     

Bootstrapping quantiles in a fixed design regression model with censored data

We consider the problem of estimating the quantiles of a distribution function in a fixed design regression model in which the observations are subject to random right censoring. The quantile estimator is defined via a conditional Kaplan-Meier type estimator for the distribution at a given design point. We establish an a.s. asymptotic representation for this quantile estimator, from which we obtain its asymptotic normality. Because a complicated estimation procedure is necessary for estimating the asymptotic bias and variance, we use a resampling procedure, which provides us, via an asymptotic representation for the bootstrapped estimator, with an alternative for the normal approximation.     

Quantile estimation in the proportional hazards model of random censorship

The asymptotic theory is given for quantile estimation in the proportional hazards model of random censorship. In this model, the tail of the censoring distribution function is some power of the tail of the survival distribution function. The quantile estimator is based on the maximum likelihood estimator for the survival time distribution, due to Abdushukurov, Cheng and Lin.     

On using percentiles to fit data by a johnson distribution

The well-known Johnson system of distributions was developed by N. L. Johnson (1949). Slifker and Shapiro (1980) presented a criterion for choosing a member from the three distributional classes (S_B,S_L, and S_v) in the Johnson system to fit a set of data. The criterion is based on the value of a quantile ratio which depends on a specified positive z value and the parameters of the distribution. In this paper, we present some properties of the quantile ratio for various distributions and for some selected z values. Some comments are made on using the criterion for selecting a Johnson distribution to fit empirical data.     

Bias reduction for high quantiles

High quantile estimation is of importance in risk management. For a heavy-tailed distribution, estimating a high quantile is done via estimating the tail index. Reducing the bias in a tail index estimator can be achieved by using either the same order or a larger order of number of the upper order statistics in comparison with the theoretical optimal one in the classical tail index estimator. For the second approach, one can either estimate all parameters simultaneously or estimate the first and second order parameters separately. Recently, the first method and the second method via external estimators for the second order parameter have been applied to reduce the bias in high quantile estimation. Theoretically, the second method obviously gives rise to a smaller order of asymptotic mean squared error than the first one. In this paper we study the second method with simultaneous estimation of all parameters for reducing bias in high quantile estimation.     

A fractional order statistic towards defining a smooth quantile function for discrete data

This work is motivated in part by a recent publication by Ma et al. (2011) who resolved the asymptotic non-normality problem of the classical sample quantiles for discrete data through defining a new mid-distribution based quantile function. This work is the motivation for defining a new and improved smooth population quantile function given discrete data. Our definition is based on the theory of fractional order statistics. The main advantage of our definition as compared to its competitors is the capability to distinguish the uth quantile across different discrete distributions over the whole interval, u∈(0,1). In addition, we define the corresponding estimator of the smooth population quantiles and demonstrate the convergence and asymptotic normal distribution of the corresponding sample quantiles. We verify our theoretical results through a Monte Carlo simulation, and illustrate the utilization of our quantile function in a Q-Q plot for discrete data.     

Using Weight Functions in Spatial Point Pattern Analysis with Application to Plant Ecology Data

This paper presents procedures for percentile estimation in the three-parameter inverse Gaussian (IG3) and the two-parameter inverse Gaussian (IG2) distributions. All procedures require first the estimation of distribution parameters and second the computation of the desired quantile at the estimated parameters. Parameter estimation is accomplished by maximum likelihood (ML) or a mixed moments (MXM) method. A Newton-Rahpson (NR) procedure is used for inverting the CDF. Simulation and asymptotic results are given for the resulting estimators. Two sets of real data are used to illustrate the procedures.     

On a log-symmetric quantile tobit model applied to female labor supply data

The study of female labor supply has been a topic of relevance in the economic literature. Generally, the data are left-censored and the classic tobit model has been extensively used in the modeling strategy. This model, however, assumes normality for the error distribution and is not recommended for data with positive skewness, heavy-tails and heteroscedasticity, as is the case of female labor supply data. Moreover, it is well-known that the quantile regression approach accounts for the influences of different quantiles in the estimated coefficients. We take all these features into account and propose a parametric quantile tobit regression model based on quantile log-symmetric distributions. The proposed method allows one to model data with positive skewness (which is not suitable for the classic tobit model), to study the influence of the quantiles of interest, and to account for heteroscedasticity. The model parameters are estimated by maximum likelihood and a Monte Carlo experiment is performed to evaluate alternative estimators. The new method is applied to two distinct female labor supply data sets. The results indicate that the log-symmetric quantile tobit model fits better the data than the classic tobit model.     

Automatic Bayesian quantile regression curve fitting

Quantile regression, including median regression, as a more completed statistical model than mean regression, is now well known with its wide spread applications. Bayesian inference on quantile regression or Bayesian quantile regression has attracted much interest recently. Most of the existing researches in Bayesian quantile regression focus on parametric quantile regression, though there are discussions on different ways of modeling the model error by a parametric distribution named asymmetric Laplace distribution or by a nonparametric alternative named scale mixture asymmetric Laplace distribution. This paper discusses Bayesian inference for nonparametric quantile regression. This general approach fits quantile regression curves using piecewise polynomial functions with an unknown number of knots at unknown locations, all treated as parameters to be inferred through reversible jump Markov chain Monte Carlo (RJMCMC) of Green (Biometrika 82:711–732, 1995). Instead of drawing samples from the posterior, we use regression quantiles to create Markov chains for the estimation of the quantile curves. We also use approximate Bayesian factor in the inference. This method extends the work in automatic Bayesian mean curve fitting to quantile regression. Numerical results show that this Bayesian quantile smoothing technique is competitive with quantile regression/smoothing splines of He and Ng (Comput. Stat. 14:315–337, 1999) and P-splines (penalized splines) of Eilers and de Menezes (Bioinformatics 21(7):1146–1153, 2005).     

Quantile Regression via the EM Algorithm

The three-parameter asymmetric Laplace distribution (ALD) has received increasing attention in the field of quantile regression due to an important feature between its location and asymmetric parameters. On the basis of the representation of the ALD as a normal-variance–mean mixture with an exponential mixing distribution, this article develops EM and generalized EM algorithms, respectively, for computing regression quantiles of linear and nonlinear regression models. It is interesting to show that the proposed EM algorithm and the MM (Majorization–Minimization) algorithm for quantile regressions are really the same in terms of computation, since the updating formula of them are the same. This provides a good example that connects the EM and MM algorithms. Simulation studies show that the EM algorithm can successfully recover the true parameters in quantile regressions.     

分位数回归技术综述

普通最小二乘回归建立了在自变量X=x下因变量Y的条件均值与X的关系的线性模型。而分位数回归（Quantile Regression）则利用自变量X和因变量y的条件分位数进行建模。与普通的均值回归相比,它能充分反映自变量X对于因变量y的分布的位置、刻度和形状的影响,有着十分广泛的应用,尤其是对于一些非常关注尾部特征的情况。文章介绍了分位数回归的概念以及分位数回归的估计、检验和拟合优度,回顾了分位数回归的发展过程以及其在一些经济研究领域中的应用,最后做了总结。     

On multivariate quantile regression analysis

This paper investigates the estimation of parameters in a multivariate quantile regression model when the investigator wants to evaluate the associated distribution function. It proposes a new directional quantile estimator with the following properties: (1) it applies to an arbitrary number of random variables; (2) it is equivalent to estimating the distribution function allowing for non-convex distribution contours; (3) it satisfies nice equivariance properties; (4) it has desirable statistical properties (i.e., consistency and asymptotic normality); and (5) its implementation involves a modest computational burden: our proposed estimator can be obtained by solving parametric linear programming problems. As such, this paper expands the range of applications of quantile estimation for multivariate regression models.     

函数型变量倾斜分位回归模型及其应用

本文考虑函数型数据的结构特征,针对两类函数型变量分位回归模型（函数型因变量对标量自变量和函数型因变量对函数型自变量）,基于函数型倾斜分位曲线的定义构建新型函数型倾斜分位回归模型。对于第二类模型,本文分别考虑样条基函数对模型系数展开和函数型主成分基函数对函数型自变量展开,得到倾斜分位回归模型的基本形式。参数估计采用成分梯度Boosting算法最小化加权非对称损失函数,提高计算效率。在理论上证明了倾斜分位回归模型的系数估计量均服从渐近正态分布。模拟和实证研究结果显示,倾斜分位回归模型比已有的逐点分位回归模型具有更好的拟合效果。根据积分均方预测误差准则,本文提出的模型有一致较好的预测能力。