bayesian - 我们可以从 JAGS 输出中得到残差吗？

Question

假设我们使用 JAGS(或 WinBUGS)拟合贝叶斯线性混合模型，输出对象是否包括模型残差？我们怎样才能找到残差？

谢谢!

Answer 1

JAGS (BUGS) 模型只是输出您让它监控的模型中的节点值。要获得残差，您需要在模型中定义残差，然后监视它们。例如

model {
 bResponse ~ dnorm(0, 5^-2)
 sResponse ~ dunif(0, 5)
 for (i in 1:length(Response)) {
   eResponse[i] <- bResponse
   Response[i] ~ dlnorm(eResponse[i], sResponse^-2)
   Residual[i] <- (log(Response[i]) - log(eResponse[i])) / sResponse
 }
} 

were Residual[i] 定义 Response 的每个 i 值的残差。请注意，上面的示例不涉及指定要监视的值。