R/unsupervised3.R

In sean-cho/toolkit: Toolkit for data manipulation and plotting.

#' Unsupervised heatmap by standard deviation
#' 
#' \code{unsupervised.3} Performs unsupervised clustering on the most
#'    variable features using heatmap.3 which allows for multiple Row/ColSideCol.
#'    
#' @param x Numeric matrix with samples in columns and features in rows.
#' @param n Number of features.
#' @param col Color for heatmap. Defaults to blue-black-yellow color scheme.
#' @param get_hm Logical vector whether to return heatmap with feature names.
#' @param main Title.
#' @param ... Passes arguments to \code{heatmap.3}
#' 
#' @details Finds the most variable features across all samples by standard
#'    deviation, clusters them, and plots it on a heatmap using the
#'    \code{heatmap.3} function.
#'    
#' @return A heatmap.3 object and plot.
#' @export

unsupervised.3 = function(x,n,col=colorRampPalette(c("blue","black","yellow"))(50),
                          get_hm = FALSE,
                          main = NULL,...){
  dumsd = apply(x,1,sd,na.rm=TRUE)
  dumrows = dim(x)[1]
  dum = x[dumsd>=quantile(dumsd,(dumrows-n)/dumrows,na.rm=TRUE),]
  if(is.null(main)){
    Main = paste("Top ", round(n)," Probes\nN = ", ncol(x), sep="")
  } else {
    Main = main
  }
  hm <- myheatmap3(dum,main=Main,col = col,...)
  f <- rev(colnames(hm$carpet))
  if(get_hm){
    return(list(features = f, hm = hm))
  } else{
    return(f)
  }
}