data='socatt.TXT';

model:'A. standard binomial count logistic regression'  regression binomial    1 /
  mmncfac=0  cfnde=10  mmnclass=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 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;
model:'B. nonparametric random intercept at level 2 (2 classes)'  regression binomial    2 /
  mmncfac=0  cfnde=10  mmnclass=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=1681905  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;
model:'C. nonparametric random intercept at level 2 (3 classes)'  regression binomial    3 /
  mmncfac=0  cfnde=10  mmnclass=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=342731  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;
model:'D. nonparametric random intercept at level 2 (4 classes)'  regression binomial    4 /
  mmncfac=0  cfnde=10  mmnclass=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=746172  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;
model:'E. nonparametric random intercept at level 2 (5 classes)'  regression binomial    5 /
  mmncfac=0  cfnde=10  mmnclass=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=2192496  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;
model:'F. nonparametric random intercept at level 2 and 3 (4 classes and 2 GClasses)'  regression binomial    4 /
  mmncfac=0  cfnde=10  mmnclass=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=1716242  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;
model:'G. parametric random intercept at level 2 (CFactor)'  regression binomial    1 /
  cfn=1  cfnde=20  cfcsp=0
  mmncfac=0  cfnde=20  mmnclass=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=95993  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
cfactor 1 <- Response;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;
model:'H. parametric random intercept at level 2 and 3 (CFactor and GCFactor)'  regression binomial    1 /
  cfn=1  cfnde=20  cfcsp=0
  mmncfac=1  cfnde=20  mmnclass=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=162382  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
gcfactor 1 <- __intercept__;
cfactor 1 <- Response;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;
model:'I. Model F + class-specific time effect'  regression binomial    4 /
  mmncfac=0  cfnde=10  mmnclass=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=2602348  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Religion;

end;
model:'J. Model F + GClass-specific religion effect'  regression binomial    4 /
  mmncfac=0  cfnde=10  mmnclass=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=378300  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__ Religion;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;
model:'K. Model F + GClasses also affecting Classes'  regression binomial    4 /
  mmncfac=0  cfnde=10  mmnclass=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=2519255  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
mlx  0 <- __intercept__;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;
model:'L. Model F + interaction between Classes and GClasses'  regression binomial    4 /
  mmncfac=0  cfnde=10  mmnclass=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=1353774  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
dependent Response;
exposure 7;
replicate Case_ID;
predictor Year Religion;
gclass  __intercept__;
mlrestrict __intercept__ : 0  1 2 3 4;
groupid District_ID;
attr Year nominal ;
attr Response binomial ;
attr Religion nominal ;
classind Year Religion;

end;