A long-term agricultural experiment had n=90 grassland plots, each 25 m×25 m, differing in biomass, soil pH, and species richness (the count of species in the whole plot). While it was well-known that species richness declines with increasing biomass, it was not known how this relationship depends on soil pH. In the experiment, there were 30 plots of "low pH", 30 of "medium pH" and 30 of "high pH". Three lines of the data are reproduced here as an aid.
Briefly explain the commands below. That is, explain the models being fitted.
Let H1,H2 and H3 denote the hypotheses represented by the three models and fits. Based on the output of the code below, what hypotheses are being tested, and which of the models seems to give the best fit to the data? Why?
![](https://cdn.mathpix.com/cropped/2022_04_28_f42726bf90cfbe3fb2f8g-14.jpg?height=364&width=561&top_left_y=877&top_left_x=220)
Finally, what is the value obtained by the following command?
![](https://cdn.mathpix.com/cropped/2022_04_28_f42726bf90cfbe3fb2f8g-14.jpg?height=128&width=824&top_left_y=1286&top_left_x=222)
>grass[c(1,31,61),] pH Biomass Species 1 high 0.46929723031 mid 0.17576272961 low 0.100847918> fit1 <- glm(Species Biomass, family = poisson) > fit2 <- glm(Species pH+ Biomass, family = poisson) > fit3 <- glm(Species pH∗ Biomass, family = poisson)