Hi,
I have two corresponded count table, they are RNA-seq count data.
Example, actual matrix much larger:
Matrix A:
3 3 6 2 0 0
15 10 0 8 0 0
0 0 0 0 0 0
0 0 0 0 0 0
3 14 2 2 2 7
Matrix B:
29.0 31.5 27 29.5 0.0 0
33.0 37.5 0 34.0 0.0 0
0.0 0.0 0 0.0 0.0 0
0.0 0.0 0 0.0 0.0 0
24.5 22.5 26 23.5 23.5 24
The second matrix is a average of two integer count table, so it has *.5 numbers.
The two matrix are corresponding with each other, that means what in the first matrix is 0, it is 0 at the second.
The problem is, I want normalize the first matrix with the second one, i.e, by element-wise division of first matrix with the second.
- the data are extremely skewed with a lot of '0's.
- the number in second matrix is general much larger than the first. After division, I end up with a matrix containing very tinny numbers and '0's.
I want to find differentially expression genes across conditions. (each 3 columns are replicates). I'm now stacked with how to properly normalize and what kind of model to fit into a statistic test.
Any help would be appreciated.
Best
1 answer
Your matrix A and B - are they environmental replicates? Integrate these two matrices separately into your analysis if so - try filtering instead. Below is a formula I use in edgeR to filter low-expressed genes across RNA-Seq count tables - then use the exact test to compare these groups. First, try using some code I from a very good Nature Protocols paper:
Anders, S. et al. (2013) Count-based differential expression analysis of RNA sequencing data using R and Bioconductor. Nat. Protoc. 8, 1765–86
library(edgeR)
group <- factor(c(rep("your_group_1", 3) , rep("your_group_2", 3)))
d <- DGEList(counts = counts, group = group, genes = genes)
k <- rowSums(cpm(d)>1)>=3
d <- d[k,]
d<-estimateCommonDisp(d,verbose=TRUE)
d<-calcNormFactors(d,method="TMM")
d<-estimateTagwiseDisp(d,verbose=TRUE)
d.est<-estimateSmoothing(d)
et=exactTest(d,pair=c(1:2))
tt=topTags(et,n=nrow(d))
rn=rownames(tt$table)
deg=rn[tt$table$FDR<0.05]
plotMeanVar(d,show.tagwise.vars=TRUE,NBline=TRUE)
plotBCV(d)
plotSmear(d,de.tags=deg)
Cheers!
•
written
by
Craig Henry •
Log in to answer this question.
-
FastQC Quality per tile and per sequence behaving strange after using Cutadapt
written by Sergio A.S. •
-
Error while running nf-core/rnaseq pipeline
written by eesha28112001 •
-
How to calculate position-specific fastq stats?
written by O.rka
-
Anyone could help me with cutadapt?
written by ALEJANDRO •
-
What is the right design matrix in a Microarray data analysis?
written by Paul Paodumai •
-
Sorting a matrix rowname based on column of another matrix
written by Duckula •
-
Tutorial: raw count data- normalization
written by ghmdsr •
-
iterating through list of tuples #python
written by KVC_bioinfo
-
Cutadapt not trimming adapters in PE mode
written by tlorin
-
Extracting Specific Reads From Vcf File
written by Bnfoguy •
You can perform gene-wise normalization in limma, edgeR and DESeq2 (this is how conditional quantile normalization works), but the bigger questions are:
Answering #2 will also have answer the appropriate type of model to use.
It actually hard to explain, my goal is quite specific. Short answer is the first count data is dependent on the second, as I'm interested in the part of variance from matrix A, which are independent of the matrix B, I would like to normalize the matrix A against B.
OK, well then doing something like the following in DESeq2 may be the easiest:
Note that
ddsis aDESeqDataSetobject. You can do something similar in edgeR and limma.I should mention that using matrix B directly as I showed may be the opposite of what you want. This way will essentially bias estimates toward values with high matrix B values. If you instead want the exact opposite of this then you'll need to modify matrix B accordingly before using it (still divide it by its mean).