Кто-нибудь когда-нибудь пытался создать нулевую модель последовательности ДНК у мыши, используя пакет gkmSVM ?Он отлично работает для человека, но для мыши - нет.Мне интересно, если кто-нибудь из вас никогда не использовал этот пакет для мышей и имел такую же проблему.
При использовании функции genNullseqs возникает ошибка:
library(gkmSVM)
library(BSgenome.Mmusculus.UCSC.mm10)
library(BSgenome.Mmusculus.UCSC.mm10.masked)
library(IRanges)
genome=BSgenome.Mmusculus.UCSC.mm10.masked
fileBedBreaks="Rep1.intersec.Rep2.cov.major2.sort.uniq.bed"
fileFastaPos="Rep1.intersec.Rep2.cov.major2.all.sort.uniq.bed.pos.fa"
fileBedNeg="Rep1.intersec.Rep2.cov.major2.all.sort.uniq.Random.gkmSVM.bed"
fileFastaNeg="Rep1.intersec.Rep2.cov.major2.all.sort.uniq.Random.gkmSVM.fa"
genNullSeqs(inputBedFN=fileBedBreaks,nMaxTrials=5,xfold=2,genome=genome,
outputPosFastaFN=fileFastaPos,outputBedFN=fileBedNeg,
outputNegFastaFN=fileFastaNeg)
Ошибка
Error in normalizeDoubleBracketSubscript(i, x, exact = exact) :
subscript "TRF" matches no name