程序如下

model
{
for (t in 1:n) {
z~dnorm(2.25,3.33);
p<-0.3*z+0.25;
v[t]<-478.55*exp(p);
}
}
list(z=c(1,2,1,1,3,6,4,2,5,7,1,1,6,0,,5,4,1,0,1,7),
n=20)