BookmarkSubscribeRSS Feed
🔒 This topic is solved and locked. Need further help from the community? Please sign in and ask a new question.
boriquez
Calcite | Level 5

I would like to analyze a multi-location complete randomized block design  where resposne variable is days to maturity (of a crop).  After running a proc severity the weibull distribution was selected.  My questions are :

 

Is there any information about using the weibull distribution in the glimmix procedure, if so which is the appropriate link function?

 

Do I need a fully saturated model with weibull distribution  as in a binomial distribution?

 

 

1 ACCEPTED SOLUTION

Accepted Solutions
Rick_SAS
SAS Super FREQ

I don't think a univariate histogram/analysis of the response variable tells you very much. The DIST= option in GENMOD and GLIMMIX specify the (conditional) error distribution, which might be very different from the distribution of the response. See the pictures in the article "Error distributions and exponential regression models."  

 

To answer your implied question, the GLIMMIX procedure does not support a Weibull distribution function for the errors. It sounds like you are treating the response variable as a continuous variable. If you think you have long-tailed errors, you might look at the gamma or lognormal distributions. (The gamma is sometimes used for waiting times.) Look at the goodness-of-fit statistics and the diagnostic residual plots to guide you in choosing a good model.

View solution in original post

5 REPLIES 5
Rick_SAS
SAS Super FREQ

You say "after running PROC SEVERITY, a Weibull distribution was selected." Did you use PROC SEVERITY to run a fixed-model regression, or did you just analyze the response variable in isolation (a univariate model)?

boriquez
Calcite | Level 5
In Isolation (Univariate model)




Rick_SAS
SAS Super FREQ

I don't think a univariate histogram/analysis of the response variable tells you very much. The DIST= option in GENMOD and GLIMMIX specify the (conditional) error distribution, which might be very different from the distribution of the response. See the pictures in the article "Error distributions and exponential regression models."  

 

To answer your implied question, the GLIMMIX procedure does not support a Weibull distribution function for the errors. It sounds like you are treating the response variable as a continuous variable. If you think you have long-tailed errors, you might look at the gamma or lognormal distributions. (The gamma is sometimes used for waiting times.) Look at the goodness-of-fit statistics and the diagnostic residual plots to guide you in choosing a good model.

boriquez
Calcite | Level 5

Thank you for taking the time  to help Rick.

 

boriquez
Calcite | Level 5

Once again thank you Rick,  for pointing me in the right direction.  I will take a look at your sugestions and explore both the lognormal and gamma distribution. 

SAS Innovate 2025: Call for Content

Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!

Submit your idea!

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
  • 5 replies
  • 1988 views
  • 2 likes
  • 2 in conversation