Valid Standard Errors for Bayesian Quantile Regression With Clustered and Independent Data.

Authors

Feng Ji, JoonHo Lee, Sophia Rabe-Hesketh

Published in

Journal of educational and behavioral statistics : a quarterly publication sponsored by the American Educational Research Association and the American Statistical Association. Volume 51. Issue 5. Pages 1055-1084. Epub Nov 06, 2025.

Abstract

Bayesian quantile regression typically uses the asymmetric Laplace distribution as working likelihood, not because it is a plausible data-generating distribution but because the corresponding maximum likelihood estimator is identical to the classical estimator by Koenker and Bassett. While point estimation is consistent, credible intervals tend to have poor frequentist coverage. We propose using infinitesimal jackknife (IJ) standard errors introduced by Giordano and Broderick, which do not require resampling and can be obtained from a single Markov chain Monte Carlo run. Simulations and applications to real data show that IJ standard errors have good frequentist properties for both independent and clustered data. We provide an R package, IJSE, that computes IJ standard errors after estimation of any model with the brms wrapper for Stan.

PMID:
42669008
Bibliographic data and abstract were imported from PubMed on 30 Aug 2026.

Read full publication at:
Please sign in to see all details.

Stats

1-terrible, 9-excellent. How would you rate this publication? Sign in in to submit your rating.

Post a comment

You need to be signed in to post comments. You can sign in here.

Comments

There are no comments yet.