Help using Base SAS procedures

Zero-inflated model using proc GLIMMIX

Reply
Contributor
Posts: 20

Zero-inflated model using proc GLIMMIX

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;

PROC Star
Posts: 554

Re: Zero-inflated model using proc GLIMMIX

Can you provide some sample data?

Contributor
Posts: 20

Re: Zero-inflated model using proc GLIMMIX

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

Contributor
Posts: 20

Re: Zero-inflated model using proc GLIMMIX

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 ... 

Ask a Question
Discussion stats
  • 3 replies
  • 225 views
  • 0 likes
  • 2 in conversation