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

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- PROC MIANALYZE procedure after PROC GLIMMIX binary outcome (NO odds ra...

Options

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

🔒 This topic is **solved** and **locked**.
Need further help from the community? Please
sign in and ask a **new** question.

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

Posted 11-29-2019 08:44 AM
(2545 views)

Hi SAS users,

I used Multiple imputation with GLIMMIX for binary outcome variable with one binary independent variable (as example). When I ran PROC MIANALYZE, the results only showed the estimate and SE but did not show 95% CL nor the significance (attach pic of the table). Also did not provides the odds ratio table.

I used the following code:

PROC GLIMMIX DATA=shis.shis_mi3 METHOD=QUAD empirical=classical;

CLASS rgn_code SEX1 (ref=FIRST);

MODEL PA (EVENT=LAST)= SEX1 /CL DIST=BINARY LINK=LOGIT SOLUTION

ODDSRATIO (DIFF=LAST LABEL);

RANDOM INTERCEPT / SUBJECT=rgn_code S CL TYPE=VC;

by imputation_;

nloptions gconv=0;

ods output ParameterEstimates=PAparm1;

COVTEST /WALD;

RUN;

/* Prepare pooled estimates for MI (excluding original unimputed data) */

data imputed;

set PAparm1;

_imputation_=imputation_;

If _imputation_>0;

/* Calculate and print pooled estimates for MI */

proc mianalyze parms = imputed;

CLASS SEX1;

modeleffects Intercept SEX1;

run;

1 ACCEPTED SOLUTION

Accepted Solutions

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

You likely have a WARNING message in your LOG regarding the between imputation variance being zero. This essentially means that the estimates for each of the imputations is identical. When the between imputation variance is zero then the number of Degrees of Freedom is undefined so you cannot get a confidence interval or p-values.

Unfortunately, there is not a good approach to take in this case. This is one of the limitations of multiple imputation in general. You might try doing a literature search to see if you can find any references that deal with this issue.

Because GLIMMIX does not report a standard error for the odds ratio directly, you would not be able to combine those either. I suspect you would run into the same issue with the between imputation variance anyway.

Normally you would combine the log (OR), that is the parameter estimates, and then exponentiate those values in a subsequent data step.

2 REPLIES 2

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

You likely have a WARNING message in your LOG regarding the between imputation variance being zero. This essentially means that the estimates for each of the imputations is identical. When the between imputation variance is zero then the number of Degrees of Freedom is undefined so you cannot get a confidence interval or p-values.

Unfortunately, there is not a good approach to take in this case. This is one of the limitations of multiple imputation in general. You might try doing a literature search to see if you can find any references that deal with this issue.

Because GLIMMIX does not report a standard error for the odds ratio directly, you would not be able to combine those either. I suspect you would run into the same issue with the between imputation variance anyway.

Normally you would combine the log (OR), that is the parameter estimates, and then exponentiate those values in a subsequent data step.

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

Well explained Rob. Thanks.

**Available on demand!**

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

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.