title:
  data from page 166 in eqs program.
  input adapted from appendix 6.3 page 166.
  figure 6.1 page 77.
  output. 

    table 6.2 (combined columns), page 79.
    note. correlation is standardized value of
      int with slope.

    also table 6.3 (cohort sequential column), computed
    based on intercept and slope values, e.g.
      age 12: .304 + 0*.394 = .304
      age 13: .304 + 1*.394 = .698
      age 14: .304 + 2*.394 = 1.092
      age 15: .304 + 3*.394 = 1.486
      age 16: .304 + 4*.394 = 1.88
data:  
  file is ch6co.dat;
  format is free;
  type is means covariance;
  ngroups = 4;
  nobservations = 248 250 231 215;
variable:  
  names are y1 y2 y3 y4 y5; 
  usevar =  y1-y5;
analysis:  
  type = general mgroup meanstructure;          
  iterations = 200;
model: 
  ! model for group 4
  int by y1@1 y2@1 y3@1 y4@1 y5@1;
  slp by y1@0 y2@1 y3@2 y4@3 y5@4;
  [y1@0 y2@0 y3@0 y4@0 y5@0];
  [int]        (m1);
  [slp]        (m2);
  int with slp (r1);
  int          (d1);
  slp          (d2);
  y1           (e1);
  y2           (e2);
  y3           (e3);
  y4           (e4);
  y5           (e5);
model g1:
  ! group 1
  int by y1@1 y2@1 y3@1 y4@0 y5@0; ! different from overall
  slp by y1@0 y2@1 y3@2 y4@0 y5@0; ! different from overall
  [y1@0 y2@0 y3@0 y4@0 y5@0];
  [int]        (m1);
  [slp]        (m2);
  int with slp (r1);
  int          (d1);
  slp          (d2);
  y1           (e1);
  y2           (e2);
  y3           (e3);
  y4@1; ! different from overall
  y5@1; ! different from overall
model g2:
  ! group 2
  int by y1@0 y2@1 y3@1 y4@1 y5@0; ! different from overall
  slp by y1@0 y2@1 y3@2 y4@3 y5@0; ! different from overall
  [y1@0 y2@0 y3@0 y4@0 y5@0];
  [int]        (m1);
  [slp]        (m2);
  int with slp (r1);
  int          (d1);
  slp          (d2);
  y1           (e1);
  y2           (e2);
  y3           (e3);
  y1@1; ! different from overall
  y5@1; ! different from overall
model g3:
  ! group 3
  int by y1@0 y2@0 y3@1 y4@1 y5@1; ! different from overall
  slp by y1@0 y2@0 y3@2 y4@3 y5@4; ! different from overall
  [y1@0 y2@0 y3@0 y4@0 y5@0];
  [int]        (m1);
  [slp]        (m2);
  int with slp (r1);
  int          (d1);
  slp          (d2);
  y1           (e1);
  y2           (e2);
  y3           (e3);
  y1@1; ! different from overall
  y2@1; ! different from overall
model g4:
  ! no changes, same as overall model
  ! model for group 4
  int by y1@1 y2@1 y3@1 y4@1 y5@1;
  slp by y1@0 y2@1 y3@2 y4@3 y5@4;
  [y1@0 y2@0 y3@0 y4@0 y5@0];
  [int]        (m1);
  [slp]        (m2);
  int with slp (r1);
  int          (d1);
  slp          (d2);
  y1           (e1);
  y2           (e2);
  y3           (e3);
  y4           (e4);
  y5           (e5);
output: 
  tech1 samp standardized;