Normalizzation indipentent data from RNAAeq experiment
2
0
Entering edit mode
aristotele_m ▴ 30
@aristotele_m-6821
Last seen 4.3 years ago
Italy

Dear All,

I have 5 RNA-Seq experiments  and I want to compare the expression values. I have the raw counts. How
can normalize using deseq2? I can't draw any design on this because there is no meaning.

What do you suggest ?
Thanks in advance for any help

 

rna normalization • 685 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 12 hours ago
United States

If you only want to compare across samples (not across genes), do this (you have to make up an empty column data):

coldata = DataFrame(row.names=1:ncol(counts))
dds = DESeqDataSetFromMatrix(counts, coldata, ~ 1)
dds = estimateSizeFactors(dds)
counts(dds, normalized=TRUE)

(edit) but this leads me to wonder: what are you trying to do with these normalized counts?

ADD COMMENT
0
Entering edit mode
aristotele_m ▴ 30
@aristotele_m-6821
Last seen 4.3 years ago
Italy

 

These samples are create after subsampling one big RNASeq experiments. I want to be sure the proportion of the expression are the same. So, the simple way is to verify if the  gene I'm studing are remain with the same ratio.

 

 

Thanks for the help!!

ADD COMMENT

Login before adding your answer.

Traffic: 169 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6