Browse code

Updated version to match Bioconductor

Joshua D. Campbell authored on 28/12/2021 19:47:07
Showing 120 changed files

... ...
@@ -1,6 +1,6 @@
1 1
 Package: celda
2 2
 Title: CEllular Latent Dirichlet Allocation
3
-Version: 1.9.3
3
+Version: 1.10.0
4 4
 Authors@R: c(person("Joshua", "Campbell", email = "camp@bu.edu",
5 5
     role = c("aut", "cre")),
6 6
     person("Shiyi", "Yang", email="syyang@bu.edu", role = c("aut")),
... ...
@@ -1,3 +1,7 @@
1
+# celda v1.10.0 (2021-12-28)
2
+
3
+* Update to match Bioconductor release version
4
+
1 5
 # celda v1.9.3 (2021-10-04)
2 6
 
3 7
 * Fixed bug in checking background matrix with decontX
... ...
@@ -18,12 +18,19 @@ if (!requireNamespace("BiocManager", quietly = TRUE))
18 18
 BiocManager::install("celda")
19 19
 ```
20 20
 
21
-To install the development version (R >= 3.6) of **celda** from GitHub using `devtools`:
21
+The latest stable version  of **celda** can be installed from GitHub using `devtools`:
22
+
22 23
 ```
23 24
 library(devtools)
24 25
 install_github("campbio/celda")
25 26
 ```
26 27
 
28
+The development version of **celda** can also be installed from GitHub using `devtools`:
29
+```
30
+library(devtools)
31
+install_github("campbio/celda@devel")
32
+```
33
+
27 34
 **NOTE** For MAC OSX users, `devtools::install_github()` requires installation of **libgit2.** This can be installed via homebrew:
28 35
 ```
29 36
 brew install libgit2
... ...
@@ -71,7 +71,7 @@
71 71
       </button>
72 72
       <span class="navbar-brand">
73 73
         <a class="navbar-link" href="index.html">celda</a>
74
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
74
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
75 75
       </span>
76 76
     </div>
77 77
 
... ...
@@ -71,7 +71,7 @@
71 71
       </button>
72 72
       <span class="navbar-brand">
73 73
         <a class="navbar-link" href="index.html">celda</a>
74
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
74
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
75 75
       </span>
76 76
     </div>
77 77
 
... ...
@@ -71,7 +71,7 @@
71 71
       </button>
72 72
       <span class="navbar-brand">
73 73
         <a class="navbar-link" href="index.html">celda</a>
74
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
74
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
75 75
       </span>
76 76
     </div>
77 77
 
... ...
@@ -31,7 +31,7 @@
31 31
       </button>
32 32
       <span class="navbar-brand">
33 33
         <a class="navbar-link" href="../../index.html">celda</a>
34
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
34
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
35 35
       </span>
36 36
     </div>
37 37
 
... ...
@@ -91,14 +91,14 @@
91 91
 
92 92
       
93 93
 
94
-      </header><script src="celda_pbmc3k_files/header-attrs-2.7/header-attrs.js"></script><script src="celda_pbmc3k_files/kePrint-0.0.1/kePrint.js"></script><link href="celda_pbmc3k_files/lightable-0.0.1/lightable.css" rel="stylesheet">
94
+      </header><script src="celda_pbmc3k_files/header-attrs-2.7/header-attrs.js"></script><script src="celda_pbmc3k_files/accessible-code-block-0.0.1/empty-anchor.js"></script><script src="celda_pbmc3k_files/kePrint-0.0.1/kePrint.js"></script><link href="celda_pbmc3k_files/lightable-0.0.1/lightable.css" rel="stylesheet">
95 95
 <div class="row">
96 96
   <div class="col-md-9 contents">
97 97
     <div class="page-header toc-ignore">
98 98
       <h1 data-toc-skip>Celda - Analysis of PBMC3K</h1>
99 99
                         <h4 class="author">Joshua Campbell, Zhe Wang</h4>
100 100
             
101
-            <h4 class="date">Compiled October 04, 2021</h4>
101
+            <h4 class="date">Compiled December 28, 2021</h4>
102 102
       
103 103
       <small class="dont-index">Source: <a href="https://github.com/campbio/celda/blob/master/vignettes/articles/celda_pbmc3k.Rmd"><code>vignettes/articles/celda_pbmc3k.Rmd</code></a></small>
104 104
       <div class="hidden name"><code>celda_pbmc3k.Rmd</code></div>
... ...
@@ -124,12 +124,12 @@
124 124
 <span class="fu"><a href="https://rdrr.io/pkg/SingleCellExperiment/man/assays.html">counts</a></span><span class="op">(</span><span class="va">sce</span><span class="op">)</span> <span class="op">&lt;-</span> <span class="fu">as</span><span class="op">(</span><span class="fu"><a href="https://rdrr.io/pkg/SingleCellExperiment/man/assays.html">counts</a></span><span class="op">(</span><span class="va">sce</span><span class="op">)</span>, <span class="st">"dgCMatrix"</span><span class="op">)</span></code></pre></div>
125 125
 <p>If you have the <a href="https://bioconductor.org/packages/release/bioc/html/singleCellTK.html" target="_blank">singleCellTK</a> package installed, then this dataset can be imported and converted with a single command:</p>
126 126
 <div class="sourceCode" id="cb2"><pre class="downlit sourceCode r">
127
-<code class="sourceCode R"><span class="kw"><a href="https://rdrr.io/r/base/library.html">library</a></span><span class="op">(</span><span class="va"><a href="https://compbiomed.github.io/sctk_docs/">singleCellTK</a></span><span class="op">)</span>
128
-<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/importExampleData.html">importExampleData</a></span><span class="op">(</span><span class="st">"pbmc3k"</span><span class="op">)</span></code></pre></div>
127
+<code class="sourceCode R"><span class="kw"><a href="https://rdrr.io/r/base/library.html">library</a></span><span class="op">(</span><span class="va"><a href="https://www.camplab.net/sctk/">singleCellTK</a></span><span class="op">)</span>
128
+<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/importExampleData.html">importExampleData</a></span><span class="op">(</span><span class="st">"pbmc3k"</span><span class="op">)</span></code></pre></div>
129 129
 <p>To get your own data into a <code>SingleCellExperiment</code> object, the <code>singleCellTK</code> package has several importing functions for different preprocessing tools including CellRanger, STARsolo, BUStools, Optimus, DropEST, SEQC, and Alevin/Salmon. For example, the following code can be used as a template to read in multiple samples processed with CellRanger:</p>
130 130
 <div class="sourceCode" id="cb3"><pre class="downlit sourceCode r">
131
-<code class="sourceCode R"><span class="kw"><a href="https://rdrr.io/r/base/library.html">library</a></span><span class="op">(</span><span class="va"><a href="https://compbiomed.github.io/sctk_docs/">singleCellTK</a></span><span class="op">)</span>
132
-<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/importCellRanger.html">importCellRanger</a></span><span class="op">(</span>sampleDirs <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html">c</a></span><span class="op">(</span><span class="st">"path/to/sample1/"</span>, <span class="st">"path/to/sample2/"</span><span class="op">)</span><span class="op">)</span></code></pre></div>
131
+<code class="sourceCode R"><span class="kw"><a href="https://rdrr.io/r/base/library.html">library</a></span><span class="op">(</span><span class="va"><a href="https://www.camplab.net/sctk/">singleCellTK</a></span><span class="op">)</span>
132
+<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/importCellRanger.html">importCellRanger</a></span><span class="op">(</span>sampleDirs <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html">c</a></span><span class="op">(</span><span class="st">"path/to/sample1/"</span>, <span class="st">"path/to/sample2/"</span><span class="op">)</span><span class="op">)</span></code></pre></div>
133 133
 <p><strong>Note:</strong> As a reminder, you can view the assays, column annotation, and row annotation stored in the SCE with the commands <code>assays(sce)</code>, <code>colData(sce)</code>, and <code>rowData(sce)</code>, respectively.</p>
134 134
 <p>Finally, we set the <code>rownames</code> of the SCE to the gene symbol:</p>
135 135
 <div class="sourceCode" id="cb4"><pre class="downlit sourceCode r">
... ...
@@ -139,33 +139,33 @@
139 139
 <h1 class="hasAnchor">
140 140
 <a href="#quality-control" class="anchor"></a>Quality Control</h1>
141 141
 <p>Quality control and filtering of cells is often needed before down-stream analyses such as dimensionality reduction and clustering. Typical filtering procedures include exclusion of poor quality cells with low numbers of counts/UMIs, estimation and removal of ambient RNA, and identification of potential doublet/multiplets. Many tools and packages are available to perform these operations and users are free to apply their tool(s) of choice as the celda clustering functions will work with any matrix stored in an SCE object. The celda package does contain a Bayesian method called <a href="https://doi.org/10.1186/s13059-020-1950-6" target="_blank">decontX</a> to estimate and remove transcript contamination in individual cells in a scRNA-seq dataset.</p>
142
-<p>To perform QC, we suggest using the <code>runCellQC</code> function in <code>singleCellTK</code> package. This is a wrapper for several methods for calculation of QC metrics, doublet detection, and estimation of ambient RNA (including decontX). Below is a quick example of how to perform standard QC before applying celda. If you have another preferred approach or your data has already been QC’ed, you can move to <a href="#featureselection">Feature selection</a> section. For this tutorial, we will only run one doublet detection algorithm and one decontamination algorithms. For a full list of algorithms that this function runs by default, see <code><a href="https://compbiomed.github.io/sctk_docs/reference/runCellQC.html">?runCellQC</a></code>. We will also quantify the percentage of mitochondrial genes in each cell as this is often used as a measure of cell viability.</p>
142
+<p>To perform QC, we suggest using the <code>runCellQC</code> function in <code>singleCellTK</code> package. This is a wrapper for several methods for calculation of QC metrics, doublet detection, and estimation of ambient RNA (including decontX). Below is a quick example of how to perform standard QC before applying celda. If you have another preferred approach or your data has already been QC’ed, you can move to <a href="#featureselection">Feature selection</a> section. For this tutorial, we will only run one doublet detection algorithm and one decontamination algorithms. For a full list of algorithms that this function runs by default, see <code><a href="https://rdrr.io/pkg/singleCellTK/man/runCellQC.html">?runCellQC</a></code>. We will also quantify the percentage of mitochondrial genes in each cell as this is often used as a measure of cell viability.</p>
143 143
 <div class="sourceCode" id="cb5"><pre class="downlit sourceCode r">
