model
{
for (i in 1:T){
for (j in 1:m){
y[i,j]~dnorm(ymean[i,j],tau)
}}
for (i in 2:T){
for(j in 1:m){
ymean[i,j]<-theta*y[i-1,j]+enta
}}
# prior
theta~dnorm(0,0.001)
enta~dnorm(0,tau)
tau~dgamma(0.001,0.001)
}
INITS
list(theta=0,tau=1,enta=1)
DATA
list(T=46,m=14400))
y[,1] y[,2] y[,3] y[,4] y[,5] y[,6] y[,7] y[,8] y[,9] y[,10] y[,11] y[,12] y[,13] y[,14] y[,15] y[,16] y[,17] y[,18] y[,19] y[,20] y[,21] y[,22] y[,23] y[,24] y[,25] y[,26] y[,27] y[,28] y[,29] y[,30] y[,31] y[,32] y[,33] y[,34] y[,35] y[,36] y[,37] y[,38] y[,39] y[,40].............y[,14400]
compile 的时候,提示:made use of undefined ymean
不知道是怎么回事!!!!
这是一个一阶自回归模型。