输入的data---sst是一个120*120的矩阵,lat和log分别是120个矢量,表示经纬度。意思是每一个矩阵元素的均值是他上下左右元素的加权均值。
for (i in 2:119){
for (j in 2: 119){
sst[i,j]~dnorm(mu[i,j],tau)
mu[i,j]<-mu0[i,j]+alfa*((mu[i-1,j]-mu0[i-1,j])+(mu[i+1,j]-mu0[i+1,j]))+beta*((mu[i,j-1]-mu0[i,j-1])+(mu[i,j+1]-mu0[i,j+1]))+errorsigma[i,j]
mu0[i,j]<-mu01 + mu02* lat[j]+ mu03*log[i]
errorsigma[i,j]~dgamma(0.001,0.001)
}
}
在compile的时候,总是提示:array index is greater than array uper bound for mu 是怎么回事啊。


雷达卡



京公网安备 11010802022788号







