Filtering Gene by mean expression
0
0
Entering edit mode
5 months ago
synat.keam ▴ 50

Dear All,

Sorry for another post. I have RNAseq data and used DESeq2 pipeline for differntially expressed genes. Now I am doing WGCNA and the author suggest to filter genes based on mean expression rather than variance. I am wondering whether anyone has ever filtered genes based on mean expression and could suggest some code chunks.

Kind Regards,

Synat,

Gene meanExpression Filter • 514 views
ADD COMMENT
1
Entering edit mode

If you need only exp values based on mean use cufflinks.

ADD REPLY
0
Entering edit mode

Hi Marcro,

I'm not sure whether DESeq2 filter genes based on mean expression.

Kashiff, I will try your suggestion. Thank for you both of you for taking time to respond.

Regards,

synat

ADD REPLY
2
Entering edit mode

Hey! DESeq2 does not filter genes with low read counts across samples by default, but you need to perform such filtering step, simple as this:

#after generating a dds object
keep <- rowSums(counts(dds)) >= 10
dds <- dds[keep,]
ADD REPLY
0
Entering edit mode

You mean filtering genes with very low read counts as also suggested in the DESeq2 pipeline?

ADD REPLY

Login before adding your answer.

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