Built-in models_Unconditional linear growth curve models

Code (Edit)

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:4){
      y[i, t] ~ dnorm(muY[i,t], Inv_Sig_e2)
      muY[i,t]<-LS[i,1]+LS[i,2]*t
    }
  }

  #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)
  Sig_e2 <- 1/Inv_Sig_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(bL=c(0), bS = c(0), Inv_Sig_e2 = 1, 
Inv_cov= structure(.Data = c(1,0,0,1),.Dim=c(2,2)))