0

我有两个床文件作为 R 中的数据框,我想将所有重叠区域相互映射(类似于最接近的床工具能够做的事情)。

床A:

chr   start   end
 2       100     500
 2       200     250
 3       275     300

床B:

chr    start    end
  2       210      265
  2       99       106
  8       275      290

卧床:

chr   start.A   end.A  start.B  end.B
 2       100     500      210      265
 2       100     500      99       106
 2       200     250      210      265

现在,我发现了这个非常相似的问题,建议使用 iRanges。使用建议的方式似乎可行,但我不知道如何将输出转换为像“BedOut”这样的数据框。

4

2 回答 2

0

这是使用该data.table软件包的解决方案。

library(data.table)

chr = c(2,2,3)
start.A = c(100, 200, 275)
end.A = c(500, 250, 300)
df_A = data.table(chr, start.A, end.A)

chr = c(2,2,8)
start.B = c(210, 99, 275)
end.B = c(265, 106, 290)
df_B = data.table(chr, start.B, end.B)

首先,内连接键上的数据表chr

df_out = df_B[df_A, on="chr", nomatch=0]

然后过滤重叠区间:

df_out = df_out[(start.A>=start.B & start.A<=end.B) | (start.B>=start.A & start.B<=end.A)]
setcolorder(df_out, c("chr", "start.A", "end.A", "start.B", "end.B"))

   chr start.A end.A start.B end.B
1:   2     100   500     210   265
2:   2     100   500      99   106
3:   2     200   250     210   265
于 2020-05-11T11:50:21.487 回答
0

另一个data.table选项使用foverlaps

setkeyv(BedA, names(BedA))
setkeyv(BedB, names(BedB))
ans <- foverlaps(BedB, BedA, nomatch=0L)
setnames(ans, c("start","end","i.start","i.end"), c("start.A","end.A","start.B","end.B"))

输出:

   chr start.A end.A start.B end.B
1:   2     100   500      99   106
2:   2     100   500     210   265
3:   2     200   250     210   265

数据:

library(data.table)
BedA <- fread("chr   start   end
2       100     500
2       200     250
3       275     300")

BedB <- fread("chr    start    end
2       210      265
2       99       106
8       275      290")
于 2020-05-12T00:58:41.920 回答