假设我们使用 JAGS(或 WinBUGS)拟合贝叶斯线性混合模型,输出对象是否包含模型残差?我们怎样才能找到残差?
谢谢!
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
。请注意,上面的示例不涉及指定要监视的值。