model

{
 for (i in 1:4)
     {y[i,1:4]~dmulti(p[i,],n[i])
      for (k in 1:4)
          {p[i,k]<-phi[i,k]/sum(phi[i,])
           phi[i,k]<-exp(eta[i,k])
           }
   
      for (k in 1:4)
          {eta[i,k]<-beta0[k]+betaa[k]*age[i]+betas[k]*smoke[i]+betaas[k]*age[i]*smoke[i]
           }
      }

 beta0[1]<-0
 betaa[1]<-0
 betas[1]<-0
 betaas[1]<-0
      
 for (k in 2:4)
     {beta0[k]~dnorm(mu0,1.0)
      betaa[k]~dnorm(mua,10.0)
      betas[k]~dnorm(mus,10.0)
      betaas[k]~dnorm(muas,20.0)
      }

 mu0~dnorm(-2,1.0)
 mua~dnorm(0,10.0)
 mus~dnorm(0,10.0)
 muas~dnorm(0,20.0)

 }