can base_level be multiple if has many groups?

[snakemake-workflows]

can base_level be multiple if has many groups?

[dlaehnemann]

The base_level: is only meant for calculating fold changes in expression. And this usually only makes sense against some control condition, like a set of untreated samples in the example. Thus, you can only have one base_level per model that you compare against. You could hack around this, by specifying the same model twice, only differing in the base_level that you specify. But you then have to double-check that your comparisons both make sense.

I hope that helps?

That said, maybe you want to do something else, e.g. do multiple group comparisons. For that, you can have any number of levels in the primary_variable: column (condition in the example). So you could for example have values genotype_X, genotype_Y, and genotype_Z in the respective column of your samples.tsv.

And if you want to understand more details about the differential expression testing, you could check out the sleuth documentation:
https://pachterlab.github.io/sleuth/walkthroughs

[whiteorchid]

Thanks a lot! Much appreciate! Best!