Bayesian analysis of quantile regression for censored dynamic panel data

This paper develops a Bayesian approach to analyzing quantile regression models for censored dynamic panel data. We employ a likelihood-based approach using the asymmetric Laplace error distribution and introduce lagged observed responses into the conditional quantile function. We also deal with the initial conditions problem in dynamic panel data models by introducing correlated random effects into the model. For posterior inference, we propose a Gibbs sampling algorithm based on a location-scale mixture representation of the asymmetric Laplace distribution. It is shown that the mixture representation provides fully tractable conditional posterior densities and considerably simplifies existing estimation procedures for quantile regression models. In addition, we explain how the proposed Gibbs sampler can be utilized for the calculation of marginal likelihood and the modal estimation. Our approach is illustrated with real data on medical expenditures.     

Analyzing return asymmetry and quantiles through stochastic volatility models using asymmetric Laplace error via uniform scale mixtures

This paper proposes a new approach to analyze stock return asymmetry and quantiles. We also present a new scale mixture of uniform (SMU) representation for the asymmetric Laplace distribution (ALD). The use of the SMU for a probability distribution is a data augmentation technique that simplifies the Gibbs sampler of the Bayesian Markov chain Monte Carlo algorithms. We consider a stochastic volatility (SV) model with an ALD error distribution. With the SMU representation, the full conditional distribution for some parameters is shown to have closed form. It is also known that the ALD can be used to obtain the coefficients of quantile regression models. This paper also considers a quantile SV model by fixing the skew parameter of the ALD at specific quantile level. Simulation study shows that the proposed methodology works well in both SV and quantile SV models using Bayesian approach. In the empirical study, we analyze index returns of the stock markets in Australia, Japan, Hong Kong, Thailand, and the UK and study the effect of S&P 500 on these returns. The results show the significant return asymmetry in some markets and the influence by S&P 500 in all markets at all quantile levels. Copyright © 2014 John Wiley & Sons, Ltd.     

Quantile regression based on counting process approach under semi-competing risks data

In this paper, we investigate the quantile regression analysis for semi-competing risks data in which a non-terminal event may be dependently censored by a terminal event. The estimation of quantile regression parameters for the non-terminal event is complicated. We cannot make inference on the non-terminal event without extra assumptions. Thus, we handle this problem by assuming that the joint distribution of the terminal event and the non-terminal event follows a parametric copula model with unspecified marginal distributions. We use the stochastic property of the martingale method to estimate the quantile regression parameters under semi-competing risks data. We also prove the large sample properties of the proposed estimator, and introduce a model diagnostic approach to check model adequacy. From simulation results, it shows that the proposed estimator performs well. For illustration, we apply our proposed approach to analyze a real data.     

Birnbaum‐Saunders quantile regression and its diagnostics with application to economic data

The Birnbaum‐Saunders (BS) distribution is a model that frequently appears in the statistical literature and has proved to be very versatile and efficient across a wide range of applications. However, despite the growing interest in the study of the BS distribution, quantile regression modeling has not been considered for this distribution. To fill this gap, we introduce a class of quantile regression models based on the BS distribution, which allows us to describe positive and asymmetric data when a quantile must be predicted using covariates. We use an approach based on a quantile parameterization to generate the model, permitting us to consider a similar framework to generalized linear models, providing wide flexibility. The methodology proposed includes a thorough study of theoretical properties and practical issues, such as maximum likelihood parameter estimation and diagnostic analytics based on local influence and residuals. The performance of the residuals is evaluated by simulations, whereas an illustrative example of income data is conducted using the methodology to show its potential for applications. The numerical results report an adequate performance of the approach to quantile regression, indicating that the BS distribution is a good modeling choice when dealing with data that have both positive support and asymmetry. The economic implications of our investigation are discussed in the final section. Hence, it can be a valuable addition to the tool kit of applied statisticians and econometricians.     

Composite hierachical linear quantile regression

Multilevel （hierarchical） modeling is a generalization of linear and generalized linear modeling in which regression coefficients are modeled through a model, whose parameters are also estimated from data. Multilevel model fails to fit well typically by the use of the EM algorithm once one of level error variance （like Cauchy distribution） tends to infinity. This paper proposes a composite multilevel to combine the nested structure of multilevel data and the robustness of the composite quantile regression, which greatly improves the efficiency and precision of the estimation. The new approach, which is based on the Gauss-Seidel iteration and takes a full advantage of the composite quantile regression and multilevel models, still works well when the error variance tends to infinity, We show that even the error distribution is normal, the MSE of the estimation of composite multilevel quantile regression models nearly equals to mean regression. When the error distribution is not normal, our method still enjoys great advantages in terms of estimation efficiency.     

Modal Regression Based on Nonparametric Quantile Estimator

