As part of my 3rd year university dissertation project i am studying whether sheep grazing has a effect on bird species diversity. This involved comparing the number of species found at 3 different habitat types and resulted in the following data (Habitats= habitat type, Totalsp= Total number of species found at the habitat):

```
Habitats Totalsp
```

Grazed grassland 9

Non-grazed grassland 18

Woodland 15

I am trying to carry out a regression on the data to determine if there was a significant difference in the number of species found at each habitat type, with no luck so far. This is the code i have used to try to do this:

fit <- lm(Totalsp ~ Habitats, data=Diss_origdataR)

summary(fit)

and this was the output:

Call:

lm(formula = Totalsp ~ Habitats, data = Diss_origdataR)

Residuals:

ALL 3 residuals are 0: no residual degrees of freedom!

Coefficients:

Estimate Std. Error t value Pr(>|t|)

(Intercept) 9 NaN NaN NaN

HabitatsNon-grazed grassland 9 NaN NaN NaN

HabitatsWoodland 6 NaN NaN NaN

Residual standard error: NaN on 0 degrees of freedom

Multiple R-squared: 1, Adjusted R-squared: NaN

F-statistic: NaN on 2 and 0 DF, p-value: NA

Does anyone have any advice on how i might fix this?