2

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

我遇到过，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)))？

4

0 回答 0