4/19/2017

Loading the data

# 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")

Visualizing one part of the data

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")

Are the differences significant?

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