===== Manuscripts =====
  - Version 1. {{:lab:projects:02robust_growth_curve_modeling_using_student-t_distribution:student_t_growth_curve_modeling_mar_03_2010.pdf}}
===== Read data analysis =====
==== Data ====
  * [[NLSY real data]]
  * [[R codes related to real data analysis]]
  * {{:lab:projects:02robust_growth_curve_modeling_using_student-t_distribution:digitspan.csv|NLSY digit span variable (???)}}
===== Simulations =====
  * [[BUGS codes for data simulation]]
  * [[R codes to process simulated data from BUGS]]
  * [[R codes to generate simulation submission scripts]]
  * [[R codes to process the coda files]]
  * [[R codes for final results process]]
  * [[R function for processing DIC]]

==== Using SAS to run simulation ====
  * [[SAS codes]]
  * [[Latent basis GCM in bugs]] 
  * Initial values for uniform prior <code>list(A=c(NA, .25, .5, .75, NA),k=4, bL=c(6), bS = c(1), Pre_e2 = 1, rho=0, 
Cov= structure(.Data = c(1,NA,NA,1),.Dim=c(2,2)))</code> 
  * [[GCM with Wishart prior|Latent basis GCM with Wishart prior]]
  * Initial values for Wishart prior <code>list(A=c(NA, .25, .5, .75, NA),k=4, bL=c(0), bS = c(0), Pre_e2 = 1,
Inv_cov= structure(.Data = c(1,0,0,1),.Dim=c(2,2)))
</code>
  * Script to run BUGS <code>modelCheck("model.txt")
modelData("data.txt")
modelCompile()
modelInits("inits.txt")
modelGenInits()
modelUpdate(10000)
samplesSet("para")
modelUpdate(50000)
samplesCoda("para","coda")
modelQuit()</code>


==== Simulation Results ====



===== Some useful scripts =====

  * zip files
<code>tar czvf Ndata-Tdata-0.1.tgz NdataR</code>
  * [[Simulate data and analyze in Mplus]]
  * [[Generate student t random numbers using auxiliary variable method in R]]

===== References =====
  * [[http://sekhon.berkeley.edu/library/MASS/html/cov.trob.html|Estimate the covariance matrix for multivariate t]]

~~DISCUSSION~~