Прогнозировать значения с помощью функции spBreg_err в пакете пространственного регистра в r

Question

Я использую функцию из spatialregбиблиотека в r для пространственного моделирования переменной с использованием нескольких ковариат. Результат из spBreg_errэто объект mcmc, который можно открыть с помощью codaупаковка.

Мой вопрос в том, как я могу предсказать значения в точках и получить вероятность для данного квантиля

Документация находится здесь: https://r-spatial.github.io/spatialreg/reference/SET_MCMC.html .

Я использую один из кодов, который по умолчанию похож на пример.

          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))   

Спасибо