I have percentage areas values from 4 treatments (25 replicates per treatment). I would like to compare these percentages.

I am supposed to use a beta regression, because my response variable is a proportion (not resulting from a count). For other variables I use GLMM and then use the glht() function to get the significance of pairwise comparisons. This does not seem to work with my betareg model. Is there any way to do that ?

Also I would like to include random effects because of nested sampling protocol. From this answer this does not seem very appropriate/easy to do this with beta regression.

Do you have any suggestion ?


Thanks to @rvl's comment I could calculate pairwise comparisons and extract group letters.

   betalive = betareg(Live ~ Crop)
   live.rg <- ref.grid(betalive)        
   live.lsm = lsmeans(live.rg, "Crop")
   cld(live.lsm, alpha =  0.05)$.group
   > [1] " 1 " " 1 " " 1 " "  2"

However I could not solve my problem of random effects yet.

  • $\begingroup$ The lsmeans package supports betareg models. Look at ?models after loading the package for more info. $\endgroup$
    – Russ Lenth
    Apr 30, 2016 at 17:19
  • $\begingroup$ Thanks @rvl, I could compute pairwise comparisons with lsmeans $\endgroup$
    – Nausi
    May 2, 2016 at 7:56


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.