During the version update of Maaslin2 to 1.16.0 for Galaxy platform, I tried passing two reference factors ( for site and diagnosis) separated by a semi-colon (as mentioned in the tutorials) for a test case where no fixed_effects were provided ( default : ALL). But I keep on getting the following error.
Job in error state.. tool_id: maaslin2, exit_code: 1, stderr: Error in relevel.factor(metadata[, i], ref = ref) :
'ref' must be an existing level
Calls: Maaslin2 -> relevel -> relevel.factor
Execution halted
If I just give reference for site, it gives error asking for reference of diagnosis factor. Please let me know how I can resolve this issue. Also, is there a way to just skip/disable reference altogether, I mean if the user does not give any reference then instead of throwing an error , it should just accept default alphabetical order as was the case in previous versions. But unfortunately it seems that it is mandatory to pass a reference. Would love any suggestion on this.
Please find attached the metadata and feature file I am testing on .
HMP2_taxonomy.tsv (487.4 KB)
HMP2_metadata.tsv (106.3 KB)