How to not include start region in analyzing ribosome profiling data?
0
0
Entering edit mode
@nicholsonangela7-13852
Last seen 7.2 years ago
UCSD

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

 

 

 

riboprofiling ribosome profiling • 841 views
ADD COMMENT

Login before adding your answer.

Traffic: 780 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