Title:  
    First step of 2-group BCH LTA
    with x1 as group variable

data:
    file = LTADistalN2K.dat;

variable:  
    NAMES = u1-u12 x1 y1 y2 x2 y3 x3 class1 class2;
    usev = u1-u12 x1;    
    missing = all(999);  
              
    classes  = c(2) c1(3) c2(3);
              
    categorical = u1-u12 x1;

    Auxiliary = y1 y2 y3;

   
analysis:   
    type = mixture; 

    starts = 40 10;

model:
 
   %overall%

 !  c2 on c1;
 !  c2 c1 on c;

   Model C:

   %c#1%
   [x1$1@15]; ! P(x1=0)=1 for c class = 1
   %c#2%
   [x1$1@-15]; ! P(x1=1)=1 for c class = 2
 

  Model C1:
  %C1#1%
  [u1$1-u6$1] (1-6);


  %C1#2%
  [u1$1-u6$1] (11-16);


  %C1#3%
  [u1$1-u3$1] (21-23);
  [u4$1-u6$1] (24-26);


  Model C2:
  %C2#1%
  [u7$1-u12$1] (1-6);
 
  %C2#2%
  [u7$1-u12$1] (11-16);
 

  %C2#3%
  [u7$1-u9$1] (21-23);
  [u10$1-u12$1] (24-26);
 

savedata:
    save = bchweights;
    file = mgbch.dat;

output:
    svalues;
   
 