BookmarkSubscribeRSS Feed
yahoo0806
Calcite | Level 5

Hi all,

I recent run into a problem when I use proc glimmix, hope someone can help.

 

the data I have is a cross over data, each subject can have treatment sequence as AB or BA, my output is a binary responder: depending on if the change from baseline (CHG) value is >= 3 or not. For example, if the CHG value is >= 3 then responder variable = Y otherwise response = N.

 

The response status (responder varaible) will be analyzed using a generalized linear mixed model for binomial data (or mixed effects logistic regression model) with the logit of the probability of response as the dependent variable, treatment and period as fixed effects and subject as random effect. The Kenward-Rogers adjustment will be used to compute the denominator degrees of freedom for the test of the fixed effect. Adaptive Gauss-Hermite quadrature will be used to approximate the marginal log likelihood (SAS procedure GLIMMIX and METHOD=QUAD).

 

Below is the model I used:

 

Proc glimmix data=model METHOD=QUAD;
class subject period treatment;
model responder = period treatment/ddfm=KENWARDROGER solution link=logit;
random subject;
lsmeans treatment / cl ilink;
run;

 

The problem I run into is SAS gives me the following error message:

ERROR: For the multinomial distribution only the following links are possible: CUMLOGIT, CUMPROBIT, CUMLOGLOG, CUMCLL, GLOGIT.

 

Do you know why? and can you please let me know how I should fix my model?

 

Many thanks in advance!

4 REPLIES 4
Reeza
Super User

I think SAS see's the "responder" variable as 'multinomial', ie it has more than two responses for some reason. Run a PROC FREQ and check the results. 

 

proc freq data=model;
table responder / missing;
run;

@yahoo0806 wrote:

Hi all,

I recent run into a problem when I use proc glimmix, hope someone can help.

 

the data I have is a cross over data, each subject can have treatment sequence as AB or BA, my output is a binary responder: depending on if the change from baseline (CHG) value is >= 3 or not. For example, if the CHG value is >= 3 then responder variable = Y otherwise response = N.

 

The response status (responder varaible) will be analyzed using a generalized linear mixed model for binomial data (or mixed effects logistic regression model) with the logit of the probability of response as the dependent variable, treatment and period as fixed effects and subject as random effect. The Kenward-Rogers adjustment will be used to compute the denominator degrees of freedom for the test of the fixed effect. Adaptive Gauss-Hermite quadrature will be used to approximate the marginal log likelihood (SAS procedure GLIMMIX and METHOD=QUAD).

 

Below is the model I used:

 

Proc glimmix data=model METHOD=QUAD;
class subject period treatment;
model responder = period treatment/ddfm=KENWARDROGER solution link=logit;
random subject;
lsmeans treatment / cl ilink;
run;

 

The problem I run into is SAS gives me the following error message:

ERROR: For the multinomial distribution only the following links are possible: CUMLOGIT, CUMPROBIT, CUMLOGLOG, CUMCLL, GLOGIT.

 

Do you know why? and can you please let me know how I should fix my model?

 

Many thanks in advance!


 

yahoo0806
Calcite | Level 5

I did the proc freq as you suggested, but I do see only two results: Y and N.

 

Reeza
Super User

Do you have a format applied to the variable?

 


@yahoo0806 wrote:

I did the proc freq as you suggested, but I do see only two results: Y and N.

 


 

yahoo0806
Calcite | Level 5

I don't think so, here is how I created the responder variable:

data model; set data;
if CHG >=3 then responder = "Y"; else responder = "N";
run;

 

 

sas-innovate-2024.png

 

Secure your spot at the must-attend AI and analytics event of 2024: SAS Innovate 2024! Get ready for a jam-packed agenda featuring workshops, super demos, breakout sessions, roundtables, inspiring keynotes and incredible networking events.

 

Register by March 1 to snag the Early Bird rate of just $695! Don't miss out on this exclusive offer. 

 

Register now!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.

Discussion stats
  • 4 replies
  • 1437 views
  • 0 likes
  • 2 in conversation