if (require(clipper) & require(ALL)){
k <- pathways("hsapiens", "kegg")
path <- convertIdentifiers(k$'Chronic myeloid leukemia', "entrez")
genes <- nodes(path)
data(ALL)
all <- as.matrix(exprs(ALL[1:length(genes),1:20]))
classes <- c(rep(1,10), rep(2,10))
rownames(all) <- genes
runClipper(path, all, classes, "mean", pathThr=0.1)
}
Run the code above in your browser using DataCamp Workspace