The weight parameters were optimized based on MA-plots such that spike-in controls show their expected fold change values. lmFit function was used to fit a linear model to probe intensities across arrays, and differential expression was calculated by empirical Bayes method using the eBayes function. Control of FDR was employed as correction for multiple testing.