# NOT RUN {
## Generate artificial data
theta = list(xi=c(3,5), psi=c(2,4), G=diag(2), nu=6)
n = 1000
p = 2
X = rmvSE(n, p, 'ST', theta, seed=NULL)
y = X$y
latentVars = list(z=X$z, v=X$v)
## CML estimates for theta.
thetaHat = cmlSE('ST', y, z=X$z, v=X$v)
# }
Run the code above in your browser using DataLab