data(GeneOverlap)
go.obj <- newGeneOverlap(hESC.ChIPSeq.list$H3K27me3,
hESC.RNASeq.list$"Exp Medium",
genome.size=gs.RNASeq)
getGenomeSize(go.obj)
v.gs <- c(12e3, 14e3, 16e3, 18e3, 20e3)
setNames(sapply(v.gs, function(g) {
setGenomeSize(go.obj) <- g
go.obj <- testGeneOverlap(go.obj)
getPval(go.obj)
}), v.gs)
Run the code above in your browser using DataLab