# Regional QA for Agilent ratio.image <- matrix( nr = max(X), nc = max(Y)) image.index.vector <- X + max(X)*(Y-1) log.standard <- apply( log.intensity, 1, mean, trim= 0.1 ) for ( ii in 1:dim(log.intensity)[2]) { png(paste(colnames(log.intensity)[ii],'png',sep='.'), width=640,height=640) ratio.image[ image.index.vector ] <- log.intensity[,ii] - log.standard ratio.image[ ratio.image > .8] <- 1 ratio.image[ ratio.image < -.8] <- -1 image( ratio.image, col= rainbow(16, end = .6)[c(16, 13:4 ,1)],xaxt='n',yaxt='n' ) dev.off() }