libname in 'e:\addir\compare\threelev\mixed\sparse';

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

%glimmix(data=in.simspar,
             procopt= noclprint covtest method=ml,
             stmts=%str(class docnum patnum;
                        model y3=foreign male disc / solution;
                        random int / subject=docnum;
                        random int / subject=patnum(docnum);
                        parms (0) (0) (1) /eqcons=3;
                        ),
             error=binomial,
             options=mql
             )
run;