On quantile quantile plots for generalized linear models


Augustin, N., Sauleau, E.-A. and Wood, S., 2012. On quantile quantile plots for generalized linear models. Computational Statistics & Data Analysis, 56 (8), pp. 2404-2409.

Related documents:

PDF (Author's version) - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
Download (23MB) | Preview

    Related URLs:


    The distributional assumption for a generalized linear model is often checked by plotting the ordered deviance residuals against the quantiles of a standard normal distribution. Such plots can be difficult to interpret, because even when the model is correct, the plot often deviates substantially from a straight line. To rectify this problem Garcia Ben and Yohai (2004, J. Comput. Graph. Stat. 13: 36-47) proposed plotting the deviance residuals against their theoretical quantiles, under the assumption that the model is correct. Such plots are closer to a straight line, when the model is correct, making them much more useful for model checking. However the quantile computation proposed in Garcia Ben and Yohai is, in general, relatively complicated to implement and computationally expensive, so that general purpose software for these plots is only available for the Poisson and binary cases in the R package {\tt robust}. As an alternative the theoretical quantiles can efficiently and simply be estimated by repeatedly simulating new response data from the fitted model and computing the corresponding residuals. This method also provides reference bands for judging the significance of departures of QQ-plots from ideal straight line form. A second alternative is to estimate the quantiles using quantiles of the response variable distribution according to the estimated model. This latter alternative generally has lower computational cost than the first, but does not yield QQ-plot reference bands. In simulations the quantiles produced by the new methods give results indistinguishable from the original Garc\'ia Ben and Yohai quantile computations, but the scaling of computational cost with sample size is much improved so that a 500 fold reduction in computation time was observed at sample size 50000. Application of the methods to generalized linear models fitted to prostate cancer incidence data suggest that they are particularly useful in large dataset cases that might otherwise be incorrectly viewed as zero-inflated. The new approaches are simple enough to implement for any exponential family distribution and for several alternative types of residual, and this has been done for all the families available for use with generalized linear models in the basic distribution of R.


    Item Type Articles
    CreatorsAugustin, N., Sauleau, E.-A. and Wood, S.
    Related URLs
    URLURL Type
    DepartmentsFaculty of Science > Mathematical Sciences
    Publisher Statementqq_gam_resub.pdf: NOTICE: this is the author’s version of a work that was accepted for publication Computational Statistics & Data Analysis. Changes resulting from the publishing process, such as peer review, editing, corrections, structural formatting, and other quality control mechanisms may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. A definitive version was subsequently published in Computational Statistics & Data Analysis, vol 56, issue 8, 2012, DOI 10.1016/j.csda.2012.01.026
    ID Code27091


    Actions (login required)

    View Item

    Document Downloads

    More statistics for this item...