topic Re: Zero-inflated model using proc GLIMMIX in SAS Procedures

Zero-inflated model using proc GLIMMIX

jackice — Wed, 02 Nov 2016 18:24:23 GMT

Hello all,

I have a data set of count data (Variable is TOTAL) that on histogram is definitly following Poisson and have many zeros values as well

When I run this programmation below it runs fine in its default normal mode but as soon as I put Dist Poisson Link=log" in the programmation then SAS cant do it ... I am guessing its because of Log(0)=err but the same model in R with family poisson runs fine ... why and what is my solution in this SAS programmation please?

proc glimmix data=FAMAIN plots=ALL;
class HiveOrientation SunPosition Box Bloc Boxid HiveID Day;
model TOTAL = HiveOrientation|SunPosition|Box /*/dist=poisson link=log solution*/;
random Bloc HiveID(Bloc) Day(HiveID Bloc);
lsmeans SunPosition /pdiff ilink lines;
run;

Re: Zero-inflated model using proc GLIMMIX

PeterClemmensen — Wed, 02 Nov 2016 18:27:36 GMT

Can you provide some sample data?

Re: Zero-inflated model using proc GLIMMIX

jackice — Wed, 02 Nov 2016 19:16:43 GMT

there u go draycut and thanks for looking into this for me

Re: Zero-inflated model using proc GLIMMIX

jackice — Wed, 02 Nov 2016 19:29:43 GMT

perhaps i have to change SAS Main procedure and include zero inflated ? or I dont have enough zero and I could manually transformed log(x+1) - count data variance increasing with count value is my case for sure ...