# Function to merge topTable output with original expression values

mergett.r

# Function to merge topTable output with original expression values
mergett <- function(tt, eset) {
    if (class(tt)!="data.frame") stop("tt should be a data.frame")
    if (class(eset)!="ExpressionSet") stop("eset should be an ExpressionSet")
    exprset <- data.frame(exprs(eset))
    exprset$ID <- rownames(exprset)
    rownames(exprset) <- NULL
    print(paste("nrow(tt):", nrow(tt)))
    print(paste("nrow(exprset):", nrow(exprset)))
    if (nrow(tt)==nrow(exprset)) {
        newtt <- merge(tt, exprset) # merge with intensities
    } else {
        stop("nrow(tt) != nrow(exprset)")
    }
    print(paste("nrow(newtt):", nrow(newtt)))
    newtt <- newtt[order(newtt$P.Val), ] #reorder by p-value
    rownames(newtt) <- NULL
    newtt
}