#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 )) )