...
|
...
|
@@ -139,7 +139,7 @@ g
|
139
|
139
|
|
140
|
140
|
```{r cpg_discretization,echo=FALSE,warning=FALSE,message=FALSE,fig.height=5, fig.width=5, eval=TRUE,fig.align='center'}
|
141
|
141
|
# CpG Discretization
|
142
|
|
-discretizedCpG <- cpgDiscretization(bs,covVec,subSample=params$nCpGs,offset=params$offset)
|
|
142
|
+discretizedCpG <- cpgDiscretization(bs,subSample=params$nCpGs,offset=params$offset,coverageVec=covVec)
|
143
|
143
|
#percentDiscarded <- discretizedCpG[3]
|
144
|
144
|
percentDiscarded <- discretizedCpG[2]
|
145
|
145
|
sampleNames <- bsseq::sampleNames(bs)
|
...
|
...
|
@@ -200,7 +200,7 @@ g
|
200
|
200
|
## Methylation Distribution
|
201
|
201
|
```{r meth_dist,echo=FALSE,warning=FALSE,message=FALSE,fig.height=7, fig.width=10,eval=TRUE,fig.align='center'}
|
202
|
202
|
# Methylation distribution for all the cells
|
203
|
|
-methylDistMatrix <- scmeth::methylationDist(bs,covVec,subSample=params$nCpGs,offset=params$offset)
|
|
203
|
+methylDistMatrix <- scmeth::methylationDist(bs,subSample=params$nCpGs,offset=params$offset,coverageVec=covVec)
|
204
|
204
|
|
205
|
205
|
g <- ggplot2::ggplot(methylDistMatrix, ggplot2::aes_string(x='Var2',
|
206
|
206
|
y='Var1',fill='value'))
|