想在winbugs里面建立一个VAR方程,可是不知道怎么设置协方差的先验分布,以及残差的表达。
model{
for(j in 1:p){
for(i in 3:n){
y[i,j]~dmnorm(mu[j],phi[i,j])
u[j]<-mean(y[i,j])
ephat[i,j]<-u[j]-mu[j]
mu[j]<-b[1,j]*y[i-1,1]+b[2,j]*y[i-2,1]+b[3,j]*y[i-1,2]+b[4,j]*y[i-2,2]+b[5,j]*y[i-1,3]+b[6,j]*y[i-2,3]+b[7,j]*y[i-1,4]+b[8,j]*y[i-2,4]+b[9,j]
}}
for(i in 1:9){b[i,1]~dnorm(8.81,0.68)
b[i,2]~dnorm(11.48,0.99)
b[i,3]~dnorm(5.19,0.53)
b[i,4]~dnorm(0.29,0.80)}
phi[1:4,1:4]~dwish(R[1:4,1:4],4)
phx[1:4,1:4]<-inverse(phi[1:4,1:4])
}
可是总是运行有问题,求高手指点。


雷达卡



京公网安备 11010802022788号







