小编use*_*480的帖子

我有三张不同基因组间隔的表.这是一个例子:

> a
   chr interval.start interval.end names
1 chr1              5           10     a
2 chr1              6           10     b
3 chr2              7           10     c
4 chr3              8           10     d

> b
   chr interval.start interval.end names
1 chr1              6           15     e
2 chr1              7           15     f
3 chr1              8           15     g

> c
   chr interval.start interval.end names
1 chr1              7           12     h
2 chr1              8           12     i
3 chr5              9           12     j
4 chr10             10          12     k
5 chr20             11          12     l

我试图在将信息转换为GRanges后找到所有表之间的公共间隔.基本上我想做像交叉(c,intersect(a,b))之类的东西.但是,因为我正在使用基因组坐标,所以我必须使用GRanges和GenomicRanges包,我不熟悉. …