Bayes Book
Chapter 6 : Simulations

Exercise 6.10

model{
 for(i in 1:n){ y[i] ~ dweib(alpha,lambda) }
 lambda ~ dgamma(1.53,26.3)
 alpha ~ dgamma(1,1)
 median <- log(2)/lambda
 S24 <- exp(-lambda*pow(24,alpha))
}
list(n=17,
      y=c(65,156,100,134,16,108,121,4,39,143,56,26,22,1,1,5,65))
list(lambda=0.05, alpha =1)