topic Re: How to specify event of outcome variable in Proc GLIMMIX in SAS Programming

How to specify event of outcome variable in Proc GLIMMIX

UcheOkoro — Thu, 11 Aug 2022 16:32:46 GMT

Please, I need help specifying the direction of the event in proc GLIMMIX for a longitudinal study. I used the (event='1') statement but when I flipped it the results were the same. It appears that the 'event' statement does not work. The outcome is binomial.

The following is my SAS code

PROC GLIMMIX DATA=COVERED.Demographics method=laplace;

CLASS  record_id wkly_contact_simsx(ref='0');

MODEL hospital_adherence(event='1') =  wkly_contact_simsx
/SOLUTION DIST=bin LINK=Logit ;

random int / subject=record_id;

lsmeans  wkly_contact_simsx/  diff e cl oddsratio ;
  run;

Thank you so much for your assistance.

Re: How to specify event of outcome variable in Proc GLIMMIX

Reeza — Thu, 11 Aug 2022 16:39:34 GMT

Can you please show the log from each run?

Does your variable

hospital_adherence

have a format attached to it?

Re: How to specify event of outcome variable in Proc GLIMMIX

UcheOkoro — Thu, 11 Aug 2022 16:47:31 GMT

Thank you for your prompt response the following are the log and results.

Thank you again

Re: How to specify event of outcome variable in Proc GLIMMIX

Reeza — Thu, 11 Aug 2022 17:02:11 GMT

Well, the warning does say it's ignoring the EVENT statement because the outcome isn't binary.

Can you show an proc freq of the hospital_adherence variable?

proc freq data=covered.demographics;
table hospital_adherence*wkly_gathering / missing;
run;

@UcheOkoro wrote:

Thank you for your prompt response the following are the log and results.

Thank you again

Re: How to specify event of outcome variable in Proc GLIMMIX

UcheOkoro — Thu, 11 Aug 2022 17:07:33 GMT

Here it is . Thank you

Re: How to specify event of outcome variable in Proc GLIMMIX

Reeza — Thu, 11 Aug 2022 17:23:26 GMT

Ok, it looks you have some missing values across your outcome and fixed effect. Try removing those and re-running your proc to see what happens.

I would expect the warning to disappear and if not that's part of the problem.

Re: How to specify event of outcome variable in Proc GLIMMIX

UcheOkoro — Thu, 11 Aug 2022 17:28:33 GMT

Thank you for your response. I deleted the missing observations but the results still remained the same

Re: How to specify event of outcome variable in Proc GLIMMIX

PaigeMiller — Thu, 11 Aug 2022 17:35:59 GMT

GLIMMIX produces output of Class Level Information and Number of Observations like this:

Please show this to us from running PROC GLIMMIX using your data set.

Re: How to specify event of outcome variable in Proc GLIMMIX

Reeza — Thu, 11 Aug 2022 17:37:39 GMT

Is there a warning in the log still?

Re: How to specify event of outcome variable in Proc GLIMMIX

UcheOkoro — Thu, 11 Aug 2022 18:20:03 GMT

Re: How to specify event of outcome variable in Proc GLIMMIX

Reeza — Thu, 11 Aug 2022 19:21:07 GMT

I don't know enough about proc glimmix to be 100% certain but I'm guess it has to do with the subject being each record id.

Is each record unique with a single value of adherence?

Re: How to specify event of outcome variable in Proc GLIMMIX

Ksharp — Fri, 12 Aug 2022 11:47:54 GMT

It should be "dist=binary " ,not "dist=bin " . "dist=bin" would make sas take it as "dist=binomial"