spss='heinen3reg.sav';

model:'Rating Scale model with sum delta(t) = 0'  regression ordinal    1 /
  cfn=1  cfnde=10  cfcsp=0
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=1  bayess2=1  bayeslat=1
  bayespoi=0
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=0  nseedboot=0
  nrand=10
  usemiss=No
  sewald=yes  dummy=no
 outsect=0x1c1f
;
dependent response;
replicate id;
predictor itemnr;
cfactor 1 <- response;
restrict __intercept__ : 1  1;
attr itemnr nominal ;
attr response ordinal ;

end;

model:'Mixture Rating Scale model with sum delta(t)=0'  regression ordinal    2 /
  cfn=1  cfnde=10  cfcsp=0
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=1  bayess2=1  bayeslat=1
  bayespoi=0
  iterem=250  iternr=50  itersv=100
  iterboot=500
  nseed=0  nseedboot=0
  nrand=20
  usemiss=No
  sewald=yes  dummy=no
 outsect=0x1c1f
;
dependent response;
replicate id;
predictor itemnr;
cfactor 1 <- response;
restrict __intercept__ : 1  1 2;
attr itemnr nominal ;
attr response ordinal ;

end;