biopipen 0.34.14__py3-none-any.whl → 0.34.15__py3-none-any.whl

biopipen/__init__.py

@@ -1 +1 @@
-__version__ = "0.34.14"
+__version__ = "0.34.15"

biopipen/scripts/scrna/MarkersFinder.R

@@ -600,6 +600,28 @@ run_case <- function(name) {
 
             if (length(allenrich_plots) > 0 && !is.null(enriches) && nrow(enriches) > 0) {
                 log$info("- Visualizing all enrichments together ...")
+                # add other metadata columns if any by mapping groupname
+                # only add the metadata columns from object if there is a single value mapped
+                metacols <- srtobj@meta.data %>% group_by(!!sym(each)) %>%
+                    summarize(across(everything(), ~ n_distinct(.) == 1), .groups = "keep") %>%
+                    select(where(~ all(. == TRUE))) %>%
+                    colnames()
+
+                if (length(metacols) > 1) {
+                    metadf <- srtobj@meta.data[, metacols, drop = FALSE]  %>%
+                        distinct(!!sym(each), .keep_all = TRUE)
+
+                    for (col in setdiff(metacols, each)) {
+                        if (col %in% colnames(enriches)) {
+                            warning("Column name conflict: {col}, adding with suffix '_meta'", immediate. = TRUE)
+                            metadf[[paste0(col, "_meta")]] <- metadf[[col]]
+                            metadf[[col]] <- NULL
+                        }
+                    }
+
+                    enriches <- left_join(enriches, metadf, by = each)
+                }
+
                 process_allenriches(enriches, allenrich_plots, name, each)
             }
         }
@@ -659,6 +681,29 @@ run_case <- function(name) {
 
         if (length(allenrich_plots) > 0) {
             log$info("- Visualizing all enrichments together ...")
+            # add other metadata columns if any by mapping groupname
+            # only add the metadata columns from object if there is a single value mapped
+            metacols <- subobj@meta.data %>% group_by(!!sym(case$group_by)) %>%
+                summarize(across(everything(), ~ n_distinct(.) == 1), .groups = "keep") %>%
+                select(where(~ all(. == TRUE))) %>%
+                colnames()
+
+            if (length(metacols) > 1) {
+                metadf <- subobj@meta.data[, metacols, drop = FALSE]  %>%
+                    distinct(!!sym(case$group_by), .keep_all = TRUE)
+
+                for (col in setdiff(metacols, case$group_by)) {
+                    if (col %in% colnames(enriches[[1]])) {
+                        warning("Column name conflict: {col}, adding with suffix '_meta'", immediate. = TRUE)
+                        metadf[[paste0(col, "_meta")]] <- metadf[[col]]
+                        metadf[[col]] <- NULL
+                    }
+                }
+
+                for (ne in names(enriches)) {
+                    enriches[[ne]] <- left_join(enriches[[ne]], metadf, by = case$group_by)
+                }
+            }
             process_allenriches(enriches, allenrich_plots, name, case$group_by)
         }
     } else {

{biopipen-0.34.14.dist-info → biopipen-0.34.15.dist-info}/METADATA

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

{biopipen-0.34.14.dist-info → biopipen-0.34.15.dist-info}/RECORD

@@ -1,4 +1,4 @@
-biopipen/__init__.py,sha256=EumHAzPCe3UipaK_9iqfQf47_8LP5jWSnXBVUZmzQTA,24
+biopipen/__init__.py,sha256=Rz3xEpDF80HTbA8ImEKPXgpJjUW3trCJqVbTVpBlx44,24
 biopipen/core/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 biopipen/core/config.py,sha256=edK5xnDhM8j27srDzsxubi934NMrglLoKrdcC8qsEPk,1069
 biopipen/core/config.toml,sha256=lZV_vbYWk6uqm19ZWJcsZCcSNqAdIfN2fOfamzxZpg4,2148
@@ -155,7 +155,7 @@ biopipen/scripts/scrna/ExprImputation-rmagic.R,sha256=ePgbMZ_3bKbeUrjsMdkdtBM_MS
 biopipen/scripts/scrna/ExprImputation-scimpute.R,sha256=MI_bYfvCDKJsuGntUxfx_-NdrssBoQgL95-DGwJVE5s,1191
 biopipen/scripts/scrna/ExprImputation.R,sha256=GcdZJpkDpq88hRQjtLZY5-byp8V43stEFm5T-pQbU6A,319
 biopipen/scripts/scrna/LoomTo10X.R,sha256=c6F0p1udsL5UOlb84-53K5BsjSDWkdFyYTt5NQmlIec,1059
-biopipen/scripts/scrna/MarkersFinder.R,sha256=nz18T_70ni9sfoFYp8weQMCiojphF61l5yHe9f_gOSM,26320
+biopipen/scripts/scrna/MarkersFinder.R,sha256=SBh3cA2WZQIlZTysAlU-U9knzO9gU0XW-yi8mk2iJG4,28588
 biopipen/scripts/scrna/MetaMarkers.R,sha256=BgYaWYEj6obwqaZaDWqNPtxb1IEEAnXAeBE0Ji9PvBA,12426
 biopipen/scripts/scrna/ModuleScoreCalculator.R,sha256=_mvo35a-wk5miUb_kMIVwvKK0b6InRa1NKtN8zznGwk,4457
 biopipen/scripts/scrna/PseudoBulkDEG.R,sha256=IuM4hl-KHZ5aaaTqZeylw4b1ZenMZaY4qobD5qxAlHs,25199
@@ -284,7 +284,7 @@ biopipen/utils/misc.py,sha256=pDZ-INWVNqHuXYvcjmu8KqNAigkh2lsHy0BxX44CPvc,4048
 biopipen/utils/reference.py,sha256=Oc6IlA1giLxymAuI7DO-IQLHQ7-DbsWzOQE86oTDfMU,5955
 biopipen/utils/reporter.py,sha256=VwLl6xyVDWnGY7NEXyqBlkW8expKJoNQ5iTyZSELf5c,4922
 biopipen/utils/vcf.py,sha256=MmMbAtLUcKPp02jUdk9TzuET2gWSeoWn7xgoOXFysK0,9393
-biopipen-0.34.14.dist-info/METADATA,sha256=uIS3ZdJr3bVM1TspZIe2NmoVavOa_XGKe4qsuRk7gdw,1027
-biopipen-0.34.14.dist-info/WHEEL,sha256=zp0Cn7JsFoX2ATtOhtaFYIiE2rmFAD4OcMhtUki8W3U,88
-biopipen-0.34.14.dist-info/entry_points.txt,sha256=BYqHGBQJxyFDNLYqgH64ycI5PYwnlqwYcCFsMvJgzAU,653
-biopipen-0.34.14.dist-info/RECORD,,
+biopipen-0.34.15.dist-info/METADATA,sha256=ArJuTPoapr8eXqqyw0QMZkzalGdJLJED0QpmttmNSA0,1027
+biopipen-0.34.15.dist-info/WHEEL,sha256=zp0Cn7JsFoX2ATtOhtaFYIiE2rmFAD4OcMhtUki8W3U,88
+biopipen-0.34.15.dist-info/entry_points.txt,sha256=BYqHGBQJxyFDNLYqgH64ycI5PYwnlqwYcCFsMvJgzAU,653
+biopipen-0.34.15.dist-info/RECORD,,