- model {
- for (i in 1:N) {
- y[i] ~ dnorm(mu[i],tau)
- mu[i] <- alpha + beta*X[i]
- }
- # Blocking of coefficients
- coef[1:2] ~ dmnorm(m0[1:2],prec[1:2,1:2])
- alpha <- coef[1]
- beta <- coef[2]
- lsigma ~ dunif(-k,k)
- sigma <- exp(lsigma)
- tau <- pow(sigma,-2)
- }


雷达卡

京公网安备 11010802022788号







