Bayesian Quantile Regression

Why this is not working, how can we hint it?..

myModel <- brms::brm(brms::bf(myFormula, quantile=seq(0.1, 0.9, by = 0.1)), data = df_train, family = asym_laplace(),chains=myChains,iter=myIter)
Caused by error:
! Expecting a single value when fixing parameter ‘quantile’.

hello @Steviey this error is occurring because you aren’t allowed to specify more than one conditional quantile. You will need a separate model for each quantile.

Models for multiple quantiles are generally supported by frequentist methods, but the asymmetric Laplace based approach to quantile regression will not allow this. The model works by setting the target quantile equal to the location of the distribution, which means that it can only handle one quantile at a time. Find the full description here (https://doi.org/10.1016/S0167-7152(01)00124-9).

There are some described approaches to Bayesian quantile modelling for multiple quantiles, but to my understanding they would not be readily compatible with Stan.

2 Likes