144
-<code class="sourceCode R"><span class="kw"><a href="https://rdrr.io/r/base/library.html">library</a></span><span class="op">(</span><span class="va"><a href="https://compbiomed.github.io/sctk_docs/">singleCellTK</a></span><span class="op">)</span>
144
+<code class="sourceCode R"><span class="kw"><a href="https://rdrr.io/r/base/library.html">library</a></span><span class="op">(</span><span class="va"><a href="https://www.camplab.net/sctk/">singleCellTK</a></span><span class="op">)</span>
145 145
 
146 146
 <span class="co"># Get list of mitochondrial genes</span>
147 147
 <span class="va">mito.genes</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/grep.html">grep</a></span><span class="op">(</span><span class="st">"^MT-"</span>, <span class="fu"><a href="https://rdrr.io/r/base/colnames.html">rownames</a></span><span class="op">(</span><span class="va">sce</span><span class="op">)</span>, value <span class="op">=</span> <span class="cn">TRUE</span><span class="op">)</span>
148 148
 
149 149
 <span class="co"># Run QC</span>
150
-<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/runCellQC.html">runCellQC</a></span><span class="op">(</span><span class="va">sce</span>, sample <span class="op">=</span> <span class="cn">NULL</span>, algorithms <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html">c</a></span><span class="op">(</span><span class="st">"QCMetrics"</span>, <span class="st">"scDblFinder"</span>, <span class="st">"decontX"</span><span class="op">)</span>, geneSetList <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html">list</a></span><span class="op">(</span>mito<span class="op">=</span><span class="va">mito.genes</span><span class="op">)</span>, geneSetListLocation <span class="op">=</span> <span class="st">"rownames"</span><span class="op">)</span></code></pre></div>
150
+<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/runCellQC.html">runCellQC</a></span><span class="op">(</span><span class="va">sce</span>, sample <span class="op">=</span> <span class="cn">NULL</span>, algorithms <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html">c</a></span><span class="op">(</span><span class="st">"QCMetrics"</span>, <span class="st">"scDblFinder"</span>, <span class="st">"decontX"</span><span class="op">)</span>, geneSetList <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/list.html">list</a></span><span class="op">(</span>mito<span class="op">=</span><span class="va">mito.genes</span><span class="op">)</span>, geneSetListLocation <span class="op">=</span> <span class="st">"rownames"</span><span class="op">)</span></code></pre></div>
151 151
 <p><strong>Note:</strong> If you have cells from multiple samples stored in the SCE object, make sure to supply the <code>sample</code> parameter as the QC tools need to be applied to cells from each sample individually.</p>
152 152
 <p>Individual sets of QC metrics can be plotted with specific functions. For example to plot distributions of total numbers of UMIs derived from <code>runPerCellQC</code>, doublet scores from <code>runScDblFinder</code>, and contamination scores from <code>runDecontX</code> (all of which were run by the <code>runCellQC</code> function), the following plotting functions can be used:</p>
153 153
 <div class="sourceCode" id="cb6"><pre class="downlit sourceCode r">
154
-<code class="sourceCode R"><span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/plotRunPerCellQCResults.html">plotRunPerCellQCResults</a></span><span class="op">(</span><span class="va">sce</span><span class="op">)</span></code></pre></div>
154
+<code class="sourceCode R"><span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/plotRunPerCellQCResults.html">plotRunPerCellQCResults</a></span><span class="op">(</span><span class="va">sce</span><span class="op">)</span></code></pre></div>
155 155
 <p><img src="celda_pbmc3k_files/figure-html/sctk_qc_plot-1.png" width="700"></p>
156 156
 <div class="sourceCode" id="cb7"><pre class="downlit sourceCode r">
157
-<code class="sourceCode R"><span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/plotScDblFinderResults.html">plotScDblFinderResults</a></span><span class="op">(</span><span class="va">sce</span>, reducedDimName <span class="op">=</span> <span class="st">"decontX_UMAP"</span><span class="op">)</span></code></pre></div>
157
+<code class="sourceCode R"><span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/plotScDblFinderResults.html">plotScDblFinderResults</a></span><span class="op">(</span><span class="va">sce</span>, reducedDimName <span class="op">=</span> <span class="st">"decontX_UMAP"</span><span class="op">)</span></code></pre></div>
158 158
 <p><img src="celda_pbmc3k_files/figure-html/sctk_qc_plot-2.png" width="700"></p>
159 159
 <div class="sourceCode" id="cb8"><pre class="downlit sourceCode r">
160
-<code class="sourceCode R"><span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/plotDecontXResults.html">plotDecontXResults</a></span><span class="op">(</span><span class="va">sce</span>, reducedDimName <span class="op">=</span> <span class="st">"decontX_UMAP"</span><span class="op">)</span></code></pre></div>
160
+<code class="sourceCode R"><span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/plotDecontXResults.html">plotDecontXResults</a></span><span class="op">(</span><span class="va">sce</span>, reducedDimName <span class="op">=</span> <span class="st">"decontX_UMAP"</span><span class="op">)</span></code></pre></div>
161 161
 <p><img src="celda_pbmc3k_files/figure-html/sctk_qc_plot-3.png" width="700"></p>
162 162
 <p>An comprehensive HTML report can be generated to visualize and explore the QC metrics in greater detail:</p>
163 163
 <div class="sourceCode" id="cb9"><pre class="downlit sourceCode r">
164
-<code class="sourceCode R"><span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/reportCellQC.html">reportCellQC</a></span><span class="op">(</span><span class="va">sce</span><span class="op">)</span></code></pre></div>
164
+<code class="sourceCode R"><span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/reportCellQC.html">reportCellQC</a></span><span class="op">(</span><span class="va">sce</span><span class="op">)</span></code></pre></div>
165 165
 <p>After examining the distributions of various QC metrics, poor quality cells will need to be removed. Typically, thresholds for QC metrics should exclude cells that are outliers of the distribution (i.e. long tails in the violin or density plots). Cells can be removed using the <code>subsetSCECols</code> function. Metrics stored in the <code>colData</code> of the SCE object can be filtered using the <code>colData</code> parameter. Here we will limit to cells with at least 600 counts and 300 genes detected:</p>
166 166
 <div class="sourceCode" id="cb10"><pre class="downlit sourceCode r">
167 167
 <code class="sourceCode R"><span class="co"># Filter SCE</span>
168
-<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/subsetSCECols.html">subsetSCECols</a></span><span class="op">(</span><span class="va">sce</span>, colData <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html">c</a></span><span class="op">(</span><span class="st">"total &gt; 600"</span>, <span class="st">"detected &gt; 300"</span><span class="op">)</span><span class="op">)</span>
168
+<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/subsetSCECols.html">subsetSCECols</a></span><span class="op">(</span><span class="va">sce</span>, colData <span class="op">=</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html">c</a></span><span class="op">(</span><span class="st">"total &gt; 600"</span>, <span class="st">"detected &gt; 300"</span><span class="op">)</span><span class="op">)</span>
169 169
 
170 170
 <span class="co"># See number of cells after filtering</span>
171 171
 <span class="fu"><a href="https://rdrr.io/r/base/nrow.html">ncol</a></span><span class="op">(</span><span class="va">sce</span><span class="op">)</span></code></pre></div>
... ...
@@ -190,9 +190,9 @@
190 190
 <p><strong>Note:</strong> The subsetted matrix is stored in the “alternative experiment” slot (<code>altExp</code>) within the SCE. This allows for a matrix with a different number of rows to be stored within the same SCE object (rather than creating two SCE objects). The celda functions described in the next several sections operate on a matrix stored in the <code>altExp</code> slot. The default name given to the alternative experiment and used in all downstream celda functions is <code>"featureSubset"</code>. If the <code>altExpName</code> parameter is changed here, then it will need to be supplied to downstream plotting functions as well. The list of alternative experiments in an SCE can be view with <code><a href="https://rdrr.io/pkg/SingleCellExperiment/man/altExps.html">altExpNames(sce)</a></code>. If you have already have an SCE with selected features or do not want to perform feature selection, then you need to set the alternative experiment directly with a command like <code>altExp(sce, "featureSubset") &lt;- assay(sce, "counts")</code>. In the future, this will be updated to be more simple by utilizing the <code>ExperimentSubset</code> package.</p>
191 191
 <p>If the number of features is still relatively large (e.g. &gt;5000), an alternative approach is to select highly variable features that can be used in the downstream clustering. The advantage of this approach is that it can greatly speed up celda and can improve with module detection among highly variable features with overall lower expression. The disadvantage of this approach is that features that do not fall into the highly variable group will not be clustered into modules. The celda package does not include methods for selection of highly variable genes (HVGs). However, the <code>singleCellTK</code> provides wrappers for methods used in <a href="https://doi.org/10.1016/j.cell.2019.05.031" target="_blank">Seurat</a> and <a href="https://bioconductor.org/packages/release/bioc/html/scran.html">Scran</a>. We recommend keeping at least 2,000-5,000 HVGs for clustering. Here is some example code of how to select the top 5,000 most variable genes and store it back in the SCE as an <code>altExp</code>:</p>
192 192
 <div class="sourceCode" id="cb14"><pre class="downlit sourceCode r">
