频率分布图和散点图SOP，R教程
一、 频率分布图
1.1 数据准备
需要将数据和代码放入同一文件夹：
ATPAN.csv文件
RTPAN.csv文件
1.2 数据格式
1.3 画图
##画左边一列2个图
dev.new()
par(mfcol=c(2,2))
otu=read.csv("ATPAN.csv",row.names=1,header=T)
B<-otu$B
h<-hist(B,freq=FALSE,breaks=10,col="red",ylim=c(0,5),xlim=c(0,1),main="Abundant bloom")
lines(density(B),col="blue",lwd=2)
box()
NB<-otu$NB
h<-hist(NB,freq=FALSE,breaks=10,col="red",xlim=c(0,1),ylim=c(0,5),main="Abundant non-bloom")
lines(density(NB),col="blue",lwd=2)
box()
##画右边一列2个图
otu=read.csv("RTPAN.csv",row.names=1,header=T)
B<-otu$B
B<-na.omit(B)
h<-hist(B,freq=FALSE,breaks=20,col="red",ylim=c(0,5),xlim=c(0,1),main="Rare bloom")
lines(density(B),col="blue",lwd=2)
box()
NB<-otu$NB
h<-hist(NB,freq=FALSE,breaks=20,col="red",xlim=c(0,1),ylim=c(0,5),main="Rare non-bloom")
lines(density(NB),col="blue",lwd=2)
box()
- 文章剩余内容<<<<