-1

我正在尝试在rtracklayer包中使用 LiftOver 函数,但出现一个错误,这是我的代码:

library(rtracklayer)
library(gwascat)
library(IRanges)

chain <- import.chain("bosTau6.hg19.all.chain")

df <- read.table(textConnection("chr  start  end
chr5  58024309  58039917"), header =T)

library(GenomicRanges)

gr <- makeGRangesFromDataFrame(df, TRUE)

humcon <- liftOver(gr, chain)

要将链文件下载并解压缩到我使用的工作文件夹中:

wget http://hgdownload.cse.ucsc.edu/goldenPath/bosTau6/vsHg19/bosTau6.hg19.all.chain.gz
gunzip gunzip bosTau6.hg19.all.chain.gz

错误如下:

Warning message:
'mseqapply' is deprecated.
Use 'as(mapply(FUN, ...), 'List')' instead.
See help("Deprecated")

一些想法为什么它不起作用?谢谢!

4

1 回答 1

1

这只是一个警告,而不是错误。

一个简单的网络搜索发现:

https://www.biostars.org/p/122005/

于 2014-12-03T13:13:07.680 回答