### Name: arrayWeights ### Title: Array Quality Weights ### Aliases: arrayWeights arrayWeightsSimple ### Keywords: regression models ### ** Examples library(sma) # Subset of data from ApoAI case study in Limma User's Guide data(MouseArray) # Avoid non-positive intensities RG <- backgroundCorrect(mouse.data, method="half") MA <- normalizeWithinArrays(RG, mouse.setup) MA <- normalizeBetweenArrays(MA, method="Aq") targets <- data.frame(Cy3=I(rep("Pool",6)),Cy5=I(c("WT","WT","WT","KO","KO","KO"))) design <- modelMatrix(targets, ref="Pool") arrayw <- arrayWeightsSimple(MA, design) fit <- lmFit(MA, design, weights=arrayw) fit2 <- contrasts.fit(fit, contrasts=c(-1,1)) fit2 <- eBayes(fit2) # Use of array weights increases the significance of the top genes topTable(fit2)