aroundPromoter function in RiboSeqR gives error
0
0
Entering edit mode
@amirmahboubi79-18020
Last seen 5.4 years ago

Hi,

I am performing RiboseqR analysis now and ran into a problem with aroundPromoter command that I don't understand and have no idea what to do. The problem is that when I run my command as below:

#Install

source("http://bioconductor.org/biocLite.R")

#Load the library

library(RiboProfiling)

#BAM file import

library(GenomicAlignments)

aln <- readGAlignments(BamFile("/Users/amir/Documents/Amir files/iSeq 100 data/Iseq data analyzed/‌bam files/BAM_4D-Suc.bam"))

#txdb object with annotations

biocLite("TxDb.Athaliana.BioMart.plantsmart28")

library(TxDb.Athaliana.BioMart.plantsmart28)

txdb <- TxDb.Athaliana.BioMart.plantsmart28

oneBinRanges <- aroundPromoter(txdb, alnGRanges, percBestExpressed=0.001)

 

I then get the error:

 

Error in if (quantCounts <= 0) { : missing value where TRUE/FALSE needed

 

Can you please help me what should I do to get rid of this problem?

 

Thanks

Amir Mahboubi

Postdoc at Umea plant science centre

Umea, Sweden

riboseqr • 906 views
ADD COMMENT
0
Entering edit mode

I have exact same problem analyzing footprints from mouse liver. Any insight would be much appreciated.

ADD REPLY
0
Entering edit mode

Hi everyone! I am also running into the same error message with the aroundPromoter function. Any progress? Has anyone been successful in contacting the researchers maintaining this package? Thank you!

ADD REPLY

Login before adding your answer.

Traffic: 669 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6