r - 将 BED 文件转换为 WIG 文件

Question

你知道有什么方法可以轻松地将 BED 文件转换为 WIG（通过 R 或其他程序）吗？

你能给我一些指导吗？

Answer 1

查看rtracklayer包，特别是以下手册页：

 ?import
 ?export

Answer 2

这是我将如何转换.bed为.wig. 正如@Paolo 暗示的那样，这是一个简单的过程：

library(rtracklayer) #bioconductor

bed_loaded <- import(con="~/Downloads/my_bed.bed.gz", format="bed") #no need to unzip .gz
# bed_loaded <- import.bed(con="~/Downloads/my_bed.bed") #if you unzip 

export.wig(object=bed_loaded, con="~/Downloads/bed2wig.wig")

请注意，您import和export两者都有方法（假发、床、大人物或 bw 等）。您可以直接使用它们而不指定方法但指定format参数。

这个GitHub 教程会很有帮助。