# load experiment
load(system.file("extdata/E-GEOD-21942.topTable.RData", package = "ROntoTools"))
fc <- top$logFC[top$adj.P.Val <= .01]
names(fc) <- top$entrez[top$adj.P.Val <= .01]
ref <- top$entrez
# load the set of pathways
kpg <- keggPathwayGraphs("hsa")
kpg <- setEdgeWeights(kpg)
kpg <- setNodeWeights(kpg, defaultWeight = 1)
# perform the pathway analysis
pDisRes <- pDis(fc, graphs = kpg, ref = ref, nboot = 100, verbose = TRUE)
# obtain summary of results
head(summary(pDisRes))
kpn <- keggPathwayNames("hsa")
head(summary(pDisRes))
head(summary(pDisRes, pathNames = kpn, totalpDis = FALSE,
pORA = FALSE, comb.pv = NULL, order.by = "pDis"))
Run the code above in your browser using DataLab