How to decide the best threshold for filtering probes in limma

How to decide the best threshold for filtering probes in limma

0

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)


limma

• 23 views

Read more here: Source link