add test case for new zero intercept models

Author: atrigila

tests/assets/zero_intercept_contrasts.yml

@@ -0,0 +1,23 @@
+contrasts:
+  - id: treatment_effect_classic_comparison
+    comparison: ["treatment", "Control", "Treated"]
+
+  - id: treatment_plus_genotype
+    formula: "~ treatment + genotype"
+    make_contrasts_str: "treatmentTreated"
+
+  - id: interaction_genotype_treatment
+    formula: "~ genotype * treatment"
+    make_contrasts_str: "genotypeWT.treatmentTreated"
+
+  - id: full_model_with_interactions
+    formula: "~ genotype * treatment * time"
+    make_contrasts_str: "genotypeWT.treatmentTreated.time"
+
+  - id: zero_intercept_model_1
+    formula: "~ 0 + genotype"
+    make_contrasts_str: "genotypeWT"
+
+  - id: zero_intercept_model_2
+    formula: "~ genotype - 1"
+    make_contrasts_str: "genotypeWT"

tests/test_rnaseq_dream.nf.test

@@ -67,4 +67,33 @@ nextflow_pipeline {
         }
     }
 
+    test("Test rnaseq dream profile - complex contrast - with zero intercept") {
+        when {
+            params {
+                outdir        = "$outputDir"
+                contrasts_yml = "$projectDir/tests/assets/zero_intercept_contrasts.yml"
+            }
+        }
+
+        then {
+            // stable_name: All files + folders in ${params.outdir}/ with a stable name
+            def stable_name = getAllFilesFromDir(params.outdir, relative: true, includeDir: true, ignore: ['pipeline_info/*.{html,json,txt}'])
+            // stable_path: All files in ${params.outdir}/ with stable content
+            def stable_path = getAllFilesFromDir(params.outdir, ignoreFile: 'tests/.nftignore')
+            assertAll(
+                { assert workflow.success},
+                { assert snapshot(
+                    // Number of successful tasks
+                    workflow.trace.succeeded().size(),
+                    // pipeline versions.yml file for multiqc from which Nextflow version is removed because we tests pipelines on multiple Nextflow versions
+                    removeNextflowVersion("$outputDir/pipeline_info/nf_core_differentialabundance_software_versions.yml"),
+                    // All stable path name, with a relative path
+                    stable_name,
+                    // All files with stable contents
+                    stable_path
+                ).match() }
+                )
+        }
+    }
+
 }

tests/test_rnaseq_dream.nf.test.snap

@@ -98,5 +98,59 @@
             "nextflow": "25.04.6"
         },
         "timestamp": "2025-07-30T16:32:49.712383975"
+    },
+    "Test rnaseq dream profile - complex contrast - with zero intercept": {
+        "content": [
+            14,
+            {
+                "CUSTOM_FILTERDIFFERENTIALTABLE": {
+                    "pandas": "1.5.2"
+                },
+                "GTF_TO_TABLE": {
+                    "atlas-gene-annotation-manipulation": "1.1.1"
+                },
+                "GUNZIP_GTF": {
+                    "gunzip": 1.1
+                },
+                "VALIDATOR": {
+                    "r-shinyngs": "2.2.4"
+                },
+                "VARIANCEPARTITION_DREAM": {
+                    "r-base": "4.3.3",
+                    "bioconductor-edger": "4.0.16",
+                    "variancePartition": "1.32.2"
+                },
+                "Workflow": {
+                    "nf-core/differentialabundance": "v1.6.0dev"
+                }
+            },
+            [
+                "pipeline_info",
+                "pipeline_info/collated_versions.yml",
+                "pipeline_info/nf_core_differentialabundance_software_versions.yml",
+                "tables",
+                "tables/annotation",
+                "tables/annotation/SRP254919.anno.tsv",
+                "tables/differential",
+                "tables/differential/full_model_with_interactions_SRP254919.dream.results.tsv",
+                "tables/differential/full_model_with_interactions_SRP254919.dream.results_filtered.tsv",
+                "tables/differential/interaction_genotype_treatment_SRP254919.dream.results.tsv",
+                "tables/differential/interaction_genotype_treatment_SRP254919.dream.results_filtered.tsv",
+                "tables/differential/treatment_plus_genotype_SRP254919.dream.results.tsv",
+                "tables/differential/treatment_plus_genotype_SRP254919.dream.results_filtered.tsv",
+                "tables/differential/zero_intercept_model_1_SRP254919.dream.results.tsv",
+                "tables/differential/zero_intercept_model_1_SRP254919.dream.results_filtered.tsv",
+                "tables/differential/zero_intercept_model_2_SRP254919.dream.results.tsv",
+                "tables/differential/zero_intercept_model_2_SRP254919.dream.results_filtered.tsv"
+            ],
+            [
+                "SRP254919.anno.tsv:md5,c1d7f21e64bd00f845ec6545c123a1fb"
+            ]
+        ],
+        "meta": {
+            "nf-test": "0.9.2",
+            "nextflow": "25.04.7"
+        },
+        "timestamp": "2025-10-01T17:12:51.093884289"
     }
 }