spss='hagdu.SAV';

model:'A. 1-cluster'  cluster 1 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=1  bayeslat=0
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=0  nseedboot=0
  nrand=10
  usemiss=No
  sewald=yes
 outsect=0x1c17
;
casewgt freq;
indicator syva_DFA syva_EIA abbott_EIA genprobe sanofi_EIA culture;
attr syva_DFA ordinal ;
attr syva_EIA ordinal ;
attr abbott_EIA ordinal ;
attr genprobe ordinal ;
attr sanofi_EIA ordinal ;
attr culture ordinal ;

end;
model:'B. 2-cluster'  cluster 2 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=1  bayeslat=0
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=2395682  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c17
;
casewgt freq;
indicator syva_DFA syva_EIA abbott_EIA genprobe sanofi_EIA culture;
attr syva_DFA ordinal ;
attr syva_EIA ordinal ;
attr abbott_EIA ordinal ;
attr genprobe ordinal ;
attr sanofi_EIA ordinal ;
attr culture ordinal ;

end;
model:'C. 2-cluster + CFactor'  cluster 2 /
  cfn=1  cfnde=20  cfcsp=0
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=1  bayeslat=0
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=2372325  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c17
;
casewgt freq;
indicator syva_DFA syva_EIA abbott_EIA genprobe sanofi_EIA culture;
attr syva_DFA ordinal ;
attr syva_EIA ordinal ;
attr abbott_EIA ordinal ;
attr genprobe ordinal ;
attr sanofi_EIA ordinal ;
attr culture ordinal ;

end;
model:'D. 2-cluster + CFactor with equal effects (=mixture Rasch)'  cluster 2 /
  cfn=1  cfnde=20  cfcsp=0
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=1  bayeslat=0
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=1695501  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c17
;
equxy -1;
casewgt freq;
indicator syva_DFA syva_EIA abbott_EIA genprobe sanofi_EIA culture;
attr syva_DFA ordinal ;
attr syva_EIA ordinal ;
attr abbott_EIA ordinal ;
attr genprobe ordinal ;
attr sanofi_EIA ordinal ;
attr culture ordinal ;

end;
model:'E. 2-cluster + CFactor with equal and class specific effects'  cluster 2 /
  cfn=1  cfnde=20  cfcsp=1
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=1  bayeslat=0
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=2207713  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c17
;
equxy -1;
casewgt freq;
indicator syva_DFA syva_EIA abbott_EIA genprobe sanofi_EIA culture;
attr syva_DFA ordinal ;
attr syva_EIA ordinal ;
attr abbott_EIA ordinal ;
attr genprobe ordinal ;
attr sanofi_EIA ordinal ;
attr culture ordinal ;

end;
model:'F. as B with culture as gold standard'  cluster 2 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=1  bayeslat=0
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=498998  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c17
;
casewgt freq;
indicator syva_DFA syva_EIA abbott_EIA genprobe sanofi_EIA;
covariate culture;
attr syva_DFA ordinal ;
attr syva_EIA ordinal ;
attr abbott_EIA ordinal ;
attr genprobe ordinal ;
attr sanofi_EIA ordinal ;
attr culture perfect ;
know  2  * 2
 0 : 1 0 :
 1 : 0 1 :
 ;

end;
model:'G. as C with culture as gold standard'  cluster 2 /
  cfn=1  cfnde=20  cfcsp=0
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=1  bayeslat=0
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=1733639  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c17
;
casewgt freq;
indicator syva_DFA syva_EIA abbott_EIA genprobe sanofi_EIA;
covariate culture;
attr syva_DFA ordinal ;
attr syva_EIA ordinal ;
attr abbott_EIA ordinal ;
attr genprobe ordinal ;
attr sanofi_EIA ordinal ;
attr culture perfect ;
know  2  * 2
 0 : 1 0 :
 1 : 0 1 :
 ;

end;
model:'H. as D with culture as gold standard'  cluster 2 /
  cfn=1  cfnde=20  cfcsp=0
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=1  bayeslat=0
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=1973479  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c17
;
equxy -1;
casewgt freq;
indicator syva_DFA syva_EIA abbott_EIA genprobe sanofi_EIA;
covariate culture;
attr syva_DFA ordinal ;
attr syva_EIA ordinal ;
attr abbott_EIA ordinal ;
attr genprobe ordinal ;
attr sanofi_EIA ordinal ;
attr culture perfect ;
know  2  * 2
 0 : 1 0 :
 1 : 0 1 :
 ;

end;
model:'I. as E with culture as gold standard'  cluster 2 /
  cfn=1  cfnde=20  cfcsp=1
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=1  bayeslat=0
  bayespoi=1
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=2183334  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c17
;
equxy -1;
casewgt freq;
indicator syva_DFA syva_EIA abbott_EIA genprobe sanofi_EIA;
covariate culture;
attr syva_DFA ordinal ;
attr syva_EIA ordinal ;
attr abbott_EIA ordinal ;
attr genprobe ordinal ;
attr sanofi_EIA ordinal ;
attr culture perfect ;
know  2  * 2
 0 : 1 0 :
 1 : 0 1 :
 ;

end;