aov(BD~Site,data=BDdata) #ANOVA for bulk density data summary(aov(BD~Site,data=BDdata)) TukeyHSD(aov(BD~Site,data=BDdata)) #postHOC for BD data aov(WC~Site,data=WCdata) #ANOVA for water content data TukeyHSD(aov(WC~Site,data=WCdata)) #postHOC for WC data kruskal.test(TTR~Site,data=TTRdata) #Kruskal wallis rank sum for time to runoff and bd data cor(CORdata$TTR,CORdata$BD) #correlation between time to runoff and bd chisq.test(RVdata$Site,RVdata$RV) #chi square test for runoff volume data