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