head(Cone11.12.3)
names(Cone11.12.3)<-c("Habitat", "Presence")
aov11.12.3<-aov(Presence~Habitat, data=Cone11.12.3)
summary(aov11.12.3)
TukeyHSD(aov11.12.3)