/*
**   Confirmatory Factor Model Simulation
**
*/


rndseed 1786897;


lambda = { 1  0,
          .8  0,
           0  1,
           0 .8 };

phi = {  1  .5,
        .5   1  };

theta = { .3  0  0  0,
           0 .3  0  0,
           0  0 .3  0,
           0  0  0 .3 };


sigma = lambda * phi * lambda' + theta;

dat = rndn(100,4) * chol(sigma);

call saved(dat,"maxfact","Y");