パッケージで LiftOver 関数を使用しようとしていrtracklayer
ますが、1 つのエラーが発生します。これが私のコードです。
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")
なぜそれが機能しないのかいくつかのアイデアはありますか? ありがとうございました!