193
-<code class="sourceCode R"><span class="kw"><a href="https://rdrr.io/r/base/library.html">library</a></span><span class="op">(</span><span class="va"><a href="https://compbiomed.github.io/sctk_docs/">singleCellTK</a></span><span class="op">)</span>
194
-<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/seuratFindHVG.html">seuratFindHVG</a></span><span class="op">(</span><span class="va">sce</span>, useAssay <span class="op">=</span> <span class="va">useAssay</span>, hvgMethod <span class="op">=</span> <span class="st">"vst"</span><span class="op">)</span>
195
-<span class="va">g</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/getTopHVG.html">getTopHVG</a></span><span class="op">(</span><span class="va">sce</span>, method <span class="op">=</span> <span class="st">"vst"</span>, n <span class="op">=</span> <span class="fl">5000</span><span class="op">)</span>
193
+<code class="sourceCode R"><span class="kw"><a href="https://rdrr.io/r/base/library.html">library</a></span><span class="op">(</span><span class="va"><a href="https://www.camplab.net/sctk/">singleCellTK</a></span><span class="op">)</span>
194
+<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/seuratFindHVG.html">seuratFindHVG</a></span><span class="op">(</span><span class="va">sce</span>, useAssay <span class="op">=</span> <span class="va">useAssay</span>, hvgMethod <span class="op">=</span> <span class="st">"vst"</span><span class="op">)</span>
195
+<span class="va">g</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/getTopHVG.html">getTopHVG</a></span><span class="op">(</span><span class="va">sce</span>, method <span class="op">=</span> <span class="st">"vst"</span>, n <span class="op">=</span> <span class="fl">5000</span><span class="op">)</span>
196 196
 <span class="fu"><a href="https://rdrr.io/pkg/SingleCellExperiment/man/altExps.html">altExp</a></span><span class="op">(</span><span class="va">sce</span>, <span class="va">altExpName</span><span class="op">)</span> <span class="op">&lt;-</span> <span class="va">sce</span><span class="op">[</span><span class="va">g</span>, <span class="op">]</span></code></pre></div>
197 197
 <p>For the rest of the analysis with the PBMC3K data, we will use the first approach where features with at least 3 counts in 3 cells were included.</p>
198 198
 </div>
... ...
@@ -203,9 +203,9 @@
203 203
 <h2 class="hasAnchor">
204 204
 <a href="#bi-clustering-with-known-numbers-of-clusters" class="anchor"></a>Bi-clustering with known numbers of clusters</h2>
205 205
 <p>As mentioned earlier, celda is discrete Bayesian model that is able to simultaneously bi-cluster features into modules and cells into cell clusters. The primary bi-clustering model can be accessed with the function <code>celda_CG</code>. This function operates on a matrix stored as an alternative experiment in the <code>altExp</code> slot. If you did not perform feature selection as recommended in the previous section and your matrix of interest is not currently located in an <code>altExp</code> slot, the following code can be used to copy a matrix in the main assay slot to the <code>altExp</code> slot:</p>
206
-<div class="sourceCode" id="cb15"><pre class="sourceCode r"><code class="sourceCode r"><span id="cb15-1"><a href="#cb15-1" aria-hidden="true" tabindex="-1"></a>useAssay <span class="ot">&lt;-</span> <span class="st">"counts"</span></span>
207
-<span id="cb15-2"><a href="#cb15-2" aria-hidden="true" tabindex="-1"></a>altExpName <span class="ot">&lt;-</span> <span class="st">"featureSubset"</span></span>
208
-<span id="cb15-3"><a href="#cb15-3" aria-hidden="true" tabindex="-1"></a><span class="fu">altExp</span>(sce, altExpName) <span class="ot">&lt;-</span> <span class="fu">assay</span>(sce, useAssay)<span class="st">`</span><span class="at">. </span></span></code></pre></div>
206
+<div class="sourceCode" id="cb15"><pre class="sourceCode r"><code class="sourceCode r"><span id="cb15-1"><a href="#cb15-1"></a>useAssay &lt;-<span class="st"> "counts"</span></span>
207
+<span id="cb15-2"><a href="#cb15-2"></a>altExpName &lt;-<span class="st"> "featureSubset"</span></span>
208
+<span id="cb15-3"><a href="#cb15-3"></a><span class="kw">altExp</span>(sce, altExpName) &lt;-<span class="st"> </span><span class="kw">assay</span>(sce, useAssay)<span class="st">`</span><span class="dt">. </span></span></code></pre></div>
209 209
 <p>The two major adjustable parameters in this model are <code>L</code>, the number of modules, and <code>K</code>, the number of cell populations. The following code bi-clusters the PBMC3K dataset into 100 modules and 15 cell populations:</p>
210 210
 <div class="sourceCode" id="cb16"><pre class="downlit sourceCode r">
211 211
 <code class="sourceCode R"><span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="../../reference/celda_CG.html">celda_CG</a></span><span class="op">(</span><span class="va">sce</span>, L <span class="op">=</span> <span class="fl">100</span>, K <span class="op">=</span> <span class="fl">15</span>, useAssay <span class="op">=</span> <span class="va">useAssay</span>, altExpName <span class="op">=</span> <span class="va">altExpName</span><span class="op">)</span></code></pre></div>
... ...
@@ -28243,12 +28243,12 @@ RCE1
28243 28243
 <span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/pkg/scuttle/man/logNormCounts.html">logNormCounts</a></span><span class="op">(</span><span class="va">sce</span>, exprs_values <span class="op">=</span> <span class="va">useAssay</span>, name <span class="op">=</span> <span class="st">"logcounts"</span><span class="op">)</span>
28244 28244
 
28245 28245
 <span class="co"># Run differential expression analysis</span>
28246
-<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/findMarkerDiffExp.html">findMarkerDiffExp</a></span><span class="op">(</span><span class="va">sce</span>, useAssay <span class="op">=</span> <span class="st">"logcounts"</span>, method <span class="op">=</span> <span class="st">"wilcox"</span>, cluster <span class="op">=</span> <span class="fu"><a href="../../reference/celdaClusters.html">celdaClusters</a></span><span class="op">(</span><span class="va">sce</span><span class="op">)</span>, minMeanExpr <span class="op">=</span> <span class="fl">0</span>, fdrThreshold <span class="op">=</span> <span class="fl">0.05</span>, log2fcThreshold <span class="op">=</span> <span class="fl">0</span>, minClustExprPerc <span class="op">=</span> <span class="fl">0</span>, maxCtrlExprPerc <span class="op">=</span> <span class="fl">1</span><span class="op">)</span></code></pre></div>
28246
+<span class="va">sce</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/findMarkerDiffExp.html">findMarkerDiffExp</a></span><span class="op">(</span><span class="va">sce</span>, useAssay <span class="op">=</span> <span class="st">"logcounts"</span>, method <span class="op">=</span> <span class="st">"wilcox"</span>, cluster <span class="op">=</span> <span class="fu"><a href="../../reference/celdaClusters.html">celdaClusters</a></span><span class="op">(</span><span class="va">sce</span><span class="op">)</span>, minMeanExpr <span class="op">=</span> <span class="fl">0</span>, fdrThreshold <span class="op">=</span> <span class="fl">0.05</span>, log2fcThreshold <span class="op">=</span> <span class="fl">0</span>, minClustExprPerc <span class="op">=</span> <span class="fl">0</span>, maxCtrlExprPerc <span class="op">=</span> <span class="fl">1</span><span class="op">)</span></code></pre></div>
28247 28247
 <p>The function <code>plotMarkerDiffExp</code> can be used to plot the results in a heatmap. The <code>topN</code> parameter will plot the top N ranked genes for each cluster.</p>
28248 28248
 <div class="sourceCode" id="cb47"><pre class="downlit sourceCode r">
28249 28249
 <code class="sourceCode R"><span class="co"># Plot differentially expressed genes that pass additional thresholds 'minClustExprPerc' and 'maxCtrlExprPerc'</span>
28250
-<span class="fu"><a href="https://compbiomed.github.io/sctk_docs/reference/plotMarkerDiffExp.html">plotMarkerDiffExp</a></span><span class="op">(</span><span class="va">sce</span>, topN <span class="op">=</span> <span class="fl">5</span>, log2fcThreshold <span class="op">=</span> <span class="fl">0</span>, rowLabel <span class="op">=</span> <span class="cn">TRUE</span>, fdrThreshold <span class="op">=</span> <span class="fl">0.05</span>, minClustExprPerc <span class="op">=</span> <span class="fl">0.6</span>, maxCtrlExprPerc <span class="op">=</span> <span class="fl">0.4</span>, minMeanExpr <span class="op">=</span> <span class="fl">0</span><span class="op">)</span></code></pre></div>
28251
-<p><img src="celda_pbmc3k_files/figure-html/diffExp_plot-1.png" width="700"></p>
28250
+<span class="fu"><a href="https://rdrr.io/pkg/singleCellTK/man/plotMarkerDiffExp.html">plotMarkerDiffExp</a></span><span class="op">(</span><span class="va">sce</span>, topN <span class="op">=</span> <span class="fl">5</span>, log2fcThreshold <span class="op">=</span> <span class="fl">0</span>, rowLabel <span class="op">=</span> <span class="cn">TRUE</span>, fdrThreshold <span class="op">=</span> <span class="fl">0.05</span>, minClustExprPerc <span class="op">=</span> <span class="fl">0.6</span>, maxCtrlExprPerc <span class="op">=</span> <span class="fl">0.4</span>, minMeanExpr <span class="op">=</span> <span class="fl">0</span><span class="op">)</span></code></pre></div>
28251
+<p><img src="celda_pbmc3k_files/figure-html/diffExp_plot-1.png" width="700"><img src="celda_pbmc3k_files/figure-html/diffExp_plot-2.png" width="700"></p>
28252 28252
 <p>Other parameters such as <code>minClustExprPerc</code> (the minimum number of cells expressing the marker gene in the cluster) and <code>maxCtrlExprPerc</code> (the maximum number of cells expression the marker gene in other clusters) can be used to control how specific each marker gene is to each cell populations. Similarly, adding a log2 fold-change cutoff (e.g. 1) can select for markers that are more strongly up-regulated in a cell population.</p>
28253 28253
 </div>
28254 28254
 <div id="violin-plots-for-marker-genes" class="section level2">
... ...
@@ -28321,9 +28321,8 @@ RCE1
28321 28321
 <div id="session-information" class="section level1">
28322 28322
 <h1 class="hasAnchor">
28323 28323
 <a href="#session-information" class="anchor"></a>Session information</h1>
