BookmarkSubscribeRSS Feed
UcheOkoro
Lapis Lazuli | Level 10

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.

 

11 REPLIES 11
Reeza
Super User

Can you please show the log from each run?

 

Does your variable

hospital_adherence

have a format attached to it?

 

 

UcheOkoro
Lapis Lazuli | Level 10

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

 

Thank you again

 

UcheOkoro_0-1660236338563.png

 

 

UcheOkoro_1-1660236373297.png

 

UcheOkoro_2-1660236400563.png

 

 

Reeza
Super User

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

 

UcheOkoro_0-1660236338563.png

 

 

UcheOkoro_1-1660236373297.png

 

UcheOkoro_2-1660236400563.png

 

 


 

UcheOkoro
Lapis Lazuli | Level 10

Here it is . Thank you

 

UcheOkoro_0-1660237636807.png

 

Reeza
Super User

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. 

 

 

UcheOkoro
Lapis Lazuli | Level 10

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

PaigeMiller
Diamond | Level 26

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

 

PaigeMiller_0-1660239297366.png

 

 

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

--
Paige Miller
Reeza
Super User

Is there a warning in the log still?

UcheOkoro
Lapis Lazuli | Level 10

UcheOkoro_1-1660241794764.png

 

UcheOkoro_2-1660241953341.png

UcheOkoro_3-1660241967800.png

UcheOkoro_4-1660241988091.png

 

 

 

Reeza
Super User

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?

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

sas-innovate-2024.png

Available on demand!

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

 

Register now!

How to Concatenate Values

Learn how use the CAT functions in SAS to join values from multiple variables into a single value.

Find more tutorials on the SAS Users YouTube channel.

Click image to register for webinarClick image to register for webinar

Classroom Training Available!

Select SAS Training centers are offering in-person courses. View upcoming courses for:

View all other training opportunities.

Discussion stats
  • 11 replies
  • 994 views
  • 8 likes
  • 4 in conversation