...
|
...
|
@@ -35,7 +35,6 @@ nSamples <- length(Biobase::sampleNames(bs))
|
35
|
35
|
```{r read_metrics,echo=FALSE,warning=FALSE,message=FALSE,fig.width=8, fig.height=7,fig.align='center'}
|
36
|
36
|
# Methylation distribution for all the cells
|
37
|
37
|
dat <- scmeth::readmetrics(bs)
|
38
|
|
-<<<<<<< HEAD
|
39
|
38
|
|
40
|
39
|
SummaryTable <- dat
|
41
|
40
|
SummaryTable$sample <- sub('\\..*$', '', SummaryTable$sample)
|
...
|
...
|
@@ -45,11 +44,8 @@ Summary_read <- summary(dat[,c('total','mapped','unmapped')])
|
45
|
44
|
|
46
|
45
|
|
47
|
46
|
|
48
|
|
-=======
|
49
|
|
-o <- order(dat$total,dat$sample)
|
50
|
|
-sampleOrder <- dat$sample[o]
|
51
|
|
-
|
52
|
|
->>>>>>> upstream/master
|
|
47
|
+
|
|
48
|
+
|
53
|
49
|
m <- reshape2::melt(dat[,c("sample","mapped","unmapped")], id.vars="sample",
|
54
|
50
|
variable.name="Mapping_status")
|
55
|
51
|
m$sample <- factor(m$sample,levels=sampleOrder)
|
...
|
...
|
@@ -130,11 +126,9 @@ g
|
130
|
126
|
```{r bs_conversion,echo=FALSE,warning=FALSE,message=FALSE,fig.height=5,fig.width=5,fig.align='center'}
|
131
|
127
|
bscDf <- scmeth::bsConversionPlot(bs)
|
132
|
128
|
|
133
|
|
-<<<<<<< HEAD
|
|
129
|
+
|
134
|
130
|
SummaryTable$Bisulfite.converion.Rate <- round(bscDf$bsc,4)
|
135
|
131
|
|
136
|
|
-=======
|
137
|
|
->>>>>>> upstream/master
|
138
|
132
|
g <- ggplot2::ggplot(bscDf, ggplot2::aes_string(x="''", y='bsc'))
|
139
|
133
|
g <- g+ggplot2::geom_boxplot()
|
140
|
134
|
g <- g+ggplot2::ylim(max(min(bscDf$bsc)-0.02,0),min(max(bscDf$bsc)+0.02,1))
|