Date:         Wed, 24 Jan 2007 11:01:54 -0500
Reply-To:     Sheetal Kamble <skamble@GMAIL.COM>
Sender:       "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU>
From:         Sheetal Kamble <skamble@GMAIL.COM>
Subject:      Re: PROC GLIMMIX - cross classified random effects - query
Comments: To: Dale McLerran <stringplayer_2@yahoo.com>
In-Reply-To:  <688115.9250.qm@web32205.mail.mud.yahoo.com>
Content-Type: text/plain; charset=ISO-8859-1; format=flowed

Thank you, Dale. I tried the following syntax-

 proc glimmix data = one;

class  female agevar racevar insurance ch_ren_fail copd cvd chf diabet sAMI
pAMI hospadm hv4 pv6 MD Hospital;

 model hem (descending) = female agevar racevar insurance ch_ren_fail copd
cvd chf diabet pAMI sAMI hospadm hv4|pv6 / dist= binary link= logit s chisq
CL ODDSRATIO;

random  hospital md_operating  ;

run;

However, the model did not converge. I am not sure what is wrong with the
syntax.

Could you or anyone on this listserve guide me with the same?

Thanks,
Shital


On 1/23/07, Dale McLerran <stringplayer_2@yahoo.com> wrote:
>
> --- Shital Kamble <skamble@GMAIL.COM> wrote:
>
> > Hi,
> >
> > I am using proc glimmix procedure in SAS 9.1 to model the interaction
> > effect of physician volume (pv6) and hospital volume (hv4) on adverse
> > dichotomous outcome "hem". The data structure is cross classified
> > (patients
> > nested within hospitals, patients nested within physicians, but
> > approximately 25% of physicians in my dataset are not nested within
> > hospitals). I am using proc glimmix syntax for cross classified
> > random
> > effects, however, my model does not converge for some reason. I would
> > like
> > to seek your advise pertaining to the same.
> >
> > Following is the proc glimmix syntax-
> >
> >
> > proc glimmix data = one;
> >
> > class hem female agevar racevar insurance ch_ren_fail copd cvd chf
> > diabet
> > sAMI pAMI hospadm hv4 pv6 MD Hospital;
> >
> > model hem (descending) = female agevar racevar insurance ch_ren_fail
> > copd
> > cvd chf diabet pAMI sAMI hospadm hv4|pv6 / dist= binary link= logit s
> > chisq
> > CL ODDSRATIO;
> >
> > random int hospital md_operating / s ;
> >
> > run;
> >
> > *Note- hem = 0,1 (dichotomous variable) ; Primary regressor -
> > interaction
> > of physician vol (pv6 = 0=low, 1=high) and hospital volume (hv4 =
> > 0=low,
> > 1=high); hospital = hospital id and MD = physician id. All other
> > independent variables are dichotomous or categorical variables.
> >
> > Please guide me peratining the correct syntax for this SAS procedure.
> >
> > Thanks,
> > Shital
> >
> > Shital Kamble
> > Ph.D. Student
> >
>
> Shital,
>
> Generally, an intercept term is specified as a random effect only
> when you also specify a SUBJECT on the RANDOM statement.  You have
> not named a SUBJECT for your RANDOM statement (and rightly so).
>
> You might take a look at Example 2: Mating Experiment with Crossed
> Random Effects in the GLIMMIX documentation.  Note that their
> RANDOM statement is
>
>    RANDOM FPOP*FNUM MPOP*MNUM;
>
> and that FPOP*FNUM exactly specifies a single female salamander -
> salamander FNUM in either the FPOP='rb' or FPOP='ws' salamander
> population.  Similarly, mpop*mnum exactly specifies a single male
> salamander.  This is, then, analogous to specifying
>
>    RANDOM HOSPITAL MD_OPERATING;
>
>
> HTH,
>
> Dale
>
>
> ---------------------------------------
> Dale McLerran
> Fred Hutchinson Cancer Research Center
> mailto: dmclerra@NO_SPAMfhcrc.org
> Ph:  (206) 667-2926
> Fax: (206) 667-5977
> ---------------------------------------
>
>
>
>
> ____________________________________________________________________________________
> Sucker-punch spam with award-winning protection.
> Try the free Yahoo! Mail Beta.
> http://advision.webevents.yahoo.com/mailbeta/features_spam.html
>