当前位置:  开发笔记 > 编程语言 > 正文

扩展两个方向的范围

如何解决《扩展两个方向的范围》经验,为你挑选了1个好方法。



1> Ven Yao..:

这很简单.你可以使用startend功能GenomicRanges.

gr <- GRanges(c('chr1','chr1'), IRanges(start=c(20, 120), +')
gr
# GRanges object with 2 ranges and 0 metadata columns:
#       seqnames     ranges strand
#              
#   [1]     chr1 [ 20,  29]      +
#   [2]     chr1 [120, 129]      +
#   -------
#   seqinfo: 1 sequence from an unspecified genome; no seqlengths

start(gr) <- start(gr) - 10
end(gr) <- end(gr) + 10
gr
# GRanges object with 2 ranges and 0 metadata columns:
#      seqnames     ranges strand
#             
#  [1]     chr1 [ 10,  39]      +
#  [2]     chr1 [110, 139]      +
#  -------
#  seqinfo: 1 sequence from an unspecified genome; no seqlengths


当然.但是由于有很多函数可以改变范围(调整大小,缩小等),或者根据现有的(侧翼,启动器)定义新的范围,我期望一个专用函数能够做到这一点相当简单,但(我认为)是非常有用的操作.
推荐阅读
jerry613
这个屌丝很懒,什么也没留下!
DevBox开发工具箱 | 专业的在线开发工具网站    京公网安备 11010802040832号  |  京ICP备19059560号-6
Copyright © 1998 - 2020 DevBox.CN. All Rights Reserved devBox.cn 开发工具箱 版权所有