R/STATegRa_fused.R

In lunacab/STATegRa: Classes and methods for multi-omics data integration

#' @export
#' @import Biobase
#' @title createOmicsExpressionSet
#' @aliases createOmicsExpressionSet,matrix-method
#'
#' @description
#' This function allow to the user to create a ExpressionSet object from a matrix representing an omics dataset. It allows to include the experimental design and annotation in the ExpressionSet object.
#'
#' @details
#' In Data matrix, samples has to be in columns and variables has to be in rows
#'
#' @param Data Omics data
#' @param pData Data associated with the samples/phenotype
#' @param pDataDescr Description of the phenotypic data
#' @param feaData Data associated with the variables/features annotation
#' @param feaDataDescr Description of the feature annotation
#' @return ExpressionSet with the data provided
#' @author Patricia Sebastian-Leon
#' @keywords datagen
#'
#' @examples
#' data(STATegRa_S3)
#' B1 <- createOmicsExpressionSet(Data=Block1.PCA,pData=ed.PCA,
#'                                pDataDescr=c("classname"))
#' B2 <- createOmicsExpressionSet(Data=Block2.PCA,pData=ed.PCA,
#'                                pDataDescr=c("classname"))
setGeneric(
    name="createOmicsExpressionSet",
    def=function(Data,pData=NULL,pDataDescr=NULL,feaData=NULL,feaDataDescr=NULL){
        standardGeneric("createOmicsExpressionSet")}
)

setMethod(
    f="createOmicsExpressionSet",
    signature=signature(
        Data="matrix"),
    definition=function(Data,pData,pDataDescr,feaData,feaDataDescr){
        ## Data
        exprs <- as.matrix(Data)
        ## Phenotypic data
        
        if(is.null(pData)){
            if (is.null(colnames(exprs))) {
                stop("If the data has no colnames, pData with names must be provided")
            }
            pData <- data.frame(Names=colnames(exprs))
            rownames(pData) <- colnames(exprs)
            pDataDescr <-c("Sample name")
        } else if (is.null(pDataDescr)) {
            pDataDescr <- colnames(pData)
        }
        pData <- as.data.frame(pData)
        if (is.null(rownames(pData)) | !all(rownames(pData)==colnames(exprs))){
            stop("pData rownames and exprs colnames must be the same")
        }
        metadata <- data.frame(labelDescription=pDataDescr)
        rownames(metadata) <- colnames(pData)
        phenoData <- new("AnnotatedDataFrame",data=pData,varMetadata=metadata)
        ## Feature data
        if(is.null(feaData)){
            if (is.null(rownames(exprs))) {
                stop("If the data has no rownames, feaData with names must be provided")
            }
            feaData <- data.frame(Name=rownames(exprs))
            rownames(feaData) <- rownames(exprs)
            feaDataDescr <- c("Feature name")
        } else if (is.null(feaDataDescr)) {
            feaDataDescr <- colnames(feaData)
        }
        fData <- as.data.frame(feaData)
        if (is.null(rownames(fData)) | !all(rownames(fData)==rownames(exprs))){
            stop("annot rownames and exprs rownames must be the same")
        }
        metadata2 <- data.frame(labelDescription=feaDataDescr)
        rownames(metadata2) <- colnames(fData)
        featureData <- new("AnnotatedDataFrame",data=fData,varMetadata=metadata2)
        ## ExpressionSet
        es <- new("ExpressionSet",exprs=exprs,phenoData=phenoData,featureData=featureData)
        return(es)
    }
)


#' @export
#' @importFrom utils browseURL
#' @title STATegRaUsersGuide
#'
#' @description Finds the location of the STATegRa User's Guide and optionally opens it.
#'
#' @param view Whether to open a browser
#' @return The path to the documentation
#' @author David Gomez-Cabrero
#'
#' @examples
#' STATegRaUsersGuide(view=FALSE)
STATegRaUsersGuide <- function(view=TRUE)
    #	Find and optionally view edgeR STATegRa Guide
    #	David Gomez-Cabrero, following Gordon Smyth example on edgeR
    #	23 May 2012.
{
    f <- system.file("doc", "STATegRa.html", package="STATegRa")
    if(view) {
        browseURL(f)
    }
    return(f)
}