deseq2 multicomp

maxplanck-ie · Apr 10, 2024 · 1280216 · 1280216
1 parent 4385081
commit 1280216
Showing 1 changed file with 5 additions and 15 deletions.
diff --git a/snakePipes/shared/rules/DESeq2.multipleComp.snakefile b/snakePipes/shared/rules/DESeq2.multipleComp.snakefile
@@ -65,20 +65,10 @@ rule DESeq2_Salmon_basic:
         fdr = fdr,
         importfunc = os.path.join(maindir, "shared", "rscripts", "DE_functions.R"),
         allele_info = 'FALSE',
-        tx2gene_file = "Annotation/genes.filtered.t2g",
+        tx2gene_file = os.path.join(outdir,"Annotation/genes.filtered.t2g"),
         rmdTemplate = os.path.join(maindir, "shared", "rscripts", "DESeq2Report.Rmd"),
-        formula = config["formula"]
+        formula = config["formula"],
+        counts_table = lambda wildcards,input: os.path.join(outdir,input.counts_table),
+        symbol_file = lambda wildcards,input: os.path.join(outdir,input.symbol_file)
     conda: CONDA_RNASEQ_ENV
-    shell:
-        "cd {params.outdir} && "
-        "Rscript {params.script} "
-        "{params.sampleSheet} " # 1
-        "../{input.counts_table} " # 2
-        "{params.fdr} " # 3
-        "../{input.symbol_file} " # 4
-        "{params.importfunc} " # 5
-        "{params.allele_info} " # 6
-        "../{input.tx2gene_file} " # 7
-        "{params.rmdTemplate} " # 8
-        "{params.formula} "
-        " > ../{log.out} 2> ../{log.err}"
+    script: "{params.script}"