
model{
  # Model specification for linear growth curve model
  for (i in 1:N){
    LS[i,1:2]~dmnorm(muLS[i,1:2], Inv_cov[1:2,1:2])
    muLS[i,1]<-bL[1]
    muLS[i,2]<-bS[1]
    for (t in 1:5){
      y[i, t] ~ dt(muY[i,t], Inv_Sig_e2, k)
      muY[i,t]<-LS[i,1]+LS[i,2]*t
    }
  }

  k~dexp(.1)

  #Priors for model parameter
  for (i in 1:1){
    bL[i] ~ dnorm(0, 1.0E-6)
    bS[i] ~ dnorm(0, 1.0E-6)
  }
  Inv_cov[1:2,1:2]~dwish(R[1:2,1:2], 2)
  R[1,1]<-1
  R[2,2]<-1
  R[2,1]<-R[1,2]
  R[1,2]<-0
  Inv_Sig_e2 ~ dgamma(.001, .001)
  para[6] <- 1/Inv_Sig_e2
  Cov[1:2,1:2]<-inverse(Inv_cov[1:2,1:2])
  para[1]<-bL[1]
  para[2]<-bS[1] 
  para[3] <- Cov[1,1]
  para[4] <- Cov[2,2]
  para[5] <- Cov[1,2]
  para[7] <- k
}