# NOT RUN {
## Set n, p and the parameter values
theta = list(xi=c(3,5), psi=c(2,4), G=diag(2), nu=6)
n = 1000
p = 2
## Generate the data
X = rmvSE(n, p, 'ST', theta, seed=NULL)
## X contains the data matrix and the vectors z and v of latent variables:
y = X$y
z = X$z
v = X$v
# }
Run the code above in your browser using DataCamp Workspace