veg<-read.csv2("veg.csv", row.names = 1)
veg<-read.csv("veg.csv", row.names = 1,
              sep=";", dec=",")


env<-read.csv2("env.csv", row.names = 1)

?read.csv

table(rownames(veg)==rownames(env))

summary(veg)
summary(env)

class(veg)
dim(veg)

# install.packages('vegan')
library(vegan)
SpRich<-specnumber(veg)
quantile(SpRich)

hist(SpRich)

boxplot(SpRich~env$hemiparasite)

pdf("scatter.pdf", 5,5)
plot(SpRich~env$productivity)
dev.off()

svg("scatter.svg", 5,5)
plot(SpRich~env$productivity)
dev.off()

png("scatter.png", 5*300,5*300, res=300)
plot(SpRich~env$productivity)
dev.off()

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

veg$Bromus.erectus[env$productivity>500]