model
{
a ~ dnorm (0, 1.0E-6) # prior for the intercept
b ~ dnorm (0, 1.0E-6) # prior for the slope
prec ~ dgamma(0.001, 0.001) # prior for the precision
for (i in 1:16)
{
mean <- a + b*Dens
CWD ~ dnorm(mean, prec)
}
}
list(Dens = c(1270, 1210, 1800, 1875, 1300, 2150,
1330, 964, 961, 1400, 1280, 976, 771, 833, 883, 956),
CWD = c(121, 41, 183, 130, 127, 134, 65, 52, 12, 46, 54,
97, 1, 4, 1, 4))
list(a=0, b=0, prec=100)
我想计算当Dens=1500时,mean 的估计值,及变动区间(在winbugs怎么完成?)
predicted_mean <- a + b*1500
predicted_obs ~ dnorm(predicted_mean, prec)
|