BIMSBbioinfo / pigx_rnaseq

Bulk RNA-seq Data Processing, Quality Control, and Downstream Analysis Pipeline
GNU General Public License v3.0
20 stars 11 forks source link

suggesting check of covariates against sample sheet column names #122

Closed smoe closed 2 years ago

smoe commented 2 years ago

Hello,

label: run_deseq2
design formula:~ sex + age + AnalysisGroup
Quitting from lines 148-192 (deseqReport.Rmd) 
Error in DESeqDataSet(se, design = design, ignoreRank) : 
  all variables in design formula must be columns in colData
Calls: runReport ... withVisible -> eval -> eval -> <Anonymous> -> DESeqDataSet

Execution halted

In my case, the case sensitivity triggered the error. I know, it is trivial, but I find that this is not so very easy to find in snakemake and the extra check would be nice so one can start the job and leave it for the night with a bit more of a confidence.

Many thanks! Steffen