#model
model{
for (i in 1:N){risk[i,1:2] ~ dmnorm(beta[i,1:2],inverse(structure(.Data=c(risk[i,3],risk[i,4],risk[i,5],risk[i,6]),.Dim=c(2,2))))}
for (i in 1:N){beta[i,1:2] ~ dmnorm(mu[1:2],pre.sigma[1:2,1:2])}
mu[1:2] ~ dmnorm(c(0,0),inverse(structure(.Data=c(1,0,0,1),.Dim=c(2,2))))
pre.sigma[1:2,1:2] ~ dwish(structure(.Data=c(1,0,0,1),.Dim=c(2,2)),2)
pre.sigma[1:2,1:2] <- inverse(sigma[1:2,1:2])
}

#data
list(N=3,risk=structure(.Data=c(0.0001164745,-0.0006360886,5.454619*10^(-8),-2.15988*10^(-6),-2.15988*10^(-6),0.02313066,
0.0002386487,-0.0001279721,1.419456*10^(-7),-3.196332*10^(-6),-3.196332*10^(-6),0.03024943,
-6.353768*10^(-5),-0.0001478020,2.418857*10^(-7),-4.84376*10^(-6),-4.84376*10^(-6),0.08127031),.Dim=c(N,6)))
)

#initial values
list(beta=structure(.Data=rep(c(1,0),N),.Dim=c(N,2)),mu=c(1,1),sigma2=structure(.Data=c(1,0,0,1),.Dim=c(2,2)))

运行代码check model时即报错expected variable name,小弟初次接触WinBUGS，请大家多多指教