最近在研究winbugs,编写了如下程序,但在 “load data”时一直显示:expected variable name.请各位大神看看。程序如下:
model
{
for (i in 1:N) {
R~dlnorm(mu,tau)
R<-p*s
p~dlnorm(mu1,tau1)
s~dlnorm(mu2,tau2)
}
mu1~dnorm(0,0.000001)
tau1~dgamma(0.01,0.001)
mu2~dnorm(0,0.000001)
tau2~dgamma(0.01,0.001)
mu~dnorm(0,0.000001)
tau~dgamma(0.01,0.001)
sigma<-1/sqrt(tau)
sigma1<-1/sqrt(tau1)
sigma2<-1/sqrt(tau2)
}
list(p=c(4.3846E-04,8.7846E-04,5.2489E-04,1.3811E-04,2.7605E-04,2.1212E-04,1.5369E-04,3.2355E-04,9.4635E-04,4.4001E-04,1.0022E-04,9.6560E-04,4.3810E-04,2.1000E-04,2.33997E-04,9.9860E-04,4.4112E-04,5.6568E-04,6.0024E-04,3.4635E-04,6.0156E-04,8.6350E-04,2.3253E-04,5.6837E-04),s=c(0.8,2.00,0.75,1.20,0.90,1.60,2.59,3.75,0.60,1.81,2.00,0.75,2.20,1.99,1.20,1.38,0.75,1.20,0.65,0.05,0.60,0.41,0.90,0.86,,N=24)
list(mu1=0,tau1=1,mu2=0,tau2=1,tau=1,mu=0)