【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: 还有如下三个问题有待于后面总结上来:

  1. 为什么颜色我标注的是红色,但是显示出来的是黑色
  2. 如何像上面的那个图一样每个条形之间有空隙
  3. 如果数据集的x轴过大,我可不可以把x轴中间切割空隙,空隙之间搞一个跨度,这样就能把所有的x轴都包含过来
个人公众号,比较懒,很少更新,可以在上面提问题:

更多精彩,请移步公众号阅读:

Sam avatar
About Sam
专注生物信息 专注转化医学