libname in 'e:\multilev\compare\threelev\mixed\rc';

%INCLUDE 'c:\macros\glmm612.sas' / nosource;

/* ddfm=bw options changes default denominator degrees of freedom, reducing  
   amount of memory needed 
   
   parms (1) (0) (1) (1) (0) (1) (1) / eqcons = 7
         lowerb=0.1,.,0.1,0.1,.,0.1 accomplishes the following
   
   The parms portion specifies the starting values for the covariance matrix    
   and dispersion parameter
   
   The eqcons=7 contstrains the seventh parameter listed in the parms statement, 
   here the dispersion parameter to be equal to the starting value
   
   The lowerb options specifies lower bounds for the covariance parameters.

*/

%glimmix(data=in.simrc,
             procopt= noclprint covtest update method=ml,
             stmts=%str(class docnum patnum;
                        model y3=foreign disc agecat / ddfm=bw solution;
                        random int agecat / subject=docnum type=un g;
                        random int disc / subject=patnum(docnum) type=un;
                        parms (1) (0) (1) (1) (0) (1) (1) / eqcons=7
                            lowerb=0.1,.,0.1,0.1,.,0.1;
                        ),
             error=binomial,
             options=pql noprev
             )
run;