Why DESeq2 filter and remove genes with high count?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Differential Gene Expression Analysis with DESeq2: Is it correct to just get a subset of genes that…
written by Savindu •My understanding is that if we were to subset a set of genes and carry out DESeq2 just for those genes, it would generate inaccurate …
-
what can cause high(> 80%) spike-ins and low read count (<100 reads) in low complexity ampseq run
written by eli_bayat •Hello, I have data output from low complexity ampseq run. The read counts are very low (<100 reads per well) and the spike-ins% is very …
-
Removal of 'low' counts before DESeq2 - what threshold?
written by dnljmrs •Hi all, I have run DESeq2 successfully in my RNAseq experiment and have observed some differentially expressed genes, so thats fine. There are 2300 genes …
-
Expression thresholds for enrichment analysis of RNA-Seq data
written by James Cook •Hi Guys I have a question about filtering genes detected in an RNA-Seq experiment for subsequent enrichment analysis. I know it's common to use a …
-
Deduplication using UMItools
written by Ati •I have some RNAseq data with a high duplication rate but the reads have UMI (Unique Molecular Identifiers). The UMI length is 5 bp. I …
-
Survival plot between low and high expression of gene
written by BiologistHi, I wanted to make a survival plot showing between low and high expression samples of a gene. I followed this [cutpoint using maxstat package][1] …
-
From read counts to genotypes, or even better, a vcf
written by Joanna.ilska •Hi all, I am fairly new to dealing with sequence data and have come across a problem which I am not entirely sure how to …
-
number of grouping affect log2 fold change in DESeq2 analysis?
written by starI have RNA-seq data for samples in 10 different groups and I would like to find deferentially expressed genes. I have used Deseq2 R package. …
-
DESeq2 with 2 factors: problem to fit the negative binomial and weird p-value distributions
written by Aurelie MLBHello, I have RNA-seq data from 27 samples to analyse. I have 2 differents treatments (Compound1 and compound 2) + control (DMSO) and I have …
-
Why Is Overdispersion Worse At Higher Coverage?
written by Stan Letovsky •<p>In <a href="http://biostar.stackexchange.com/questions/17360/is-rnaseq-overdispersion-due-to-sample-prep-bottlenecks">this post</a> ZK noted that Poisson is a better approximation at low coverage than at high, something I have also observed. Does anyone …
Outliers are removed, on both sides. Are those genes highly expressed in all samples or just a few?
If removed genes have high counts in all samples what is your explanation? @WouterDeCoster
How high is "high coverage"? DESeq2 performs a few different kinds of filtering. At the gene level, though, it only filters for maximal power. Regardless of how high the counts are, presumably that's why the gene got filtered. You can disable that (look for "independent filtering" in the documentation) if it's a problem.