biopipen 0.23.7__py3-none-any.whl → 0.23.8__py3-none-any.whl

biopipen/__init__.py

@@ -1 +1 @@
-__version__ = "0.23.7"
+__version__ = "0.23.8"

biopipen/ns/scrna.py

@@ -1463,6 +1463,7 @@ class ScFGSEA(Proc):
         ident-1: The first group of cells to compare
         ident-2: The second group of cells to compare, if not provided, the rest of the cells that are not `NA`s in `group-by` column are used for `ident-2`.
         each: The column name in metadata to separate the cells into different subsets to do the analysis.
+        subset: An expression to subset the cells.
         section: The section name for the report. Worked only when `each` is not specified. Otherwise, the section name will be constructed from `each` and its value.
             This allows different cases to be put into the same section in the report.
         gmtfile: The pathways in GMT format, with the gene names/ids in the same format as the seurat object.
@@ -1513,6 +1514,7 @@ class ScFGSEA(Proc):
         "ident-1": None,
         "ident-2": None,
         "each": None,
+        "subset": None,
         "section": "DEFAULT",
         "gmtfile": "",
         "method": "s2n",

biopipen/scripts/scrna/ScFGSEA.R

@@ -14,6 +14,7 @@ group.by <- {{envs["group-by"] | r}}  # nolint
 ident.1 <- {{envs["ident-1"] | r}}  # nolint
 ident.2 <- {{envs["ident-2"] | r}}  # nolint
 each <- {{envs.each | r}}  # nolint
+subset <- {{envs.subset | r}}  # nolint
 section <- {{envs.section | r}}  # nolint
 gmtfile <- {{envs.gmtfile | r}}  # nolint
 method <- {{envs.method | r}}  # nolint
