Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Problem with my glimmix model

Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 08-23-2018 02:41 PM
(1436 views)

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

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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!

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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;

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.

** **

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.