菜鸟一个,想请教大家一个问题,用openbugs具体如何实现量化退火。
模型如下:
model
{
for( i in 1 : N ) {
y[i] ~ dnorm(mu,tau) }
mu ~ dnorm(0,0.01)
tau ~ dgamma(0.01,0.01)
sigma.squared <- 1 / tau
sigma <- sqrt(sigma.squared)
}
data list(N = 10, y = c(-1.76, 0.38, 1.23, -0.67, -0.47, -1.36, 1.41, -0.07, -1.23, 2.35))
intial list(mu = 1, tau = 2)
intial list(mu = 2, tau = 4)