## Load data
library('derfinder')
## Annotate regions
annotatedRegions <- annotateRegions(regions = genomeRegions$regions,
genomicState = genomicState$fullGenome, minoverlap = 1)
## Make venn diagram
venn <- vennRegions(annotatedRegions)
## Add title and choose text color
venn2 <- vennRegions(annotatedRegions, main = 'Venn diagram', counts.col =
'blue')
## Subset to only significant regions, so you don't have to annotate them
## again
venn3 <- vennRegions(annotatedRegions, subsetIndex =
genomeRegions$regions$significant == 'TRUE', main = 'Significant only')
Run the code above in your browser using DataLab