R/normScaleMarray.R

In maigesPack: Functions to handle cDNA microarray data, including several methods of data analysis

## Function to scale the microarray data (using marray function)
##
## Parameters: obj -> object of type maigesRaw to be normalized
##             ... -> additional parameters for function maNormScale
##
##
## Gustavo H. Esteves
## 15/05/07
##
##


normScaleMarray <- function(obj=NULL, ...) {
    
    
    ## converting the obj to an object of class limma
    if(class(obj) == "maigesRaw")
        toNorm <- as(obj, "marrayRaw")
    else if(class(obj) == "maiges")
        toNorm <- as(obj, "marrayNorm")
    else
        stop("The 'obj' object isn't of class 'maigesRaw' or 'maiges'.")
    
    
    
    ## Defining a new object for normalized data
    norm <- new("maiges")
    norm@BadSpots <- obj@BadSpots
    norm@UseSpots <- obj@UseSpots
    norm@GeneGrps <- obj@GeneGrps
    norm@Paths <- obj@Paths
    norm@Layout <- obj@Layout
    norm@Glabels <- obj@Glabels
    norm@Slabels <- obj@Slabels
    norm@Date <- date()
    
    ## Picking R and packages version information
    tmp <- sessionInfo()
    vInfo <- list()
    vInfo$R.version <- tmp$R.version$version.string
    vInfo$BasePacks <- tmp$basePkgs
    tmp1 <- NULL
    for (i in 1:length(tmp$otherPkgs))
        tmp1 <- c(tmp1, paste(tmp$otherPkgs[[i]]$Package, "version",
        tmp$otherPkgs[[i]]$Version))
    
    vInfo$AddPacks <- tmp1
    
    norm@V.info <- vInfo
    
    
    ## Merging UseSpots and BadSpots
    for(i in 1:dim(obj@UseSpots)[2])
        obj@UseSpots[, i] <- obj@UseSpots[, i] & !obj@BadSpots
    
    
    ## Geting samples where Ref is labelled with ch2 to do the (-1) M multiplier
    idxMult <- tolower(getLabels(obj, "Ref")) == "red"
    
    
    ## Catching the method of normalization
    add <- list(...)
    if(is.null(add$norm))
        type <- "globalMAD"
    else
        type <- add$norm
    
    
    
    ## Normalization by the limma method
    idxTmp <- apply(obj@UseSpots, 1, sum) == dim(obj)[2]
    tmp <- maNormScale(toNorm, subset=idxTmp)
    tmp1 <- as.matrix(maM(tmp))
    if (sum(idxMult) > 0)
        tmp1[, idxMult] <- (-1)*tmp1[, idxMult]
    norm@W <- tmp1
    norm@A <- as.matrix(maA(tmp))
    norm@Notes <- paste(obj@Notes, ", scale adjusted by ", type,
    " (from marray)", sep="")
    
    
    ## Returning the object
    return(norm)
    
}