https://www.researchgate.net/post/How-to-use-Generalized-Linear-Models-to-compare-means-in-R