Question: How to not include start region in analyzing ribosome profiling data?
gravatar for nicholsonangela7
7 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 7 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: 315 users visited in the last hour