Browse code

Changed visualization for methylation distribution visualization

Divyagash authored on 20/09/2018 02:01:30
Showing 5 changed files

... ...
@@ -21,8 +21,6 @@ Imports:
21 21
   GenomeInfoDb,
22 22
   Biostrings,
23 23
   DT,
24
-  scales,
25
-  viridis,
26 24
   HDF5Array (>= 1.7.5)
27 25
 Suggests:
28 26
   BSgenome.Mmusculus.UCSC.mm10,
... ...
@@ -28,10 +28,8 @@ importFrom(annotatr,build_annotations)
28 28
 importFrom(annotatr,builtin_genomes)
29 29
 importFrom(annotatr,summarize_annotations)
30 30
 importFrom(bsseq,getCoverage)
31
-importFrom(scales,comma)
32 31
 importFrom(stats,aggregate)
33 32
 importFrom(stats,dbinom)
34 33
 importFrom(stats,sd)
35 34
 importFrom(utils,read.csv)
36 35
 importFrom(utils,read.table)
37
-importFrom(viridis,scale_fill_viridis)
... ...
@@ -25,8 +25,6 @@
25 25
 #'mbiasDirectory=system.file("extdata",package='scmeth')
26 26
 #'outDir <- system.file(package='scmeth')
27 27
 #'report(bs,outDir,Hsapiens,'hg38',mbiasDir=mbiasDirectory,small=TRUE)
28
-#'@importFrom scales comma
29
-#'@importFrom viridis scale_fill_viridis
30 28
 #'@importFrom HDF5Array loadHDF5SummarizedExperiment
31 29
 #'@import knitr
32 30
 #'@import DT
... ...
@@ -259,12 +259,11 @@ g
259 259
 methylDistMatrix <- scmeth::methylationDist(bs,subSample=params$nCpGs,offset=params$offset,coverageVec=covVec)
260 260
 
261 261
 g <- ggplot2::ggplot(methylDistMatrix, ggplot2::aes_string(x='Var2',
262
-                                                        y='Var1',fill='value'))
263
-g <- g + ggplot2::geom_tile(color="white", size=0.1)
264
-g <- g + viridis::scale_fill_viridis(name="Fraction of \n CpGs observed", label=scales::comma)
262
+                                                        y='value',group='Var1'))
263
+g <- g + ggplot2::geom_point()
264
+g <- g + ggplot2::geom_line()
265 265
 #g <- g + ggplot2::coord_equal() # Adds same aspect ratio
266 266
 g <- g + ggplot2::labs(x=NULL, y=NULL, title="Methylation Distribution")
267
-g <- g + ggthemes::theme_tufte(base_family="Helvetica")
268 267
 g <- g + ggplot2::theme(plot.title=ggplot2::element_text(hjust=0))
269 268
 g <- g + ggplot2::theme(axis.text.y=ggplot2::element_blank())
270 269
 # get rid of tick marks
... ...
@@ -70,7 +70,6 @@ Code chunk below shows how it can be loaded.
70 70
 </p>
71 71
 
72 72
 ```{r, eval=FALSE}
73
-# Works in devtools::check but not in bioconductor
74 73
 directory <- system.file("extdata","bismark_data",package='scmeth')
75 74
 bsObject <- HDF5Array::loadHDF5SummarizedExperiment(directory)
76 75
 ```