Question: How to not include start region in analyzing ribosome profiling data?
gravatar for nicholsonangela7
13 months ago by
nicholsonangela70 wrote:

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 13 months ago by nicholsonangela70
Please log in to add an answer.


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