data='sids.dat';

model:'A. 1-class'  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=0x1d17
;
dependent '#DEATHS';
exposure POPSIZE;
attr '#DEATHS' poisson ;

end;
model:'B. 2-class'  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=424820  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent '#DEATHS';
exposure POPSIZE;
attr '#DEATHS' poisson ;

end;
model:'C. 3-class'  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=239468  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent '#DEATHS';
exposure POPSIZE;
attr '#DEATHS' poisson ;

end;
model:'D. 4-class (=NPML estimator)'  regression poisson    4 /
  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=424820  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent '#DEATHS';
exposure POPSIZE;
attr '#DEATHS' poisson ;

end;