Modal regression based on nonparametric quantile estimator is given. Unlike the traditional mean and median regression, modal regression uses mode but not mean or median to represent the center of a conditional distribution, which helps the model to be more robust for outliers, asymmetric or heavy-taileddistribution. Most of solutions for modal regression are based on kernel estimation of density. This paper studies a new solution for modal regression by means of nonparametric quantile estimator. This method builds on the fact that the distribution function is the inverse of the quantile function, then the flexibility of nonparametric quantile estimator is utilized to improve the estimation of modal function. The simulations and application show that the new model outperforms the modal regression model via linear quantile function estimation.     

Pyramid Quantile Regression

We describe a Bayesian model for simultaneous linear quantile regression at several specified quantile levels. More specifically, we propose to model the conditional distributions by using random probability measures, known as quantile pyramids, introduced by Hjort and Walker. Unlike many existing approaches, this framework allows us to specify meaningful priors on the conditional distributions, while retaining the flexibility afforded by the nonparametric error distribution formulation. Simulation studies demonstrate the flexibility of the proposed approach in estimating diverse scenarios, generally outperforming other competitive methods. We also provide conditions for posterior consistency. The method is particularly promising for modeling the extremal quantiles. Applications to extreme value analysis and in higher dimensions are also explored through data examples. Supplemental material for this article is available online.     

Quantile regression for robust bank efficiency score estimation

We discuss quantile regression techniques as a robust and easy to implement alternative for estimating Farell technical efficiency scores. The quantile regression approach estimates the production process for benchmark banks located at top conditional quantiles. Monte Carlo simulations reveal that even when generating data according to the assumptions of the stochastic frontier model (SFA), efficiency estimates obtained from quantile regressions resemble SFA-efficiency estimates. We apply the SFA and the quantile regression approach to German bank data for three banking groups, commercial banks, savings banks and cooperative banks to estimate efficiency scores based on a simple value added function and a multiple-input–multiple-output cost function. The results reveal that the efficient (benchmark) banks have production and cost elasticities which differ considerably from elasticities obtained from conditional mean functions and stochastic frontier functions.     

Logistic响应分布分位数的鞍点近似置信区间

论文基于响应数据，应用鞍点近似方法，给出构造Logistic响应分布分位数的近似置信区间的方法. 论文还对这种置信区间进行了模拟，并将该方法应用于QD8电雷管. 模拟和实例结果表明，当样本量较小时，该方法能够较好地推断Logistic响应分布的分位数     

A goodness-of-fit statistic for Pareto-type behaviour

The fit of a statistical model can be visually assessed by inspection of a quantile–quantile or QQ plot. For the strict Pareto distribution, since log-transformed Pareto random variables are exponentially distributed, it is natural to consider an exponential quantile plot based on the log-transformed data. In case the data originate from a Pareto-type distribution, the Pareto quantile plot will be linear but only in some of the largest observations. In this paper we modify the Jackson statistic, originally proposed as a goodness-of-fit statistic for testing exponentiality, in such a way that it measures the linearity of the k largest observations on the Pareto quantile plot. Further, by taking the second-order tail behaviour of a Pareto-type model into account we construct a bias-corrected Jackson statistic. For both statistics the limiting distribution is derived. Next to these asymptotic results we also evaluate the small sample behaviour on the basis of a simulation study. The method is illustrated on two practical case studies.     

Bankruptcy prediction of small and medium enterprises using a flexible binary generalized extreme value model

We introduce a binary regression accounting-based model for bankruptcy prediction of small and medium enterprises (SMEs). The main advantage of the model lies in its predictive performance in identifying defaulted SMEs. Another advantage, which is especially relevant for banks, is that the relationship between the accounting characteristics of SMEs and response is not assumed a priori (eg, linear, quadratic or cubic) and can be determined from the data. The proposed approach uses the quantile function of the generalized extreme value distribution as link function as well as smooth functions of accounting characteristics to flexibly model covariate effects. Therefore, the usual assumptions in scoring models of symmetric link function and linear or pre-specified covariate-response relationships are relaxed. Out-of-sample and out-of-time validation on Italian data shows that our proposal outperforms the commonly used (logistic) scoring model for different default horizons.     

左截断相依数据下条件分位数的双核局部线性估计

本文在左截断相依数据下,利用局部线性估计的方法,先提出了条件分布函数的双核估计;然后利用该估计导出了条件分位数的双核局部线性估计,并建立了这些估计的渐近正态性结果;最后,通过模拟显示该估计在偏移和边界点调节上要比一般的核估计更好.     

A generalized Newton algorithm for quantile regression models

This paper formulates the quadratic penalty function for the dual problem of the linear programming associated with the \(L_1\) constrained linear quantile regression model. We prove that the solution of the original linear programming can be obtained by minimizing the quadratic penalty function, with the formulas derived. The obtained quadratic penalty function has no constraint, thus could be minimized efficiently by a generalized Newton algorithm with Armijo step size. The resulting algorithm is easy to implement, without requiring any sophisticated optimization package other than a linear equation solver. The proposed approach can be generalized to the quantile regression model in reproducing kernel Hilbert space with slight modification. Extensive experiments on simulated data and real-world data show that, the proposed Newton quantile regression algorithms can achieve performance comparable to state-of-the-art.     

