自己编写了一段程序,在运行导入数据时,一直报错,“expected avriable 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(1.45347E-05,4.3846E-05,6.33331E-05,8.2211E-05,5.01096E-05,6.02634E-05,
3.69578E-05,0.000173611,2.20814E-05,0.000526567,0.000313453,0.000285528,0.000367107,1.61614E-05,1.46316E-05,2.90693E-05,8.28202E-05,3.89742E-05,0.000113268,8.3516E-06,6.02634E-05,9.44477E-05,8.68056E-05,2.20814E-05,0.000975904,0.000572861,0.000611845,0.001019742,8.08068E-05,0.000190211,7.26734E-06,9.74355E-06,4.87177E-06,5.48073E-06,3.76646E-06,1.23193E-05,5.61672E-05,6.48522E-05,4.07897E-05,0.000163159,1.45347E-05,3.89742E-05,2.43589E-05,5.48073E-06,6.2637E-06,2.63652E-05,1.23193E-05,0.001625,2.10627E-05,3.24261E-05,8.15794E-05,3.23227E-05,2.92633E-05,4.10642E-06,3.24261E-05,4.07897E-05,8.15794E-05,7.26734E-06,4.87177E-06,1.82691E-06,1.40418E-05,4.32348E-05,4.07897E-05,4.87177E-06,4.87177E-06,1.64422E-05,4.1758E-06,7.53293E-06,4.21254E-05,5.40435E-05,0.000163159,8.46238E-05,1.61614E-05),
s=c(0.6,3.2,1.046153846,1.493333333,1.075,0.94375,0.4,10,0.05,0.628,1.220689655,0.485714286,2.022222222,0.2,0.2,0.2,0.323529412,1.075,0.57983871,3.0625,0.2125,1.093478261,0.05,0.2,0.660431655,0.869811321,0.903333333,0.93,0.28,0.296153846,0.2,0.4,6.0,0.533333333,0.2,0.416666667,0.35,2.266666667,0.2,0.7,0.3,0.20625,0.2,0.15,0.416666667,0.407142857,0.15,0.284615385,0.2,0.2,0.525,0.5,0.425,0.2,0.15,0.4,0.4,0.8,0.05,0.2,0.6,0.2,0.8,0.8,0.2,0.411111111,0.6,0.3,0.15,0.4,0.55,0.2,0.6),N=73)
list(mu1 = 0, tau1 = 1,mu2 = 0,tau2=1, tau = 1,mu=0)