Predict values with spBreg_err function in spatialreg package in r

Question

I am using spBreg_err function from spatialreg library in r to modelling spatially a variable using several covariates. The outcame from spBreg_err is a mcmc object that is possible to open with coda package.

My question is how can I to predict the values in the points and obtain the probability for a given quantile

The documentation is here: https://r-spatial.github.io/spatialreg/reference/SET_MCMC.html

I am using the one of the codes that have by default like example.

    library(spdep)
    library(spatialreg)
    library(coda)
    data(oldcol, package="spdep")
    lw <- spdep::nb2listw(COL.nb, style="W")
    set.seed(1)
    COL.err.Bayes <- spBreg_err(CRIME ~ INC + HOVAL, data=COL.OLD, listw=lw)
    print(summary(COL.err.Bayes))
    print(raftery.diag(COL.err.Bayes, r=0.01))   

Thanks