Can loess.normalization be used for normalizing RNASeq performed with two ERCC mixes? · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Can loess.normalization be used for normalizing RNASeq performed with two ERCC mixes?

Hi,

I am working on RNAseq data where my colleague has used two ERCC mixes , Mix 1 for tests and Mix 2 for WT. https://www.thermofisher.com/order/catalog/product/4456740 So, I have 92 ERCC transcripts with different expected fold change. I am not sure how to provide the information for two mixes when I normalize the data using loess.normalization. I know I should provide a subset of rows (ERCC transcript) and perform normalization to the whole matrix in R as shown below.

result <- format(loess.normalize(mat,mat[15337:15427,1:4], epsilon = 10^-2, maxit = 1, log.it =TRUE, verbose = TRUE, span = 2/3, family.loess ="symmetric"),digits=2,nsmall=2)

where, mat is the matrix of FPKM values for WT and three different tests, mat[15337:15427,1:4] are the rows corresponding to ERCC FPKM.

I have two questions: 1) Does loess.normalize normalizes other genes of the matrix column wise corresponding to the FPKM values of ERCC for that column or it also takes care of fold change 2)How to provide fold change information, if it doesn't consider fold change

Thanks for any help.

rna-seq normalization ercc

• 3,478 views • link

updated 8.6 years ago by Kevin Blighe

• written 9.6 years ago by manjovial •

Not sure about the answer, but maybe you can find a hint at https://www.bioconductor.org/packages/erccdashboard ?

• 0 views • link

written 8.6 years ago by Charles Plessy

0 answers

No answers yet.

Log in to answer this question.

More posts like this

Tutorial: Bulk RNA-seq: Different Expression Analysis with DEseq2
written 2.6 years ago by Julia Ma •

Content taken verbatim from: https://omicverse.readthedocs.io/en/latest/Tutorials-bulk/t_deseq2/ An important task of bulk rna-seq analysis is the different expression , which we can perform with omicverse. For different …
RiP RLE normalisation using spike-in peaks in DiffBind for ChIP-seq
written 3.7 years ago by spg •

Hi, I have two conceptual/technical questions about spike-in normalisation of ChIP-seq data. We did spike-in for PolII ChIP-seq from a mouse cell line with ActiveMotif …
Tutorial: Unsupervised clustering on gene expression data
written 4.8 years ago by Hamid Ghaedi

Clustering is a data mining method to identify unknown possible groups of items solely based on intrinsic features and no external variables. Basically, clustering includes …
Deseq2 Single-Cell analysis with two variables
written 7.3 years ago by michele.tebaldi.92 •

Hello everyone. I'm performing a single-cell RNA-seq on a MusMusculus dataset following the DESeq2 instructions in https://bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html In my experiment there are two variables of …
RNA-Seq time series analysis using a DESeq2 spline approach yields far too many significant genes
written 7.6 years ago by stu111538

Hi, I performed RNA-Seq time series analysis on my set of 12 patients under therapy. Since I only have 5 time points (and for some …
Plot heatmap of differentially expressed genes identified by DESeq2 a R coding problem
written 9.5 years ago by Wet&DryImmunology

I was trying to a differential gene expression analysis by using Deseq2 [1] with samples like this [2], and would like to plot a heatmap …
Multi-factor time series design model in DESeq2 analysis
written 11.0 years ago by prppmet •

I have the following data: a). 3 strains (wild type & two mutants); b). 2 conditions: untreated & treated; c). 3 time points (min): 0, …