Quantile process for left truncated and right censored data

In this paper, we consider the product-limit quantile estimator of an unknown quantile function when the data are subject to random left truncation and right censorship. This is a parallel problem to the estimation of the unknown distribution function by the product-limit estimator under the same model. Simultaneous strong Gaussian approximations of the product-limit process and product-limit quantile process are constructed with rate . A functional law of the iterated logarithm for the maximal deviation of the estimator from the estimand is derived from the construction. Work partially supported by NSC Grant 89-2118-M-259-011.     

Multi-stage nested classification credibility quantile regression model

In insurance (or in finance) practice, in a regression setting, there are cases where the error distribution is not normal and other cases where the set of data is contaminated due to outlier events. In such cases the classical credibility regression models lead to an unsatisfactory behavior of credibility estimators, and it is more appropriate to use quantile regression instead of the ordinary least squares estimation. However, these quantile credibility models cannot perform effectively when the set of data has nested (hierarchical) structure. This paper develops credibility models for regression quantiles with nested classification as an alternative to Norberg’s (1986) approach of random coefficient regression model with multi-stage nested classification. This paper illustrates two types of applications, one with insurance data and one with Fama/French financial data.     

Bayesian quantile regression model for claim count data

Quantile regression model estimates the relationship between the quantile of a response distribution and the regression parameters, and has been developed for linear models with continuous responses. In this paper, we apply Bayesian quantile regression model for the Malaysian motor insurance claim count data to study the effects of change in the estimates of regression parameters (or the rating factors) on the magnitude of the response variable (or the claim count). We also compare the results of quantile regression models from the Bayesian and frequentist approaches and the results of mean regression models from the Poisson and negative binomial. Comparison from Poisson and Bayesian quantile regression models shows that the effects of vehicle year decrease as the quantile increases, suggesting that the rating factor has lower risk for higher claim counts. On the other hand, the effects of vehicle type increase as the quantile increases, indicating that the rating factor has higher risk for higher claim counts.     

Statistical inference for quantiles in the frequency domain

For second-order stationary processes, the spectral distribution function is uniquely determined by the autocovariance function of the process. We define the quantiles of the spectral distribution function in frequency domain. The estimation of quantiles for second-order stationary processes is considered by minimizing the so-called check function. The quantile estimator is shown to be asymptotically normal. We also consider a hypothesis testing for quantiles in frequency domain and propose a test statistic associated with our quantile estimator, which asymptotically converges to standard normal under the null hypothesis. The finite sample performance of the quantile estimator is shown in our numerical studies.     

A new weighted Gompertz distribution with applications to reliability data

A new weighted version of the Gompertz distribution is introduced. It is noted that the model represents a mixture of classical Gompertz and second upper record value of Gompertz densities, and using a certain transformation it gives a new version of the two-parameter Lindley distribution. The model can be also regarded as a dual member of the log-Lindley-X family. Various properties of the model are obtained, including hazard rate function, moments, moment generating function, quantile function, skewness, kurtosis, conditional moments, mean deviations, some types of entropy, mean residual lifetime and stochastic orderings. Estimation of the model parameters is justified by the method of maximum likelihood. Two real data sets are used to assess the performance of the model among some classical and recent distributions based on some evaluation goodness-of-fit statistics. As a result, the variance-covariance matrix and the confidence interval of the parameters, and some theoretical measures have been calculated for such data for the proposed model with discussions.     

Quantile Processes in the Presence of Auxiliary Information

We employ the empirical likelihood method to propose a modified quantile process under a nonparametric model in which we have some auxiliary information about the population distribution. Furthermore, we propose a modified bootstrap method for estimating the sampling distribution of the modified quantile process. To explore the asymptotic behavior of the modified quantile process and to justify the bootstrapping of this process, we establish the weak convergence of the modified quantile process to a Gaussian process and the almost-sure weak convergence of the modified bootstrapped quantile process to the same Gaussian process. These results are demonstrated to be applicable, in the presence of auxiliary information, to the construction of asymptotic bootstrap confidence bands for the quantile function. Moreover, we consider estimating the population semi-interquartile range on the basis of the modified quantile process. Results from a simulation study assessing the finite-sample performance of the proposed semi-interquartile range estimator are included.     

A generalization of the half-normal distribution

In this paper we introduce an extension of the half-normal distribution in order to model a great variety of non-negative data. Its hazard rate function can be decreasing or increasing, depending on its parameters. Some properties of this new distribution are presented. For example, we give a general expression for the moments and a stochastic representation. Also, the cumulative distribution function, the hazard rate function, the survival function and the quantile function can be easily evaluated. Maximum likelihood estimators can be computed by using numerical procedures. Finally, a real-life dataset has been presented to illustrate its applicability.