想做一个潜变量的分析。模型如下:其中x1是潜变量，观测不到的。
y1,y2 观测的到，y1是内生变量(endogenous V), 不能直接放在(2)等式中，只能用e1放在其中，这样b2反应的就是把x1控制之后的y1对y2的影响了

y1=a0+a1*x1+e1                      (1)

y2=b0+b1*x1+b2*e1+e2          (2)

我编的程序如下，请大牛指教一下可不可以：
model
{  
for(i in 1: N) {
x1[i] ~ dnorm(-1, z)

y2 ~ dnorm(y2mu[i], tau.y2)
y1~ dnorm (y1mu[i], tau.y1)

y2mu[i] <- b0+b1*x1[i]+b2*e1[i]
y1mu[i] <- a0+a1*x1[i]+e1[i]  
}


tau.y1 ~ dgamma(0.5, 0.5)
tau.y2 ~ dgamma(0.5, 0.5)

z ~ dgamma(0.5, 0.5)
a0~dnorm(0.0, 0.01)
a1~dnorm(0.0, 0.01)

b0~dnorm(0.0, 0.01)
b1~dnorm(0.0, 0.01)
b2~dnorm(0.0, 0.01)


}