Browse code

Simplified code for plotting.

Robert Ivánek authored on 17/11/2020 11:04:06
Showing 1 changed files

... ...
@@ -106,7 +106,10 @@ hasBiomartConnection <- (!is(try(download.file("http://www.biomart.org", tempfil
106 106
                          !is(try(biomaRt::listMarts(), silent=TRUE), "try-error"))
107 107
 options(timeout=oto)
108 108
 
109
-## Uncommenting this helps when the UCSC server has a hickup but still lets you connect:
109
+## ## clear BiomaRt cache
110
+## biomaRt::biomartCacheClear()
111
+
112
+## ## Uncommenting this helps when the UCSC server has a hickup but still lets you connect:
110 113
 ## hasUcscConnection <- !is(try(rtracklayer::browserSession(), silent=TRUE), "try-error") && 
111 114
 ##   !is(try(IdeogramTrack(genome="hg19", chromosome=7), silent=TRUE), "try-error")
112 115
 ```
... ...
@@ -942,22 +945,21 @@ set is plotted for the horizon type for the sake of clarity.]
942 945
 ```{r typeGroupedPlots, fig.width=7.5, fig.height=6, echo=FALSE, results='hide'}
943 946
 pushViewport(viewport(layout=grid.layout(nrow=9, ncol=1)))
944 947
 i <- 1
945
-for(t in c("a", "s", "confint", "smooth", "histogram", "boxplot", "heatmap", "horizon")) {
946
-    pushViewport(viewport(layout.pos.col=((i-1)%%1)+1, layout.pos.row=((i-1)%/%1)+1))
947
-    if(t != "horizon") {
948
-        names(dTrack) <- t
949
-        plotTracks(dTrack, type=t, add=TRUE, cex.title=0.8, groups=rep(1:2, each=3), margin=0.5)
950
-    } else {
951
-        plotTracks(dtHoriz[c(1,8),], type="horizon", add=TRUE, cex.title=0.8, margin=0.5, showAxis=FALSE, horizon.origin=0.3,
952
-                   groups=1:2)
953
-    }
954
-    i <- i+1
955
-    popViewport(1)
948
+for(t in c("a", "s", "confint", "smooth", "histogram", "boxplot", "heatmap", "horizon", "hor. type")) {
949
+  pushViewport(viewport(layout.pos.col=((i-1)%%1)+1, layout.pos.row=((i-1)%/%1)+1))
950
+  names(dTrack) <- t
951
+  if(!t %in% c("horizon", "hor. type")) {
952
+    plotTracks(dTrack, type=t, add=TRUE, cex.title=0.8, groups=rep(1:2, each=3), margin=0.5)
953
+  } else if (t == "horizon") {
954
+    plotTracks(dtHoriz[c(1,8),], type="horizon", add=TRUE, cex.title=0.8, margin=0.5, showAxis=FALSE, horizon.origin=0.3,
955
+               groups=1:2)
956
+  } else {
957
+    plotTracks(dTrack, type="histogram", stackedBars=FALSE, add=TRUE, cex.title=0.8, groups=rep(1:2, each=3), margin=0.5)
958
+  }
959
+  i <- i+1
960
+  popViewport(1)
956 961
 }
957
-pushViewport(viewport(layout.pos.col=((i-1)%%1)+1, layout.pos.row=((i-1)%/%1)+1))
958
-names(dTrack) <- "hor. hist."
959
-plotTracks(dTrack, type="histogram", stackedBars=FALSE, add=TRUE, cex.title=0.8, groups=rep(1:2, each=3), margin=0.5)
960
-popViewport(2)
962
+popViewport(1)
961 963
 names(dTrack) <- "uniform"
962 964
 ```
963 965