Hi  I have a data set with outcome as proportions.   When i saw the distribution of the outcome, it looked like Poisson distribution and hence did a GLMM with Poisson distribution.   However, the reviewer for the manuscript did not accept this and has asked me to do GLMM with beta regression.   The problem is, the model wont run. PROC mixed and PROC GLMM (Poisson) runs ok. Not the beta regression.   I am adding the code i am using here. Any ideas will be welcome.   proc glimmix data=Cyril_nov29; class Treatment Day Rep; model freqsowcont= Treatment Day Treatment*Day/ddfm=kr dist=poisson link=log; Random intercept/subject=Day(Treatment*Rep); lsmeans Treatment Day Treatment*Day/pdiff lines ilink; run;   This one runs ok.   proc glimmix data=Cyril_nov29 order=formatted exphessian; class Treatment Day Rep; model freqsowcont= Treatment Day Treatment*Day/dist=beta solution link=logit s; Random intercept/subject=Day(Treatment*Rep); lsmeans Treatment Day Treatment*Day/pdiff lines ilink; run;   This doesnt run properly. Why?   I have attached data too as attachment, if anybody want to try it. Please use the first excel sheet (group) ... View more