Question: Filtering 2-fold differentially expressed genes
0
gravatar for vinayjrao
7 weeks ago by
vinayjrao140
JNCASR, India
vinayjrao140 wrote:

Hi,

I am working with a microarray dataset, which is in the format -

Control      Control2      Treatment      Treatment2

The dataset has been log transformed. I want to know if there a way to filter the data such that I get genes whose fold difference is greater than 2 fold in Treatment as compared to Control.

For this purpose, I have already calculated the rowMeans() of each replicate to give

Control      Treatment

Thank you.

microarray R • 119 views
ADD COMMENTlink modified 7 weeks ago by Nicolas Rosewick8.0k • written 7 weeks ago by vinayjrao140
1

It is better to do DE analysis first, then after that if necessary you can include a Fold change filter. Try to use limma for your analysis, it has a very good manual.

ADD REPLYlink written 7 weeks ago by Benn7.5k
2
gravatar for Nicolas Rosewick
7 weeks ago by
Belgium, Brussels
Nicolas Rosewick8.0k wrote:

You should use limma to compute gene expression statistics (p-value, foldchange) and use this value to filter out your data.

https://bioconductor.org/packages/release/bioc/html/limma.html

Otherwise you can do something like this :

# control_samples # vector of control sample names
# treatment_samples # vector of treatment sample names
# M # gene expression matrix

fc <- rowMeans(M[,colnames(M) %in% control_samples]) / rowMeans(M[,colanames(M) %in% treatment_samples])

M.filtered <- M[ fc > 2,]
ADD COMMENTlink written 7 weeks ago by Nicolas Rosewick8.0k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1705 users visited in the last hour