【1.9】直方图(hist)
上面的这个图是我要模仿的对象,捣弄了一段时间,大致弄了一一副这个图
>mm<-subset(d,length<10000)
#因为数据集d的x轴太大了,所以我就只好取前面的一部分了
> hist(mm$length,breaks=1000,xlim=c(0,9000),main="Frequency of Length",xlab="Length(bp)",ylab="Scaffold Number",xaxt="n")
#xaxt=”n”为关闭x轴的坐标轴
>axis= axis(1,at=seq(0,10000,500))
#添加上坐标轴:1代表横坐标轴,2代表中纵坐标,0,10000代表坐标起始和终止位点,500代表间隔
labels = TRUE,表示在每列头上表注数目
freq=FALSE 纵坐标可以显示概率,而不是数目
部分参数: >help(“hist”)
hist(x, ...)
## Default S3 method:
hist(x, breaks = "Sturges",
freq = NULL, probability = !freq,
include.lowest = TRUE, right = TRUE,
density = NULL, angle = 45, col = NULL, border = NULL,
main = paste("Histogram of" , xname),
xlim = range(breaks), ylim = NULL,
xlab = xname, ylab,
axes = TRUE, plot = TRUE, labels = FALSE,
nclass = NULL, warn.unused = TRUE, ...)
以某一个数为x坐标,以对应的另一个数的数值的总和为纵坐标。
> library(plotrix)
> weighted.hist(mm$gc,mm$leth,breaks=0:50)
> weighted.hist(mm$gc,mm$length,breaks=0:100) #这个breaks是亮点啊,灵活处理的
ps: 还有如下三个问题有待于后面总结上来:
- 为什么颜色我标注的是红色,但是显示出来的是黑色
- 如何像上面的那个图一样每个条形之间有空隙
- 如果数据集的x轴过大,我可不可以把x轴中间切割空隙,空隙之间搞一个跨度,这样就能把所有的x轴都包含过来
这里是一个广告位,,感兴趣的都可以发邮件聊聊:tiehan@sina.cn
个人公众号,比较懒,很少更新,可以在上面提问题,如果回复不及时,可发邮件给我: tiehan@sina.cn
个人公众号,比较懒,很少更新,可以在上面提问题,如果回复不及时,可发邮件给我: tiehan@sina.cn