TITLE:

  DATA:
  FILE IS "simlist.dat";
  TYPE = MONTECARLO;

  VARIABLE:
  names are modav x1 x2 zyg pat;
  USEVARIABLES ARE modav x1 x2 zyg;
  CONSTRAINT IS modav;
  MISSING ARE ALL (999);

  grouping is zyg (1=mz 2=dz);

  ANALYSIS:
  MODEL = NOCOVARIANCES;

  MODEL:
  !ace components for delta.
  T2A1 by x1*(l2a);
  T2A2 by x2*(l2a);

  T2C by x1*(l2c);
  T2C by x2*(l2c);

  T2E1 by x1*(l2e);
  T2E2 by x2*(l2e);

  x1@0;
  x2@0;

  !COVARY OUT SES
  x1 x2 ON modAV*(s);

  !T2BS1 T2BS2 ON RACE*(r);

  !set residual/variances at 1
  T2A1@1;
  T2A2@1;

  T2C@1;

  T2E1@1;
  T2E2@1;

  !set means at 0
  [T2A1@0];
  [T2A2@0];

  [T2C@0];

  [T2E1@0];
  [T2E2@0];

  !make intercepts equal
  [x1*](miq);
  [x2*](miq);


  !differentiate mzs and dzs
  MODEL MZ:
  T2A1 with T2A2@1;

  MODEL DZ:
  T2A1 with T2A2@.5;
  [T2A1@0];
  [T2A2@0];


  MODEL CONSTRAINT:
  NEW (l2a0*.327 l2a1*.05 !vary startvalues to match pop values
  l2c0*.678 l2c1*-.015
  l2e0*-.422 l2e1*.034);
  !l2a2 l2c2 l2e2);

  l2a=l2a0 +l2a1*modAV;

  l2c=l2c0 +l2c1*modAV;

  l2e=l2e0 +l2e1*modAV;


  OUTPUT:
  STANDARDIZED SAMPSTAT PATTERNS;

  SAVEDATA: 
  RESULTS ARE read_simresults.dat;