28324
-<details><summary>
28325
-sessionInfo()
28326
-</summary><pre><code>## R version 4.0.4 (2021-02-15)
28324
+<details><p><summary>sessionInfo()</summary></p>
28325
+<pre><code>## R version 4.0.4 (2021-02-15)
28327 28326
 ## Platform: x86_64-apple-darwin17.0 (64-bit)
28328 28327
 ## Running under: macOS Big Sur 10.16
28329 28328
 ## 
... ...
@@ -28341,7 +28340,7 @@ sessionInfo()
28341 28340
 ## other attached packages:
28342 28341
 ##  [1] scater_1.18.6               kableExtra_1.3.4           
28343 28342
 ##  [3] knitr_1.31                  ggplot2_3.3.5              
28344
-##  [5] celda_1.9.3                 singleCellTK_2.2.0         
28343
+##  [5] celda_1.10.0                singleCellTK_2.4.1         
28345 28344
 ##  [7] TENxPBMCData_1.8.0          HDF5Array_1.18.1           
28346 28345
 ##  [9] rhdf5_2.34.0                DelayedArray_0.16.2        
28347 28346
 ## [11] Matrix_1.3-2                SingleCellExperiment_1.12.0
... ...
@@ -28388,7 +28387,7 @@ sessionInfo()
28388 28387
 ##  [67] zlibbioc_1.36.0               purrr_0.3.4                  
28389 28388
 ##  [69] RCurl_1.98-1.2                dbscan_1.1-6                 
28390 28389
 ##  [71] GetoptLong_1.0.5              viridis_0.5.1                
28391
-##  [73] cowplot_1.1.1                 cluster_2.1.1                
28390
+##  [73] cowplot_1.1.1                 cluster_2.1.0                
28392 28391
 ##  [75] ggrepel_0.9.1                 fs_1.5.0                     
28393 28392
 ##  [77] magrittr_2.0.1                data.table_1.14.0            
28394 28393
 ##  [79] RSpectra_0.16-0               magick_2.7.0                 
28395 28394
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/cell_clusters-1.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/cell_clusters-1.png differ
28396 28395
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/cell_identities-1.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/cell_identities-1.png differ
28397 28396
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/cell_split_perplexity-1.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/cell_split_perplexity-1.png differ
28398 28397
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/cell_split_perplexity-2.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/cell_split_perplexity-2.png differ
28399 28398
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/diffExp_plot-1.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/diffExp_plot-1.png differ
28400 28399
new file mode 100644
28401 28400
Binary files /dev/null and b/docs/articles/articles/celda_pbmc3k_files/figure-html/diffExp_plot-2.png differ
28402 28401
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/module_split_perplexity-1.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/module_split_perplexity-1.png differ
28403 28402
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/module_split_rpc-1.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/module_split_rpc-1.png differ
28404 28403
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/sctk_qc_plot-1.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/sctk_qc_plot-1.png differ
28405 28404
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/sctk_qc_plot-2.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/sctk_qc_plot-2.png differ
28406 28405
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/sctk_qc_plot-3.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/sctk_qc_plot-3.png differ
28407 28406
Binary files a/docs/articles/articles/celda_pbmc3k_files/figure-html/violin-1.png and b/docs/articles/articles/celda_pbmc3k_files/figure-html/violin-1.png differ
... ...
@@ -31,7 +31,7 @@
31 31
       </button>
32 32
       <span class="navbar-brand">
33 33
         <a class="navbar-link" href="../../index.html">celda</a>
34
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
34
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
35 35
       </span>
36 36
     </div>
37 37
 
... ...
@@ -91,13 +91,13 @@
91 91
 
92 92
       
93 93
 
94
-      </header><script src="decontX_pbmc4k_files/header-attrs-2.7/header-attrs.js"></script><div class="row">
94
+      </header><script src="decontX_pbmc4k_files/header-attrs-2.7/header-attrs.js"></script><script src="decontX_pbmc4k_files/accessible-code-block-0.0.1/empty-anchor.js"></script><div class="row">
95 95
   <div class="col-md-9 contents">
96 96
     <div class="page-header toc-ignore">
97 97
       <h1 data-toc-skip>DecontX - Decontamination of PBMC4K</h1>
98 98
                         <h4 class="author">Joshua Campbell, Shiyi Yang, Zhe Wang, Yuan Yin</h4>
99 99
             
100
-            <h4 class="date">Compiled October 04, 2021</h4>
100
+            <h4 class="date">Compiled December 28, 2021</h4>
101 101
       
102 102
       <small class="dont-index">Source: <a href="https://github.com/campbio/celda/blob/master/vignettes/articles/decontX_pbmc4k.Rmd"><code>vignettes/articles/decontX_pbmc4k.Rmd</code></a></small>
103 103
       <div class="hidden name"><code>decontX_pbmc4k.Rmd</code></div>
... ...
@@ -305,7 +305,7 @@
305 305
 ## [11] GenomicRanges_1.42.0        GenomeInfoDb_1.26.4        
306 306
 ## [13] IRanges_2.24.1              S4Vectors_0.28.1           
307 307
 ## [15] BiocGenerics_0.36.0         MatrixGenerics_1.2.1       
308
-## [17] matrixStats_0.58.0          celda_1.9.3                
308
+## [17] matrixStats_0.58.0          celda_1.10.0               
309 309
 ## 
310 310
 ## loaded via a namespace (and not attached):
311 311
 ##   [1] AnnotationHub_2.22.0          BiocFileCache_1.14.0         
... ...
@@ -31,7 +31,7 @@
31 31
       </button>
32 32
       <span class="navbar-brand">
33 33
         <a class="navbar-link" href="../../index.html">celda</a>
34
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
34
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
35 35
       </span>
36 36
     </div>
37 37
 
... ...
@@ -91,7 +91,7 @@
91 91
 
92 92
       
93 93
 
94
-      </header><script src="installation_files/header-attrs-2.7/header-attrs.js"></script><div class="row">
94
+      </header><script src="installation_files/header-attrs-2.7/header-attrs.js"></script><script src="installation_files/accessible-code-block-0.0.1/empty-anchor.js"></script><div class="row">
95 95
   <div class="col-md-9 contents">
96 96
     <div class="page-header toc-ignore">
97 97
       <h1 data-toc-skip></h1>
... ...
@@ -120,12 +120,12 @@
120 120
 <p>The latest stable version of <strong>celda</strong> can be installed from GitHub using <code>devtools</code>:</p>
121 121
 <div class="sourceCode" id="cb2"><pre class="downlit sourceCode r">
122 122
 <code class="sourceCode R"><span class="kw"><a href="https://rdrr.io/r/base/library.html">library</a></span><span class="op">(</span><span class="va"><a href="https://devtools.r-lib.org/">devtools</a></span><span class="op">)</span>
123
-<span class="fu"><a href="https://devtools.r-lib.org//reference/remote-reexports.html">install_github</a></span><span class="op">(</span><span class="st">"campbio/celda"</span><span class="op">)</span></code></pre></div>
123
+<span class="fu"><a href="https://devtools.r-lib.org/reference/remote-reexports.html">install_github</a></span><span class="op">(</span><span class="st">"campbio/celda"</span><span class="op">)</span></code></pre></div>
124 124
 <p>The development version of <strong>celda</strong> can also be installed from GitHub using <code>devtools</code>:</p>
125 125
 <div class="sourceCode" id="cb3"><pre class="downlit sourceCode r">
126 126
 <code class="sourceCode R"><span class="kw"><a href="https://rdrr.io/r/base/library.html">library</a></span><span class="op">(</span><span class="va"><a href="https://devtools.r-lib.org/">devtools</a></span><span class="op">)</span>
127
-<span class="fu"><a href="https://devtools.r-lib.org//reference/remote-reexports.html">install_github</a></span><span class="op">(</span><span class="st">"campbio/celda@devel"</span><span class="op">)</span></code></pre></div>
128
-<p><strong>NOTE</strong> For MAC OSX users, <code><a href="https://devtools.r-lib.org//reference/remote-reexports.html">devtools::install_github()</a></code> requires installation of <strong>libgit2.</strong> This can be installed via homebrew:</p>
127
+<span class="fu"><a href="https://devtools.r-lib.org/reference/remote-reexports.html">install_github</a></span><span class="op">(</span><span class="st">"campbio/celda@devel"</span><span class="op">)</span></code></pre></div>
128
+<p><strong>NOTE</strong> For MAC OSX users, <code><a href="https://devtools.r-lib.org/reference/remote-reexports.html">devtools::install_github()</a></code> requires installation of <strong>libgit2.</strong> This can be installed via homebrew:</p>
129 129
 <pre><code>brew install libgit2</code></pre>
130 130
 </div>
131 131
 <div id="troubleshooting" class="section level1">
... ...
@@ -31,7 +31,7 @@
31 31
       </button>
32 32
       <span class="navbar-brand">
33 33
         <a class="navbar-link" href="../index.html">celda</a>
34
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
34
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
35 35
       </span>
36 36
     </div>
37 37
 
... ...
@@ -91,7 +91,7 @@
91 91
 
92 92
       
93 93
 
94
-      </header><script src="celda_files/header-attrs-2.7/header-attrs.js"></script><div class="row">
94
+      </header><script src="celda_files/header-attrs-2.7/header-attrs.js"></script><script src="celda_files/accessible-code-block-0.0.1/empty-anchor.js"></script><div class="row">
95 95
   <div class="col-md-9 contents">
96 96
     <div class="page-header toc-ignore">
97 97
       <h1 data-toc-skip>Analysis of single-cell genomic data with celda</h1>
... ...
@@ -107,7 +107,7 @@
107 107
             <address class="author_afil">
108 108
       Boston University School of Medicine<br><h4 class="author">Yusuke Koga</h4>
109 109
             <address class="author_afil">
110
-      Boston University School of Medicine<br><h4 class="date">2021-10-04</h4>
110
+      Boston University School of Medicine<br><h4 class="date">2021-12-28</h4>
111 111
       
112 112
       <small class="dont-index">Source: <a href="https://github.com/campbio/celda/blob/master/vignettes/celda.Rmd"><code>vignettes/celda.Rmd</code></a></small>
113 113
       <div class="hidden name"><code>celda.Rmd</code></div>
