model smoking



{ for (i in 1:36)
      
      {deaths[i]~dbin(p[i],n[i])

       logit(p[i])<-alpha[group[i]]+beta[group[i]]*(age[i]-60)

       }


 for (g in 1:4)

     {alpha[g]~dnorm(alpha0,0.2)

      beta[g]~dnorm(beta0,10)

      }


 alpha0~dnorm(-2,0.01)

 beta0~dnorm(0.0,8)

 }