I am using an approach where I fit the model using all the probe sets and then subsequently at the eBayes step, set a filter and then on the eBayes object use "plotSA" to see if there is a monotonous pattern between the residual standard deviation and log average expression. Would this be a reasonable approach? Below is the code I am using:
no2_lm=lmFit(bryois_no2,no2_design_mat)
no2_lmfit_eb=eBayes(no2_lm[no2_lm$Amean > 8,],trend = TRUE,robust = TRUE)
plotSA(no2_lmfit_eb,main="Plot of residual standard deviation vs Avg exp for fitted model using NO2")
summary(decideTests(no2_lmfit_eb))
no2.table <- topTable(no2_lmfit_eb, sort.by = "F", n = 8793)
•
written
by
das2000sidd •
1 answer
Yes, it is a reasonable approach, if you are analysing a data type where a decreasing variance trend is expected.
Your reference to "probe-sets" suggests that you might be analyzing Affymetrix microarray data. If you are analyzing Affymetrix microarray data with RMA normalization, then a quadratic trend (increasing then decreasing) is expected. In that case, the procedure you propose will work fine but will be conservative in the sense that you will filter some probe-sets that might be significantly DE.
Log in to answer this question.
-
Bias in FDR distribution
written by Daniel •
-
Flitering Genes in RNA-Seq causes more significant FDR Adjusted P-Val DEGs?
written by Saran •
-
Can we get the results for Student t test and Welch t test by using limma?.
written by Minushi •
-
A check over my code using Limma DE, starting from normalized data
written by AlexStar •
-
Plotting residuals for limma lmfit model
written by ritz •
-
Error with gsms
written by smrutimayipanda •
-
R analysis for read counts using edgeR
written by Liftedkris •
-
Same Values In Limma/Microarray Differential Expression
written by Wan Fahmi •
-
Single-color Agilent array analyzing in R
written by Leite
-
Differential Gene Expression Analysis using data_RNA_Seq_v2_expression_median RSEM.Normalized
written by gokce.ouz
