Sample code for thresher

thresher-demo

#' Load the package
library(Thresher)

#'  Generate random data
set.seed(3928270)
ranData <- matrix(rnorm(100*12), ncol=12)
colnames(ranData) <- paste("G", 1:12, sep='')


#' Perform dim reduction
thresh1 <- Thresher(ranData)
#' Outlier removal and clustering
reap1 <- Reaper(thresh1)


#' Noise features
colnames(ranData)[!reap1@keep]

#' Plot remaining features projected onto PC space
plot(reap1)

#' Optimal number of groups 
reap1@nGroups

#' Heat map of data using only the good features
heat(reap1)