-1
我試圖使用LiftOver功能rtracklayer
包,但我得到一個錯誤,這是我的代碼:LiftOver在R(錯誤)
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")
一些想法,爲什麼它不工作?謝謝!