回复 第1楼 的 tingerzouxian：
白话说，下面波浪符号前的参数都要给出initial estimate。
<br />
lamada ~ dbeta(2,1)<br />
alpha0 ~ dnorm( 0.0,1)<br />
alpha1 ~ dnorm( 0.0,1)<br />
alpha2 ~ dnorm( 0.0,1)<br />
alpha3 ~ dnorm( 0.0,1)<br />
lata ~ dnorm( 0.0,2)<br />
isigmma2 ~ dgamma(4,2)<br />
muz ~ dnorm( -0.01,20)<br />
tauz ~ dgamma(4,2)<br />
lamada ~ dbeta(2,1)<br />
sigmma<-1/sqrt(isigmma2)<br />
sigmmaz<-1/sqrt(tauz)<br />
</p>
我没看懂的是：模型是个微分方程，可是WingBugs code给出的似乎是个analytical solution equation.
<br />
for( i in 1 : N ) {mu[i] <- r[i] + alpha0 + alpha1 * r[i] + alpha2 * pow(r[i],2) + alpha3 / r[i] + za[i] * ja[i]<br />
}<br />
</p>
