First official release

DESCRIPTION

@@ -1,7 +1,7 @@
 Package: numbat
-Title: Accurate reconstruction of copy number profiles in scRNAseq data using population haplotype phasing
+Title: Haplotype-aware CNV anlaysis from scRNA-seq
 URL: https://github.com/kharchenkolab/numbat, https://kharchenkolab.github.io/numbat
-Version: 0.1.3
+Version: 1.0.0
 Authors@R: c(person("Teng","Gao", email="[email protected]", role=c("cre", "aut")), person("Hirak", "Sarkar", email="[email protected]", role="aut"), person("Evan", "Biederstedt", email="[email protected]", role="aut"), person("Peter", "Kharchenko", email = "[email protected]", role = "aut"))
 Description: Numbat is a computational method that infers copy number variations (CNVs) in cancer scRNA-Seq data and reconstructs the tumor phylogeny. The R package integrates signals from gene expression, allelic ratio, and population haplotype structures to accurately infer allele-specific CNVs in single cells and reconstruct their lineage relationship. Numbat can be used to: 1. detect allele-specific copy number variations from single-cells; 2. differentiate tumor versus normal cells in the tumor microenvironment; 3. infer the clonal architecture and evolutionary history of profiled tumors. Numbat does not require tumor/normal-paired DNA or genotype data, but operates solely on the donor scRNA-data data (for example, 10x Cell Ranger output).
 License: MIT

R/main.R

@@ -235,7 +235,7 @@ run_numbat = function(
             p = plot_bulks(bulk_subtrees, min_LLR = min_LLR, use_pos = TRUE)
             ggsave(
                 glue('{out_dir}/bulk_subtrees_{i}.png'), p, 
-                width = 12, height = 2*length(unique(bulk_subtrees$sample)), dpi = 200
+                width = 12, height = 2*length(unique(bulk_subtrees$sample)), dpi = 250
             )
         }
 
@@ -253,7 +253,6 @@ run_numbat = function(
         bulk_subtrees = retest_bulks(
                 bulk_subtrees,
                 segs_consensus, 
-                # segs_loh = segs_loh,
                 diploid_chroms = diploid_chroms, 
                 min_LLR = min_LLR,
                 ncores = ncores
@@ -304,7 +303,7 @@ run_numbat = function(
             p = plot_bulks(bulk_clones, min_LLR = min_LLR, use_pos = TRUE)
             ggsave(
                 glue('{out_dir}/bulk_clones_{i}.png'), p, 
-                width = 12, height = 2*length(unique(bulk_clones$sample)), dpi = 200
+                width = 12, height = 2*length(unique(bulk_clones$sample)), dpi = 250
             )
         }
 
@@ -474,7 +473,7 @@ run_numbat = function(
                     clone_bar = TRUE
                 )
 
-                ggsave(glue('{out_dir}/panel_{i}.png'), panel, width = 8, height = 3.5, dpi = 200)
+                ggsave(glue('{out_dir}/panel_{i}.png'), panel, width = 8, height = 3.5, dpi = 250)
 
             },
             error = function(e) { 
@@ -556,7 +555,7 @@ run_numbat = function(
         p = plot_bulks(bulk_clones, min_LLR = min_LLR, use_pos = TRUE)
         ggsave(
             glue('{out_dir}/bulk_clones_final.png'), p, 
-            width = 12, height = 2*length(unique(bulk_clones$sample)), dpi = 200
+            width = 12, height = 2*length(unique(bulk_clones$sample)), dpi = 250
         )
     }