turn on suggestions

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

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Proc Glimmix model resulting in all zero standard ...

Topic Options

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

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

04-28-2017 10:52 AM

Hi.

First of all, I do not consider myself being a stats-guru and my "extended" SAS and stats-knowledge is mostly self-tought. My question deals with the SAS procedure GLIMMIX and, more precisely, the generated output table.

My research question is related to animal sciences and with my statistical analysis I am trying to identify how different feeds and ingredients may interact with one another affect animal performance (my response variable).

I thought of analyzing my data using a full factorial GLMM. I am having 2 class variables and I came up with the following model:

ods graphics on;

proc glimmix

Data=Data2 plots =pearsonpanel;

class RefDiet TestIngr;

model ADCprotI=RefDiet TestIngr RefDiet*TestIngr / dist=beta

link=logit solution;

lsmeans RefDiet*TestIngr / adjust=bon ilink oddsratio;

run;

ods graphics off;

The beta distribution was chosen according to the residual panel, as it seems to provide a better fit and be more appropriate than a gaussian-shaped distribution (response variables are expressed as proportions with values between 0 and 1). With one specific response variable I am now getting standard errors of 0 for all class variables, resulting in infinite t values. As a result, all class variables turn out to likely have a highly significant effect on my response variable. However, something seems wrong here, but I am currently stumped... Does anyone has a good suggestion on how to overcome this issue or am I doing something completely wrong here? Please let me know if you need more information + thanks a lot!

The result for e.g. fixed effects looks like this (lines should be labelled with: RefDiet TestIngr RefDiet*TestIngr)

Type III Tests of Fixed Effects Effect

Num DF Den DF F Value Pr > F

1 | 16 | Infty | <.0001 |

3 | 16 | Infty | <.0001 |

3 | 16 | Infty | <.0001 |

And for the parameter estimates (just the upper lines, the shown line names are: RefDiet 1 RefDiet 2 TestIngr 10 TestIngr 20 TestIngr 30...)

Parameter Estimates Effect RefDiet TestIngr Estimate Standard

Error DF t Value Pr > |t|

Error DF t Value Pr > |t|

2.1766 | 0 | 16 | Infty | <.0001 |

2.4569 | 0 | 16 | Infty | <.0001 |

-0.4995 | 0 | 16 | -Infty | <.0001 |

-0.1410 | 0 | 16 | -Infty | <.0001 |

-0.01848 | 0 | 16 | -Infty | <.0001 |

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

04-28-2017 11:03 AM

You get infinity for a t-test value when there is zero variability in the residuals. In other words, the model has fit the data exactly.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

04-28-2017 03:01 PM

Thank you very much for your quick reply, Paige. In other words: would you consider it legitimate to report this outcome "as is" (standard-error = 0 and F-value infinity)? Of course this makes running a multiple comparison test a bit difficult (impossible)...

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

04-28-2017 03:25 PM

I'm sure there are some real-world situations where this is a legitimate outcome, but there are also some real-world situations where this is not a legitimate outcome.

The zero error situation could arise from mis-specifiying (actually over-specifying) the model (i.e. statistical analysis error), it could also arise from degenerate data (perhaps some error in data collection or data transmission or error in entering the data into the database or computer coding error).

It could be caused by a many different things.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

04-28-2017 11:47 PM

Why not use DIST=BINOMIAL since your link function is LOGIT ?

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

04-29-2017 01:11 PM

Your model is likely overspecified (too many terms). Also note (very important): the beta distribution is defined for 0 < y < 1. Thus, 0s and 1s are not allowed. All proportions equal to 0 or 1 become missing in the analysis. If you don't have 0s and 1s, then you are fine, but otherwise you are not using all your data. Be careful.

If your proportions are for a continuous variable, then binomial is not appropriate.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

05-02-2017 02:34 PM

Apologies for the late reply and thanks for the notification. Indeed, binomial would not be appropriate and, likely, I should go with a beta (proportions are for a continuous variable). The model stated above seemed to work for a bunch of variables, but is failing for the one listed above...