R/juxta.R

juxta <-
function(chrname="chr22",...)
  {
  require(ERBS)
  data(HepG2)
  data(GM12878)
  require(ggbio)
  require(GenomicRanges)
  ap1 = autoplot(GenomicRanges::subset(HepG2,seqnames==chrname))
  ap2 = autoplot(GenomicRanges::subset(GM12878,seqnames==chrname))
  tracks(HepG2=ap1,Bcell=ap2,...) 
 }
melsiddieg/epigenoviz documentation built on May 22, 2019, 6:50 p.m.