yury123
模型VARMA模型如下
[img]file:///C:/Users/liqi/AppData/Local/Temp/IICQ$%7DDC%25WUUF[9(C[X%25I%252.jpg[/img] 我的程序如下
model
{
for(i in 2:n)
{for(j in 1:m)
{y[i,j] ~dmnorm (mu[i,j],tauC)
mu[i,j]<-b[i,j]-theta*equals(j,1)*b[i-1,j]+phi*equals(j,2)*y[i-1,j]}
b[i,1:m]~dnorm(0,tau[1:m,1:m])
tau[1:m,1:m]~dwish(R[1:m,1:m],m)
R[1:m,1:m]~dgamma(0.1,0.1)}
phi~dbeta(1,1)
theta~dbeta(1,1)
tauC~dgamma(0.1,0.1)
}
请问我这样写对么？
check model 无误的
