Оценка псевдо-R2 для модели бета-регрессии Зоиба

Question

У меня есть данные и код ниже, и буду признателен за любой вклад в pseudo-Rsquared или Bayesian version of R-squared можно оценить для модели бета-регрессии с zoib,

library(zoib)
data("GasolineYield", package = "zoib")
re.md <- zoib(yield ~ temp | 1 | 1, data=GasolineYield, 
joint = FALSE, random=1, EUID=GasolineYield$batch, 
zero.inflation = FALSE, one.inflation = FALSE, 
n.iter=3200, n.thin=15, n.burn=200)
sample2 <- re.md$coeff
summary(sample2)

Пакет vignette и авторы статьи не обсуждают, как pseudo-Rsquared можно получить в zoib,

Подобные пакеты, такие как betareg включите это в сводку модели, как показано ниже:

library(betareg)
data("GasolineYield", package = "betareg")
md <- betareg(yield ~ temp, data = GasolineYield, link = "cloglog")
summary(md)

Поэтому я был бы признателен за любые идеи о подходах, чтобы получить это для zoib модель.