rnadeseq
rnadeseq copied to clipboard
Contrast that is not included in design throws error
Description of the bug
When running rnadeseq with a contrast list that contains a factor not described in the design file, the following error is thrown:
Error in cleanContrast(object, contrast, expanded = isExpanded, listValues = listValues, : as 16H is the reference level, was expecting condition_time_1H_vs_16H to be present in 'resultsNames(object)'
This example error resulted from the test dataset with the design and contrast list attached to this issue. I should add a check that controls for all the contrast options, not just list, if the respective contrasts are also provided as design elements, and if not, throws a clearer error to explain to users what is the problem
Command used and terminal output
No response
Relevant files
System information
No response