BookmarkSubscribeRSS Feed
superbibi
Obsidian | Level 7

Hi SAS users,

 

I have a data set with binary repeated measure outcome.

 

I want to know whether the incidence of GSRSINCI6 does not change over visits. Also, I want to know the difference (p-value) of each visit compared to visit 2. 

 

My question is whether there is a way that I can adjust the p-value for multiple between visit comparison in proc glimmix.

 

My code and data set have been uploaded.

 

Thanks. 

 

PROC GLIMMIX DATA = GSRS;
CLASS SUBNUM VISIT (REF = '2') ;
MODEL GSRSINCI6(EVENT ='1') = VISIT/DIST= BINARY SOLUTION CHISQ ODDSRATIO ;
RANDOM _RESIDUAL_ /SUBJECT = SUBNUM TYPE = VC ;
RUN;

3 REPLIES 3
PGStats
Opal | Level 21

I guess you may want:

 

PROC GLIMMIX DATA = sasforum.bin;
CLASS ID VISIT (REF = '2') ;
MODEL GSRSINCI6(EVENT ='1') = VISIT/DIST= BINARY SOLUTION CHISQ ODDSRATIO ;
RANDOM _RESIDUAL_ /SUBJECT = ID TYPE = VC ;
lsmeans visit / diff=control('2') adjust=dunnett;
RUN;
PG
superbibi
Obsidian | Level 7
Thank you. But since the data is binary, is it appropriate to use lsmeans?
PGStats
Opal | Level 21

Oops, I'm sorry. I guess not.

 

You can always adjust multiple p-values (or a subset of them) with proc multtest after exporting with ods ouitput.

PG

Ready to join fellow brilliant minds for the SAS Hackathon?

Build your skills. Make connections. Enjoy creative freedom. Maybe change the world. Registration is now open through August 30th. Visit the SAS Hackathon homepage.

Register today!
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
  • 3 replies
  • 647 views
  • 0 likes
  • 2 in conversation