... ...
@@ -394,7 +394,7 @@
394 394
 ##  [5] GenomeInfoDb_1.26.4         IRanges_2.24.1             
395 395
 ##  [7] S4Vectors_0.28.1            BiocGenerics_0.36.0        
396 396
 ##  [9] MatrixGenerics_1.2.1        matrixStats_0.58.0         
397
-## [11] celda_1.9.3                 BiocStyle_2.18.1           
397
+## [11] celda_1.10.0                BiocStyle_2.18.1           
398 398
 ## 
399 399
 ## loaded via a namespace (and not attached):
400 400
 ##   [1] bitops_1.0-6               fs_1.5.0                  
... ...
@@ -446,7 +446,7 @@
446 446
 ##  [93] xfun_0.22                  assertive.types_0.0-3     
447 447
 ##  [95] RcppEigen_0.3.3.9.1        ragg_1.1.3                
448 448
 ##  [97] tibble_3.1.0               iterators_1.0.13          
449
-##  [99] memoise_2.0.0              cluster_2.1.1             
449
+##  [99] memoise_2.0.0              cluster_2.1.0             
450 450
 ## [101] ellipsis_0.3.1</code></pre>
451 451
 </div>
452 452
   </div>
453 453
Binary files a/docs/articles/celda_files/figure-html/module_split_rpc-1.png and b/docs/articles/celda_files/figure-html/module_split_rpc-1.png differ
454 454
Binary files a/docs/articles/celda_files/figure-html/plot_grid_search-1.png and b/docs/articles/celda_files/figure-html/plot_grid_search-1.png differ
455 455
Binary files a/docs/articles/celda_files/figure-html/rpc_cell-1.png and b/docs/articles/celda_files/figure-html/rpc_cell-1.png differ
456 456
Binary files a/docs/articles/celda_files/figure-html/rpc_cell-2.png and b/docs/articles/celda_files/figure-html/rpc_cell-2.png differ
457 457
Binary files a/docs/articles/celda_files/figure-html/unnamed-chunk-3-1.png and b/docs/articles/celda_files/figure-html/unnamed-chunk-3-1.png differ
... ...
@@ -31,7 +31,7 @@
31 31
       </button>
32 32
       <span class="navbar-brand">
33 33
         <a class="navbar-link" href="../index.html">celda</a>
34
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
34
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
35 35
       </span>
36 36
     </div>
37 37
 
... ...
@@ -91,7 +91,7 @@
91 91
 
92 92
       
93 93
 
94
-      </header><script src="decontX_files/header-attrs-2.7/header-attrs.js"></script><div class="row">
94
+      </header><script src="decontX_files/header-attrs-2.7/header-attrs.js"></script><script src="decontX_files/accessible-code-block-0.0.1/empty-anchor.js"></script><div class="row">
95 95
   <div class="col-md-9 contents">
96 96
     <div class="page-header toc-ignore">
97 97
       <h1 data-toc-skip>Decontamination of ambient RNA in single-cell genomic data with DecontX</h1>
... ...
@@ -106,7 +106,7 @@
106 106
       Boston University School of Medicine<br><a class="author_email" href="mailto:#"></a><a href="mailto:camp@bu.edu" class="email">camp@bu.edu</a>
107 107
       </address>
108 108
                   
109
-            <h4 class="date">2021-10-04</h4>
109
+            <h4 class="date">2021-12-28</h4>
110 110
       
111 111
       <small class="dont-index">Source: <a href="https://github.com/campbio/celda/blob/master/vignettes/decontX.Rmd"><code>vignettes/decontX.Rmd</code></a></small>
112 112
       <div class="hidden name"><code>decontX.Rmd</code></div>
... ...
@@ -317,7 +317,7 @@
317 317
 ## [11] GenomicRanges_1.42.0        GenomeInfoDb_1.26.4        
318 318
 ## [13] IRanges_2.24.1              S4Vectors_0.28.1           
319 319
 ## [15] BiocGenerics_0.36.0         MatrixGenerics_1.2.1       
320
-## [17] matrixStats_0.58.0          celda_1.9.3                
320
+## [17] matrixStats_0.58.0          celda_1.10.0               
321 321
 ## [19] BiocStyle_2.18.1           
322 322
 ## 
323 323
 ## loaded via a namespace (and not attached):
324 324
new file mode 100644
... ...
@@ -0,0 +1,15 @@
1
+// Hide empty <a> tag within highlighted CodeBlock for screen reader accessibility (see https://github.com/jgm/pandoc/issues/6352#issuecomment-626106786) -->
2
+// v0.0.1
3
+// Written by JooYoung Seo (jooyoung@psu.edu) and Atsushi Yasumoto on June 1st, 2020.
4
+
5
+document.addEventListener('DOMContentLoaded', function() {
6
+  const codeList = document.getElementsByClassName("sourceCode");
7
+  for (var i = 0; i < codeList.length; i++) {
8
+    var linkList = codeList[i].getElementsByTagName('a');
9
+    for (var j = 0; j < linkList.length; j++) {
10
+      if (linkList[j].innerHTML === "") {
11
+        linkList[j].setAttribute('aria-hidden', 'true');
12
+      }
13
+    }
14
+  }
15
+});
... ...
@@ -71,7 +71,7 @@
71 71
       </button>
72 72
       <span class="navbar-brand">
73 73
         <a class="navbar-link" href="../index.html">celda</a>
74
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
74
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
75 75
       </span>
76 76
     </div>
77 77
 
... ...
@@ -71,7 +71,7 @@
71 71
       </button>
72 72
       <span class="navbar-brand">
73 73
         <a class="navbar-link" href="index.html">celda</a>
74
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
74
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
75 75
       </span>
76 76
     </div>
77 77
 
... ...
@@ -37,7 +37,7 @@
37 37
       </button>
38 38
       <span class="navbar-brand">
39 39
         <a class="navbar-link" href="index.html">celda</a>
40
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
40
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
41 41
       </span>
42 42
     </div>
43 43
 
... ...
@@ -112,10 +112,13 @@
112 112
 <pre><code>if (!requireNamespace("BiocManager", quietly = TRUE))
113 113
     install.packages("BiocManager")
114 114
 BiocManager::install("celda")</code></pre>
115
-<p>To install the development version (R &gt;= 3.6) of <strong>celda</strong> from GitHub using <code>devtools</code>:</p>
115
+<p>The latest stable version of <strong>celda</strong> can be installed from GitHub using <code>devtools</code>:</p>
116 116
 <pre><code><a href="https://devtools.r-lib.org/">library(devtools)
117 117
 install_github("campbio/celda")</a></code></pre>
118
-<p><strong>NOTE</strong> For MAC OSX users, <code><a href="https://devtools.r-lib.org//reference/remote-reexports.html">devtools::install_github()</a></code> requires installation of <strong>libgit2.</strong> This can be installed via homebrew:</p>
118
+<p>The development version of <strong>celda</strong> can also be installed from GitHub using <code>devtools</code>:</p>
119
+<pre><code><a href="https://devtools.r-lib.org/">library(devtools)
120
+install_github("campbio/celda@devel")</a></code></pre>
121
+<p><strong>NOTE</strong> For MAC OSX users, <code><a href="https://devtools.r-lib.org/reference/remote-reexports.html">devtools::install_github()</a></code> requires installation of <strong>libgit2.</strong> This can be installed via homebrew:</p>
119 122
 <pre><code>brew install libgit2</code></pre>
120 123
 <p>Also, if you receive installation errors when Rcpp is being installed and compiled, try following the steps outlined here to solve the issue:</p>
121 124
 <p><a href="https://thecoatlessprofessor.com/programming/cpp/r-compiler-tools-for-rcpp-on-macos/" class="uri">https://thecoatlessprofessor.com/programming/cpp/r-compiler-tools-for-rcpp-on-macos/</a></p>
... ...
@@ -71,7 +71,7 @@
71 71
       </button>
72 72
       <span class="navbar-brand">
73 73
         <a class="navbar-link" href="../index.html">celda</a>
74
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
74
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
75 75
       </span>
76 76
     </div>
77 77
 
... ...
@@ -138,6 +138,13 @@
138 138
       <small>Source: <a href='https://github.com/campbio/celda/blob/master/NEWS.md'><code>NEWS.md</code></a></small>
139 139
     </div>
140 140
 
141
+    <div id="celda-v1100-2021-12-28" class="section level1">
142
+<h1 class="page-header" data-toc-text="1.10.0">
143
+<a href="#celda-v1100-2021-12-28" class="anchor"></a>celda v1.10.0 (2021-12-28)</h1>
144
+<ul>
145
+<li>Update to match Bioconductor release version</li>
146
+</ul>
147
+</div>
141 148
     <div id="celda-v193-2021-10-04" class="section level1">
142 149
 <h1 class="page-header" data-toc-text="1.9.3">
143 150
 <a href="#celda-v193-2021-10-04" class="anchor"></a>celda v1.9.3 (2021-10-04)</h1>
... ...
@@ -1,4 +1,4 @@
1
-pandoc: 2.11.4
1
+pandoc: 2.9.2.1
2 2
 pkgdown: 1.6.1
3 3
 pkgdown_sha: ~
4 4
 articles:
... ...
@@ -7,5 +7,5 @@ articles:
7 7
   articles/installation: installation.html
8 8
   celda: celda.html
9 9
   decontX: decontX.html
10
-last_built: 2021-10-04T16:07Z
10
+last_built: 2021-12-28T19:10Z
11 11
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -77,7 +77,7 @@
77 77
       </button>
78 78
       <span class="navbar-brand">
79 79
         <a class="navbar-link" href="../index.html">celda</a>
80
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
80
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
81 81
       </span>
82 82
     </div>
83 83
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -79,7 +79,7 @@
79 79
       </button>
80 80
       <span class="navbar-brand">
81 81
         <a class="navbar-link" href="../index.html">celda</a>
82
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
82
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
83 83
       </span>
84 84
     </div>
85 85
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -75,7 +75,7 @@
75 75
       </button>
76 76
       <span class="navbar-brand">
77 77
         <a class="navbar-link" href="../index.html">celda</a>
78
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
78
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
79 79
       </span>
