RLE plot of DESeq2 normalisation / extract matrix of normalised counts from DESeqDataSet
1
0
Entering edit mode
f_pardow • 0
@felicitaspardow-22887
Last seen 22 months ago

Hi,

I am comparing different normalisation methods for my count table (I'm doing RNAseq analysis). One of them is DESeq2, the other one is RUVs. One way to compare them is to look at the unwanted variation that is still in the dataset after normalisation. This can be visualised using RLE plots.

RLE plots can be made using the function plotRLE. The main argument for this function is the object to be plotted. This can be a numeric matrix or a SeqExpressionSet object containing the gene expression.

For the RUVs normalisation making an RLE plot is simple as the normalised object is a SeqExpressionSet. However, for the DESeqDataSet, it is more complicated. Could somebody tell me how to extract the count matrix with normalised counts from the DESeqDataSet?

For reference, my code to plot the RLE plot for the RUVs normalisation is the following:

library(RUVSeq)
set <- newSeqExpressionSet(counts = as.matrix(counts), phenoData = sample_info)
setUQ <- betweenLaneNormalization(set, which="upper")
genes <- rownames(counts)
replicates <- makeGroups(control) # control is a column of sample_info and therefore information from the SeqExpressionSet
setRUVs9 <- RUVs(setUQ, genes, k = 9, replicates)
plotRLE(setRUVs9)


Thanks a lot in advance. Felicitas

deseq2 RLE • 971 views
0
Entering edit mode
@mikelove
Last seen 1 day ago
United States

Take a look at our help pages. There is a function counts with an argument normalized.