From 4873e5597a3165c396250ef501b6742694333eb5 Mon Sep 17 00:00:00 2001
From: jamespeapen <39574560+jamespeapen@users.noreply.github.com>
Date: Mon, 2 Dec 2024 17:10:51 +0000
Subject: [PATCH] =?UTF-8?q?Deploying=20to=20gh-pages=20from=20@=20huishenl?=
 =?UTF-8?q?ab/iscream@cc91128e02f4824a3de376e040dbe4615f8e22e9=20?=
 =?UTF-8?q?=F0=9F=9A=80?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 dev/LICENSE-text.html                     |  2 +-
 dev/LICENSE.html                          |  2 +-
 dev/authors.html                          | 10 +++++++++-
 dev/index.html                            |  4 +++-
 dev/pkgdown.yml                           |  2 +-
 dev/reference/Cpp_query_all.html          |  2 +-
 dev/reference/Cpp_query_chroms.html       |  2 +-
 dev/reference/Cpp_query_interval.html     |  2 +-
 dev/reference/Cpp_set_log_level.html      |  2 +-
 dev/reference/Cpp_summarize_regions.html  |  2 +-
 dev/reference/check_files_exist.html      |  2 +-
 dev/reference/check_thread_count.html     |  2 +-
 dev/reference/dot-get_threads.html        |  2 +-
 dev/reference/get_granges_string.html     |  2 +-
 dev/reference/get_log_level.html          |  2 +-
 dev/reference/get_omp_threads.html        |  2 +-
 dev/reference/get_threads.html            |  2 +-
 dev/reference/index.html                  |  2 +-
 dev/reference/query_all.html              | 10 +++++-----
 dev/reference/query_chroms.html           |  2 +-
 dev/reference/scan_tabix.html             |  2 +-
 dev/reference/set_log_level.html          |  2 +-
 dev/reference/set_threads.html            |  2 +-
 dev/reference/summarize_regions.html      | 14 +++++++-------
 dev/reference/tabix.html                  |  2 +-
 dev/reference/validate_log_level.html     |  2 +-
 dev/reference/verify_aligner_or_stop.html |  2 +-
 dev/reference/verify_files_or_stop.html   |  2 +-
 dev/reference/verify_filetype.html        |  2 +-
 dev/search.json                           |  2 +-
 30 files changed, 50 insertions(+), 40 deletions(-)

diff --git a/dev/LICENSE-text.html b/dev/LICENSE-text.html
index d0d2b0f..a12f60c 100644
--- a/dev/LICENSE-text.html
+++ b/dev/LICENSE-text.html
@@ -38,7 +38,7 @@
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/LICENSE.html b/dev/LICENSE.html
index f2dca38..351d18c 100644
--- a/dev/LICENSE.html
+++ b/dev/LICENSE.html
@@ -42,7 +42,7 @@
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/authors.html b/dev/authors.html
index 56e92bb..b3274df 100644
--- a/dev/authors.html
+++ b/dev/authors.html
@@ -36,6 +36,14 @@ <h2>Authors</h2>
           <p><strong>James Eapen</strong>. Author, maintainer. <a href="https://orcid.org/0000-0001-6016-3598" target="orcid.widget" aria-label="ORCID" class="external-link"><span class="fab fa-orcid orcid" aria-hidden="true"></span></a>
           </p>
         </li>
+        <li>
+          <p><strong>Jacob Morrison</strong>. Contributor. <a href="https://orcid.org/0000-0001-8592-4744" target="orcid.widget" aria-label="ORCID" class="external-link"><span class="fab fa-orcid orcid" aria-hidden="true"></span></a>
+          </p>
+        </li>
+        <li>
+          <p><strong>Hui Shen</strong>. Funder. <a href="https://orcid.org/0000-0001-9767-4084" target="orcid.widget" aria-label="ORCID" class="external-link"><span class="fab fa-orcid orcid" aria-hidden="true"></span></a>
+          </p>
+        </li>
       </ul></div>
 
     <div class="section level2">
@@ -60,7 +68,7 @@ <h2 id="citation">Citation</h2>
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/index.html b/dev/index.html
index 6215a55..d8bee94 100644
--- a/dev/index.html
+++ b/dev/index.html
@@ -142,6 +142,8 @@ <h2 data-toc-skip>Citation</h2>
 <h2 data-toc-skip>Developers</h2>
 <ul class="list-unstyled">
 <li>James Eapen <br><small class="roles"> Author, maintainer </small> <a href="https://orcid.org/0000-0001-6016-3598" target="orcid.widget" aria-label="ORCID" class="external-link"><span class="fab fa-orcid orcid" aria-hidden="true"></span></a> </li>
+<li>Hui Shen <br><small class="roles"> Funder </small> <a href="https://orcid.org/0000-0001-9767-4084" target="orcid.widget" aria-label="ORCID" class="external-link"><span class="fab fa-orcid orcid" aria-hidden="true"></span></a> </li>
+<li><a href="authors.html">More about authors...</a></li>
 </ul>
 </div>
 
@@ -158,7 +160,7 @@ <h2 data-toc-skip>Dev status</h2>
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/pkgdown.yml b/dev/pkgdown.yml
index 14d5ad6..e18b2db 100644
--- a/dev/pkgdown.yml
+++ b/dev/pkgdown.yml
@@ -2,7 +2,7 @@ pandoc: 3.1.11
 pkgdown: 2.1.1
 pkgdown_sha: ~
 articles: {}
-last_built: 2024-12-02T16:55Z
+last_built: 2024-12-02T17:10Z
 urls:
   reference: https://huishenlab.github.io/iscream/reference
   article: https://huishenlab.github.io/iscream/articles
diff --git a/dev/reference/Cpp_query_all.html b/dev/reference/Cpp_query_all.html
index e00f999..4e15b41 100644
--- a/dev/reference/Cpp_query_all.html
+++ b/dev/reference/Cpp_query_all.html
@@ -71,7 +71,7 @@ <h2 id="arguments">Arguments<a class="anchor" aria-label="anchor" href="#argumen
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/Cpp_query_chroms.html b/dev/reference/Cpp_query_chroms.html
index eeb36ce..a188ede 100644
--- a/dev/reference/Cpp_query_chroms.html
+++ b/dev/reference/Cpp_query_chroms.html
@@ -58,7 +58,7 @@ <h2 id="value">Value<a class="anchor" aria-label="anchor" href="#value"></a></h2
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/Cpp_query_interval.html b/dev/reference/Cpp_query_interval.html
index e2555f7..b808f63 100644
--- a/dev/reference/Cpp_query_interval.html
+++ b/dev/reference/Cpp_query_interval.html
@@ -59,7 +59,7 @@ <h2 id="arguments">Arguments<a class="anchor" aria-label="anchor" href="#argumen
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/Cpp_set_log_level.html b/dev/reference/Cpp_set_log_level.html
index 3ee51a3..30c0de5 100644
--- a/dev/reference/Cpp_set_log_level.html
+++ b/dev/reference/Cpp_set_log_level.html
@@ -64,7 +64,7 @@ <h2 id="value">Value<a class="anchor" aria-label="anchor" href="#value"></a></h2
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/Cpp_summarize_regions.html b/dev/reference/Cpp_summarize_regions.html
index 5167d9c..18993d2 100644
--- a/dev/reference/Cpp_summarize_regions.html
+++ b/dev/reference/Cpp_summarize_regions.html
@@ -99,7 +99,7 @@ <h2 id="details">Details<a class="anchor" aria-label="anchor" href="#details"></
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/check_files_exist.html b/dev/reference/check_files_exist.html
index dcd2093..5566ea4 100644
--- a/dev/reference/check_files_exist.html
+++ b/dev/reference/check_files_exist.html
@@ -63,7 +63,7 @@ <h2 id="value">Value<a class="anchor" aria-label="anchor" href="#value"></a></h2
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/check_thread_count.html b/dev/reference/check_thread_count.html
index 8e42c9b..3ca4bba 100644
--- a/dev/reference/check_thread_count.html
+++ b/dev/reference/check_thread_count.html
@@ -72,7 +72,7 @@ <h2 id="value">Value<a class="anchor" aria-label="anchor" href="#value"></a></h2
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/dot-get_threads.html b/dev/reference/dot-get_threads.html
index 3d46848..557553b 100644
--- a/dev/reference/dot-get_threads.html
+++ b/dev/reference/dot-get_threads.html
@@ -45,7 +45,7 @@ <h2 id="ref-usage">Usage<a class="anchor" aria-label="anchor" href="#ref-usage">
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/get_granges_string.html b/dev/reference/get_granges_string.html
index 21b7ab6..3810776 100644
--- a/dev/reference/get_granges_string.html
+++ b/dev/reference/get_granges_string.html
@@ -69,7 +69,7 @@ <h2 id="details">Details<a class="anchor" aria-label="anchor" href="#details"></
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/get_log_level.html b/dev/reference/get_log_level.html
index 289082f..b60f4d7 100644
--- a/dev/reference/get_log_level.html
+++ b/dev/reference/get_log_level.html
@@ -53,7 +53,7 @@ <h2 id="value">Value<a class="anchor" aria-label="anchor" href="#value"></a></h2
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/get_omp_threads.html b/dev/reference/get_omp_threads.html
index 244a65f..5949c50 100644
--- a/dev/reference/get_omp_threads.html
+++ b/dev/reference/get_omp_threads.html
@@ -76,7 +76,7 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/get_threads.html b/dev/reference/get_threads.html
index cbb5f7a..992fe53 100644
--- a/dev/reference/get_threads.html
+++ b/dev/reference/get_threads.html
@@ -72,7 +72,7 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/index.html b/dev/reference/index.html
index 0999b1a..df678f0 100644
--- a/dev/reference/index.html
+++ b/dev/reference/index.html
@@ -145,7 +145,7 @@ <h3 id="logging">Logging<a class="anchor" aria-label="anchor" href="#logging"></
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/query_all.html b/dev/reference/query_all.html
index 37ae10e..4fe0d74 100644
--- a/dev/reference/query_all.html
+++ b/dev/reference/query_all.html
@@ -157,11 +157,11 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 <span class="r-in"><span><span class="co"># make a vector of regions</span></span></span>
 <span class="r-in"><span><span class="va">regions</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"chr1:1-6"</span>, <span class="st">"chr1:7-10"</span>, <span class="st">"chr1:11-14"</span><span class="op">)</span></span></span>
 <span class="r-in"><span><span class="fu">query_all</span><span class="op">(</span><span class="va">bedfiles</span>, <span class="va">regions</span><span class="op">)</span></span></span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [16:55:47.927007] [iscream::query_all] [info] Querying 3 regions from 4 bedfiles</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [17:10:40.479319] [iscream::query_all] [info] Querying 3 regions from 4 bedfiles</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> </span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [16:55:47.927390] [iscream::query_all] [info] Creating metadata vectors</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [16:55:47.927440] [iscream::query_all] [info] nrows 10000 - 9993 extra rows allocated with 0 resizes</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [16:55:47.927446] [iscream::query_all] [info] Creating dense matrix</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [17:10:40.479756] [iscream::query_all] [info] Creating metadata vectors</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [17:10:40.479803] [iscream::query_all] [info] nrows 10000 - 9993 extra rows allocated with 0 resizes</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [17:10:40.479807] [iscream::query_all] [info] Creating dense matrix</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> $M</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span>         a b c d</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> chr1:2  1 0 0 1</span>