80 80
     </div>
81 81
 
... ...
@@ -75,7 +75,7 @@
75 75
       </button>
76 76
       <span class="navbar-brand">
77 77
         <a class="navbar-link" href="../index.html">celda</a>
78
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
78
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
79 79
       </span>
80 80
     </div>
81 81
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -74,7 +74,7 @@
74 74
       </button>
75 75
       <span class="navbar-brand">
76 76
         <a class="navbar-link" href="../index.html">celda</a>
77
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
77
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
78 78
       </span>
79 79
     </div>
80 80
 
... ...
@@ -79,7 +79,7 @@
79 79
       </button>
80 80
       <span class="navbar-brand">
81 81
         <a class="navbar-link" href="../index.html">celda</a>
82
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
82
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
83 83
       </span>
84 84
     </div>
85 85
 
... ...
@@ -79,7 +79,7 @@
79 79
       </button>
80 80
       <span class="navbar-brand">
81 81
         <a class="navbar-link" href="../index.html">celda</a>
82
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
82
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
83 83
       </span>
84 84
     </div>
85 85
 
... ...
@@ -78,7 +78,7 @@
78 78
       </button>
79 79
       <span class="navbar-brand">
80 80
         <a class="navbar-link" href="../index.html">celda</a>
81
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
81
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
82 82
       </span>
83 83
     </div>
84 84
 
... ...
@@ -357,7 +357,7 @@ is only used when `zInitialize = 'predfined'`. Default NULL.</p></td>
357 357
     K <span class='op'>=</span> <span class='va'>celdaCSim</span><span class='op'>$</span><span class='va'>K</span>,
358 358
     sampleLabel <span class='op'>=</span> <span class='va'>celdaCSim</span><span class='op'>$</span><span class='va'>sampleLabel</span>,
359 359
     nchains <span class='op'>=</span> <span class='fl'>1</span><span class='op'>)</span>
360
-</div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Starting Celda_C: Clustering cells.</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:37 2021 .. Initializing 'z' in chain 1 with 'split' </span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:37 2021 .... Completed iteration: 1 | logLik: -1282027.27277705</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:37 2021 .... Completed iteration: 2 | logLik: -1282027.27277705</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:37 2021 .. Finished chain 1</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Completed Celda_C. Total time: 0.137078 secs</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div></pre>
360
+</div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Starting Celda_C: Clustering cells.</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:10:56 2021 .. Initializing 'z' in chain 1 with 'split' </span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:10:57 2021 .... Completed iteration: 1 | logLik: -1282027.27277705</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:10:57 2021 .... Completed iteration: 2 | logLik: -1282027.27277705</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:10:57 2021 .. Finished chain 1</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Completed Celda_C. Total time: 0.123713 secs</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div></pre>
361 361
   </div>
362 362
   <div class="col-md-3 hidden-xs hidden-sm" id="pkgdown-sidebar">
363 363
     <nav id="toc" data-toggle="toc" class="sticky-top">
... ...
@@ -78,7 +78,7 @@
78 78
       </button>
79 79
       <span class="navbar-brand">
80 80
         <a class="navbar-link" href="../index.html">celda</a>
81
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
81
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
82 82
       </span>
83 83
     </div>
84 84
 
... ...
@@ -402,7 +402,7 @@ is only used when `zInitialize = 'predfined'`. Default NULL.</p></td>
402 402
     L <span class='op'>=</span> <span class='va'>celdaCGSim</span><span class='op'>$</span><span class='va'>L</span>,
403 403
     sampleLabel <span class='op'>=</span> <span class='va'>celdaCGSim</span><span class='op'>$</span><span class='va'>sampleLabel</span>,
404 404
     nchains <span class='op'>=</span> <span class='fl'>1</span><span class='op'>)</span>
405
-</div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Starting Celda_CG: Clustering cells and genes.</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:38 2021 .. Initializing 'z' in chain 1 with 'split' </span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:38 2021 .. Initializing 'y' in chain 1 with 'split' </span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 1 | logLik: -1215542.98684529</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 2 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 3 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 4 | logLik: -1215542.98684529</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 5 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 6 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 7 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 8 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 9 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Determining if any gene clusters should be split.</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... No additional splitting was performed.</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Determining if any cell clusters should be split.</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... No additional splitting was performed.</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 10 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Determining if any cell clusters should be split.</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... No additional splitting was performed.</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .... Completed iteration: 11 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .. Finished chain 1</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Completed Celda_CG. Total time: 3.250677 secs</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div></pre>
405
+</div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Starting Celda_CG: Clustering cells and genes.</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:10:57 2021 .. Initializing 'z' in chain 1 with 'split' </span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:10:57 2021 .. Initializing 'y' in chain 1 with 'split' </span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 1 | logLik: -1215542.98684529</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 2 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 3 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 4 | logLik: -1215542.98684529</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 5 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 6 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 7 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 8 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 9 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Determining if any gene clusters should be split.</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... No additional splitting was performed.</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Determining if any cell clusters should be split.</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... No additional splitting was performed.</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 10 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Determining if any cell clusters should be split.</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... No additional splitting was performed.</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .... Completed iteration: 11 | logLik: -1215541.0958389</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .. Finished chain 1</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Completed Celda_CG. Total time: 3.006528 secs</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div></pre>
406 406
   </div>
407 407
   <div class="col-md-3 hidden-xs hidden-sm" id="pkgdown-sidebar">
408 408
     <nav id="toc" data-toggle="toc" class="sticky-top">
... ...
@@ -78,7 +78,7 @@
78 78
       </button>
79 79
       <span class="navbar-brand">
80 80
         <a class="navbar-link" href="../index.html">celda</a>
81
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
81
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
82 82
       </span>
83 83
     </div>
84 84
 
... ...
@@ -341,7 +341,7 @@ Default NULL.</p></td>
341 341
     <h2 class="hasAnchor" id="examples"><a class="anchor" href="#examples"></a>Examples</h2>
342 342
     <pre class="examples"><div class='input'><span class='fu'><a href='https://rdrr.io/r/utils/data.html'>data</a></span><span class='op'>(</span><span class='va'>celdaGSim</span><span class='op'>)</span>
343 343
 <span class='va'>sce</span> <span class='op'>&lt;-</span> <span class='fu'>celda_G</span><span class='op'>(</span><span class='va'>celdaGSim</span><span class='op'>$</span><span class='va'>counts</span>, L <span class='op'>=</span> <span class='va'>celdaGSim</span><span class='op'>$</span><span class='va'>L</span>, nchains <span class='op'>=</span> <span class='fl'>1</span><span class='op'>)</span>
344
-</div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Starting Celda_G: Clustering genes.</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:41 2021 .. Initializing 'y' in chain 1 with 'split' </span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:43 2021 .... Completed iteration: 1 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:43 2021 .... Completed iteration: 2 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:43 2021 .... Completed iteration: 3 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:43 2021 .... Completed iteration: 4 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:43 2021 .... Completed iteration: 5 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:43 2021 .... Completed iteration: 6 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:43 2021 .... Completed iteration: 7 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:44 2021 .... Completed iteration: 8 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:44 2021 .... Completed iteration: 9 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:44 2021 .... Determining if any gene clusters should be split.</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:44 2021 .... No additional splitting was performed.</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:44 2021 .... Completed iteration: 10 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:44 2021 .... Completed iteration: 11 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:44 2021 .. Finished chain 1</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Completed Celda_G. Total time: 2.44341 secs</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div></pre>
344
+</div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Starting Celda_G: Clustering genes.</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:00 2021 .. Initializing 'y' in chain 1 with 'split' </span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:02 2021 .... Completed iteration: 1 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:02 2021 .... Completed iteration: 2 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:02 2021 .... Completed iteration: 3 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:02 2021 .... Completed iteration: 4 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:02 2021 .... Completed iteration: 5 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:02 2021 .... Completed iteration: 6 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:02 2021 .... Completed iteration: 7 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:02 2021 .... Completed iteration: 8 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:02 2021 .... Completed iteration: 9 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:02 2021 .... Determining if any gene clusters should be split.</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:03 2021 .... No additional splitting was performed.</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:03 2021 .... Completed iteration: 10 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:03 2021 .... Completed iteration: 11 | logLik: -290669.046132139</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:03 2021 .. Finished chain 1</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Completed Celda_G. Total time: 2.186054 secs</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div></pre>
345 345
   </div>
346 346
   <div class="col-md-3 hidden-xs hidden-sm" id="pkgdown-sidebar">
347 347
     <nav id="toc" data-toggle="toc" class="sticky-top">
... ...
@@ -76,7 +76,7 @@
76 76
       </button>
77 77
       <span class="navbar-brand">
78 78
         <a class="navbar-link" href="../index.html">celda</a>
79
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
79
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
80 80
       </span>
81 81
     </div>
82 82
 
... ...
@@ -77,7 +77,7 @@
77 77
       </button>
78 78
       <span class="navbar-brand">
79 79
         <a class="navbar-link" href="../index.html">celda</a>
80
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
80
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
81 81
       </span>
82 82
     </div>
83 83
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -73,7 +73,7 @@ in single cell genomic datasets." />
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -387,7 +387,7 @@ cell cluster labels will be stored in
387 387
 #&gt; <span class='message'></span>
388 388
 #&gt; <span class='message'>    anyMissing, rowMedians</span></div><div class='input'><span class='va'>sce</span> <span class='op'>&lt;-</span> <span class='fu'><a href='https://rdrr.io/pkg/SingleCellExperiment/man/SingleCellExperiment.html'>SingleCellExperiment</a></span><span class='op'>(</span><span class='fu'><a href='https://rdrr.io/r/base/list.html'>list</a></span><span class='op'>(</span>counts <span class='op'>=</span> <span class='va'>s</span><span class='op'>$</span><span class='va'>observedCounts</span><span class='op'>)</span><span class='op'>)</span>
389 389
 <span class='va'>sce</span> <span class='op'>&lt;-</span> <span class='fu'>decontX</span><span class='op'>(</span><span class='va'>sce</span><span class='op'>)</span>
