model{
  for (i in 1:N)
    {
		TTA[i] ~ dgamma(shape, rate)
	  }
  ttastar ~ dgamma(shape, rate) 
  ht <- step(120-ttastar)
	shape ~ dgamma(0.1, 0.1)
	rate ~ dgamma(0.1, 0.1)
}