- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
Hello everyone,
This is my final survival model for a bird breeding colony, and I would take the residuals. I know that in proc mixed I must write "/outp=" just behind the model, but I do not know what should I use in proc glimmix.
proc glimmix data=Survival INFOCRIT=PQ noreml;
class treartment attempt year nest ring;
model Survival = treatment Year attempt Treatment*attempt size/ dist= binomial link=logit solution cl;
random Year/ subject=ring(nest) type=Cs residual;
lsmeans treatment;
run;
Thanks very much for your help.
Best wishes,
Jaime
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
The documentation for PROC GLIMMIX is very helpful. It says you should add an OUTPUT statement
Paige Miller
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
What are you going to do with the residuals from a binomial distribution? Recall that there is no "error term" outside of the model formulation for a generalized linear mixed model that employs a binomial distribution, as there would be for models employing a normal or lognormal distribution. Original scale "residuals" will consist of a set of negative ones, zeroes, and positive ones. These are good for calculating concordances or ROC curves, if that is where you are going.
Steve Denham
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
Thank you very much!
In case of not getting the model residuals, I would not integrating the effect of "size". Then... Would it be correct to use only raw data to examine the results graphically without taking into account other significant covariates?
Best whishes,
Jaime