I do not see any clear patterns. There is lots of overlap between all the neighbors.I set my notches to false because they were messing looking with nothched. I can put them back on though. I thing presenting all the colonization data together would be best. Thoughts?
rm(list=ls())
colon <- read.csv("F:/Thesis/colon.csv")
boxplot(amf.all ~ neighbor, data = colon, notch=F ,col = "seagreen4",xlab="Neighbor",ylab="Colonization proportion", main= "Total colonization")
boxplot(amf.all ~ neighbor, data = colon,subset = site == "w", col = "seagreen4",xlab="Neighbor",ylab="Colonization proportion", main="Whaley Draw")
boxplot(amf.all ~ neighbor, data = colon,subset = site == "2", col = "seagreen4",xlab="Neighbor",ylab="Colonization proportion", main="Site 2")
boxplot(am.ves ~ neighbor, data = colon, col = "seagreen4",xlab="Neighbor",ylab="Vesicle proportion")
boxplot(am.ves ~ neighbor, data = colon,subset = site == "2", col = "seagreen4",xlab="Neighbor",ylab="Vesicle proportion", main="Site 2")
boxplot(am.ves ~ neighbor, data = colon,subset = site == "w", col = "seagreen4",xlab="Neighbor",ylab="Vesicle proportion", main="Whaley Draw")