# Load in the species matrix and the environmental metadata together otu_matrix_w_enviro <- read.csv("data/2011-02-06_OTU_counts_and_metadata_QC_NBZ.csv")
4/19/2017
# Load in the species matrix and the environmental metadata together otu_matrix_w_enviro <- read.csv("data/2011-02-06_OTU_counts_and_metadata_QC_NBZ.csv")
boxplot(OTU_0001 ~ approx_annual_rainfall_mm, data = otu_matrix_w_enviro, xlab = "Mean annual rainfall (mm)", ylab = "Number of sequences (abundance)", main = "OTU0001 abundance and mean annual rainfall")
kruskal.test(OTU_0001 ~ approx_annual_rainfall_mm, data = otu_matrix_w_enviro)
## ## Kruskal-Wallis rank sum test ## ## data: OTU_0001 by approx_annual_rainfall_mm ## Kruskal-Wallis chi-squared = 80.677, df = 6, p-value = 2.589e-15