@@ -203,7 +203,7 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/query_chroms.html b/dev/reference/query_chroms.html
index 64043a4..e077cb4 100644
--- a/dev/reference/query_chroms.html
+++ b/dev/reference/query_chroms.html
@@ -63,7 +63,7 @@ <h2 id="value">Value<a class="anchor" aria-label="anchor" href="#value"></a></h2
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/scan_tabix.html b/dev/reference/scan_tabix.html
index 9848ef0..b1c14eb 100644
--- a/dev/reference/scan_tabix.html
+++ b/dev/reference/scan_tabix.html
@@ -59,7 +59,7 @@ <h2 id="arguments">Arguments<a class="anchor" aria-label="anchor" href="#argumen
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/set_log_level.html b/dev/reference/set_log_level.html
index 39c4dae..cc2fbc5 100644
--- a/dev/reference/set_log_level.html
+++ b/dev/reference/set_log_level.html
@@ -68,7 +68,7 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/set_threads.html b/dev/reference/set_threads.html
index abba89e..77681bd 100644
--- a/dev/reference/set_threads.html
+++ b/dev/reference/set_threads.html
@@ -102,7 +102,7 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/summarize_regions.html b/dev/reference/summarize_regions.html
index 62ed5a1..0e83e3c 100644
--- a/dev/reference/summarize_regions.html
+++ b/dev/reference/summarize_regions.html
@@ -166,8 +166,8 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 <span class="r-in"><span><span class="co"># make a vector of regions</span></span></span>
 <span class="r-in"><span><span class="va">regions</span> <span class="op">&lt;-</span> <span class="fu"><a href="https://rdrr.io/r/base/c.html" class="external-link">c</a></span><span class="op">(</span><span class="st">"chr1:1-6"</span>, <span class="st">"chr1:7-10"</span>, <span class="st">"chr1:11-14"</span><span class="op">)</span></span></span>
 <span class="r-in"><span><span class="fu">summarize_regions</span><span class="op">(</span><span class="va">bedfiles</span>, <span class="va">regions</span><span class="op">)</span></span></span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [16:55:55.914749] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [16:55:55.914764] [iscream::summarize_regions] [info] using sum, mean, median, stddev, variance, min, max, range, cpg_count</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [17:10:48.581726] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [17:10:48.581740] [iscream::summarize_regions] [info] using sum, mean, median, stddev, variance, min, max, range, cpg_count</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span>       Feature Sample coverage.sum M.sum coverage.mean    M.mean coverage.median</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> 1    chr1:1-6      a            4     2      1.333333 0.6666667             1.0</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> 2   chr1:7-10      a            3     1      1.500000 0.5000000             1.5</span>
@@ -208,8 +208,8 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 <span class="r-out co"><span class="r-pr">#&gt;</span> 11     0            2     1              1       1         2</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> 12     1            1     1              0       0         1</span>
 <span class="r-in"><span><span class="fu">summarize_regions</span><span class="op">(</span><span class="va">bedfiles</span>, <span class="va">regions</span>, mval <span class="op">=</span> <span class="cn">FALSE</span><span class="op">)</span></span></span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [16:55:55.932774] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [16:55:55.932789] [iscream::summarize_regions] [info] using sum, mean, median, stddev, variance, min, max, range, cpg_count</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [17:10:48.599579] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [17:10:48.599591] [iscream::summarize_regions] [info] using sum, mean, median, stddev, variance, min, max, range, cpg_count</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span>       Feature Sample coverage.sum beta.sum coverage.mean beta.mean</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> 1    chr1:1-6      a            4      2.0      1.333333 0.6666667</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> 2   chr1:7-10      a            3      0.5      1.500000 0.2500000</span>
@@ -263,8 +263,8 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 <span class="r-out co"><span class="r-pr">#&gt;</span> 11        0.5         2</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> 12        0.0         1</span>
 <span class="r-in"><span><span class="fu">summarize_regions</span><span class="op">(</span><span class="va">bedfiles</span>, <span class="va">regions</span>, fun <span class="op">=</span> <span class="st">"sum"</span><span class="op">)</span></span></span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [16:55:55.949919] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles</span>
-<span class="r-out co"><span class="r-pr">#&gt;</span> [16:55:55.949933] [iscream::summarize_regions] [info] using sum</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [17:10:48.616708] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles</span>
+<span class="r-out co"><span class="r-pr">#&gt;</span> [17:10:48.616717] [iscream::summarize_regions] [info] using sum</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span>       Feature Sample coverage.sum M.sum</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> 1    chr1:1-6      a            4     2</span>
 <span class="r-out co"><span class="r-pr">#&gt;</span> 2   chr1:7-10      a            3     1</span>
@@ -285,7 +285,7 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/tabix.html b/dev/reference/tabix.html
index aa9d5d2..6f667ea 100644
--- a/dev/reference/tabix.html
+++ b/dev/reference/tabix.html
@@ -107,7 +107,7 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/validate_log_level.html b/dev/reference/validate_log_level.html
index 356a48d..032168b 100644
--- a/dev/reference/validate_log_level.html
+++ b/dev/reference/validate_log_level.html
@@ -71,7 +71,7 @@ <h2 id="ref-examples">Examples<a class="anchor" aria-label="anchor" href="#ref-e
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/verify_aligner_or_stop.html b/dev/reference/verify_aligner_or_stop.html
index a695390..6375f7f 100644
--- a/dev/reference/verify_aligner_or_stop.html
+++ b/dev/reference/verify_aligner_or_stop.html
@@ -58,7 +58,7 @@ <h2 id="value">Value<a class="anchor" aria-label="anchor" href="#value"></a></h2
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/verify_files_or_stop.html b/dev/reference/verify_files_or_stop.html
index 1b50e5d..17ddaca 100644
--- a/dev/reference/verify_files_or_stop.html
+++ b/dev/reference/verify_files_or_stop.html
@@ -63,7 +63,7 @@ <h2 id="value">Value<a class="anchor" aria-label="anchor" href="#value"></a></h2
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/reference/verify_filetype.html b/dev/reference/verify_filetype.html
index a00ec0d..43aeb0c 100644
--- a/dev/reference/verify_filetype.html
+++ b/dev/reference/verify_filetype.html
@@ -67,7 +67,7 @@ <h2 id="value">Value<a class="anchor" aria-label="anchor" href="#value"></a></h2
 
 
     <footer><div class="pkgdown-footer-left">
-  <p>Developed by James Eapen.</p>
+  <p>Developed by James Eapen, Hui Shen.</p>
 </div>
 
 <div class="pkgdown-footer-right">
