if (require(clipper) & require(ALL)){
k <- pathways("hsapiens", "kegg")
paths <- convertIdentifiers(k[1:5], "entrez")
genes <- unlist(lapply(paths, nodes))
data(ALL)
all <- as.matrix(exprs(ALL[1:length(genes),1:20]))
classes <- c(rep(1,10), rep(2,10))
rownames(all) <- genes
runClipperMulti(paths, all, classes, "mean", pathThr=0.1)
}
Run the code above in your browser using DataLab