Search
Question: How to not include start region in analyzing ribosome profiling data?
0
gravatar for nicholsonangela7
11 weeks ago by
UCSD
nicholsonangela70 wrote:

Hi,
I am new to ribosome profiling but am hoping I can use the RiboProfiling package for analyzing my data. I want to customize where I am looking for protected fragments though because I would like to not include the start region (the first 50 nucleotides) for all transcripts as it can sometimes result in a biased comparison between genes. Is there a way I can use this package and specify that I would like to not count reads that fall on the first 50 nucleotides of every transcript? I tried using the listshiftvalue parameter but the output I got in the data.frame is the same as when I didn't add 50 there so I think this must be the wrong place. 
covData <- riboSeqFromBAM(listInputBam, txdb = txdb, listShiftValue=c(50))

Thank you! 

Angela Nicholson

 

 

 

ADD COMMENTlink written 11 weeks ago by nicholsonangela70
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 104 users visited in the last hour