data='empltran.dat';

model:'A. Competing risks multinomial logit model'  regression nominal    1 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=1  bayess2=1  bayeslat=1
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=0  nseedboot=0
  nrand=10
  usemiss=No
  sewald=yes  dummy=first
 outsect=0x1d17
;
repwgt REPLICATION_WEIGHT;
dependent OUTCOME;
replicate ID;
predictor COHORT2 COHORT3 FEMALE;
attr COHORT2 ordinal ;
attr COHORT3 ordinal ;
attr FEMALE ordinal ;
attr OUTCOME nominal ;
classind COHORT2 COHORT3 FEMALE;

end;
model:'B. 2-class version of A, with class-specific intercept'  regression nominal    2 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=1  bayess2=1  bayeslat=1
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=1382786  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes  dummy=first
 outsect=0x1d17
;
repwgt REPLICATION_WEIGHT;
dependent OUTCOME;
replicate ID;
predictor COHORT2 COHORT3 FEMALE;
attr COHORT2 ordinal ;
attr COHORT3 ordinal ;
attr FEMALE ordinal ;
attr OUTCOME nominal ;
classind COHORT2 COHORT3 FEMALE;

end;
model:'C. 3-class version of A, with class-specific intercept'  regression nominal    3 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=1  bayess2=1  bayeslat=1
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=704720  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes  dummy=first
 outsect=0x1d17
;
repwgt REPLICATION_WEIGHT;
dependent OUTCOME;
replicate ID;
predictor COHORT2 COHORT3 FEMALE;
attr COHORT2 ordinal ;
attr COHORT3 ordinal ;
attr FEMALE ordinal ;
attr OUTCOME nominal ;
classind COHORT2 COHORT3 FEMALE;

end;
model:'D. Model A + random intercept terms (using CFactor)'  regression nominal    1 /
  cfn=1  cfnde=10  cfcsp=0
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=1  bayess2=1  bayeslat=1
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=2208176  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes  dummy=first
 outsect=0x1d17
;
repwgt REPLICATION_WEIGHT;
dependent OUTCOME;
replicate ID;
predictor COHORT2 COHORT3 FEMALE;
cfactor 1 <- OUTCOME;
attr COHORT2 ordinal ;
attr COHORT3 ordinal ;
attr FEMALE ordinal ;
attr OUTCOME nominal ;
classind COHORT2 COHORT3 FEMALE;

end;