Given this observation, we opt to fit a model for each species using a Gaussian family model. Details about model family [here](https://www.rdocumentation.org/packages/stats/versions/3.6.2/topics/family).
Given this observation, we opt to fit a model for each species using a Gaussian family model (default). Details about model family [here](https://www.rdocumentation.org/packages/stats/versions/3.6.2/topics/family).