Comments (3)
Hi Freek,
I was investigating the effects of rRNA on final sequencing results. One of the first questions I had: Does >(in iRAP) rRNA abundace influence the final quantification? Meaning: are there still rRNA reads present > when the final (TPM) normalization is performed?
The short answer to your two questions is: it will depend on the sequencing protocol (e.g., are the rRNA depleted?) and/or the analysis protocol (e.g., the quantification can be made for all biotypes - including rRNA - or only for protein coding, in the first case, depending on the protocol, the rRNA may affect the TPMs while in the second it does not).
Note that you can remove the rRNA genes from the quantification matrix and then recompute the TPMs manually (using the irap_raw2metric script).
I do not find anything about rRNA in the original publication of iRAP but I do find it referenced in this iRAP script: https://github.com/nunofonseca/irap/blob/master/aux/R/irap_utils.R (line 2015, 2024 and 2037). Is there a way in iRAP to quantify rRNA abundance?
Each BAM file may get a companion .bam.stats file with the number of alignments per biotype (including rRNA).
Cheers.
from irap.
Hi! As I mentioned before, in preparation of the coming release (0.9.0) I'm going through all open issues and trying to address them. Please feel free to reopen the issue if necessary. Cheers.
from irap.
Ok, very nice thank you.
By the way, I just recently indeed filtered out all biotypes except protein coding and re-normalized, this indeed produces more stable results (regardless of silva/bbduk based rRNA removal which removes half the unmaped reads typically).
Keep up the great work!
from irap.
Related Issues (20)
- aux/R/irap_de.R::process.cmdline.args mis-casts --min (opt$min_count) as a string, should be integer HOT 1
- sample name do not match
- fastq_pre_barcodes.c:34:17: fatal error: bam.h: No such file or directory
- Downloads to a submodule?
- iRAP pipeline for EMBL-EBI Expression Atlas HOT 1
- Configuration file for strand specific protocol
- mapper options in config file is ignored HOT 3
- irap_raw2metric no html file, not gzipped HOT 7
- our_prefix error when analysis name is same as the raw read file name without suffix HOT 2
- failed to load GTF file HOT 3
- RSEM strandness HOT 12
- TPM values rounded to 2 decimals HOT 3
- Docker link in wiki is wrong HOT 1
- Error running Docker image HOT 2
- iRAP's Kallisto call seems to provide read length at fragment size parameter HOT 1
- sh: 1: set: Illegal option -o pipefail error in "irap_gtf_to_fasta". HOT 4
- iRAP does not recognise stage0 completion for large genomes with HISAT2 HOT 5
- docker image for v1.0.1 HOT 2
- docker image 1.0.1 not working HOT 1
- No rule to make target HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from irap.