小编Min*_*hau的帖子

我试图拟合一个多项式逻辑回归模型，其rjags结果是一个具有 3 个级别的分类（名义）变量（结果），解释变量是年龄（连续）和组（具有 3 个级别的分类）。在此过程中，我想获得Age和Group的后验均值和基于分位数的 95% 区域。

我不太擅长，for loop我认为这就是我为模型编写的代码无法正常工作的原因。

我的 beta 先验遵循正态分布，\xce\xb2j \xe2\x88\xbc Normal(0,100) for j \xe2\x88\x88 {0, 1, 2}。

可重现的 R 代码

library(rjags)\n\nset.seed(1)\ndata <- data.frame(Age = round(runif(119, min = 1, max = 18)),\n                   Group = c(rep("pink", 20), rep("blue", 18), rep("yellow", 81)), \n                   Outcome = c(rep("A", 45), rep("B", 19), rep("C", 55)))\n\nX <- as.matrix(data[,c("Age", "Group")]) \nJ <- ncol(X)\nN <- nrow(X)\n\n## Step …