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)
}

list(N=19, TTA=c(224, 188, 133,  98, 113,  91,  97, 159, 124, 128, 267, 316, 279, 210, 422, 102,  75, 413, 177))
list(shape=.5, rate=.5, ttastar=100)