@@ -43,6 +44,7 @@ expand_cases <- function() {
                 ident.1 = ident.1,
                 ident.2 = ident.2,
                 each = each,
+                subset = subset,
                 section = section,
                 gmtfile = gmtfile,
                 method = method,
@@ -63,6 +65,7 @@ expand_cases <- function() {
                 ident.1 = ident.1,
                 ident.2 = ident.2,
                 each = each,
+                subset = subset,
                 section = section,
                 gmtfile = gmtfile,
                 method = method,
@@ -136,6 +139,9 @@ do_case <- function(name, case) {
     # prepare expression matrix
     log_info("  Preparing expression matrix...")
     sobj <- srtobj %>% filter(!is.na(!!sym(case$group.by)))
+    if (!is.null(case$subset)) {
+        sobj <- sobj %>% filter(!!!parse_exprs(case$subset))
+    }
     if (!is.null(case$ident.2)) {
         sobj <- sobj %>% filter(!!sym(case$group.by) %in% c(case$ident.1, case$ident.2))
     }

biopipen/scripts/scrna/SeuratPreparing.R

@@ -99,8 +99,8 @@ load_sample = function(sample) {
     }
     obj <- CreateSeuratObject(exprs, project=sample)
     # filter the cells that don't have any gene expressions
-    cell_exprs = colSums(obj@assays$RNA)
-    obj = subset(obj, cells = names(cell_exprs[cell_exprs > 0]))
+    # cell_exprs = colSums(obj@assays$RNA)
+    # obj = subset(obj, cells = names(cell_exprs[cell_exprs > 0]))
     obj = RenameCells(obj, add.cell.id = sample)
     # Attach meta data
     for (mname in names(mdata)) {
@@ -128,13 +128,7 @@ log_info("Reading samples individually ...")
 obj_list = lapply(samples, load_sample)
 
 log_info("Merging samples ...")
-if (length(obj_list) >= 2) {
-    y = c()
-    for (i in 2:length(obj_list)) y = c(y, obj_list[[i]])
-    sobj = merge(obj_list[[1]], y)
-} else {
-    sobj = obj_list[[1]]
-}
+sobj = Reduce(merge, obj_list)
 
 log_info("Adding metadata for QC ...")
 sobj$percent.mt = PercentageFeatureSet(sobj, pattern = "^MT-")
@@ -297,28 +291,36 @@ add_report(
     h1 = "Filters and QC"
 )
 
+.formatArgs <- function(args) {
+    paste(capture.output(str(args)), collapse = ", ")
+}
+
 log_info("Performing transformation/scaling ...")
 # Not joined yet
 # sobj[["RNA"]] <- split(sobj[["RNA"]], f = sobj$Sample)
 if (envs$use_sct) {
     log_info("- Running SCTransform ...")
     SCTransformArgs <- envs$SCTransform
+    log_info("  SCTransform: {.formatArgs(SCTransformArgs)}")
     SCTransformArgs$object <- sobj
     sobj <- do_call(SCTransform, SCTransformArgs)
     # Default is to use the SCT assay
 } else {
     log_info("- Running NormalizeData ...")
     NormalizeDataArgs <- envs$NormalizeData
+    log_info("  NormalizeData: {.formatArgs(NormalizeDataArgs)}")
     NormalizeDataArgs$object <- sobj
     sobj <- do_call(NormalizeData, NormalizeDataArgs)
 
     log_info("- Running FindVariableFeatures ...")
     FindVariableFeaturesArgs <- envs$FindVariableFeatures
+    log_info("  FindVariableFeatures: {.formatArgs(FindVariableFeaturesArgs)}")
     FindVariableFeaturesArgs$object <- sobj
     sobj <- do_call(FindVariableFeatures, FindVariableFeaturesArgs)
 
     log_info("- Running ScaleData ...")
     ScaleDataArgs <- envs$ScaleData
+    log_info("  ScaleData: {.formatArgs(ScaleDataArgs)}")
     ScaleDataArgs$object <- sobj
     sobj <- do_call(ScaleData, ScaleDataArgs)
 }
@@ -326,13 +328,13 @@ if (envs$use_sct) {
 log_info("- Running RunPCA ...")
 RunPCAArgs <- envs$RunPCA
 RunPCAArgs$npcs <- if (is.null(RunPCAArgs$npcs)) { 50 } else { min(RunPCAArgs$npcs, ncol(sobj) - 1) }
+log_info("  RunPCA: {.formatArgs(RunPCAArgs)}")
 RunPCAArgs$object <- sobj
 sobj <- do_call(RunPCA, RunPCAArgs)
 
 if (!envs$no_integration) {
     log_info("- Running IntegrateLayers ...")
     IntegrateLayersArgs <- envs$IntegrateLayers
-    IntegrateLayersArgs$object <- sobj
     method <- IntegrateLayersArgs$method
     if (!is.null(IntegrateLayersArgs$reference) && is.character(IntegrateLayersArgs$reference)) {
         log_info("  Using reference samples: {paste(IntegrateLayersArgs$reference, collapse = ', ')}")
@@ -359,6 +361,8 @@ if (!envs$no_integration) {
     if (is.null(IntegrateLayersArgs$new.reduction)) {
         IntegrateLayersArgs$new.reduction <- new_reductions[[method]]
     }
+    log_info("  IntegrateLayers: {.formatArgs(IntegrateLayersArgs)}")
+    IntegrateLayersArgs$object <- sobj
     sobj <- do_call(IntegrateLayers, IntegrateLayersArgs)
     # Save it for dimension reduction plots
     sobj@misc$integrated_new_reduction <- IntegrateLayersArgs$new.reduction

{biopipen-0.23.7.dist-info → biopipen-0.23.8.dist-info}/METADATA

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

{biopipen-0.23.7.dist-info → biopipen-0.23.8.dist-info}/RECORD

@@ -1,4 +1,4 @@
-biopipen/__init__.py,sha256=cC_HEsX6d06QmpjNLfRezzJvL501PPwiQYXWDM2_2UI,23
+biopipen/__init__.py,sha256=gx1o8BSQCB5X7Zn8JZ7cxhUiCNJydk8QQUY_fY1zEC8,23
 biopipen/core/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 biopipen/core/config.py,sha256=edK5xnDhM8j27srDzsxubi934NMrglLoKrdcC8qsEPk,1069
 biopipen/core/config.toml,sha256=Rn7Cta7WsMtmQkKGC4h9d5dU_STaIVBgR8UliiGgL6o,1757
@@ -20,7 +20,7 @@ biopipen/ns/gsea.py,sha256=EsNRAPYsagaV2KYgr4Jv0KCnZGqayM209v4yOGGTIOI,7423
 biopipen/ns/misc.py,sha256=fzn0pXvdghMkQhu-e3MMapPNMyO6IAJbtTzVU3GbFa0,3246
 biopipen/ns/plot.py,sha256=yguxmErUOH-hOM10JfuI_sXw2p49XF8yGR_gXfbd5yQ,4066
 biopipen/ns/rnaseq.py,sha256=l4vFeRasGhkexopGTM_VfSyIFewOxg-9L5niFzhWUNA,565
-biopipen/ns/scrna.py,sha256=42cM6n7rNy8sze9Lhl90RNkpxWT5w6LKPDrGsAK_Y7U,95808
+biopipen/ns/scrna.py,sha256=aV5HNM2TJW2LRAhE5ekEm4114rIeHXBio_nwI-bSPpQ,95883
 biopipen/ns/scrna_metabolic_landscape.py,sha256=9s1NvH3aMaNDXyfwy9TdzGcSP_lIW4JqhLgknNZcIKE,28313
 biopipen/ns/tcgamaf.py,sha256=AFbUJIxiMSvsVY3RcHgjRFuMnNh2DG3Mr5slLNEyz6o,1455
 biopipen/ns/tcr.py,sha256=IcP1uD8U9XD6UbOgdjA_Lk5PK6r4R84Gi7511uvXoy8,84411
@@ -124,7 +124,7 @@ biopipen/scripts/scrna/MetaMarkers.R,sha256=kgjk65EmewZ1uh8AoiENDLmNFe4lzmYvssJe
 biopipen/scripts/scrna/ModuleScoreCalculator.R,sha256=JSHd-_-KiFqW8avCGxgU4T-C5BtDr2u0kwIvEu2lFIg,4188
 biopipen/scripts/scrna/RadarPlots.R,sha256=iR4JKtO2b3hGfqv_KAI7BR9tq02EAYfKeqp7tzAicKs,14808
 biopipen/scripts/scrna/SCImpute.R,sha256=dSJOHhmJ3x_72LBRXT72dbCti5oiB85CJ-OjWtqONbk,2958
-biopipen/scripts/scrna/ScFGSEA.R,sha256=2nJN1AAe55SyR-jwVj5ds2pyt32fTIn548zjvEB4Jt0,7332
+biopipen/scripts/scrna/ScFGSEA.R,sha256=E6Rx-0TjplP_nptDCE-LqMyipsOFaMU3hM7FjjFFFDY,7538
 biopipen/scripts/scrna/SeuratClusterStats-dimplots.R,sha256=pZKv1SnSNEGXDeE0_2VYp0GAikYitohW2FR5YGKjs8Q,2351
 biopipen/scripts/scrna/SeuratClusterStats-features.R,sha256=SaKTJloP1fttRXZQeb2ApX0ej7al13wOoEYkthSk13k,15489
 biopipen/scripts/scrna/SeuratClusterStats-hists.R,sha256=GTZfs1yOkuoMUM1Zb19i_My8B8b1Qtve8je55pU_w-g,5054
@@ -136,7 +136,7 @@ biopipen/scripts/scrna/SeuratFilter.R,sha256=BrYK0MLdaTtQvInMaQsmOt7oH_hlks0M1zy
 biopipen/scripts/scrna/SeuratLoading.R,sha256=ekWKnHIqtQb3kHVQiVymAHXXqiUxs6KKefjZKjaykmk,900
 biopipen/scripts/scrna/SeuratMap2Ref.R,sha256=tisYmoSaCX8Kl8y6euuuUroWdDsJ2NGI27J5AWr9Niw,4392
 biopipen/scripts/scrna/SeuratMetadataMutater.R,sha256=Pp4GsF3hZ6ZC2vroC3LSBmVa4B1p2L3hbh981yaAIeQ,1093
-biopipen/scripts/scrna/SeuratPreparing.R,sha256=3X2-yCjNdYVMFNozRU38RxhNBOL-lPo9Za1qUT7V5wI,12227
+biopipen/scripts/scrna/SeuratPreparing.R,sha256=plg8hn6lvVC630iXohTDiuszOuISQM44eQI5kWS3CP8,12565
 biopipen/scripts/scrna/SeuratSplit.R,sha256=vdK11V39_Uo_NaOh76QWCtxObGaEr5Ynxqq0hTiSvsU,754
 biopipen/scripts/scrna/SeuratSubClustering.R,sha256=6b1J98YYK2gwJ_qkpuvYIdlj6uWxJA8IpUXvaCjgC_U,6334
 biopipen/scripts/scrna/SeuratSubset.R,sha256=yVA11NVE2FSSw-DhxQcJRapns0tNNHdyDYi5epO6SKM,1776
@@ -218,7 +218,7 @@ biopipen/utils/reference.py,sha256=6bPSwQa-GiDfr7xLR9a5T64Ey40y24yn3QfQ5wDFZkU,4
 biopipen/utils/rnaseq.R,sha256=Ro2B2dG-Z2oVaT5tkwp9RHBz4dp_RF-JcizlM5GYXFs,1298
 biopipen/utils/single_cell.R,sha256=bKduqOQjSC8BtZJuwfUShR49omoEMbB57n3Gi6dYlqA,4147
 biopipen/utils/vcf.py,sha256=ajXs0M_QghEctlvUlSRjWQIABVF02wPdYd-0LP4mIsU,9377
-biopipen-0.23.7.dist-info/METADATA,sha256=3y-jNYKvsvX_BQuJEnxQ1nhoCvl0L3KwfLI3ll-fHbI,886
-biopipen-0.23.7.dist-info/WHEEL,sha256=FMvqSimYX_P7y0a7UY-_Mc83r5zkBZsCYPm7Lr0Bsq4,88
-biopipen-0.23.7.dist-info/entry_points.txt,sha256=16Apdku3RFwghe1nb0JR7eVo4IzLae6hCWjU1VxYUn0,525
-biopipen-0.23.7.dist-info/RECORD,,
+biopipen-0.23.8.dist-info/METADATA,sha256=NDuFljfwMBKGx9__4fsjm2ZN7yAuW49L04kzmbxblmc,886
+biopipen-0.23.8.dist-info/WHEEL,sha256=FMvqSimYX_P7y0a7UY-_Mc83r5zkBZsCYPm7Lr0Bsq4,88
+biopipen-0.23.8.dist-info/entry_points.txt,sha256=16Apdku3RFwghe1nb0JR7eVo4IzLae6hCWjU1VxYUn0,525
+biopipen-0.23.8.dist-info/RECORD,,