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
 }
} 

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