大家好，我刚开始用winbugs做贝叶斯，很多地方都不懂，下面我写了一个面板数据的ordered probit模型，大家帮我看看哪里不对，老是不对报错。
model{
for (i in 1:N){for (t in 1:T){
fVie[i,t]~dcat(p[i,t,])
p[i,t,1]<-  phi(r-mu[i,t])
p[i,t,2]<-  phi(r+d-mu[i,t])        -phi(r-mu[i,t])
p[i,t,3]<-  phi(r+d+e-mu[i,t])      -phi(r+d-mu[i,t])
p[i,t,4]<-  phi(r+d+e+f-mu[i,t])    -phi(r+d+e-mu[i,t])
p[i,t,5]<-  phi(r+d+e+f+g-mu[i,t])  -phi(r+d+e+f-mu[i,t])
p[i,t,6]<-  phi(r+d+e+f+g+h-mu[i,t])-phi(r+d+e+f+g-mu[i,t])
p[i,t,7]<-1-phi(r+d+e+f+g+h-mu[i,t])
mu[i,t]<-alpha+b1*marriage[i,t]+b2*ad[i,t]+b3*kid[i,t]+b4*heal[i,t]+b5*dlRev[i,t]+b6*my[i,t]
}}


for (i in 1:N){alpha~dnorm(alphabar,phi)}

alphabar~dnorm(0,0.0001)
b1~dnorm(0,0.0001)
b2~dnorm(0,0.0001)
b3~dnorm(0,0.0001)
b4~dnorm(0,0.0001)
b5~dnorm(0,0.0001)
b6~dnorm(0,0.0001)
phi~dgamma(0.1,0.1)
r~dnorm(0,0.001)
d~dgamma(0.1,0.1)
e~dgamma(0.1,0.1)
f~dgamma(0.1,0.1)
g~dgamma(0.1,0.1)
h~dgamma(0.1,0.1)
}