set.seed(123)
## artifical sets -- affy probes grouped by protein family
library(hgu95av2.db)
map <- select(hgu95av2.db, keys(hgu95av2.db), "PFAM")
sets <- Filter(function(x) length(x) >= 10, split(map$PROBEID, map$PFAM))
universe <- unlist(sets, use.names=FALSE)
siggenes <- sample(universe, length(universe) / 20) ## simulate
sigsets <- Map(function(x, y) x[x %in% y], sets, MoreArgs=list(y=siggenes))
result <- hyperg(sets, sigsets, universe)
head(result)
Run the code above in your browser using DataLab