plot函数

plot(x,y,xlab=,ylab=)
plot的第一个参数显示在水平轴上，第二个参数显示在纵轴上
习惯上垂直轴绘制反应变量，水平轴绘制解释变量

plot(x=Veg$ BARESOIL,y=Veg$ R,xlab=“Exposed soil”,ylab=“Species richness”,main=“Scatter plot”,xlim=c(0,45),ylim=(4,19))

xlab和ylab用于标签
选项main应用于标题
xlim和ylim指定坐标上下限
xlim=c(min(Veg$ BARESOIL),max(Veg$ BARESOIL))
如果有缺失值，则：
xlim=c(min(Veg$ BARESOIL,na.rm=TRUE),max(Veg$ BARESOIL,na.rm=TRUE))

缺省情况下，plot函数使用空心圆，plot函数里的pch选项可以指定绘图字符

col表示颜色，以下代码可以看到8种颜色所对应的数字
x<-1:8
plot(x,col=x)

Veg$ Time2 <-Veg$ Time #生成和BARESOIL与R相同长度的新向量
Veg$ Time2[Veg$ Time <=1974] <- 15 #1958-1974用实心方块表示
Veg$ Time2[Veg$ Time >1974] <- 16 #1981-2002用实心圆点表示
Veg$ Col2 <- Veg$ Time #生成和BARESOIL与R相同长度的新向量
Veg$ Col2[Veg$ Time <=1974] <- 1 #黑色
Veg$ Col2[Veg$ Time >1974] <- 2 #红色
plot(x=Veg$ BARESOIL,y=Veg$ R,xlab=“Exposed soil”,
ylab=“Species richess”,main=“Scatter plot”,
xlim=c(0,45),ylim=c(4,19),pch=Veg$ Time2, col=Veg$Col2)

cex选项参数用于改变符号尺寸大小

Veg$ Cex2 <- Veg$ Time
Veg$ Cex2[Veg$ Time == 2002] <- 2 #大的实心点
Veg$ Cex2[Veg$ Time != 2002] <- 1 #小的实心点
plot(x=Veg$ BARESOIL,y=Veg$ R,xlab=“Exposed soil”,
ylab=“Species richess”,main=“Scatter plot”,
xlim=c(0,45),ylim=c(4,19),cex=Veg$Cex2,pch=16,col=3)

改变符号尺寸也可以用cex=1.5*Veg$ Cex2或cex= Veg$ Cex2/2

添加平滑线

由于lines命令的第一顺序是按照顺序连接的，所以上图的线条不好看。

order函数用来确定BARESOIL中元素的顺序，并允许在lines命令里把值由小到大重新排列
lwd和lty分别表示线的宽度和线的类型

总结