spss='dmft.SAV';

model:'A. Poisson distribution'  regression poisson   1 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=0  bayeslat=0
  bayespoi=0
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=0  nseedboot=0
  nrand=10
  usemiss=No
  sewald=yes
 outsect=0x1c1f
;
dependent dmft;
exposure 1;
attr dmft poisson ;

end;
model:'B. 2-class Poisson mixture (=NPML)'  regression poisson   2 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=0  bayeslat=0
  bayespoi=0
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=23106  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c1f
;
dependent dmft;
exposure 1;
attr dmft poisson ;

end;
model:'C. 3-class Poisson mixture'  regression poisson   3 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=0  bayeslat=0
  bayespoi=0
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=502238  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c1f
;
dependent dmft;
exposure 1;
attr dmft poisson ;

end;
model:'D. zeroinflated Poisson distribution'  regression poisson zeroinfl  1 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=0  bayeslat=0
  bayespoi=0
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=2588363  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c1f
;
dependent dmft;
exposure 1;
attr dmft poisson ;

end;
model:'E. normal random effect Poisson distribution'  regression poisson   1 /
  cfn=1  cfnde=10  cfcsp=0
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=0  bayeslat=0
  bayespoi=0
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=318560  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c1f
;
dependent dmft;
exposure 1;
cfactor 1 <- dmft;
attr dmft poisson ;

end;
model:'F. Poisson regression model'  regression poisson   1 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=0  bayeslat=0
  bayespoi=0
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=156462  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c1f
;
dependent dmft;
exposure 1;
predictor sex ethnic treatmnt;
attr dmft poisson ;
attr sex nominal ;
attr ethnic nominal ;
attr treatmnt nominal ;
classind sex ethnic treatmnt;

end;
model:'G. Poisson regression model with 2-class random intercept'  regression poisson   2 /
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=0  bayeslat=0
  bayespoi=0
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=647977  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c1f
;
dependent dmft;
exposure 1;
predictor sex ethnic treatmnt;
attr dmft poisson ;
attr sex nominal ;
attr ethnic nominal ;
attr treatmnt nominal ;
classind sex ethnic treatmnt;

end;
model:'H. Poisson regression model with normal random intercept'  regression poisson   1 /
  cfn=1  cfnde=10  cfcsp=0
  toler=1e-008  tolem=0.01  tolran=1e-005  bayes=0  bayess2=0  bayeslat=0
  bayespoi=0
  iterem=250  iternr=50  itersv=50
  iterboot=500
  nseed=2050829  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1c1f
;
dependent dmft;
exposure 1;
predictor sex ethnic treatmnt;
cfactor 1 <- dmft;
attr dmft poisson ;
attr sex nominal ;
attr ethnic nominal ;
attr treatmnt nominal ;
classind sex ethnic treatmnt;

end;