2015年9月25日

#plot，其中对颜色按数据分析

Veg <- read.table(file = "Vegetation2.txt",header = T)

names(Veg)

?par

Veg$Time2 <- Veg$Time

Veg$Time2 [Veg$Time <= 1974] <-15

Veg$Time2 [Veg$Time >1974] <-16

Veg$Col2 <- Veg$Time

Veg$Col2 [Veg$Time <= 1974] <- 1

Veg$Col2 [Veg$Time >1974] <- 2

plot(x=Veg$BARESOIL,y=Veg$R,xlab = "Exposed soil",

     ylab = "Species richness",xlim=c(0,45),ylim=c(4,19),

     main="Scatter plot",pch=Veg$Time2,col=Veg$Col2,cex=1.0)