biopipen 0.34.13__tar.gz → 0.34.14__tar.gz

{biopipen-0.34.13 → biopipen-0.34.14}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: biopipen
-Version: 0.34.13
+Version: 0.34.14
 Summary: Bioinformatics processes/pipelines that can be run from `pipen run`
 License: MIT
 Author: pwwang

biopipen-0.34.14/biopipen/__init__.py

@@ -0,0 +1 @@
+__version__ = "0.34.14"

{biopipen-0.34.13 → biopipen-0.34.14}/biopipen/scripts/scrna/CellTypeAnnotation-celltypist.R

@@ -111,7 +111,8 @@ if (file.exists(celltypist_outfile) &&
         command <- paste(command, "-v")
     }
     log$info("Running celltypist:")
-    print("- {command}")
+    # print("- {command}")
+    log$debug("  {command}")
     rc <- system(command)
     if (rc != 0) {
         stop("Failed to run celltypist. Check the job.stderr file to see the error message.")

{biopipen-0.34.13 → biopipen-0.34.14}/biopipen/scripts/scrna/MarkersFinder.R

@@ -607,8 +607,14 @@ run_case <- function(name) {
         return(invisible())
     }
 
-    case$object <- if (is.null(subset_)) srtobj else filter(srtobj, !!parse_expr(subset_))
+    # Let RunSeuratDEAnalysis handle the subset
+    case$subset <- subset_
+    case$object <- srtobj
     markers <- do_call(RunSeuratDEAnalysis, case)
+    case$object <- NULL  # Release memory
+    gc()
+
+    subobj <- if (is.null(subset_)) srtobj else filter(srtobj, !!parse_expr(subset_))
 
     if (is.null(case$ident_1)) {
         all_idents <- unique(as.character(markers[[case$group_by]]))
@@ -621,7 +627,7 @@ run_case <- function(name) {
 
             attr(ident_markers, "ident_1") <- ident
             enrich <- process_markers(ident_markers, info = info, case = list(
-                object = case$object,
+                object = subobj,
                 dbs = dbs,
                 group_by = case$group_by,
                 sigmarkers = sigmarkers,
@@ -638,7 +644,7 @@ run_case <- function(name) {
             log$info("- Visualizing all markers together ...")
             process_allmarkers(
                 markers,
-                object = case$object,
+                object = subobj,
                 comparison_by = case$group_by,
                 plotcases = allmarker_plots,
                 casename = name,
@@ -658,7 +664,7 @@ run_case <- function(name) {
     } else {
         info <- case_info(name, outdir, create = TRUE)
         enrich <- process_markers(markers, info = info, case = list(
-            object = case$object,
+            object = subobj,
             dbs = dbs,
             group_by = case$group_by,
             sigmarkers = sigmarkers,

{biopipen-0.34.13 → biopipen-0.34.14}/pyproject.toml

@@ -1,6 +1,6 @@
 [tool.poetry]
 name = "biopipen"
-version = "0.34.13"
+version = "0.34.14"
 description = "Bioinformatics processes/pipelines that can be run from `pipen run`"
 authors = ["pwwang <pwwang@pwwang.com>"]
 license = "MIT"

{biopipen-0.34.13 → biopipen-0.34.14}/setup.py

@@ -86,7 +86,7 @@ entry_points = \
 
 setup_kwargs = {
     'name': 'biopipen',
-    'version': '0.34.13',
+    'version': '0.34.14',
     'description': 'Bioinformatics processes/pipelines that can be run from `pipen run`',
     'long_description': 'None',
     'author': 'pwwang',

biopipen-0.34.13/biopipen/__init__.py

@@ -1 +0,0 @@
-__version__ = "0.34.13"