data='jobchange.dat';

model:'A. piecewise exponential model with 3 time categories'  regression poisson 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
;
casewgt FREQUENCY;
dependent Y_VARIABLE;
exposure EXPOSURE;
replicate ID;
predictor TIME_CATEGORY FIRMSIZE DUMMY_GOVERNEMENT;
attr FIRMSIZE ordinal ;
attr DUMMY_GOVERNEMENT ordinal ;
attr TIME_CATEGORY nominal ;
attr Y_VARIABLE poisson ;
classind FIRMSIZE DUMMY_GOVERNEMENT TIME_CATEGORY;

end;

model:'B. 2-class version of A, with class-specfic intercept'  regression poisson    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=30908  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
casewgt FREQUENCY;
dependent Y_VARIABLE;
exposure EXPOSURE;
replicate ID;
predictor TIME_CATEGORY FIRMSIZE DUMMY_GOVERNEMENT;
attr FIRMSIZE ordinal ;
attr DUMMY_GOVERNEMENT ordinal ;
attr TIME_CATEGORY nominal ;
attr Y_VARIABLE poisson ;
classind FIRMSIZE DUMMY_GOVERNEMENT TIME_CATEGORY;

model:'C. 3-class version of A, with class-specfic intercept'  regression poisson   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=343021  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
casewgt FREQUENCY;
dependent Y_VARIABLE;
exposure EXPOSURE;
replicate ID;
predictor TIME_CATEGORY FIRMSIZE DUMMY_GOVERNEMENT;
attr FIRMSIZE ordinal ;
attr DUMMY_GOVERNEMENT ordinal ;
attr TIME_CATEGORY nominal ;
attr Y_VARIABLE poisson ;
classind FIRMSIZE DUMMY_GOVERNEMENT TIME_CATEGORY;

end;

model:'D. Model A + random intercept (CFactor)'  regression poisson    1 /

  cfn=1  cfnde=10  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=2044344  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
casewgt FREQUENCY;
dependent Y_VARIABLE;
exposure EXPOSURE;
replicate ID;
predictor TIME_CATEGORY FIRMSIZE DUMMY_GOVERNEMENT;
cfactor 1 <- Y_VARIABLE;
attr FIRMSIZE ordinal ;
attr DUMMY_GOVERNEMENT ordinal ;
attr TIME_CATEGORY nominal ;
attr Y_VARIABLE poisson ;
classind FIRMSIZE DUMMY_GOVERNEMENT TIME_CATEGORY;

end;
model:'D. Model B, but with FirmSize as covariate'  regression poisson    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=2768034  nseedboot=0
  nrand=0
  usemiss=No
  sewald=yes
 outsect=0x1d17
;
casewgt FREQUENCY;
dependent Y_VARIABLE;
exposure EXPOSURE;
replicate ID;
predictor TIME_CATEGORY;
covariate FIRMSIZE DUMMY_GOVERNEMENT;
attr FIRMSIZE ordinal ;
attr DUMMY_GOVERNEMENT ordinal ;
attr TIME_CATEGORY nominal ;
attr Y_VARIABLE poisson ;
classind FIRMSIZE DUMMY_GOVERNEMENT TIME_CATEGORY;

end;