Some genes are filtered after DESeq2 normalization
1
0
Entering edit mode
yura.song • 0
@yurasong-20744
Last seen 4.4 years ago
Belgium

I am using DESeq2 for finding differential expressed genes.

Originally, there were around 49,000 genes in rawcount of expression, however, when I performed these steps, the number of genes were almost half of rawcount.

ddsk14 <- estimateSizeFactors(ddsk14) sizeFactors(dds_k14)

K14CTRL1 K14CTRL2 K14Exp1 K14Exp2 1.0163136 1.0524673 0.9809502 0.9752045

View(counts(dds_k14, normalized = TRUE))

I tried to understand why the genes are filtered out on normalization step, it was hard to grab with the paper and vignettes. Could you explain to me a little bit easier way about this?

Thanks in advance!

deseq2 normalization • 384 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 21 hours ago
United States

Genes are not filtered by normalization.

Take a look at:

dim(counts(dds))

or

dim(counts(dds, normalized=TRUE))

ADD COMMENT

Login before adding your answer.

Traffic: 705 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