Hi,
The QQ Plot suggests that the studentized residuals of the response variables are not normally distributed.
When I run it with a poisson distributions (or others) it is outputting completely unfeasible estimates. It should be around 0.20 (which would be 20%) and it is estimating negative numbers.
I thought could be scaling issues (multiplied by 100) but it is not.
I have tried different estimation methods too (LaPlace, Quadrature). Not sure if it could be the linkfunction?
Any thoughts?
Thanks,
Marcio
You have to be more specific. YOu have to put ilink in the lsmeans statement and in the plots opiton (on the proc statement) to see the inverse link means. Even with ilink option on lsmeans, you first get the logs, then on the right the inverse link (data scale means). Your example code does not give ilink anywhere.
Without your data or code its almost impossible to say
Thanks Reeza,
Please let me know if this helps. Note this is only one of the codes I tried. Thanks!
proc glimmix data=sows method=quad plots=residualpanel(conditional marginal) ;
CLASS week WeighGr90 ;
MODEL CV = WeighGr90 / dist=poisson ;
*random _residual_ / group=vargroup ;
random intercept / subject=week ;
lsmeans WeighGr90/pdiff;
output out=igausout pred=p resid=r STUDENT=std;
run;
The default link function for the Poissonis log. So, I think all your graphics are on the log linear predictor scale. You can add the ilink option to the plots statement to get the inverse line (original scale).
Thanks lvm, but the problem are the lsmeans estimates, even with the ilink option on the lsmeans statement it is retuning the negative values. Any thoughts?
You have to be more specific. YOu have to put ilink in the lsmeans statement and in the plots opiton (on the proc statement) to see the inverse link means. Even with ilink option on lsmeans, you first get the logs, then on the right the inverse link (data scale means). Your example code does not give ilink anywhere.
lvm,
Using the ilink on the lsmeans statement was actually giving the right information, but you are right I was only looking on the left side of the table (so used to it). It is giving me on the right side like you said, the information that I needed. Thanks a lot! I really appreciate it! Problem solved.
Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.
Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.
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.