Entering edit mode
5.6 years ago
shinobee
▴
60
One paper was mentioning removing genes with residual batch effects. What does that mean?
Is it something removing genes higher than 0.05 q values after batch correction?
pValuesComBat = f.pvalue(combat_data,mod,mod0)
qValuesComBat = p.adjust(pValuesComBat,method="BH")
combat_data = combat_data[qValuesComBat<0.05,]