|
...
|
...
|
@@ -41,23 +41,23 @@
|
|
41
|
41
|
#'
|
|
42
|
42
|
#' set.seed(1000)
|
|
43
|
43
|
#'
|
|
44
|
|
-#' #Example 1: create 20 sets of non-overlapping random regions and plot them all. Add a coverage plot on top.
|
|
45
|
|
-#' kp <- plotKaryotype("hg19", plot.type=1, chromosomes=c("chr1", "chr2"))
|
|
46
|
|
-#'
|
|
47
|
|
-#' all.regs <- GRanges()
|
|
48
|
|
-#'
|
|
49
|
|
-#' nreps <- 20
|
|
50
|
|
-#' for(i in 1:nreps) {
|
|
51
|
|
-#' regs <- createRandomRegions(nregions = 100, length.mean = 10000000, length.sd = 1000000,
|
|
52
|
|
-#' non.overlapping = TRUE, genome = "hg19", mask=NA)
|
|
53
|
|
-#' all.regs <- c(all.regs, regs)
|
|
54
|
|
-#' kpPlotRegions(kp, regs, r0 = (i-1)*(0.8/nreps), r1 = (i)*(0.8/nreps), col="#AAAAAA")
|
|
55
|
|
-#' }
|
|
56
|
|
-#'
|
|
57
|
|
-#' kpPlotCoverage(kp, all.regs, ymax = 20, r0=0.8, r1=1, col="#CCCCFF")
|
|
58
|
|
-#' kpAxis(kp, ymin = 0, ymax= 20, numticks = 2, r0 = 0.8, r1=1)
|
|
59
|
|
-#'
|
|
60
|
|
-#'
|
|
|
44
|
+# #Example 1: create 20 sets of non-overlapping random regions and plot them all. Add a coverage plot on top.
|
|
|
45
|
+# kp <- plotKaryotype("hg19", plot.type=1, chromosomes=c("chr1", "chr2"))
|
|
|
46
|
+#
|
|
|
47
|
+# all.regs <- GRanges()
|
|
|
48
|
+#
|
|
|
49
|
+# nreps <- 20
|
|
|
50
|
+# for(i in 1:nreps) {
|
|
|
51
|
+# regs <- createRandomRegions(nregions = 100, length.mean = 10000000, length.sd = 1000000,
|
|
|
52
|
+# non.overlapping = TRUE, genome = "hg19", mask=NA)
|
|
|
53
|
+# all.regs <- c(all.regs, regs)
|
|
|
54
|
+# kpPlotRegions(kp, regs, r0 = (i-1)*(0.8/nreps), r1 = (i)*(0.8/nreps), col="#AAAAAA")
|
|
|
55
|
+# }
|
|
|
56
|
+#
|
|
|
57
|
+# kpPlotCoverage(kp, all.regs, ymax = 20, r0=0.8, r1=1, col="#CCCCFF")
|
|
|
58
|
+# kpAxis(kp, ymin = 0, ymax= 20, numticks = 2, r0 = 0.8, r1=1)
|
|
|
59
|
+#
|
|
|
60
|
+#
|
|
61
|
61
|
#' #Example 2: Do the same with a single bigger set of possibly overlapping regions
|
|
62
|
62
|
#'
|
|
63
|
63
|
#' kp <- plotKaryotype("hg19", plot.type=1, chromosomes=c("chr1", "chr2"))
|
