Fitting Distribution Of The Gene Expression Levels Of Rna Seq Data

I am trying to fit a distribution of the gene expression (DESeq Normalized counts) of each sample from a RNA Seq data. And I came across this paper where they describe that the expression levels of genes have a Bimodal distribution. .

I was wondering would this be true for the count data as well or do we expect any specific distribution for Counts that will have a change of skew that may arise due to the length factor that comes in RPKM?

Can some one help me here on how to test this for count data?

Thanks

• 5,817 views • link

written 12.6 years ago by k.nirmalraman 114

This post in ResearchGate might be helpful

• 198 views • link

written 12.6 years ago by Sudeep 168

1 answer

Many groups model the raw counts in RNA-seq data as a negative binomial.

• 42 views • link

written 12.6 years ago by Sean Davis 2.7K

I thought this applies for the distribution of say Gene-i across samples. I understood it that way from the DE Seq Paper

Equation 1:

http://www.w3.org/1998/Math/MathML " name="gb-2010-11-10-r106-i1"> K_{i j} ~ NB (μ_{i j}, σ_{i j}^{2}),

• 222 views • link

written 12.6 years ago by k.nirmalraman 114

More posts like this

Calculating differential expression from RPKM
written 10.6 years ago by Mukund 0

I have a dataset (downloaded from GEO) that has expression levels in terms of RPKM. I understand that edgeR, DeSeq etc, use raw counts and …
RPKM values in normal cells
written 11.1 years ago by akantusz 0

Hello, I would like to compare gene expression level data from COSMIC to expression level data of normal cells. I've found the TCGA level 3 …
Remove Batch Effect From RNAseq with SVAseq and Combat
written 11.5 years ago by cankutcubuk 19

Hi, I need your help and comments to remove batch effect from RNA-Seq data. In total I have almost 1000 RNA-Seq data of different cancer …
Differential Gene Expression Between Control And Treated Samples, Rpkm
written 12.5 years ago by HNK 15

Hi Everyone, I am working on a set of samples (generated using Illumina) for differential gene expression between control and treated samples. I performed pre …
Representing Genome Wide Regulatory Elements
written 13.2 years ago by ancient_learner 68

Hi I want to represent the presence of certain important regulatory elements for each chromosome separately. I came across a particular paper (figure 2) and …
Scatter Plot For Correlations With Heatdensity
written 13.3 years ago by k.nirmalraman 114

Hi All, I am trying to show some correlation between two samples and would like to do a scatter plot for the same. I …
Deseq Infinite In Logfc And Na For P Value
written 13.5 years ago by k.nirmalraman 114

I am performing DE analysis using DESeq package. The single end illumina reads are bowtie aligned. I used the raw counts for DE as required …
Nomenclature: Genes With Reads In An Rna-Seq Experiment?
written 13.6 years ago by 14134125465346445 361

I would like to know what is the best way to call the set of genes that show "at least one read" in an RNA-seq …
Draw Heatmap Or Do Pca Analysis With Raw Read Counts?
written 13.7 years ago by camelbbs 71

Hi, I want to ask a question about viewing RNAseq data with raw read counts. After I get the raw reads counts from HTseq-count or …
What Makes One Probability Distribution Better For Rna-Seq Than Another?
written 15.6 years ago by Jeremy Leipzig 2.3K

Many experts claim that the negative binomial distribution is better than the Poisson distribution for modeling discrete RNA-Seq data. Is this because <ul> <li>The …