390
-</div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Starting DecontX</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:46 2021 .. Analyzing all cells</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:46 2021 .... Converting to sparse matrix</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:46 2021 .... Generating UMAP and estimating cell types</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:49 2021 .... Estimating contamination</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:49 2021 ...... Completed iteration: 9 | converge: 0.0009154</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:07:49 2021 .. Calculating final decontaminated matrix</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Completed DecontX. Total time: 2.786641 secs</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='input'>
390
+</div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Starting DecontX</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:05 2021 .. Analyzing all cells</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:05 2021 .... Converting to sparse matrix</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:05 2021 .... Generating UMAP and estimating cell types</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:08 2021 .... Estimating contamination</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:08 2021 ...... Completed iteration: 9 | converge: 0.0009154</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:11:08 2021 .. Calculating final decontaminated matrix</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='output co'>#&gt; <span class='message'>Completed DecontX. Total time: 2.749794 secs</span></div><div class='output co'>#&gt; <span class='message'>--------------------------------------------------</span></div><div class='input'>
391 391
 <span class='co'># Plot contamination on UMAP</span>
392 392
 <span class='fu'><a href='plotDecontXContamination.html'>plotDecontXContamination</a></span><span class='op'>(</span><span class='va'>sce</span><span class='op'>)</span>
393 393
 </div><div class='img'><img src='decontX-1.png' alt='' width='700' height='433' /></div><div class='input'>
... ...
@@ -73,7 +73,7 @@ a SingleCellExperiment object." />
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -75,7 +75,7 @@
75 75
       </button>
76 76
       <span class="navbar-brand">
77 77
         <a class="navbar-link" href="../index.html">celda</a>
78
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
78
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
79 79
       </span>
80 80
     </div>
81 81
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -74,7 +74,7 @@
74 74
       </button>
75 75
       <span class="navbar-brand">
76 76
         <a class="navbar-link" href="../index.html">celda</a>
77
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
77
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
78 78
       </span>
79 79
     </div>
80 80
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -75,7 +75,7 @@
75 75
       </button>
76 76
       <span class="navbar-brand">
77 77
         <a class="navbar-link" href="../index.html">celda</a>
78
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
78
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
79 79
       </span>
80 80
     </div>
81 81
 
... ...
@@ -71,7 +71,7 @@
71 71
       </button>
72 72
       <span class="navbar-brand">
73 73
         <a class="navbar-link" href="../index.html">celda</a>
74
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
74
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
75 75
       </span>
76 76
     </div>
77 77
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -78,7 +78,7 @@
78 78
       </button>
79 79
       <span class="navbar-brand">
80 80
         <a class="navbar-link" href="../index.html">celda</a>
81
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
81
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
82 82
       </span>
83 83
     </div>
84 84
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -74,7 +74,7 @@
74 74
       </button>
75 75
       <span class="navbar-brand">
76 76
         <a class="navbar-link" href="../index.html">celda</a>
77
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
77
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
78 78
       </span>
79 79
     </div>
80 80
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -72,7 +72,7 @@
72 72
       </button>
73 73
       <span class="navbar-brand">
74 74
         <a class="navbar-link" href="../index.html">celda</a>
75
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
75
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
76 76
       </span>
77 77
     </div>
78 78
 
... ...
@@ -73,7 +73,7 @@ cells colored by the estimated percentation of contamation." />
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -75,7 +75,7 @@ decontamination with decontX." />
75 75
       </button>
76 76
       <span class="navbar-brand">
77 77
         <a class="navbar-link" href="../index.html">celda</a>
78
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
78
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
79 79
       </span>
80 80
     </div>
81 81
 
... ...
@@ -76,7 +76,7 @@ decontamination with decontX." />
76 76
       </button>
77 77
       <span class="navbar-brand">
78 78
         <a class="navbar-link" href="../index.html">celda</a>
79
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
79
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
80 80
       </span>
81 81
     </div>
82 82
 
... ...
@@ -78,7 +78,7 @@
78 78
       </button>
79 79
       <span class="navbar-brand">
80 80
         <a class="navbar-link" href="../index.html">celda</a>
81
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
81
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
82 82
       </span>
83 83
     </div>
84 84
 
... ...
@@ -75,7 +75,7 @@
75 75
       </button>
76 76
       <span class="navbar-brand">
77 77
         <a class="navbar-link" href="../index.html">celda</a>
78
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
78
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
79 79
       </span>
80 80
     </div>
81 81
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -77,7 +77,7 @@
77 77
       </button>
78 78
       <span class="navbar-brand">
79 79
         <a class="navbar-link" href="../index.html">celda</a>
80
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
80
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
81 81
       </span>
82 82
     </div>
83 83
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -74,7 +74,7 @@
74 74
       </button>
75 75
       <span class="navbar-brand">
76 76
         <a class="navbar-link" href="../index.html">celda</a>
77
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
77
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
78 78
       </span>
79 79
     </div>
80 80
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -73,7 +73,7 @@
73 73
       </button>
74 74
       <span class="navbar-brand">
75 75
         <a class="navbar-link" href="../index.html">celda</a>
76
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
76
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
77 77
       </span>
78 78
     </div>
79 79
 
... ...
@@ -81,7 +81,7 @@
81 81
       </button>
82 82
       <span class="navbar-brand">
83 83
         <a class="navbar-link" href="../index.html">celda</a>
84
-        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.9.3</span>
84
+        <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">1.10.0</span>
85 85
       </span>
86 86
     </div>
87 87
 
... ...
@@ -355,36 +355,36 @@ If FALSE, then perplexity can be calculated later with
355 355
 <span class='co'>## Create models that range from K = 3 to K = 7 by recursively splitting</span>
356 356
 <span class='co'>## cell populations into two to produce \link{celda_C} cell clustering models</span>
357 357
 <span class='va'>sce</span> <span class='op'>&lt;-</span> <span class='fu'>recursiveSplitCell</span><span class='op'>(</span><span class='va'>sceCeldaCG</span>, initialK <span class='op'>=</span> <span class='fl'>3</span>, maxK <span class='op'>=</span> <span class='fl'>7</span><span class='op'>)</span>
358
-</div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Starting recursive cell population splitting.</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:56 2021 .. Initializing with 3 populations</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:56 2021 .. Current cell population 4 | logLik: -1225755.01101897</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:56 2021 .. Current cell population 5 | logLik: -1213677.60126784</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:56 2021 .. Current cell population 6 | logLik: -1213903.59449854</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:56 2021 .. Current cell population 7 | logLik: -1214081.54311397</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:56 2021 .. Calculating perplexity</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Completed recursive cell population splitting. Total time: 0.3744612 secs</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='input'>
358
+</div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Starting recursive cell population splitting.</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:14 2021 .. Initializing with 3 populations</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:14 2021 .. Current cell population 4 | logLik: -1225755.01101897</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:14 2021 .. Current cell population 5 | logLik: -1213677.60126784</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:14 2021 .. Current cell population 6 | logLik: -1213903.59449854</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:14 2021 .. Current cell population 7 | logLik: -1214081.54311397</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:14 2021 .. Calculating perplexity</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Completed recursive cell population splitting. Total time: 0.3816781 secs</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='input'>
359 359
 <span class='co'>## Alternatively, first identify features modules using</span>
360 360
 <span class='co'>## \link{recursiveSplitModule}</span>
361 361
 <span class='va'>moduleSplit</span> <span class='op'>&lt;-</span> <span class='fu'><a href='recursiveSplitModule.html'>recursiveSplitModule</a></span><span class='op'>(</span><span class='va'>sceCeldaCG</span>, initialL <span class='op'>=</span> <span class='fl'>3</span>, maxL <span class='op'>=</span> <span class='fl'>15</span><span class='op'>)</span>
362
-</div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Starting recursive module splitting.</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:56 2021 .. Collapsing to 100 temporary cell populations</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:57 2021 .. Initializing with 3 modules</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:57 2021 .. Created module 4 | logLik: -1241379.90928455</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:57 2021 .. Created module 5 | logLik: -1235212.7977535</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Created module 6 | logLik: -1232789.9817561</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Created module 7 | logLik: -1227246.66090571</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Created module 8 | logLik: -1223898.757694</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Created module 9 | logLik: -1221848.26936098</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Created module 10 | logLik: -1220147.96681948</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Created module 11 | logLik: -1220818.37022325</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Created module 12 | logLik: -1221489.07685946</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Created module 13 | logLik: -1222032.53497571</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Created module 14 | logLik: -1222712.17543857</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Created module 15 | logLik: -1223268.97596756</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:58 2021 .. Calculating perplexity</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Completed recursive module splitting. Total time: 2.250609 secs</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='input'><span class='fu'><a href='plotGridSearchPerplexity.html'>plotGridSearchPerplexity</a></span><span class='op'>(</span><span class='va'>moduleSplit</span><span class='op'>)</span>
362
+</div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Starting recursive module splitting.</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:15 2021 .. Collapsing to 100 temporary cell populations</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:15 2021 .. Initializing with 3 modules</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 4 | logLik: -1241379.90928455</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 5 | logLik: -1235212.7977535</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 6 | logLik: -1232789.9817561</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 7 | logLik: -1227246.66090571</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 8 | logLik: -1223898.757694</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 9 | logLik: -1221848.26936098</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 10 | logLik: -1220147.96681948</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 11 | logLik: -1220818.37022325</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 12 | logLik: -1221489.07685946</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 13 | logLik: -1222032.53497571</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 14 | logLik: -1222712.17543857</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Created module 15 | logLik: -1223268.97596756</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:16 2021 .. Calculating perplexity</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Completed recursive module splitting. Total time: 2.269934 secs</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='input'><span class='fu'><a href='plotGridSearchPerplexity.html'>plotGridSearchPerplexity</a></span><span class='op'>(</span><span class='va'>moduleSplit</span><span class='op'>)</span>
363 363
 </div><div class='img'><img src='recursiveSplitCell-1.png' alt='' width='700' height='433' /></div><div class='input'><span class='va'>moduleSplitSelect</span> <span class='op'>&lt;-</span> <span class='fu'><a href='subsetCeldaList.html'>subsetCeldaList</a></span><span class='op'>(</span><span class='va'>moduleSplit</span>, <span class='fu'><a href='https://rdrr.io/r/base/list.html'>list</a></span><span class='op'>(</span>L <span class='op'>=</span> <span class='fl'>10</span><span class='op'>)</span><span class='op'>)</span>
