R：用于非线性混合效应模型的贝叶斯包

Question

我正在寻找非线性混合效果模型的贝叶斯并行，特别是那些在R中使用nlme包的模型。

我遇到blme，但这似乎只适用于线性混合效果模型。 brms在这种情况下是否合适？我尝试使用函数nlme编写一些类似于brm构造的代码。

library(nlme)
model <- nlme(height ~ exp(beta1*age + 1),
              data = Loblolly,
              fixed = list(beta1 ~ 1),
              random = list(Seed = pdDiag(list(beta1 ~ 1))),
              start = list(fixed = c(beta1 = 3)))

library(brms)
bayesian_model <- brm(bf(height ~ exp(beta1*age + 1), beta1 ~ 1, nl = TRUE), 
                      data = Loblolly,
                      prior = c(prior(normal(0, 1), nlpar = beta1)))

我能够达到这一点，但我究竟如何为beta1指定随机效果？我如何指定像random = list(Seed = pdDiag(list(beta1 ~ 1)))那样的对角线方差结构？