spss='abortion.sav';

model:'A. binomial count growth model'  regression binomial    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
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind year religion;

end;
model:'B. Model A + 2-class random intercept'  regression binomial    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=1676357  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind year religion;

end;
model:'C. Model A + 3-class random intercept'  regression binomial    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=1040358  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind year religion;

end;
model:'D. Model A + 4-class random intercept'  regression binomial    4 /
  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=391953  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind year religion;

end;
model:'E. Model A + 5-class random intercept'  regression binomial    5 /
  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=1994012  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind year religion;

end;
model:'F. Model A + 2-class random intercept and slope'  regression binomial    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=672669  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind religion;

end;
model:'G. Model A + 3-class random intercept and slope'  regression binomial    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=449907  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind religion;

end;
model:'H. Model A + 4-class random intercept and slope'  regression binomial    4 /
  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=317994  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind religion;

end;
model:'I. Model A + 5-class random intercept and slope'  regression binomial    5 /
  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=2676600  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind religion;

end;
model:'G. Model A + parametric random intercept '  regression binomial    1 /
  cfn=1  cfnde=20  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=398465  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
cfactor 1 <- response;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind year religion;

end;

model:'G. Model A + parametric random intercept and slope'  regression binomial    1 /
  cfn=2  cfnde=20  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=398465  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent response;
exposure exposure;
replicate id;
predictor year religion;
cfactor 1 <- response;
cfactor 2 <- response year;
attr year nominal ;
attr religion nominal ;
attr response binomial ;
classind year religion;

end;