#WinBUGS codes generated by BAUW: http://bauw.psychstat.org
#USE WITH CAUTION!

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] ~ dnorm(muY[i,t], Inv_Sig_e2[i])
      muY[i,t]<-LS[i,1]+LS[i,2]*t
    }
    w[i]~dgamma(a, a)
    Inv_Sig_e2[i]<-w[i]*Pre_e2
  }

  #Priors for model parameter

  #k~dexp(0.1)
  k~dunif(0,100)
  a<-k/2
 

  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
  Pre_e2 ~ dgamma(.001, .001)
  Sig_e2 <- 1/Pre_e2
  Cov[1:2,1:2]<-inverse(Inv_cov[1:2,1:2])
  Sig_L2 <- Cov[1,1]
  Sig_S2 <- Cov[2,2]
  Cov_LS <- Cov[1,2]
  rho_LS <- Cov[1,2]/sqrt(Cov[1,1]*Cov[2,2])
}

# The (naive) starting values for model parameters.
list(k=4, bL=c(0), bS = c(0), Pre_e2 = 1, Inv_cov= structure(.Data = c(1,0,0,1),.Dim=c(2,2)))

# Please put your data here. Remeber you can convert your data by BAUW--Convert Data. :)

list( N = 375 , y = structure(.Data = c( NA, 6.5, 5, 6.1, NA, 4.8, 4.4, 5, NA, 
6.2, 6, 3.7, 4.2, 5.2, 5.3, 6.1, 8.5, 6.9, 7.1, 7.1, NA, 3.8, 7.2, 6.9, 7.6, 
5.9, 7.9, 6.1, 5.7, 7, 5.2, 7.2, 7.9, 6.9, NA, 7.3, 7.9, 8.2, 7.3, 9, 4.1, 6.7, 
NA, NA, 7, 6.9, 7.8, 7.2, 7, 7.4, 5.1, 4.5, 5, 7, 6.9, 6.3, 7.2, 6.9, 7.4, NA, 
8, 8.5, 8.9, 9.2, 9.7, 7.8, 8.1, 9.6, 9.7, 8.1, 7.4, 5, NA, 6.2, 3.8, 8, 6.5, 
8.5, 8.6, 8.6, 6.3, 7.2, 6.4, NA, 6.9, 6.8, 8.6, 7.3, 9, 9.5, 8.2, 9.2, 9.5, 
9.7, 9.8, 6.1, 5.8, 7.9, 7.2, 7.2, 7.6, 5.5, 6.9, 8.2, 5.8, 8, 8.1, 8.7, 8.5, 
9.2, 5, 4.8, 6.1, 7.3, 5, 6.4, NA, 5.7, 6.5, 5.5, 9.3, 9.3, 9.5, 9.7, 10, 5, 
5.8, 7.7, 7.6, 5.8, 4.9, 5.4, 7.2, 6.5, 6.4, 5.8, 6.2, 8.1, 8.1, 8.4, 6, 7, 
7.2, 8.2, 8.3, 6.7, 9, 8.4, 7.2, 8.7, 6.8, 8, 8.5, 9, 8.9, 6.6, 5.9, 6, 7, 7.1, 
7.5, 7.6, 8.1, 8.5, 7.9, 9, 8.9, 9.7, 9.5, 10, 7.3, 7.6, 8.1, 7.9, 8.5, 4.4, 5, 
NA, 5.5, 5.3, 3.1, NA, 2.8, NA, NA, 6, 5.6, 6.7, 6.5, 7.2, 7.9, 6.8, 9.8, 8.8, 
7.3, 9.2, 8, 9.1, NA, 7.8, 3.5, 6.2, 7, 7.4, 7.4, 3.9, 5, 7.1, NA, 7.5, 6.1, 
8.8, 6.3, 8.3, 7.9, 7, 6.3, 6.4, 7, 9.1, 6.9, 5.3, 6.5, 6.6, NA, 6, 6.7, 7.5, 
7, 8, 2.4, 4.8, 4.2, 5.2, 5.9, 4, 5, 6.6, 6.7, 5.9, NA, 2.9, 4.9, 4.9, 4.5, 
5.5, 6.3, 7.9, 6.9, 6.8, 6, 4.7, 5.2, 4.9, NA, 5.4, 6.4, 7.8, 6.8, 7.1, 6.1, 
7.5, 7.7, 9.3, NA, 5.8, 4.7, 6.1, 5, 5.3, 4, 5.6, 5.1, 6.2, 5.7, 8.5, 7.4, 10, 
9.2, 9.5, 6, 5.7, 8.3, 7.7, 8.2, 5.8, NA, 4.9, 5.5, 8, 6, 8.8, 8.9, 8.9, NA, 
7.5, 7.3, 7.4, 8.2, 8.9, 5.4, 4.5, 4.3, 5.7, 5, 4.4, NA, 4.6, 7.3, 7.3, 3, 5.7, 
3, 5.9, 4, 7, NA, 9.4, NA, 9.7, 5.1, 5.9, 4.9, 2.5, 6.1, 8.4, 8.9, 9.5, 9.3, 
9.5, 5.8, 6.7, NA, 8.3, 8.3, 6.1, 6, 6.5, 8.9, 9.1, 5.6, 5.1, 6.1, 6, 5.5, 6.5, 
6.8, 7.2, 6.6, 6.2, 6.6, 7.3, NA, 6.9, NA, 5.8, 5.8, 6, 7.4, 7.4, 3.9, 5.3, NA, 
NA, 6.5, 5.9, 7.6, 5, NA, 5.9, 6, NA, 6, NA, 5.7, 6.2, 7.4, 7.2, NA, 6.8, NA, 
6.6, 4.9, 6.6, NA, 8.5, 8.8, 9.1, 9.7, 9.6, 5.8, 4.9, 5.4, 3.9, 5, 5, 6.2, 5, 
7.4, 6.4, 7.8, 9.3, NA, 9.3, 9.2, 5.8, 7.5, 6.1, 5.2, 6.8, 6.3, 5.9, 8.3, 6.7, 
5.9, 6.7, 8.3, 9, 9.2, 9.3, 7.9, 7.8, 8.6, 9.4, 9.6, 4.9, 5.8, 6.6, 5.6, 6.8, 
5, 7, 6.4, 6.8, 7.4, 4.8, 6, 5.6, 5.2, 5.5, 6.7, 7.4, NA, NA, NA, 9.5, 9.6, 
9.5, 9.9, 9.8, NA, 8.6, 8.8, 9.3, NA, 6.6, 7.1, 7.2, NA, NA, 7.4, NA, 8.8, 8.9, 
NA, 5.5, 4.2, 2.7, 4.1, 5.7, 6.3, 5, 5.2, 6.9, 6.7, NA, 5.9, 7.5, 6.5, 6.2, 
7.6, 8.3, 8, 8.6, 7.7, 6.7, 6.4, 6.3, 9.2, 7.2, 6.9, 7, 7.3, 7.6, 7.1, 5.6, NA, 
4.3, 6.8, 8, 7.1, 5.5, 7.2, 6.6, 6.8, 6, 3.7, 8.8, 7.8, 8.5, 6.9, 6.6, 6.1, 
7.3, 5.6, 7.1, 8.5, 7.4, 7, 9.6, 6.2, 6.5, NA, 2.2, 6.3, 5.3, 6.5, 6.7, 7.3, 
7.1, NA, 4.1, 3.4, NA, 4.9, NA, 7.2, 6.6, 7.3, 7.5, 5.1, 6.3, 6, 6.7, NA, 5.6, 
6.8, 7.5, NA, 9, 5.2, 6.2, 7, 6, 7.1, 1.6, 2, 2.3, 2.6, 3.3, 6, 5.4, 6.2, 7.7, 
6.8, 8.8, 9.6, 8.9, NA, 9.8, 5.6, 6.5, 7.5, 8.3, 9.4, 8.5, 8.9, 9.2, NA, NA, 
6.6, 7, 6.6, 6.9, 6.8, 7.2, 7.3, 7.4, 6, 5.4, 6.8, 8.2, 7.9, 7.7, 9.1, 6, 6.7, 
6.6, 8.1, 8.2, 5.5, 4.7, 3.4, 6.5, 6, 8, 9, 8.9, 6.4, 8.8, 6.7, 8.7, 9.9, 8.7, 
10, 7.9, 8.8, 7.3, 9.2, 9, 8.1, 7.8, 8.2, 8.7, 8.9, 6.5, 8.5, 7.5, 8.8, 7.3, 
2.8, 0.2, 5.4, 6.1, 6.6, 8.8, 8.7, 9, 9.7, 9.6, 7.5, 7.6, 8.3, 8.7, 8.2, 6.5, 
6.3, 6.9, 6.3, NA, 7.5, 7.3, 7.9, 8.1, 7.5, 7, 7.2, 8.5, 8.9, 9.3, 5.8, 4.6, 
5.3, 7, 6.7, 8.3, 8.7, 8.6, 9.7, 8.7, 6.6, 6.9, 6.8, 8.9, 9.2, 6.5, 5.9, 6.6, 
7.3, NA, 3.7, 3.6, NA, 4.3, 4.5, 7.8, 8.3, 6.3, 7.4, 8.8, 4.7, 5.5, 5.8, 6.7, 
7.6, 9.7, 5, 4.5, 10, NA, 5.2, 5.6, 5.6, 10, 7.6, 4.4, 2.3, 3.9, 10, 5.6, 4.6, 
NA, 4.9, NA, 5.6, 5.5, 6.6, 6.6, 6.8, 7, 4.7, 5.2, 6.3, 6.6, 6.9, 4.9, 6.3, 
5.9, 6.2, 6.2, 5.2, 5, 5.6, 6.2, 5.9, 5, 4.7, 5.4, 6.3, 6.3, 5.6, 6.9, 5, 7.5, 
6.8, 7.3, 8.6, 7.5, 7.5, 8.2, 5.3, 5.7, 5.6, 6.9, 7.1, 7, 7.2, 7, 7.2, 7.2, 
5.7, 5.5, 5.5, 5.3, 5.7, 5.1, 5.8, 6.4, 7, 6.3, 5.3, 5.7, 6, 7, 7.1, 5, 7, 7, 
7.2, 8.1, 5.6, 5.5, 6.5, 7, 6.8, 4.9, 4.7, 8.9, 6.9, 6.9, 6.5, 7.2, 6.5, 7.6, 
7.1, 7.5, 7.8, 9.2, 9.4, 9.4, 6.8, 7.1, 5.6, 7.5, 6.7, 5.5, 6.1, 6, 6.4, 6.2, 
6.3, 7.1, 6.6, 6.7, 0.3, 5.6, 5, 0, 5.9, 8.2, 6.2, 6.7, 6.3, 5.7, 0.3, 6.2, 
5.3, 7.3, 6.5, 7.5, 5.6, 5.3, 5.1, NA, 5.7, 9.2, 7, 7.1, 7.5, 7.1, NA, 5.9, 
6.5, 6.4, 4.9, 7.1, 6.9, 6.9, 7.5, 9.3, 8.6, 8.2, 9, 9.1, 9.4, 8.6, 8.7, 9.2, 
9.4, 9.3, NA, 3.6, 4.1, 4.6, NA, 5.8, 5.4, 6.5, 6.4, 6.2, 6.9, 7.1, 8.7, 8.1, 
9.4, 8.4, 3.5, 2.9, 2.6, 9.6, 5, 7, 6.8, 8, 10, 7.6, 5.9, 7.3, 8.9, NA, 7.5, 9, 
9.2, 9.7, 9.9, 4.3, 4.3, 7.1, 4.6, 4.4, 4, 4.4, 4, 4, 4.1, 4.8, 4.8, 8, 8.4, 
8.2, 7, 7.3, 9.1, 7.1, NA, 4.6, 4.5, 4.9, 6.3, 6.5, 3.8, 5.8, 6.6, 6.5, 6.2, 
4.9, 4.6, 6.6, 6.5, NA, 7.4, NA, 7, 8.5, 6.8, 4.6, 7.5, 7.8, 7.4, 6.8, 7, 7.1, 
8.7, 9.1, 9, 6.5, 4, 6, 7, 7.2, 6.1, 5.4, 5.6, 5.5, 6.5, 4.8, 4.3, 4.4, 8, 6.1, 
3.5, 5.5, 5.7, 3.4, 5.4, 7.3, 7.2, 8.4, 8.4, 8.7, NA, 3.9, 6.2, 3.3, 3.8, 6.3, 
9.3, 8.1, 8.9, 7.2, 7.3, 7.9, 7.1, 7.5, 9.2, 6.3, 8.4, 7.5, 6.4, 6.6, NA, 6.6, 
5.5, 6.6, 5, 6.6, 6.8, 6.9, 6.8, 7.1, 6.1, 5.7, 3.8, 4.3, 7.1, 6.3, 5.3, 9, 
8.5, 8.7, NA, 5.4, NA, 7, 5.1, 6.6, 7, 6, 6.9, 4.7, 7.3, 9.4, 9.2, 9.5, 9.6, 
6.1, 5, 4.9, 5, NA, 9.1, NA, 9.4, NA, 9.9, 8.3, 7.6, 7.2, NA, 7.6, 5.1, 5.7, 
NA, NA, 6, 5, 6, 6.3, 5, 5, 6.4, 7.5, 8.8, NA, 9.1, 2.4, NA, NA, 4.5, NA, NA, 
9.1, 9.3, 9.6, 9.5, 5.1, 4.9, 4.9, 9.9, 7.2, 6.2, 8.9, 9.1, 9, 9.2, 5.6, 6.5, 
6.9, 0.4, 7.1, 6.3, 6, 6.9, 6.4, 6.8, 4.9, 6.7, 7.7, 9, 7, 6.8, 7.1, 7.1, 7.9, 
8.3, 5.9, 6.8, 6.4, 6.2, 6.3, 5.5, 5.9, 6.8, 7.1, 5, 8, 8.1, 7.6, 9.9, 8.8, 
4.8, 4.6, 4.5, 6.6, 6.1, 6, 6.8, 6.9, 7, 7.3, 5.7, 6.4, 6.8, 6.8, 6.7, 6.3, 
7.3, 7, 5.7, 6.3, 5.1, 6, 5.3, 5.1, 4.7, 6, 7.4, 6.9, 8.3, 6.9, 5.3, 5.8, 6.1, 
6.1, 6.3, 6, 6.2, 5.7, 6, 7.3, 6.6, 7.5, 7.5, NA, 7.3, 5.3, 3.8, 3.6, NA, 4, 
4.7, 3.7, 9.9, 9.8, 7.1, 6, 6, 9.4, NA, 6.2, 6.9, NA, 7.1, 8.7, 9.3, 4.8, 2.3, 
2.6, 4.1, 5.2, 5.4, 3.8, 5.3, NA, 5.9, 5.5, 5.9, 5.7, 7.4, 7.6, 6.5, 8, 6.8, 
7.2, 8.3, 3.8, 5.4, 4.4, 5.2, 5.4, 3.8, 2.8, 5.4, 6.6, 6.7, 1.7, 4, 0.8, 2.2, 
2.5, 2.5, 2.6, 7, 1.8, 6.3, NA, 7.3, 6.2, 8.5, 6.3, 4.8, 4.4, 5.9, 5.3, 6.1, 
5.8, 7, 6.9, 6.1, 6, 8, 8.8, 8.7, 9.3, 9.3, 6.4, 6.6, 8.6, 8.6, 8.2, 2.3, 0.1, 
NA, 5.5, 5.2, 6.5, 6.8, 6.4, 6.9, NA, 8.5, 7.3, 5.8, 6.4, 6, 5.5, 6.7, 6.9, 
7.7, 7.3, 4.5, 5.2, 5.6, 5.2, 7.5, 3.1, 1.9, 1.9, 5, 2.8, 2, 2.3, 3.8, 3.8, 5, 
5.6, 6.8, 7.2, 7, 6.7, 6.4, 7.8, 8, 7.4, NA, 6, 5.4, 6.6, 6.7, NA, 4.4, 3.8, 
4.7, 6.5, NA, 6, 6.1, 6.3, 6, 6.8, 5.4, 6.2, 6.9, 5.1, 6.5, 4.9, 4.7, 6.8, 6.2, 
7.2, 4.4, 4.9, 4.7, 4.7, 5.2, 6.3, 7.3, 7.2, 8.8, 8.7, 6.5, 7.7, 9.3, NA, 7.1, 
9.2, 9.3, 9.7, 9.7, 10, 6, 7.5, 8.3, 6.9, 6.8, 5.9, 6.8, 7.3, 7.5, 8.1, 8.3, 
9.6, 8.8, 9.3, 9.6, 2, 6.1, 5.8, 2.1, 4, 7.6, 8.1, 8.5, 6.7, 6.9, 7.3, 4.9, 
5.6, 4.2, 5.3, NA, 4.6, 5.9, 5.7, NA, 6.9, 8.2, NA, 7.2, 7.1, 6.5, 6.9, NA, 
7.9, 8.3, 4.2, 4.4, 4.7, 5, 5.9, 6.1, 6.8, 7.2, 7.3, 7.6, 8.7, 9, 9.2, 9.7, 9, 
6.7, 7.4, 6.9, 7.7, 9.5, 9, 9.2, 9, 9.1, 10, 6.5, 5.9, 5.1, NA, 7.2, NA, 5.9, 
6.1, 10, 6.2, 8.5, 9, 8.9, 9, 8.8, 7, 6.9, 6.4, 7.8, 8.7, 5.9, 7.5, 6.9, 7, 
7.6, 7.1, 6.5, 9.1, 9.4, NA, 7.6, 8, 7.3, 7.5, 7.8, 6.9, NA, 8, 7.1, 9.6, 4.8, 
5.7, 6, 5.9, 6.7, 5.4, 4.1, 5.9, 6.8, 3.8, 5.8, 6.7, 9.4, 9.5, 9.3, 5.4, 4.8, 
5.8, 7, 7, 6.4, 5.4, 5.2, 3.5, 5.1, 5.3, 7.2, 6.8, 6.6, 7.9, 6, 7, 6.7, 6.5, 
7.2, 3.1, NA, 4, 4.6, 5.2, 5.9, 6.1, 6.9, NA, NA, 5.9, 5, 10, 6.1, 9.7, 6.7, 
7.6, 7.4, 7.5, 8.3, 6.2, 9.2, NA, NA, 5.5, 5.2, 6.6, 9.2, 7.6, 6.4, 4.2, 0.6, 
4.8, 7.5, 0.1, 8.1, 7.3, NA, NA, 7.2, 3.4, 3.6, 4.4, 5.4, 4, 6.6, 7, NA, NA, 
9.1, 8.6, 5.6, 7, NA, 4.6, 6.6, 7.2, 7.9, NA, 7.9, 4.9, 7.3, 7, 7.8, 6.8, 6.2, 
5, 6.1, 8.1, 6.5, 9.3, 9, 9.7, 9.8, 10, 2.2, 2, NA, 2.7, 6.8, 7, 7, 7.9, 8.3, 
8.5, 8.1, NA, NA, 9.1, 9.2, 5.8, 7.3, 7.3, 9.7, 7.3, NA, 0.5, 3.2, 3.9, 4.4, 
9.5, 9.8, 9.7, 9.7, 9.5, 5.4, 6, 6, 7, 7.1, 7.2, 7, NA, NA, 7.2, NA, 4.9, 6.1, 
5.6, 5.4, NA, 8.2, 7.2, 9.5, 6.5, NA, 7.1, 8.7, 6.7, 7.6, 6.8, 5.6, 8.8, 8, NA, 
5.5, 8.8, 7.6, 9.4, 8.9, 7.8, 7.7, 8.8, 8.8, 8.7, 8.6, 9.1, 9.1, 8.8, 9.6, 6.2, 
6.9, NA, 4.9, NA, 3.1, 4.5, 5, 4.1, 5.8, 4.6, 6.4, NA, 7.3, 7, 6.1, 6.7, 6.4, 
6.5, 5.8, 5.9, 5.6, 5.9, 7.5, 5.6, NA, 6.6, 5.7, 7.1, NA, 2.9, 4.4, 4.4, 4.2, 
6.3, 4.9, 5.2, 5.4, 6.2, 5.9, 4.3, 7.8, 7, 9, 6.9, 8.3, 6.8, 7.6, 7.1, 6.9, 
6.3, 6.5, 7.9, 7.2, 6.5, 5.5, 5.9, 5.9, 5.8, 5.9, 5.1, 6.9, 5.2, 6, 3.8, 2.8, 
6.1, NA, 5.9, NA, 5.8, 5.6, 6.7, 6.6, 6.7, 5.4, 4.7, 5.8, 5.9, 5.8, NA, 4, 5.5, 
6.8, 7, 4.6, 6.1, 5.4, 7, 7.3, 4.4, 5.6, 6.6, 5.4, 6.7, 6.3, 6.9, 6.6, 6.8, 
7.1, 4.6, 5.4, 6.4, 4.9, 6.8, NA, 6.9, 7.8, 5.8, 6.4, NA, 5.7, 6, 7.4, 5.1, 
6.8, 3.7, 5.1, 5.6, 5.6, 7.4, 6.4, 7.8, 8, 7.3, 7, 6.5, 6.5, 5.9, 6.6, 8.3, 
7.3, 9.1, 9.3, 9.2, 7.4, 7, 5.3, 6.4, 6.4, 5.9, 5.6, 8, 5.7, 5.4, 7.2, NA, NA, 
4.9, 5, 7.8, 9.1, 8.9, 8.7, 8.2, 3.1, 0.1, 4.7, 4.9, 5.7, 8.4, 7.7, 6.5, 8.7, 
6.4, 6.3, 0, 6.4, 7.1, 7.3, 4.2, 5, 5.7, 5.7, 6.2, 3.1, 3.7, 5.7, 4.1, 3.9, 
5.4, 7.1, 8.8, 8.8, 8.7, 5.4, NA, 5.3, 4.6, 7.1, 6, 6.2, 6.2, 6.1, 8.7, 8.1, 
6.1, 7.2, 8, 8.8, 8.8, 9, 7, 9.2, 9.3, 2.4, 3, 2.7, 2.7, 2, 6.3, 5.8, 6.4, 5.1, 
6.4, 7.9, 8, 8.2, 8.7, 8.6, 4.6, 4.8, 6.3, 6.8, 6.5, 6.6, 6.9, 7.6, 6.9, 6.5, 
6.8, 6.9, 6.8, 7, 8.7, 5.9, 6, 6.4, 6.1, 7.2, 7.3, 8.6, 8.8, 9, NA, 6.9, 6.5, 
6.9, 7.4, NA), .Dim=c( 375, 5 )) )