model{
  for(i in 1:n){
    q[i] <- 1/n
  }
  tau ~ dcat(q[])
  for(i in 1:2){
    lambda[i] ~ dgamma(alpha,beta)
  }
  for(i in 1:n){
      mu[i] <- lambda[1]+step(i-tau-0.1)*(lambda[2]-lambda[1])
      y[i] ~ dpois(mu[i])
  }  
}