回复 第1楼 的 realme：
我也刚学WinBugs，共同讨论一下.
1. 下面这段和你最开始说的模型不符合:
L2<-(pow(10, x1)*2541.79+pow(10, x2)*3816.4+pow(10, x3)*160.95)/1000*(1-exp(-(pow(k, 2)*t)))/(pow(k, 2)*t)
最开始说的模型是:
y=(a1*x1+a2*x2+a3*x3)*(1-exp(-k*t))/(k*t)
2. 还有,前验分布中说到的参数e1, e2, e3在model中是x1, x2, x3, 这不符合.
3. 尝试着改了一下程序如下
model
{
for (i in 1:I)
{
L1 <- log(L)
L1 ~ dnorm(mu, tau)
mu <- log(L2)
L2 <- (2541.79*x1 + 3816.4*x2 + 160.95*x3) * (1 - exp(-k*t[i])) / (k*t[i])
}
log.x1 ~ dnorm(1.4181, 3.286)
log(x1) <- log.x1
log.x2 ~ dnorm(0.8009, 4.1297)
log(x2) <- log.x2
log.x3 ~ dnorm(1.4122, 3.6488)
log.x3 <- log(x3)
sqrt.k ~ dnorm(0.5288, 19.8958)
sqrt.k <- sqrt(k)
sigma ~ dgamma(0.001,0.001)
tau <- 1/(sigma*sigma)
}