364 364
 
365 365
 <span class='co'>## Then use module labels for initialization in \link{recursiveSplitCell} to</span>
366 366
 <span class='co'>## produce \link{celda_CG} bi-clustering models</span>
367 367
 <span class='va'>cellSplit</span> <span class='op'>&lt;-</span> <span class='fu'>recursiveSplitCell</span><span class='op'>(</span><span class='va'>sceCeldaCG</span>,
368 368
   initialK <span class='op'>=</span> <span class='fl'>3</span>, maxK <span class='op'>=</span> <span class='fl'>7</span>, yInit <span class='op'>=</span> <span class='fu'><a href='celdaModules.html'>celdaModules</a></span><span class='op'>(</span><span class='va'>moduleSplitSelect</span><span class='op'>)</span><span class='op'>)</span>
369
-</div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Starting recursive cell population splitting.</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:59 2021 .. Collapsing to 10 modules</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:08:59 2021 .. Initializing with 3 populations</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:00 2021 .. Current cell population 4 | logLik: -1225286.49558716</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:00 2021 .. Current cell population 5 | logLik: -1212955.15575681</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:00 2021 .. Current cell population 6 | logLik: -1212982.74290613</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:00 2021 .. Current cell population 7 | logLik: -1213005.40337891</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:00 2021 .. Calculating perplexity</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Completed recursive cell population splitting. Total time: 1.471492 secs</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='input'><span class='fu'><a href='plotGridSearchPerplexity.html'>plotGridSearchPerplexity</a></span><span class='op'>(</span><span class='va'>cellSplit</span><span class='op'>)</span>
369
+</div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Starting recursive cell population splitting.</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:17 2021 .. Collapsing to 10 modules</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:17 2021 .. Initializing with 3 populations</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:18 2021 .. Current cell population 4 | logLik: -1225286.49558716</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:18 2021 .. Current cell population 5 | logLik: -1212955.15575681</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:18 2021 .. Current cell population 6 | logLik: -1212982.74290613</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:18 2021 .. Current cell population 7 | logLik: -1213005.40337891</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:18 2021 .. Calculating perplexity</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Completed recursive cell population splitting. Total time: 1.154469 secs</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='input'><span class='fu'><a href='plotGridSearchPerplexity.html'>plotGridSearchPerplexity</a></span><span class='op'>(</span><span class='va'>cellSplit</span><span class='op'>)</span>
370 370
 </div><div class='img'><img src='recursiveSplitCell-2.png' alt='' width='700' height='433' /></div><div class='input'><span class='va'>sce</span> <span class='op'>&lt;-</span> <span class='fu'><a href='subsetCeldaList.html'>subsetCeldaList</a></span><span class='op'>(</span><span class='va'>cellSplit</span>, <span class='fu'><a href='https://rdrr.io/r/base/list.html'>list</a></span><span class='op'>(</span>K <span class='op'>=</span> <span class='fl'>5</span>, L <span class='op'>=</span> <span class='fl'>10</span><span class='op'>)</span><span class='op'>)</span>
371 371
 <span class='fu'><a href='https://rdrr.io/r/utils/data.html'>data</a></span><span class='op'>(</span><span class='va'>celdaCGSim</span>, <span class='va'>celdaCSim</span><span class='op'>)</span>
372 372
 <span class='co'>## Create models that range from K = 3 to K = 7 by recursively splitting</span>
373 373
 <span class='co'>## cell populations into two to produce \link{celda_C} cell clustering models</span>
374 374
 <span class='va'>sce</span> <span class='op'>&lt;-</span> <span class='fu'>recursiveSplitCell</span><span class='op'>(</span><span class='va'>celdaCSim</span><span class='op'>$</span><span class='va'>counts</span>, initialK <span class='op'>=</span> <span class='fl'>3</span>, maxK <span class='op'>=</span> <span class='fl'>7</span><span class='op'>)</span>
375
-</div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Starting recursive cell population splitting.</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:01 2021 .. Initializing with 3 populations</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:01 2021 .. Current cell population 4 | logLik: -1341630.1679001</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:01 2021 .. Current cell population 5 | logLik: -1327506.91718317</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:01 2021 .. Current cell population 6 | logLik: -1315227.54586167</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:01 2021 .. Current cell population 7 | logLik: -1304393.65802293</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:01 2021 .. Calculating perplexity</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Completed recursive cell population splitting. Total time: 0.3327959 secs</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='input'>
375
+</div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Starting recursive cell population splitting.</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:19 2021 .. Initializing with 3 populations</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:19 2021 .. Current cell population 4 | logLik: -1341630.1679001</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:19 2021 .. Current cell population 5 | logLik: -1327506.91718317</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:19 2021 .. Current cell population 6 | logLik: -1315227.54586167</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:19 2021 .. Current cell population 7 | logLik: -1304393.65802293</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:19 2021 .. Calculating perplexity</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Completed recursive cell population splitting. Total time: 0.334146 secs</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='input'>
376 376
 <span class='co'>## Alternatively, first identify features modules using</span>
377 377
 <span class='co'>## \link{recursiveSplitModule}</span>
378 378
 <span class='va'>moduleSplit</span> <span class='op'>&lt;-</span> <span class='fu'><a href='recursiveSplitModule.html'>recursiveSplitModule</a></span><span class='op'>(</span><span class='va'>celdaCGSim</span><span class='op'>$</span><span class='va'>counts</span>,
379 379
   initialL <span class='op'>=</span> <span class='fl'>3</span>, maxL <span class='op'>=</span> <span class='fl'>15</span><span class='op'>)</span>
380
-</div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Starting recursive module splitting.</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:01 2021 .. Collapsing to 100 temporary cell populations</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:02 2021 .. Initializing with 3 modules</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:02 2021 .. Created module 4 | logLik: -1243396.62348886</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:02 2021 .. Created module 5 | logLik: -1237610.11790137</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:02 2021 .. Created module 6 | logLik: -1232128.87013396</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:02 2021 .. Created module 7 | logLik: -1227611.8250329</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:02 2021 .. Created module 8 | logLik: -1225618.06184004</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:02 2021 .. Created module 9 | logLik: -1223967.77531912</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:02 2021 .. Created module 10 | logLik: -1222801.11395987</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:03 2021 .. Created module 11 | logLik: -1223402.66903597</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:03 2021 .. Created module 12 | logLik: -1224026.19892208</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:03 2021 .. Created module 13 | logLik: -1224675.63005464</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:03 2021 .. Created module 14 | logLik: -1225317.91966369</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:03 2021 .. Created module 15 | logLik: -1225971.50555157</span></div><div class='output co'>#&gt; <span class='message'>Mon Oct  4 12:09:03 2021 .. Calculating perplexity</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Completed recursive module splitting. Total time: 2.125219 secs</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='input'><span class='fu'><a href='plotGridSearchPerplexity.html'>plotGridSearchPerplexity</a></span><span class='op'>(</span><span class='va'>moduleSplit</span><span class='op'>)</span>
380
+</div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Starting recursive module splitting.</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:19 2021 .. Collapsing to 100 temporary cell populations</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:20 2021 .. Initializing with 3 modules</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:20 2021 .. Created module 4 | logLik: -1243396.62348886</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:20 2021 .. Created module 5 | logLik: -1237610.11790137</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:20 2021 .. Created module 6 | logLik: -1232128.87013396</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:20 2021 .. Created module 7 | logLik: -1227611.8250329</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:20 2021 .. Created module 8 | logLik: -1225618.06184004</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:21 2021 .. Created module 9 | logLik: -1223967.77531912</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:21 2021 .. Created module 10 | logLik: -1222801.11395987</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:21 2021 .. Created module 11 | logLik: -1223402.66903597</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:21 2021 .. Created module 12 | logLik: -1224026.19892208</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:21 2021 .. Created module 13 | logLik: -1224675.63005464</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:21 2021 .. Created module 14 | logLik: -1225317.91966369</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:21 2021 .. Created module 15 | logLik: -1225971.50555157</span></div><div class='output co'>#&gt; <span class='message'>Tue Dec 28 14:12:21 2021 .. Calculating perplexity</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='output co'>#&gt; <span class='message'>Completed recursive module splitting. Total time: 2.180047 secs</span></div><div class='output co'>#&gt; <span class='message'>==================================================</span></div><div class='input'><span class='fu'><a href='plotGridSearchPerplexity.html'>plotGridSearchPerplexity</a></span><span class='op'>(</span><span class='va'>moduleSplit</span><span class='op'>)</span>
381 381
 </div><div class='img'><img src='recursiveSplitCell-3.png' alt='' width='700' height='433' /></div><div class='input'><span class='va'>moduleSplitSelect</span> <span class='op'>&lt;-</span> <span class='fu'><a href='subsetCeldaList.html'>subsetCeldaList</a></span><span class='op'>(</span><span class='va'>moduleSplit</span>, <span class='fu'><a href='https://rdrr.io/r/base/list.html'>list</a></span><span class='op'>(</span>L <span class='op'>=</span> <span class='fl'>10</span><span class='op'>)</span><span class='op'>)</span>
382 382
 
383 383
 <span class='co'>## Then use module labels for initialization in \link{recursiveSplitCell} to</span>
384 384
 <span class='co'>## produce \link{celda_CG} bi-clustering models</span>
385 385
 <span class='va'>cellSplit</span> <span class='op'>&lt;-</span> <span class='fu'>recursiveSplitCell</span><span class='op'>(</span><span class='va'>celdaCGSim</span><span class='op'>$</span><span class='va'>counts</span>,
386 386
   initialK <span class='op'>=</span> <span class='fl'>3</span>, maxK <span class='op'>=</span> <span class='fl'>7</span>, yInit <span class='op'>=</span> <span class='fu'><a href='celdaModules.html'>celdaModules</a></span><span class='op'>(</span><span class='va'>moduleSplitSelect</span><span class='op'>)</span><span class='op'>)</span>