diff --git a/dev/search.json b/dev/search.json
index bc65e5f..d43fd7d 100644
--- a/dev/search.json
+++ b/dev/search.json
@@ -1 +1 @@
-[{"path":"https://huishenlab.github.io/iscream/dev/LICENSE.html","id":null,"dir":"","previous_headings":"","what":"MIT License","title":"MIT License","text":"Copyright (c) 2023 iscream authors Permission hereby granted, free charge, person obtaining copy software associated documentation files (“Software”), deal Software without restriction, including without limitation rights use, copy, modify, merge, publish, distribute, sublicense, /sell copies Software, permit persons Software furnished , subject following conditions: copyright notice permission notice shall included copies substantial portions Software. SOFTWARE PROVIDED “”, WITHOUT WARRANTY KIND, EXPRESS IMPLIED, INCLUDING LIMITED WARRANTIES MERCHANTABILITY, FITNESS PARTICULAR PURPOSE NONINFRINGEMENT. EVENT SHALL AUTHORS COPYRIGHT HOLDERS LIABLE CLAIM, DAMAGES LIABILITY, WHETHER ACTION CONTRACT, TORT OTHERWISE, ARISING , CONNECTION SOFTWARE USE DEALINGS SOFTWARE.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/authors.html","id":null,"dir":"","previous_headings":"","what":"Authors","title":"Authors and Citation","text":"James Eapen. Author, maintainer.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/authors.html","id":"citation","dir":"","previous_headings":"","what":"Citation","title":"Authors and Citation","text":"Eapen J (2024). iscream: Integrating Single-Cell Results Exploring Analyzing Methylation. R package version 0.0.0.9000, https://github.com/huishenlab/iscream/, https://huishenlab.github.io/iscream/.","code":"@Manual{,   title = {iscream: Integrating Single-Cell Results for Exploring and Analyzing Methylation},   author = {James Eapen},   year = {2024},   note = {R package version 0.0.0.9000, https://github.com/huishenlab/iscream/},   url = {https://huishenlab.github.io/iscream/}, }"},{"path":"https://huishenlab.github.io/iscream/dev/index.html","id":"iscream-","dir":"","previous_headings":"","what":"iscream","title":"Integrating Single-Cell Results for Exploring and Analyzing Methylation","text":"Integrating Single-Cell Results Exploring Analyzing Methylation Analysis visualization Whole Genome Bisulfite Sequencing (WGBS) data requires reading aligned sequencing data formats existing packages like BSseq scMET can analyze. Getting data -disk formats like bedfiles matrix methylation values can difficult , nearly 30 million CpGs, WGBS data can quite large. iscream aims efficiently read alinged (sc)WGBS data formats can used packages. iscream uses htslib query genomic regions make matrices BSSeq aggregate methylated reads scMET.","code":""},{"path":[]},{"path":[]},{"path":"https://huishenlab.github.io/iscream/dev/index.html","id":"htslib","dir":"","previous_headings":"Installation > System dependencies","what":"htslib","title":"Integrating Single-Cell Results for Exploring and Analyzing Methylation","text":"iscream depends htslib header files. may installed package manager: ubuntu/debian: libhts-dev fedora/RHEL: htslib-devel brew: htslib nixpkgs: htslib may also found among HPC modules - make sure PKG_CONFIG_PATH environment variable includes pkgconfig location installation htslib. can verify htslib development libraries installed pkg-config:","code":"pkg-config --cflags --libs htslib"},{"path":"https://huishenlab.github.io/iscream/dev/index.html","id":"gcc","dir":"","previous_headings":"Installation > System dependencies","what":"GCC","title":"Integrating Single-Cell Results for Exploring and Analyzing Methylation","text":"GNU GCC must installed OpenMP support. usually installed default Linux systems, may need manually installed MacOS use iscream multiple threads1.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/index.html","id":"github","dir":"","previous_headings":"Installation","what":"GitHub","title":"Integrating Single-Cell Results for Exploring and Analyzing Methylation","text":"can install development version Github cloning repo running can also use R devtools package: pak:","code":"git clone https://github.com/huishenlab/iscream R CMD INSTALL iscream devtools::install_github(\"huishenlab/iscream\") pak::pkg_install(\"huishenlab/iscream\")"},{"path":"https://huishenlab.github.io/iscream/dev/index.html","id":"usage","dir":"","previous_headings":"Installation","what":"Usage","title":"Integrating Single-Cell Results for Exploring and Analyzing Methylation","text":"user guide available package website. Bug reports may submitted GitHub issues.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_all.html","id":null,"dir":"Reference","previous_headings":"","what":"Make a Bsseq object — Cpp_query_all","title":"Make a Bsseq object — Cpp_query_all","text":"Make Bsseq object","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_all.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Make a Bsseq object — Cpp_query_all","text":"","code":"Cpp_query_all(bedfiles, regions, bismark, merged, sparse, prealloc, nthreads)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_all.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Make a Bsseq object — Cpp_query_all","text":"bedfiles vector bedfiles regions vector regions bismark Whether input bismark coverage file prealloc number rows initialize matrices nthreads Set number threads use overriding \"iscream.threads\" option. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_chroms.html","id":null,"dir":"Reference","previous_headings":"","what":"Query the chromosomes or seqnames from a vector of files — Cpp_query_chroms","title":"Query the chromosomes or seqnames from a vector of files — Cpp_query_chroms","text":"Query chromosomes seqnames vector files","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_chroms.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Query the chromosomes or seqnames from a vector of files — Cpp_query_chroms","text":"","code":"Cpp_query_chroms(bedfile_vec, nthreads)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_chroms.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Query the chromosomes or seqnames from a vector of files — Cpp_query_chroms","text":"bedfile_vec vector bedfile paths","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_chroms.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Query the chromosomes or seqnames from a vector of files — Cpp_query_chroms","text":"vector seqnames","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_interval.html","id":null,"dir":"Reference","previous_headings":"","what":"Get reads from a single genomic region from one tabixed bed file to return as CharacterVector — Cpp_query_interval","title":"Get reads from a single genomic region from one tabixed bed file to return as CharacterVector — Cpp_query_interval","text":"Get reads single genomic region one tabixed bed file return CharacterVector","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_interval.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get reads from a single genomic region from one tabixed bed file to return as CharacterVector — Cpp_query_interval","text":"","code":"Cpp_query_interval(bedfile, regions)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_interval.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Get reads from a single genomic region from one tabixed bed file to return as CharacterVector — Cpp_query_interval","text":"bedfile name bed file - must corresponding tabix file name .tbi extension regions vector region strings form \"chr:start-end\"","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_set_log_level.html","id":null,"dir":"Reference","previous_headings":"","what":"spdlog Logging Lever Setter — Cpp_set_log_level","title":"spdlog Logging Lever Setter — Cpp_set_log_level","text":"helper function turn logging level given string current logging level","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_set_log_level.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"spdlog Logging Lever Setter — Cpp_set_log_level","text":"","code":"Cpp_set_log_level(name)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_set_log_level.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"spdlog Logging Lever Setter — Cpp_set_log_level","text":"name string logging level. Value understood , decreasing verbosity ‘trace’, ‘debug’, ‘info’, ‘warning’, ‘error’, ‘critical’, ‘’. Unrecognised names equivalent ‘’.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_set_log_level.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"spdlog Logging Lever Setter — Cpp_set_log_level","text":"Nothing returned.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_summarize_regions.html","id":null,"dir":"Reference","previous_headings":"","what":"Apply a function over CpGs within features — Cpp_summarize_regions","title":"Apply a function over CpGs within features — Cpp_summarize_regions","text":"function called summarize_regions() since sanity checks C++ side.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_summarize_regions.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Apply a function over CpGs within features — Cpp_summarize_regions","text":"","code":"Cpp_summarize_regions(   bedfiles,   regions,   fun_vec,   mval,   bismark,   region_rownames = FALSE,   nthreads = 1L )"},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_summarize_regions.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Apply a function over CpGs within features — Cpp_summarize_regions","text":"bedfiles vector bedfile paths regions vector genomic regions fun_vec Vector armadillo-supported stats functions apply CpGs ' regions: \"sum\", \"mean\", \"median\", \"stddev\", \"variance\" \"cpg_count\", \"min\",\"max\", \"range\". mval Calculates M values TRUE, use beta values FALSE bismark input bismark column format instead BISCUIT region_rownames Whether set rownames regions strings. necessary regions vector unnamed. names, \"Feature\" column set names rownames set regions string nthreads Number cores use. See details.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_summarize_regions.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Apply a function over CpGs within features — Cpp_summarize_regions","text":"optimal number threads depends number bedfiles, set half available OpenMP cores. See ?get_threads details. can manaully set set_threads().","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_files_exist.html","id":null,"dir":"Reference","previous_headings":"","what":"Check that files exist — check_files_exist","title":"Check that files exist — check_files_exist","text":"Check files exist","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_files_exist.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Check that files exist — check_files_exist","text":"","code":"check_files_exist(files_vec, error_file_prefix = \"Bedfile\")"},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_files_exist.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Check that files exist — check_files_exist","text":"files_vec vector file paths error_file_prefix Error message prefix 'Bedfile' vs 'Tabix file'","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_files_exist.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Check that files exist — check_files_exist","text":"TRUE input bedfiles associated tabix index file. FALSE ","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_thread_count.html","id":null,"dir":"Reference","previous_headings":"","what":"Check that the required threads are available — check_thread_count","title":"Check that the required threads are available — check_thread_count","text":"Check required threads available","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_thread_count.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Check that the required threads are available — check_thread_count","text":"","code":"check_thread_count(   n_threads,   avail_threads = availableCores(),   opt_set = FALSE )"},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_thread_count.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Check that the required threads are available — check_thread_count","text":"n_threads number threads check availability avail_threads number threads available system. Defaults parallelly::availableCores() opt_set Whether iscream.threads options set","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_thread_count.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Check that the required threads are available — check_thread_count","text":"n_threads requested number threads available stops ","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/dot-get_threads.html","id":null,"dir":"Reference","previous_headings":"","what":"Get the number of threads to use in functions from param or option — .get_threads","title":"Get the number of threads to use in functions from param or option — .get_threads","text":"Get number threads use functions param option","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/dot-get_threads.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get the number of threads to use in functions from param or option — .get_threads","text":"","code":".get_threads(nthreads)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_granges_string.html","id":null,"dir":"Reference","previous_headings":"","what":"GRanges to region strings — get_granges_string","title":"GRanges to region strings — get_granges_string","text":"Convert GRanges object vector strings","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_granges_string.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"GRanges to region strings — get_granges_string","text":"","code":"get_granges_string(gr, sep = c(\":\", \"-\"))"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_granges_string.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"GRanges to region strings — get_granges_string","text":"gr GRanges object sep Vector separators use genomic string. First element used separate chromosome coordinates, second separator used separate start end coordinates.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_granges_string.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"GRanges to region strings — get_granges_string","text":"character vector","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_granges_string.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"GRanges to region strings — get_granges_string","text":"function adapted stuart-lab/signac.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_log_level.html","id":null,"dir":"Reference","previous_headings":"","what":"Get the current log level — get_log_level","title":"Get the current log level — get_log_level","text":"Can handle spdlogs levels, iscream functions supports \"info\" \"debug\"","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_log_level.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get the current log level — get_log_level","text":"","code":"get_log_level()"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_log_level.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Get the current log level — get_log_level","text":"current logging level string","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_omp_threads.html","id":null,"dir":"Reference","previous_headings":"","what":"Get the number of available threads from OpenMP. — get_omp_threads","title":"Get the number of available threads from OpenMP. — get_omp_threads","text":"queries number available threads usign OpenMP, reliably provide accurate available thread count. get reliable count accounts environment variables HPC schedulers, use get_threads()`. function pulled github.com/rdatatable/data.table","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_omp_threads.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get the number of available threads from OpenMP. — get_omp_threads","text":"","code":"get_omp_threads(verbose)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_omp_threads.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Get the number of available threads from OpenMP. — get_omp_threads","text":"verbose Whether verbose available omp threads","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_omp_threads.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Get the number of available threads from OpenMP. — get_omp_threads","text":"","code":"get_omp_threads(verbose = TRUE) #>   OpenMP version (_OPENMP)       201511 #>   omp_get_num_procs()            4 #>   omp_get_thread_limit()         2147483647 #>   omp_get_max_threads()          4 #> [1] 4"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_threads.html","id":null,"dir":"Reference","previous_headings":"","what":"Get the number of available threads — get_threads","title":"Get the number of available threads — get_threads","text":"Gets number threads iscream currently set use, whether \"iscream.threads\" option set many threads available use. set number threads use set_threads() set iscream.threads option ~/.Rprofile. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_threads.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get the number of available threads — get_threads","text":"","code":"get_threads()"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_threads.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Get the number of available threads — get_threads","text":"named vector: use_threads = number threads iscream use opt_set = whether option set user avail_threads = number available threads reported parallelly::availableCores","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_threads.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Get the number of available threads — get_threads","text":"","code":"get_threads() #>   use_threads       opt_set avail_threads  #>             1             1             4"},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_all.html","id":null,"dir":"Reference","previous_headings":"","what":"Query all CpGs from input genomic regions — query_all","title":"Query all CpGs from input genomic regions — query_all","text":"Queries provided regions produces M Coverage matrices genomic positions. Parallelized across files using threads \"iscream.threads\" option. output query_all may used create BSseq object: .call(BSseq, query_all(...)).","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_all.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Query all CpGs from input genomic regions — query_all","text":"","code":"query_all(   bedfiles,   regions,   aligner = \"biscuit\",   merged = TRUE,   sparse = FALSE,   prealloc = 10000,   nthreads = NULL )"},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_all.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Query all CpGs from input genomic regions — query_all","text":"bedfiles vector bedfile paths regions vector genomic regions strings aligner aligner used produce BED files - one \"biscuit\", \"bismark\", \"bsbolt\". merged Whether input strands merged/collapsed sparse Whether return M coverage matrices sparse matrices (\"dgCMatrix\"). Set TRUE scWGBS data prealloc number rows initialize matrices . number methyltion loci approximately known, can reduce runtime fewer resizes need made. nthreads Set number threads use overriding \"iscream.threads\" option. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_all.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Query all CpGs from input genomic regions — query_all","text":"named list M coverage matrices character vector chromosomes numeric vector corresponding CpG base positions character vector input sample names","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_all.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Query all CpGs from input genomic regions — query_all","text":"","code":"bedfiles <- system.file(\"extdata\", package = \"iscream\") |>   list.files(pattern = \"[a|b|c|d].bed.gz$\", full.names = TRUE) # examine the bedfiles colnames <- c(\"chr\", \"start\", \"end\", \"beta\", \"coverage\") lapply(bedfiles, function(i) knitr::kable(read.table(i, col.names = colnames))) #> [[1]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|  1.0|        1| #> |chr1 |     3|   4|  1.0|        1| #> |chr1 |     5|   6|  0.0|        2| #> |chr1 |     7|   8|  0.0|        1| #> |chr1 |     9|  10|  0.5|        2| #> |chr1 |    11|  12|  1.0|        2| #> |chr1 |    13|  14|  1.0|        3| #>  #> [[2]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|    0|        2| #> |chr1 |     5|   6|    1|        2| #> |chr1 |     7|   8|    1|        1| #> |chr1 |    11|  12|    0|        2| #> |chr1 |    13|  14|    1|        1| #>  #> [[3]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     3|   4|    1|        2| #> |chr1 |     7|   8|    0|        2| #> |chr1 |     9|  10|    1|        1| #>  #> [[4]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|  1.0|        1| #> |chr1 |     3|   4|  1.0|        2| #> |chr1 |     7|   8|  0.0|        1| #> |chr1 |     9|  10|  0.5|        2| #> |chr1 |    13|  14|  1.0|        1| #>   # make a vector of regions regions <- c(\"chr1:1-6\", \"chr1:7-10\", \"chr1:11-14\") query_all(bedfiles, regions) #> [16:55:47.927007] [iscream::query_all] [info] Querying 3 regions from 4 bedfiles #>  #> [16:55:47.927390] [iscream::query_all] [info] Creating metadata vectors #> [16:55:47.927440] [iscream::query_all] [info] nrows 10000 - 9993 extra rows allocated with 0 resizes #> [16:55:47.927446] [iscream::query_all] [info] Creating dense matrix #> $M #>         a b c d #> chr1:2  1 0 0 1 #> chr1:4  1 0 2 2 #> chr1:6  0 2 0 0 #> chr1:8  0 1 0 0 #> chr1:10 1 0 1 1 #> chr1:12 2 0 0 0 #> chr1:14 3 1 0 1 #>  #> $Cov #>         a b c d #> chr1:2  1 2 0 1 #> chr1:4  1 0 2 2 #> chr1:6  2 2 0 0 #> chr1:8  1 1 2 1 #> chr1:10 2 0 1 2 #> chr1:12 2 2 0 0 #> chr1:14 3 1 0 1 #>  #> $pos #> [1]  1  3  5  7  9 11 13 #>  #> $chr #> [1] \"chr1\" \"chr1\" \"chr1\" \"chr1\" \"chr1\" \"chr1\" \"chr1\" #>  #> $sampleNames #> [1] \"a\" \"b\" \"c\" \"d\" #>  # for BSseq object run if (FALSE) { # \\dontrun{ library(bsseq) do.call(BSseq, query_all(bedfiles, regions)) } # }"},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_chroms.html","id":null,"dir":"Reference","previous_headings":"","what":"Query the chromosomes or seqnames from a vector of files — query_chroms","title":"Query the chromosomes or seqnames from a vector of files — query_chroms","text":"Query chromosomes seqnames vector files","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_chroms.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Query the chromosomes or seqnames from a vector of files — query_chroms","text":"","code":"query_chroms(bedfiles, nthreads = NULL)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_chroms.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Query the chromosomes or seqnames from a vector of files — query_chroms","text":"bedfiles vector bedfile paths nthreads Set number threads use overriding \"iscream.threads\" option. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_chroms.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Query the chromosomes or seqnames from a vector of files — query_chroms","text":"vector seqnames","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/scan_tabix.html","id":null,"dir":"Reference","previous_headings":"","what":"Get namde list of reads from a single genomic region from one tabixed bed file — scan_tabix","title":"Get namde list of reads from a single genomic region from one tabixed bed file — scan_tabix","text":"Get namde list reads single genomic region one tabixed bed file","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/scan_tabix.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get namde list of reads from a single genomic region from one tabixed bed file — scan_tabix","text":"","code":"scan_tabix(bedfile, regions)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/scan_tabix.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Get namde list of reads from a single genomic region from one tabixed bed file — scan_tabix","text":"bedfile name bed file - must corresponding tabix file name .tbi extension regions vector region strings form \"chr:start-end\"","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_log_level.html","id":null,"dir":"Reference","previous_headings":"","what":"Set logging level — set_log_level","title":"Set logging level — set_log_level","text":"Set logging level","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_log_level.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Set logging level — set_log_level","text":"","code":"set_log_level(level = \"info\")"},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_log_level.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Set logging level — set_log_level","text":"level logging verbosity level use \"info\": default gives provides basic information number files regions used function \"debug\": verbose row allocations, many CpGs found region, filename parsing etc. mode used one thread R output messages multiple threads without crashing.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_log_level.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Set logging level — set_log_level","text":"None; sets log level provided level","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_log_level.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Set logging level — set_log_level","text":"","code":"set_log_level(\"info\")"},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":null,"dir":"Reference","previous_headings":"","what":"Set the number of available threads — set_threads","title":"Set the number of available threads — set_threads","text":"Sets \"iscream.threads\" option n_threads. see many threads available see ?get_threads().","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Set the number of available threads — set_threads","text":"","code":"set_threads(n_threads)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Set the number of available threads — set_threads","text":"n_threads number threads use","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Set the number of available threads — set_threads","text":"NULL. Sets iscream.threads option requested number threads available","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Set the number of available threads — set_threads","text":"iscream uses OpenMP parallelize certain functions. can use many threads available system varying degrees performance improvements. get_threads() function uses parallelly::availableCores() report number available threads. Although OpenMP can detect number available cores, high preformance computers (HPCs) resource allocating job schedulers like SLURM, OpenMP may detect available threads across HPC limit cores allocated scheduler. system administrator set limits, may result job taking resources jobs. limits, trying use threads available reduce iscream's performance. Job schedulers typically environment variable (e.g. SLURM_CPUS_ON_NODE SLURM) gives actual number available cores. , hyperthreaded systems, count may double available processors. Using hyperthreading guarantee performance improvement - may better set number threads half reported number. parallelly::availableCores() takes HPC scheduler/CRAN/Bioconductor limits account reporting number available threads may reliably report hyperthreading ('system' 'nproc'). set number threads without call set_threads() every session, put   .Rprofile See help('Rprofile') information startup options. Functions currently using multithreading: tabix() query_chroms() query_all() summarize_regions()","code":"options(iscream.threads = [n_threads])"},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Set the number of available threads — set_threads","text":"","code":"(ncores <- parallelly::availableCores()) #> system  #>      4  if (FALSE) { # \\dontrun{ set_threads(ncores) } # }"},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":null,"dir":"Reference","previous_headings":"","what":"Summarize CpGs methylation information over genomic regions — summarize_regions","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"Run summarizing functions CpGs bedfiles across genomic regions. Parallelized across files using threads \"iscream.threads\" option.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"","code":"summarize_regions(   bedfiles,   regions,   fun = \"all\",   aligner = \"biscuit\",   mval = TRUE,   set_region_rownames = FALSE,   nthreads = NULL )"},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"bedfiles vector bedfile paths regions vector genomic regions strings. named vector provided, names used feature column instead genomic regions string fun Function apply region. See details. aligner aligner used produce BED files - one \"biscuit\", \"bismark\", \"bsbolt\". mval Whether calculate M value (coverage \\(\\times \\beta\\)) use beta value applying function. set_region_rownames Use region strings returned data frame's rownames. Can useful named vector regions want rownames feature names. nthreads Set number threads use overriding \"iscream.threads\" option. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"data.frame","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"Supported fun arguments given . functions, setting mval = FALSE use beta values instead M value: Sum: \"sum\" Mean: \"mean\" Median: \"median\" Standard deviation: \"stddev\" Variance: \"variance\" Minimum: \"min\" Maximum: \"max\" Range: \"range\" . CpGs region: \"cpg_count\" summarizing computations backed Armadillo library. See https://arma.sourceforge.net/docs.html#stats_fns futher details supported functions","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"","code":"bedfiles <- system.file(\"extdata\", package = \"iscream\") |>   list.files(pattern = \"[a|b|c|d].bed.gz$\", full.names = TRUE) # examine the bedfiles colnames <- c(\"chr\", \"start\", \"end\", \"beta\", \"coverage\") lapply(bedfiles, function(i) knitr::kable(read.table(i, col.names = colnames))) #> [[1]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|  1.0|        1| #> |chr1 |     3|   4|  1.0|        1| #> |chr1 |     5|   6|  0.0|        2| #> |chr1 |     7|   8|  0.0|        1| #> |chr1 |     9|  10|  0.5|        2| #> |chr1 |    11|  12|  1.0|        2| #> |chr1 |    13|  14|  1.0|        3| #>  #> [[2]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|    0|        2| #> |chr1 |     5|   6|    1|        2| #> |chr1 |     7|   8|    1|        1| #> |chr1 |    11|  12|    0|        2| #> |chr1 |    13|  14|    1|        1| #>  #> [[3]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     3|   4|    1|        2| #> |chr1 |     7|   8|    0|        2| #> |chr1 |     9|  10|    1|        1| #>  #> [[4]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|  1.0|        1| #> |chr1 |     3|   4|  1.0|        2| #> |chr1 |     7|   8|  0.0|        1| #> |chr1 |     9|  10|  0.5|        2| #> |chr1 |    13|  14|  1.0|        1| #>   # make a vector of regions regions <- c(\"chr1:1-6\", \"chr1:7-10\", \"chr1:11-14\") summarize_regions(bedfiles, regions) #> [16:55:55.914749] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles #> [16:55:55.914764] [iscream::summarize_regions] [info] using sum, mean, median, stddev, variance, min, max, range, cpg_count #>       Feature Sample coverage.sum M.sum coverage.mean    M.mean coverage.median #> 1    chr1:1-6      a            4     2      1.333333 0.6666667             1.0 #> 2   chr1:7-10      a            3     1      1.500000 0.5000000             1.5 #> 3  chr1:11-14      a            5     5      2.500000 2.5000000             2.5 #> 4    chr1:1-6      b            4     2      2.000000 1.0000000             2.0 #> 5   chr1:7-10      b            1     1      1.000000 1.0000000             1.0 #> 6  chr1:11-14      b            3     1      1.500000 0.5000000             1.5 #> 7    chr1:1-6      c            2     2      2.000000 2.0000000             2.0 #> 8   chr1:7-10      c            3     1      1.500000 0.5000000             1.5 #> 9  chr1:11-14      c           NA    NA            NA        NA              NA #> 10   chr1:1-6      d            3     3      1.500000 1.5000000             1.5 #> 11  chr1:7-10      d            3     1      1.500000 0.5000000             1.5 #> 12 chr1:11-14      d            1     1      1.000000 1.0000000             1.0 #>    M.median coverage.stddev  M.stddev coverage.variance M.variance coverage.min #> 1       1.0       0.5773503 0.5773503         0.3333333  0.3333333            1 #> 2       0.5       0.7071068 0.7071068         0.5000000  0.5000000            1 #> 3       2.5       0.7071068 0.7071068         0.5000000  0.5000000            2 #> 4       1.0       0.0000000 1.4142136         0.0000000  2.0000000            2 #> 5       1.0       0.0000000 0.0000000         0.0000000  0.0000000            1 #> 6       0.5       0.7071068 0.7071068         0.5000000  0.5000000            1 #> 7       2.0       0.0000000 0.0000000         0.0000000  0.0000000            2 #> 8       0.5       0.7071068 0.7071068         0.5000000  0.5000000            1 #> 9        NA              NA        NA                NA         NA           NA #> 10      1.5       0.7071068 0.7071068         0.5000000  0.5000000            1 #> 11      0.5       0.7071068 0.7071068         0.5000000  0.5000000            1 #> 12      1.0       0.0000000 0.0000000         0.0000000  0.0000000            1 #>    M.min coverage.max M.max coverage.range M.range cpg_count #> 1      0            2     1              1       1         3 #> 2      0            2     1              1       1         2 #> 3      2            3     3              1       1         2 #> 4      0            2     2              0       2         2 #> 5      1            1     1              0       0         1 #> 6      0            2     1              1       1         2 #> 7      2            2     2              0       0         1 #> 8      0            2     1              1       1         2 #> 9     NA           NA    NA             NA      NA        NA #> 10     1            2     2              1       1         2 #> 11     0            2     1              1       1         2 #> 12     1            1     1              0       0         1 summarize_regions(bedfiles, regions, mval = FALSE) #> [16:55:55.932774] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles #> [16:55:55.932789] [iscream::summarize_regions] [info] using sum, mean, median, stddev, variance, min, max, range, cpg_count #>       Feature Sample coverage.sum beta.sum coverage.mean beta.mean #> 1    chr1:1-6      a            4      2.0      1.333333 0.6666667 #> 2   chr1:7-10      a            3      0.5      1.500000 0.2500000 #> 3  chr1:11-14      a            5      2.0      2.500000 1.0000000 #> 4    chr1:1-6      b            4      1.0      2.000000 0.5000000 #> 5   chr1:7-10      b            1      1.0      1.000000 1.0000000 #> 6  chr1:11-14      b            3      1.0      1.500000 0.5000000 #> 7    chr1:1-6      c            2      1.0      2.000000 1.0000000 #> 8   chr1:7-10      c            3      1.0      1.500000 0.5000000 #> 9  chr1:11-14      c           NA       NA            NA        NA #> 10   chr1:1-6      d            3      2.0      1.500000 1.0000000 #> 11  chr1:7-10      d            3      0.5      1.500000 0.2500000 #> 12 chr1:11-14      d            1      1.0      1.000000 1.0000000 #>    coverage.median beta.median coverage.stddev beta.stddev coverage.variance #> 1              1.0        1.00       0.5773503   0.5773503         0.3333333 #> 2              1.5        0.25       0.7071068   0.3535534         0.5000000 #> 3              2.5        1.00       0.7071068   0.0000000         0.5000000 #> 4              2.0        0.50       0.0000000   0.7071068         0.0000000 #> 5              1.0        1.00       0.0000000   0.0000000         0.0000000 #> 6              1.5        0.50       0.7071068   0.7071068         0.5000000 #> 7              2.0        1.00       0.0000000   0.0000000         0.0000000 #> 8              1.5        0.50       0.7071068   0.7071068         0.5000000 #> 9               NA          NA              NA          NA                NA #> 10             1.5        1.00       0.7071068   0.0000000         0.5000000 #> 11             1.5        0.25       0.7071068   0.3535534         0.5000000 #> 12             1.0        1.00       0.0000000   0.0000000         0.0000000 #>    beta.variance coverage.min beta.min coverage.max beta.max coverage.range #> 1      0.3333333            1        0            2      1.0              1 #> 2      0.1250000            1        0            2      0.5              1 #> 3      0.0000000            2        1            3      1.0              1 #> 4      0.5000000            2        0            2      1.0              0 #> 5      0.0000000            1        1            1      1.0              0 #> 6      0.5000000            1        0            2      1.0              1 #> 7      0.0000000            2        1            2      1.0              0 #> 8      0.5000000            1        0            2      1.0              1 #> 9             NA           NA       NA           NA       NA             NA #> 10     0.0000000            1        1            2      1.0              1 #> 11     0.1250000            1        0            2      0.5              1 #> 12     0.0000000            1        1            1      1.0              0 #>    beta.range cpg_count #> 1         1.0         3 #> 2         0.5         2 #> 3         0.0         2 #> 4         1.0         2 #> 5         0.0         1 #> 6         1.0         2 #> 7         0.0         1 #> 8         1.0         2 #> 9          NA        NA #> 10        0.0         2 #> 11        0.5         2 #> 12        0.0         1 summarize_regions(bedfiles, regions, fun = \"sum\") #> [16:55:55.949919] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles #> [16:55:55.949933] [iscream::summarize_regions] [info] using sum #>       Feature Sample coverage.sum M.sum #> 1    chr1:1-6      a            4     2 #> 2   chr1:7-10      a            3     1 #> 3  chr1:11-14      a            5     5 #> 4    chr1:1-6      b            4     2 #> 5   chr1:7-10      b            1     1 #> 6  chr1:11-14      b            3     1 #> 7    chr1:1-6      c            2     2 #> 8   chr1:7-10      c            3     1 #> 9  chr1:11-14      c           NA    NA #> 10   chr1:1-6      d            3     3 #> 11  chr1:7-10      d            3     1 #> 12 chr1:11-14      d            1     1"},{"path":"https://huishenlab.github.io/iscream/dev/reference/tabix.html","id":null,"dir":"Reference","previous_headings":"","what":"Query lines from a tabixed bedfile — tabix","title":"Query lines from a tabixed bedfile — tabix","text":"Query lines tabixed bedfile","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/tabix.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Query lines from a tabixed bedfile — tabix","text":"","code":"tabix(   bedfiles,   regions,   aligner = \"biscuit\",   colnames = NULL,   raw = FALSE,   nthreads = NULL )"},{"path":"https://huishenlab.github.io/iscream/dev/reference/tabix.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Query lines from a tabixed bedfile — tabix","text":"bedfiles bedfiles queried regions vector genomic region strings aligner aligner used produce BED files - one \"biscuit\", \"bismark\", \"bsbolt\". set result data.table's column names based argument. colnames vector column names result data.table. Set bedfile supported aligners general bedfile. raw Set true give named list raw strings regions style Rsamtools::scanTabix instead data.table nthreads Set number threads use overriding \"iscream.threads\" option. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/tabix.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Query lines from a tabixed bedfile — tabix","text":"data.table","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/tabix.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Query lines from a tabixed bedfile — tabix","text":"","code":"bedfiles <- system.file(\"extdata\", package = \"iscream\") |>   list.files(pattern = \"[a|b|c|d].bed.gz$\", full.names = TRUE) regions <- c(\"chr1:1-6\", \"chr1:7-10\", \"chr1:11-14\") tabix(bedfiles[1], regions, colnames = c(\"chr\", \"start\", \"end\", \"beta\", \"coverage\")) #>       chr start   end  beta coverage #>    <char> <num> <num> <num>    <num> #> 1:   chr1     1     2   1.0        1 #> 2:   chr1     3     4   1.0        1 #> 3:   chr1     5     6   0.0        2 #> 4:   chr1     7     8   0.0        1 #> 5:   chr1     9    10   0.5        2 #> 6:   chr1    11    12   1.0        2 #> 7:   chr1    13    14   1.0        3"},{"path":"https://huishenlab.github.io/iscream/dev/reference/validate_log_level.html","id":null,"dir":"Reference","previous_headings":"","what":"Validate provided logging level — validate_log_level","title":"Validate provided logging level — validate_log_level","text":"\"info\" \"debug\" currently supported, \"debug\" supported using 1 thread","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/validate_log_level.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Validate provided logging level — validate_log_level","text":"","code":"validate_log_level(level = get_log_level(), n_threads)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/validate_log_level.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Validate provided logging level — validate_log_level","text":"level logging level validate n_threads number threads next iscream function call use","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/validate_log_level.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Validate provided logging level — validate_log_level","text":"None; sets log level provide level","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/validate_log_level.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Validate provided logging level — validate_log_level","text":"","code":"set_log_level(\"info\")"},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_aligner_or_stop.html","id":null,"dir":"Reference","previous_headings":"","what":"Validate provided aligner — verify_aligner_or_stop","title":"Validate provided aligner — verify_aligner_or_stop","text":"\"biscuit\", \"bismark\", \"bsbolt\" currently supported","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_aligner_or_stop.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Validate provided aligner — verify_aligner_or_stop","text":"","code":"verify_aligner_or_stop(aligner)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_aligner_or_stop.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Validate provided aligner — verify_aligner_or_stop","text":"aligner input alinger","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_aligner_or_stop.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Validate provided aligner — verify_aligner_or_stop","text":"true; quits input among supported_aligners","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_files_or_stop.html","id":null,"dir":"Reference","previous_headings":"","what":"Verify that bedfiles are tabixed — verify_files_or_stop","title":"Verify that bedfiles are tabixed — verify_files_or_stop","text":"Verify bedfiles tabixed","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_files_or_stop.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Verify that bedfiles are tabixed — verify_files_or_stop","text":"","code":"verify_files_or_stop(bedfiles, verify_tabix = TRUE)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_files_or_stop.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Verify that bedfiles are tabixed — verify_files_or_stop","text":"bedfiles vector bedfile paths verify_tabix Whether verify presence tabix files","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_files_or_stop.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Verify that bedfiles are tabixed — verify_files_or_stop","text":"TRUE input bedfiles associated tabix index file. FALSE ","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_filetype.html","id":null,"dir":"Reference","previous_headings":"","what":"Verify that the input bedfiles are of the type specified by the input aligner — verify_filetype","title":"Verify that the input bedfiles are of the type specified by the input aligner — verify_filetype","text":"Verify input bedfiles type specified input aligner","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_filetype.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Verify that the input bedfiles are of the type specified by the input aligner — verify_filetype","text":"","code":"verify_filetype(bedfiles, aligner, stop_on_error = FALSE)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_filetype.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Verify that the input bedfiles are of the type specified by the input aligner — verify_filetype","text":"bedfiles vector bedfile paths aligner aligner chosen stop_on_error Whether warn stop aligner-filename mismatch","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_filetype.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Verify that the input bedfiles are of the type specified by the input aligner — verify_filetype","text":"TRUE input bedfiles associated tabix index file. FALSE ","code":""}]
+[{"path":"https://huishenlab.github.io/iscream/dev/LICENSE.html","id":null,"dir":"","previous_headings":"","what":"MIT License","title":"MIT License","text":"Copyright (c) 2023 iscream authors Permission hereby granted, free charge, person obtaining copy software associated documentation files (“Software”), deal Software without restriction, including without limitation rights use, copy, modify, merge, publish, distribute, sublicense, /sell copies Software, permit persons Software furnished , subject following conditions: copyright notice permission notice shall included copies substantial portions Software. SOFTWARE PROVIDED “”, WITHOUT WARRANTY KIND, EXPRESS IMPLIED, INCLUDING LIMITED WARRANTIES MERCHANTABILITY, FITNESS PARTICULAR PURPOSE NONINFRINGEMENT. EVENT SHALL AUTHORS COPYRIGHT HOLDERS LIABLE CLAIM, DAMAGES LIABILITY, WHETHER ACTION CONTRACT, TORT OTHERWISE, ARISING , CONNECTION SOFTWARE USE DEALINGS SOFTWARE.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/authors.html","id":null,"dir":"","previous_headings":"","what":"Authors","title":"Authors and Citation","text":"James Eapen. Author, maintainer. Jacob Morrison. Contributor. Hui Shen. Funder.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/authors.html","id":"citation","dir":"","previous_headings":"","what":"Citation","title":"Authors and Citation","text":"Eapen J (2024). iscream: Integrating Single-Cell Results Exploring Analyzing Methylation. R package version 0.0.0.9000, https://github.com/huishenlab/iscream/, https://huishenlab.github.io/iscream/.","code":"@Manual{,   title = {iscream: Integrating Single-Cell Results for Exploring and Analyzing Methylation},   author = {James Eapen},   year = {2024},   note = {R package version 0.0.0.9000, https://github.com/huishenlab/iscream/},   url = {https://huishenlab.github.io/iscream/}, }"},{"path":"https://huishenlab.github.io/iscream/dev/index.html","id":"iscream-","dir":"","previous_headings":"","what":"iscream","title":"Integrating Single-Cell Results for Exploring and Analyzing Methylation","text":"Integrating Single-Cell Results Exploring Analyzing Methylation Analysis visualization Whole Genome Bisulfite Sequencing (WGBS) data requires reading aligned sequencing data formats existing packages like BSseq scMET can analyze. Getting data -disk formats like bedfiles matrix methylation values can difficult , nearly 30 million CpGs, WGBS data can quite large. iscream aims efficiently read alinged (sc)WGBS data formats can used packages. iscream uses htslib query genomic regions make matrices BSSeq aggregate methylated reads scMET.","code":""},{"path":[]},{"path":[]},{"path":"https://huishenlab.github.io/iscream/dev/index.html","id":"htslib","dir":"","previous_headings":"Installation > System dependencies","what":"htslib","title":"Integrating Single-Cell Results for Exploring and Analyzing Methylation","text":"iscream depends htslib header files. may installed package manager: ubuntu/debian: libhts-dev fedora/RHEL: htslib-devel brew: htslib nixpkgs: htslib may also found among HPC modules - make sure PKG_CONFIG_PATH environment variable includes pkgconfig location installation htslib. can verify htslib development libraries installed pkg-config:","code":"pkg-config --cflags --libs htslib"},{"path":"https://huishenlab.github.io/iscream/dev/index.html","id":"gcc","dir":"","previous_headings":"Installation > System dependencies","what":"GCC","title":"Integrating Single-Cell Results for Exploring and Analyzing Methylation","text":"GNU GCC must installed OpenMP support. usually installed default Linux systems, may need manually installed MacOS use iscream multiple threads1.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/index.html","id":"github","dir":"","previous_headings":"Installation","what":"GitHub","title":"Integrating Single-Cell Results for Exploring and Analyzing Methylation","text":"can install development version Github cloning repo running can also use R devtools package: pak:","code":"git clone https://github.com/huishenlab/iscream R CMD INSTALL iscream devtools::install_github(\"huishenlab/iscream\") pak::pkg_install(\"huishenlab/iscream\")"},{"path":"https://huishenlab.github.io/iscream/dev/index.html","id":"usage","dir":"","previous_headings":"Installation","what":"Usage","title":"Integrating Single-Cell Results for Exploring and Analyzing Methylation","text":"user guide available package website. Bug reports may submitted GitHub issues.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_all.html","id":null,"dir":"Reference","previous_headings":"","what":"Make a Bsseq object — Cpp_query_all","title":"Make a Bsseq object — Cpp_query_all","text":"Make Bsseq object","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_all.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Make a Bsseq object — Cpp_query_all","text":"","code":"Cpp_query_all(bedfiles, regions, bismark, merged, sparse, prealloc, nthreads)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_all.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Make a Bsseq object — Cpp_query_all","text":"bedfiles vector bedfiles regions vector regions bismark Whether input bismark coverage file prealloc number rows initialize matrices nthreads Set number threads use overriding \"iscream.threads\" option. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_chroms.html","id":null,"dir":"Reference","previous_headings":"","what":"Query the chromosomes or seqnames from a vector of files — Cpp_query_chroms","title":"Query the chromosomes or seqnames from a vector of files — Cpp_query_chroms","text":"Query chromosomes seqnames vector files","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_chroms.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Query the chromosomes or seqnames from a vector of files — Cpp_query_chroms","text":"","code":"Cpp_query_chroms(bedfile_vec, nthreads)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_chroms.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Query the chromosomes or seqnames from a vector of files — Cpp_query_chroms","text":"bedfile_vec vector bedfile paths","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_chroms.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Query the chromosomes or seqnames from a vector of files — Cpp_query_chroms","text":"vector seqnames","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_interval.html","id":null,"dir":"Reference","previous_headings":"","what":"Get reads from a single genomic region from one tabixed bed file to return as CharacterVector — Cpp_query_interval","title":"Get reads from a single genomic region from one tabixed bed file to return as CharacterVector — Cpp_query_interval","text":"Get reads single genomic region one tabixed bed file return CharacterVector","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_interval.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get reads from a single genomic region from one tabixed bed file to return as CharacterVector — Cpp_query_interval","text":"","code":"Cpp_query_interval(bedfile, regions)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_query_interval.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Get reads from a single genomic region from one tabixed bed file to return as CharacterVector — Cpp_query_interval","text":"bedfile name bed file - must corresponding tabix file name .tbi extension regions vector region strings form \"chr:start-end\"","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_set_log_level.html","id":null,"dir":"Reference","previous_headings":"","what":"spdlog Logging Lever Setter — Cpp_set_log_level","title":"spdlog Logging Lever Setter — Cpp_set_log_level","text":"helper function turn logging level given string current logging level","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_set_log_level.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"spdlog Logging Lever Setter — Cpp_set_log_level","text":"","code":"Cpp_set_log_level(name)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_set_log_level.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"spdlog Logging Lever Setter — Cpp_set_log_level","text":"name string logging level. Value understood , decreasing verbosity ‘trace’, ‘debug’, ‘info’, ‘warning’, ‘error’, ‘critical’, ‘’. Unrecognised names equivalent ‘’.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_set_log_level.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"spdlog Logging Lever Setter — Cpp_set_log_level","text":"Nothing returned.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_summarize_regions.html","id":null,"dir":"Reference","previous_headings":"","what":"Apply a function over CpGs within features — Cpp_summarize_regions","title":"Apply a function over CpGs within features — Cpp_summarize_regions","text":"function called summarize_regions() since sanity checks C++ side.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_summarize_regions.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Apply a function over CpGs within features — Cpp_summarize_regions","text":"","code":"Cpp_summarize_regions(   bedfiles,   regions,   fun_vec,   mval,   bismark,   region_rownames = FALSE,   nthreads = 1L )"},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_summarize_regions.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Apply a function over CpGs within features — Cpp_summarize_regions","text":"bedfiles vector bedfile paths regions vector genomic regions fun_vec Vector armadillo-supported stats functions apply CpGs ' regions: \"sum\", \"mean\", \"median\", \"stddev\", \"variance\" \"cpg_count\", \"min\",\"max\", \"range\". mval Calculates M values TRUE, use beta values FALSE bismark input bismark column format instead BISCUIT region_rownames Whether set rownames regions strings. necessary regions vector unnamed. names, \"Feature\" column set names rownames set regions string nthreads Number cores use. See details.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/Cpp_summarize_regions.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Apply a function over CpGs within features — Cpp_summarize_regions","text":"optimal number threads depends number bedfiles, set half available OpenMP cores. See ?get_threads details. can manaully set set_threads().","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_files_exist.html","id":null,"dir":"Reference","previous_headings":"","what":"Check that files exist — check_files_exist","title":"Check that files exist — check_files_exist","text":"Check files exist","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_files_exist.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Check that files exist — check_files_exist","text":"","code":"check_files_exist(files_vec, error_file_prefix = \"Bedfile\")"},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_files_exist.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Check that files exist — check_files_exist","text":"files_vec vector file paths error_file_prefix Error message prefix 'Bedfile' vs 'Tabix file'","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_files_exist.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Check that files exist — check_files_exist","text":"TRUE input bedfiles associated tabix index file. FALSE ","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_thread_count.html","id":null,"dir":"Reference","previous_headings":"","what":"Check that the required threads are available — check_thread_count","title":"Check that the required threads are available — check_thread_count","text":"Check required threads available","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_thread_count.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Check that the required threads are available — check_thread_count","text":"","code":"check_thread_count(   n_threads,   avail_threads = availableCores(),   opt_set = FALSE )"},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_thread_count.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Check that the required threads are available — check_thread_count","text":"n_threads number threads check availability avail_threads number threads available system. Defaults parallelly::availableCores() opt_set Whether iscream.threads options set","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/check_thread_count.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Check that the required threads are available — check_thread_count","text":"n_threads requested number threads available stops ","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/dot-get_threads.html","id":null,"dir":"Reference","previous_headings":"","what":"Get the number of threads to use in functions from param or option — .get_threads","title":"Get the number of threads to use in functions from param or option — .get_threads","text":"Get number threads use functions param option","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/dot-get_threads.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get the number of threads to use in functions from param or option — .get_threads","text":"","code":".get_threads(nthreads)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_granges_string.html","id":null,"dir":"Reference","previous_headings":"","what":"GRanges to region strings — get_granges_string","title":"GRanges to region strings — get_granges_string","text":"Convert GRanges object vector strings","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_granges_string.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"GRanges to region strings — get_granges_string","text":"","code":"get_granges_string(gr, sep = c(\":\", \"-\"))"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_granges_string.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"GRanges to region strings — get_granges_string","text":"gr GRanges object sep Vector separators use genomic string. First element used separate chromosome coordinates, second separator used separate start end coordinates.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_granges_string.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"GRanges to region strings — get_granges_string","text":"character vector","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_granges_string.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"GRanges to region strings — get_granges_string","text":"function adapted stuart-lab/signac.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_log_level.html","id":null,"dir":"Reference","previous_headings":"","what":"Get the current log level — get_log_level","title":"Get the current log level — get_log_level","text":"Can handle spdlogs levels, iscream functions supports \"info\" \"debug\"","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_log_level.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get the current log level — get_log_level","text":"","code":"get_log_level()"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_log_level.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Get the current log level — get_log_level","text":"current logging level string","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_omp_threads.html","id":null,"dir":"Reference","previous_headings":"","what":"Get the number of available threads from OpenMP. — get_omp_threads","title":"Get the number of available threads from OpenMP. — get_omp_threads","text":"queries number available threads usign OpenMP, reliably provide accurate available thread count. get reliable count accounts environment variables HPC schedulers, use get_threads()`. function pulled github.com/rdatatable/data.table","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_omp_threads.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get the number of available threads from OpenMP. — get_omp_threads","text":"","code":"get_omp_threads(verbose)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_omp_threads.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Get the number of available threads from OpenMP. — get_omp_threads","text":"verbose Whether verbose available omp threads","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_omp_threads.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Get the number of available threads from OpenMP. — get_omp_threads","text":"","code":"get_omp_threads(verbose = TRUE) #>   OpenMP version (_OPENMP)       201511 #>   omp_get_num_procs()            4 #>   omp_get_thread_limit()         2147483647 #>   omp_get_max_threads()          4 #> [1] 4"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_threads.html","id":null,"dir":"Reference","previous_headings":"","what":"Get the number of available threads — get_threads","title":"Get the number of available threads — get_threads","text":"Gets number threads iscream currently set use, whether \"iscream.threads\" option set many threads available use. set number threads use set_threads() set iscream.threads option ~/.Rprofile. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_threads.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get the number of available threads — get_threads","text":"","code":"get_threads()"},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_threads.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Get the number of available threads — get_threads","text":"named vector: use_threads = number threads iscream use opt_set = whether option set user avail_threads = number available threads reported parallelly::availableCores","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/get_threads.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Get the number of available threads — get_threads","text":"","code":"get_threads() #>   use_threads       opt_set avail_threads  #>             1             1             4"},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_all.html","id":null,"dir":"Reference","previous_headings":"","what":"Query all CpGs from input genomic regions — query_all","title":"Query all CpGs from input genomic regions — query_all","text":"Queries provided regions produces M Coverage matrices genomic positions. Parallelized across files using threads \"iscream.threads\" option. output query_all may used create BSseq object: .call(BSseq, query_all(...)).","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_all.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Query all CpGs from input genomic regions — query_all","text":"","code":"query_all(   bedfiles,   regions,   aligner = \"biscuit\",   merged = TRUE,   sparse = FALSE,   prealloc = 10000,   nthreads = NULL )"},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_all.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Query all CpGs from input genomic regions — query_all","text":"bedfiles vector bedfile paths regions vector genomic regions strings aligner aligner used produce BED files - one \"biscuit\", \"bismark\", \"bsbolt\". merged Whether input strands merged/collapsed sparse Whether return M coverage matrices sparse matrices (\"dgCMatrix\"). Set TRUE scWGBS data prealloc number rows initialize matrices . number methyltion loci approximately known, can reduce runtime fewer resizes need made. nthreads Set number threads use overriding \"iscream.threads\" option. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_all.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Query all CpGs from input genomic regions — query_all","text":"named list M coverage matrices character vector chromosomes numeric vector corresponding CpG base positions character vector input sample names","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_all.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Query all CpGs from input genomic regions — query_all","text":"","code":"bedfiles <- system.file(\"extdata\", package = \"iscream\") |>   list.files(pattern = \"[a|b|c|d].bed.gz$\", full.names = TRUE) # examine the bedfiles colnames <- c(\"chr\", \"start\", \"end\", \"beta\", \"coverage\") lapply(bedfiles, function(i) knitr::kable(read.table(i, col.names = colnames))) #> [[1]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|  1.0|        1| #> |chr1 |     3|   4|  1.0|        1| #> |chr1 |     5|   6|  0.0|        2| #> |chr1 |     7|   8|  0.0|        1| #> |chr1 |     9|  10|  0.5|        2| #> |chr1 |    11|  12|  1.0|        2| #> |chr1 |    13|  14|  1.0|        3| #>  #> [[2]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|    0|        2| #> |chr1 |     5|   6|    1|        2| #> |chr1 |     7|   8|    1|        1| #> |chr1 |    11|  12|    0|        2| #> |chr1 |    13|  14|    1|        1| #>  #> [[3]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     3|   4|    1|        2| #> |chr1 |     7|   8|    0|        2| #> |chr1 |     9|  10|    1|        1| #>  #> [[4]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|  1.0|        1| #> |chr1 |     3|   4|  1.0|        2| #> |chr1 |     7|   8|  0.0|        1| #> |chr1 |     9|  10|  0.5|        2| #> |chr1 |    13|  14|  1.0|        1| #>   # make a vector of regions regions <- c(\"chr1:1-6\", \"chr1:7-10\", \"chr1:11-14\") query_all(bedfiles, regions) #> [17:10:40.479319] [iscream::query_all] [info] Querying 3 regions from 4 bedfiles #>  #> [17:10:40.479756] [iscream::query_all] [info] Creating metadata vectors #> [17:10:40.479803] [iscream::query_all] [info] nrows 10000 - 9993 extra rows allocated with 0 resizes #> [17:10:40.479807] [iscream::query_all] [info] Creating dense matrix #> $M #>         a b c d #> chr1:2  1 0 0 1 #> chr1:4  1 0 2 2 #> chr1:6  0 2 0 0 #> chr1:8  0 1 0 0 #> chr1:10 1 0 1 1 #> chr1:12 2 0 0 0 #> chr1:14 3 1 0 1 #>  #> $Cov #>         a b c d #> chr1:2  1 2 0 1 #> chr1:4  1 0 2 2 #> chr1:6  2 2 0 0 #> chr1:8  1 1 2 1 #> chr1:10 2 0 1 2 #> chr1:12 2 2 0 0 #> chr1:14 3 1 0 1 #>  #> $pos #> [1]  1  3  5  7  9 11 13 #>  #> $chr #> [1] \"chr1\" \"chr1\" \"chr1\" \"chr1\" \"chr1\" \"chr1\" \"chr1\" #>  #> $sampleNames #> [1] \"a\" \"b\" \"c\" \"d\" #>  # for BSseq object run if (FALSE) { # \\dontrun{ library(bsseq) do.call(BSseq, query_all(bedfiles, regions)) } # }"},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_chroms.html","id":null,"dir":"Reference","previous_headings":"","what":"Query the chromosomes or seqnames from a vector of files — query_chroms","title":"Query the chromosomes or seqnames from a vector of files — query_chroms","text":"Query chromosomes seqnames vector files","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_chroms.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Query the chromosomes or seqnames from a vector of files — query_chroms","text":"","code":"query_chroms(bedfiles, nthreads = NULL)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_chroms.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Query the chromosomes or seqnames from a vector of files — query_chroms","text":"bedfiles vector bedfile paths nthreads Set number threads use overriding \"iscream.threads\" option. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/query_chroms.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Query the chromosomes or seqnames from a vector of files — query_chroms","text":"vector seqnames","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/scan_tabix.html","id":null,"dir":"Reference","previous_headings":"","what":"Get namde list of reads from a single genomic region from one tabixed bed file — scan_tabix","title":"Get namde list of reads from a single genomic region from one tabixed bed file — scan_tabix","text":"Get namde list reads single genomic region one tabixed bed file","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/scan_tabix.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Get namde list of reads from a single genomic region from one tabixed bed file — scan_tabix","text":"","code":"scan_tabix(bedfile, regions)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/scan_tabix.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Get namde list of reads from a single genomic region from one tabixed bed file — scan_tabix","text":"bedfile name bed file - must corresponding tabix file name .tbi extension regions vector region strings form \"chr:start-end\"","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_log_level.html","id":null,"dir":"Reference","previous_headings":"","what":"Set logging level — set_log_level","title":"Set logging level — set_log_level","text":"Set logging level","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_log_level.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Set logging level — set_log_level","text":"","code":"set_log_level(level = \"info\")"},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_log_level.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Set logging level — set_log_level","text":"level logging verbosity level use \"info\": default gives provides basic information number files regions used function \"debug\": verbose row allocations, many CpGs found region, filename parsing etc. mode used one thread R output messages multiple threads without crashing.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_log_level.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Set logging level — set_log_level","text":"None; sets log level provided level","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_log_level.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Set logging level — set_log_level","text":"","code":"set_log_level(\"info\")"},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":null,"dir":"Reference","previous_headings":"","what":"Set the number of available threads — set_threads","title":"Set the number of available threads — set_threads","text":"Sets \"iscream.threads\" option n_threads. see many threads available see ?get_threads().","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Set the number of available threads — set_threads","text":"","code":"set_threads(n_threads)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Set the number of available threads — set_threads","text":"n_threads number threads use","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Set the number of available threads — set_threads","text":"NULL. Sets iscream.threads option requested number threads available","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Set the number of available threads — set_threads","text":"iscream uses OpenMP parallelize certain functions. can use many threads available system varying degrees performance improvements. get_threads() function uses parallelly::availableCores() report number available threads. Although OpenMP can detect number available cores, high preformance computers (HPCs) resource allocating job schedulers like SLURM, OpenMP may detect available threads across HPC limit cores allocated scheduler. system administrator set limits, may result job taking resources jobs. limits, trying use threads available reduce iscream's performance. Job schedulers typically environment variable (e.g. SLURM_CPUS_ON_NODE SLURM) gives actual number available cores. , hyperthreaded systems, count may double available processors. Using hyperthreading guarantee performance improvement - may better set number threads half reported number. parallelly::availableCores() takes HPC scheduler/CRAN/Bioconductor limits account reporting number available threads may reliably report hyperthreading ('system' 'nproc'). set number threads without call set_threads() every session, put   .Rprofile See help('Rprofile') information startup options. Functions currently using multithreading: tabix() query_chroms() query_all() summarize_regions()","code":"options(iscream.threads = [n_threads])"},{"path":"https://huishenlab.github.io/iscream/dev/reference/set_threads.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Set the number of available threads — set_threads","text":"","code":"(ncores <- parallelly::availableCores()) #> system  #>      4  if (FALSE) { # \\dontrun{ set_threads(ncores) } # }"},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":null,"dir":"Reference","previous_headings":"","what":"Summarize CpGs methylation information over genomic regions — summarize_regions","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"Run summarizing functions CpGs bedfiles across genomic regions. Parallelized across files using threads \"iscream.threads\" option.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"","code":"summarize_regions(   bedfiles,   regions,   fun = \"all\",   aligner = \"biscuit\",   mval = TRUE,   set_region_rownames = FALSE,   nthreads = NULL )"},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"bedfiles vector bedfile paths regions vector genomic regions strings. named vector provided, names used feature column instead genomic regions string fun Function apply region. See details. aligner aligner used produce BED files - one \"biscuit\", \"bismark\", \"bsbolt\". mval Whether calculate M value (coverage \\(\\times \\beta\\)) use beta value applying function. set_region_rownames Use region strings returned data frame's rownames. Can useful named vector regions want rownames feature names. nthreads Set number threads use overriding \"iscream.threads\" option. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"data.frame","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":"details","dir":"Reference","previous_headings":"","what":"Details","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"Supported fun arguments given . functions, setting mval = FALSE use beta values instead M value: Sum: \"sum\" Mean: \"mean\" Median: \"median\" Standard deviation: \"stddev\" Variance: \"variance\" Minimum: \"min\" Maximum: \"max\" Range: \"range\" . CpGs region: \"cpg_count\" summarizing computations backed Armadillo library. See https://arma.sourceforge.net/docs.html#stats_fns futher details supported functions","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/summarize_regions.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Summarize CpGs methylation information over genomic regions — summarize_regions","text":"","code":"bedfiles <- system.file(\"extdata\", package = \"iscream\") |>   list.files(pattern = \"[a|b|c|d].bed.gz$\", full.names = TRUE) # examine the bedfiles colnames <- c(\"chr\", \"start\", \"end\", \"beta\", \"coverage\") lapply(bedfiles, function(i) knitr::kable(read.table(i, col.names = colnames))) #> [[1]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|  1.0|        1| #> |chr1 |     3|   4|  1.0|        1| #> |chr1 |     5|   6|  0.0|        2| #> |chr1 |     7|   8|  0.0|        1| #> |chr1 |     9|  10|  0.5|        2| #> |chr1 |    11|  12|  1.0|        2| #> |chr1 |    13|  14|  1.0|        3| #>  #> [[2]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|    0|        2| #> |chr1 |     5|   6|    1|        2| #> |chr1 |     7|   8|    1|        1| #> |chr1 |    11|  12|    0|        2| #> |chr1 |    13|  14|    1|        1| #>  #> [[3]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     3|   4|    1|        2| #> |chr1 |     7|   8|    0|        2| #> |chr1 |     9|  10|    1|        1| #>  #> [[4]] #>  #>  #> |chr  | start| end| beta| coverage| #> |:----|-----:|---:|----:|--------:| #> |chr1 |     1|   2|  1.0|        1| #> |chr1 |     3|   4|  1.0|        2| #> |chr1 |     7|   8|  0.0|        1| #> |chr1 |     9|  10|  0.5|        2| #> |chr1 |    13|  14|  1.0|        1| #>   # make a vector of regions regions <- c(\"chr1:1-6\", \"chr1:7-10\", \"chr1:11-14\") summarize_regions(bedfiles, regions) #> [17:10:48.581726] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles #> [17:10:48.581740] [iscream::summarize_regions] [info] using sum, mean, median, stddev, variance, min, max, range, cpg_count #>       Feature Sample coverage.sum M.sum coverage.mean    M.mean coverage.median #> 1    chr1:1-6      a            4     2      1.333333 0.6666667             1.0 #> 2   chr1:7-10      a            3     1      1.500000 0.5000000             1.5 #> 3  chr1:11-14      a            5     5      2.500000 2.5000000             2.5 #> 4    chr1:1-6      b            4     2      2.000000 1.0000000             2.0 #> 5   chr1:7-10      b            1     1      1.000000 1.0000000             1.0 #> 6  chr1:11-14      b            3     1      1.500000 0.5000000             1.5 #> 7    chr1:1-6      c            2     2      2.000000 2.0000000             2.0 #> 8   chr1:7-10      c            3     1      1.500000 0.5000000             1.5 #> 9  chr1:11-14      c           NA    NA            NA        NA              NA #> 10   chr1:1-6      d            3     3      1.500000 1.5000000             1.5 #> 11  chr1:7-10      d            3     1      1.500000 0.5000000             1.5 #> 12 chr1:11-14      d            1     1      1.000000 1.0000000             1.0 #>    M.median coverage.stddev  M.stddev coverage.variance M.variance coverage.min #> 1       1.0       0.5773503 0.5773503         0.3333333  0.3333333            1 #> 2       0.5       0.7071068 0.7071068         0.5000000  0.5000000            1 #> 3       2.5       0.7071068 0.7071068         0.5000000  0.5000000            2 #> 4       1.0       0.0000000 1.4142136         0.0000000  2.0000000            2 #> 5       1.0       0.0000000 0.0000000         0.0000000  0.0000000            1 #> 6       0.5       0.7071068 0.7071068         0.5000000  0.5000000            1 #> 7       2.0       0.0000000 0.0000000         0.0000000  0.0000000            2 #> 8       0.5       0.7071068 0.7071068         0.5000000  0.5000000            1 #> 9        NA              NA        NA                NA         NA           NA #> 10      1.5       0.7071068 0.7071068         0.5000000  0.5000000            1 #> 11      0.5       0.7071068 0.7071068         0.5000000  0.5000000            1 #> 12      1.0       0.0000000 0.0000000         0.0000000  0.0000000            1 #>    M.min coverage.max M.max coverage.range M.range cpg_count #> 1      0            2     1              1       1         3 #> 2      0            2     1              1       1         2 #> 3      2            3     3              1       1         2 #> 4      0            2     2              0       2         2 #> 5      1            1     1              0       0         1 #> 6      0            2     1              1       1         2 #> 7      2            2     2              0       0         1 #> 8      0            2     1              1       1         2 #> 9     NA           NA    NA             NA      NA        NA #> 10     1            2     2              1       1         2 #> 11     0            2     1              1       1         2 #> 12     1            1     1              0       0         1 summarize_regions(bedfiles, regions, mval = FALSE) #> [17:10:48.599579] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles #> [17:10:48.599591] [iscream::summarize_regions] [info] using sum, mean, median, stddev, variance, min, max, range, cpg_count #>       Feature Sample coverage.sum beta.sum coverage.mean beta.mean #> 1    chr1:1-6      a            4      2.0      1.333333 0.6666667 #> 2   chr1:7-10      a            3      0.5      1.500000 0.2500000 #> 3  chr1:11-14      a            5      2.0      2.500000 1.0000000 #> 4    chr1:1-6      b            4      1.0      2.000000 0.5000000 #> 5   chr1:7-10      b            1      1.0      1.000000 1.0000000 #> 6  chr1:11-14      b            3      1.0      1.500000 0.5000000 #> 7    chr1:1-6      c            2      1.0      2.000000 1.0000000 #> 8   chr1:7-10      c            3      1.0      1.500000 0.5000000 #> 9  chr1:11-14      c           NA       NA            NA        NA #> 10   chr1:1-6      d            3      2.0      1.500000 1.0000000 #> 11  chr1:7-10      d            3      0.5      1.500000 0.2500000 #> 12 chr1:11-14      d            1      1.0      1.000000 1.0000000 #>    coverage.median beta.median coverage.stddev beta.stddev coverage.variance #> 1              1.0        1.00       0.5773503   0.5773503         0.3333333 #> 2              1.5        0.25       0.7071068   0.3535534         0.5000000 #> 3              2.5        1.00       0.7071068   0.0000000         0.5000000 #> 4              2.0        0.50       0.0000000   0.7071068         0.0000000 #> 5              1.0        1.00       0.0000000   0.0000000         0.0000000 #> 6              1.5        0.50       0.7071068   0.7071068         0.5000000 #> 7              2.0        1.00       0.0000000   0.0000000         0.0000000 #> 8              1.5        0.50       0.7071068   0.7071068         0.5000000 #> 9               NA          NA              NA          NA                NA #> 10             1.5        1.00       0.7071068   0.0000000         0.5000000 #> 11             1.5        0.25       0.7071068   0.3535534         0.5000000 #> 12             1.0        1.00       0.0000000   0.0000000         0.0000000 #>    beta.variance coverage.min beta.min coverage.max beta.max coverage.range #> 1      0.3333333            1        0            2      1.0              1 #> 2      0.1250000            1        0            2      0.5              1 #> 3      0.0000000            2        1            3      1.0              1 #> 4      0.5000000            2        0            2      1.0              0 #> 5      0.0000000            1        1            1      1.0              0 #> 6      0.5000000            1        0            2      1.0              1 #> 7      0.0000000            2        1            2      1.0              0 #> 8      0.5000000            1        0            2      1.0              1 #> 9             NA           NA       NA           NA       NA             NA #> 10     0.0000000            1        1            2      1.0              1 #> 11     0.1250000            1        0            2      0.5              1 #> 12     0.0000000            1        1            1      1.0              0 #>    beta.range cpg_count #> 1         1.0         3 #> 2         0.5         2 #> 3         0.0         2 #> 4         1.0         2 #> 5         0.0         1 #> 6         1.0         2 #> 7         0.0         1 #> 8         1.0         2 #> 9          NA        NA #> 10        0.0         2 #> 11        0.5         2 #> 12        0.0         1 summarize_regions(bedfiles, regions, fun = \"sum\") #> [17:10:48.616708] [iscream::summarize_regions] [info] Summarizing 3 regions from 4 bedfiles #> [17:10:48.616717] [iscream::summarize_regions] [info] using sum #>       Feature Sample coverage.sum M.sum #> 1    chr1:1-6      a            4     2 #> 2   chr1:7-10      a            3     1 #> 3  chr1:11-14      a            5     5 #> 4    chr1:1-6      b            4     2 #> 5   chr1:7-10      b            1     1 #> 6  chr1:11-14      b            3     1 #> 7    chr1:1-6      c            2     2 #> 8   chr1:7-10      c            3     1 #> 9  chr1:11-14      c           NA    NA #> 10   chr1:1-6      d            3     3 #> 11  chr1:7-10      d            3     1 #> 12 chr1:11-14      d            1     1"},{"path":"https://huishenlab.github.io/iscream/dev/reference/tabix.html","id":null,"dir":"Reference","previous_headings":"","what":"Query lines from a tabixed bedfile — tabix","title":"Query lines from a tabixed bedfile — tabix","text":"Query lines tabixed bedfile","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/tabix.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Query lines from a tabixed bedfile — tabix","text":"","code":"tabix(   bedfiles,   regions,   aligner = \"biscuit\",   colnames = NULL,   raw = FALSE,   nthreads = NULL )"},{"path":"https://huishenlab.github.io/iscream/dev/reference/tabix.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Query lines from a tabixed bedfile — tabix","text":"bedfiles bedfiles queried regions vector genomic region strings aligner aligner used produce BED files - one \"biscuit\", \"bismark\", \"bsbolt\". set result data.table's column names based argument. colnames vector column names result data.table. Set bedfile supported aligners general bedfile. raw Set true give named list raw strings regions style Rsamtools::scanTabix instead data.table nthreads Set number threads use overriding \"iscream.threads\" option. See ?set_threads information.","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/tabix.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Query lines from a tabixed bedfile — tabix","text":"data.table","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/tabix.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Query lines from a tabixed bedfile — tabix","text":"","code":"bedfiles <- system.file(\"extdata\", package = \"iscream\") |>   list.files(pattern = \"[a|b|c|d].bed.gz$\", full.names = TRUE) regions <- c(\"chr1:1-6\", \"chr1:7-10\", \"chr1:11-14\") tabix(bedfiles[1], regions, colnames = c(\"chr\", \"start\", \"end\", \"beta\", \"coverage\")) #>       chr start   end  beta coverage #>    <char> <num> <num> <num>    <num> #> 1:   chr1     1     2   1.0        1 #> 2:   chr1     3     4   1.0        1 #> 3:   chr1     5     6   0.0        2 #> 4:   chr1     7     8   0.0        1 #> 5:   chr1     9    10   0.5        2 #> 6:   chr1    11    12   1.0        2 #> 7:   chr1    13    14   1.0        3"},{"path":"https://huishenlab.github.io/iscream/dev/reference/validate_log_level.html","id":null,"dir":"Reference","previous_headings":"","what":"Validate provided logging level — validate_log_level","title":"Validate provided logging level — validate_log_level","text":"\"info\" \"debug\" currently supported, \"debug\" supported using 1 thread","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/validate_log_level.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Validate provided logging level — validate_log_level","text":"","code":"validate_log_level(level = get_log_level(), n_threads)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/validate_log_level.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Validate provided logging level — validate_log_level","text":"level logging level validate n_threads number threads next iscream function call use","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/validate_log_level.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Validate provided logging level — validate_log_level","text":"None; sets log level provide level","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/validate_log_level.html","id":"ref-examples","dir":"Reference","previous_headings":"","what":"Examples","title":"Validate provided logging level — validate_log_level","text":"","code":"set_log_level(\"info\")"},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_aligner_or_stop.html","id":null,"dir":"Reference","previous_headings":"","what":"Validate provided aligner — verify_aligner_or_stop","title":"Validate provided aligner — verify_aligner_or_stop","text":"\"biscuit\", \"bismark\", \"bsbolt\" currently supported","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_aligner_or_stop.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Validate provided aligner — verify_aligner_or_stop","text":"","code":"verify_aligner_or_stop(aligner)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_aligner_or_stop.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Validate provided aligner — verify_aligner_or_stop","text":"aligner input alinger","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_aligner_or_stop.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Validate provided aligner — verify_aligner_or_stop","text":"true; quits input among supported_aligners","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_files_or_stop.html","id":null,"dir":"Reference","previous_headings":"","what":"Verify that bedfiles are tabixed — verify_files_or_stop","title":"Verify that bedfiles are tabixed — verify_files_or_stop","text":"Verify bedfiles tabixed","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_files_or_stop.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Verify that bedfiles are tabixed — verify_files_or_stop","text":"","code":"verify_files_or_stop(bedfiles, verify_tabix = TRUE)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_files_or_stop.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Verify that bedfiles are tabixed — verify_files_or_stop","text":"bedfiles vector bedfile paths verify_tabix Whether verify presence tabix files","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_files_or_stop.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Verify that bedfiles are tabixed — verify_files_or_stop","text":"TRUE input bedfiles associated tabix index file. FALSE ","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_filetype.html","id":null,"dir":"Reference","previous_headings":"","what":"Verify that the input bedfiles are of the type specified by the input aligner — verify_filetype","title":"Verify that the input bedfiles are of the type specified by the input aligner — verify_filetype","text":"Verify input bedfiles type specified input aligner","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_filetype.html","id":"ref-usage","dir":"Reference","previous_headings":"","what":"Usage","title":"Verify that the input bedfiles are of the type specified by the input aligner — verify_filetype","text":"","code":"verify_filetype(bedfiles, aligner, stop_on_error = FALSE)"},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_filetype.html","id":"arguments","dir":"Reference","previous_headings":"","what":"Arguments","title":"Verify that the input bedfiles are of the type specified by the input aligner — verify_filetype","text":"bedfiles vector bedfile paths aligner aligner chosen stop_on_error Whether warn stop aligner-filename mismatch","code":""},{"path":"https://huishenlab.github.io/iscream/dev/reference/verify_filetype.html","id":"value","dir":"Reference","previous_headings":"","what":"Value","title":"Verify that the input bedfiles are of the type specified by the input aligner — verify_filetype","text":"TRUE input bedfiles associated tabix index file. FALSE ","code":""}]