bk.env<-read.delim2("env.csv", sep=";", row.names = 1)
bk.spe<-read.delim2("veg.csv", sep=";", row.names = 1)

summary(bk.env)
summary(bk.spe)

table(rownames(bk.env)==rownames(bk.spe))

class(bk.env)
# [1] "data.frame"

dim(bk.env)
# [1] 26 14
dim(bk.spe)
# [1]  26 231

install.packages("vegan")
library(vegan)

?specnumber
SpRich<-specnumber(bk.spe)

quantile(SpRich)
# 0%  25%  50%  75% 100% 
# 33.0 41.0 45.0 55.5 69.0

hist(SpRich)

boxplot(SpRich~bk.env$hemiparasite)

plot(SpRich~bk.env$productivity)

install.packages("psych")
library(psych)
pdf("PairPanels.pdf", 10, 10)
pairs.panels(bk.env[,c(1, 6:14)])
dev.off()

png("PairPanels.png", 1000, 1000)
pairs.panels(bk.env[,c(1, 6:14)])
dev.off()

svg("PairPanels.svg", 10, 10)
pairs.panels(bk.env[,c(1, 6:14)])
dev.off()

bk.rmi.spe<-bk.spe[bk.env$hemiparasite=="RMI",]
bk.rmi.env<-bk.env[bk.env$hemiparasite=="RMI",]

bk.spe$Bromus.erectus[bk.env$productivity>500]
# [1]  0 10  0 10  3 20