data1 = read.csv(file.choose(), header = TRUE)
###Use files named: Goby communities for ANOVA.csv
### AND Fish communities for ANOVA.csv

data1

############## Abundance ####################
boxplot(data1$Abundance ~ data1$Habitat)

abundance.aov = aov(data1$Abundance ~ data1$Habitat)
summary(abundance.aov)



########################## Richness #####################
boxplot(data1$Richness ~ data1$Habitat)

richness.aov = aov(data1$Richness ~ data1$Habitat)
summary(richness.aov)



################ Diversity #####################
boxplot(data1$Diversity ~ data1$Habitat)

diversity.aov = aov(data1$Diversity ~ data1$Habitat)
summary(diversity.aov)



################ Rugosity ###################
boxplot(data1$Rugosity ~ data1$Habitat)

rugosity.aov = aov(data1$Rugosity ~ data1$Habitat)
summary(rugosity.aov)