add: implement tool to add rank names to phyloseq object #6625
Conversation
| <output name="output" ftype="phyloseq"> | ||
| <!-- Check if the output contains the first character --> | ||
| <assert_contents> | ||
| <has_text text="B"/> |
There was a problem hiding this comment.
Can you please be a bit more strict in this test?
change version and add profile
| print(tax_table(physeq)) | ||
|
|
||
| # Zuordnung der Rangnamen | ||
| rank_names <- c("Kingdom", "Phylum", "Class", "Order", "Family", "Genus", "Species") # Anpassen je nach Bedarf |
There was a problem hiding this comment.
There could be more ranks, e.g.: super kingdom, subspecies, strain.
Could you make this more generic, e.g.: the use provides a comma separated list of ranks.
You could make "Kingdom", "Phylum", "Class", "Order", "Family", "Genus", "Species" the default.
If the number of ranks does not match the column numbers of the tax_table the tool fails.
| # Zuordnung der Rangnamen | ||
| rank_names <- c("Kingdom", "Phylum", "Class", "Order", "Family", "Genus", "Species") # Anpassen je nach Bedarf | ||
|
|
||
| # Füge eine leere Spalte für Species hinzu, falls sie fehlt |
There was a problem hiding this comment.
I think not needed, if the logic I proposed before is used
| print(tax_table(physeq)) | ||
|
|
||
| # Extrahiere das erste Zeichen aus dem ersten Eintrag des tax_table (z.B. Kingdom) | ||
| first_char <- substr(tax_table(physeq)[1, 1], 1, 1) |
| <tests> | ||
| <test> | ||
| <param name="input" value="output.phyloseq" ftype="phyloseq"/> | ||
| <output name="output" ftype="phyloseq"> |
There was a problem hiding this comment.
you could compare the new phylosq object to expected output, like:
<output name="ampvis" value="AalborgWWTPs-subset_samples.rds" ftype="ampvis2" compare="sim_size"/>
There was a problem hiding this comment.
Better just use an has_size assertion and remove value="AalborgWWTPs-subset_samples.rds".
| cat("Current tax_table:\n") | ||
| print(tax_table(physeq)) | ||
|
|
||
| # Zuordnung der Rangnamen |
There was a problem hiding this comment.
Guess English comments would be better.
| <tests> | ||
| <test> | ||
| <param name="input" value="output.phyloseq" ftype="phyloseq"/> | ||
| <output name="output" ftype="phyloseq"> |
There was a problem hiding this comment.
Better just use an has_size assertion and remove value="AalborgWWTPs-subset_samples.rds".
FOR CONTRIBUTOR: