# NOT RUN {
require(hgu133plus2.db)
## Two class unpaired comparison
## load data
data(marty)
## filtering data
marty <- expFilter(marty, threshold=3.5, graph=FALSE)
##Class label 0/1
marty.type.num <- ifelse(marty.type.cl=="Her2+",0,1)
## run sam analysis
gsaOUT <- runGSA(marty, marty.type.num ,
gmtfile="./c2.kegg.v2.5.symbols.gmt", chip="hgu133plus2")
# }
Run the code above in your browser using DataLab