model {
for (i in 1:6) {
z[i] <- 1
z[i] ~ dbern(p[i])
p[i] <- exp(-0.5*(y[i] - mu)*(y[i] - mu))/phi(8 - mu)
}
mu ~ dunif(0, 100)
}
Data:
list(y = c(6,6,6,7,7,7))
Inits:
list(mu = 6.5)
node mean sd MC error 2.5% median 97.5% start sample
mu 6.737 0.4965 0.00498 5.819 6.72 7.75 1001 10000