Support for "rank ~ factor" mode

cmmr / rbiom

Interact with Biological Observation Matrix files.

https://cmmr.github.io/rbiom/

Other

12 stars 0 forks source link

Support for "rank ~ factor" mode #5

Closed dansmith01 closed 11 months ago

dansmith01 commented 2 years ago

Add ".taxa" to the facet array, so it'll look like the image below, but with stat bars between each pair.

plot(
 rarefy(hmp50, 1000), "Phylum", "bd", p.top = 5, taxa = 0.02,
 color.by = "Body Site", facet.by = ".taxa",
 trans = "sqrt", scales = "free", d.alpha = 0.2 )

snip

plot(rarefy(hmp50, 1000), Phylum ~ `Body Site`)

Error in plot.BIOM(biom, Phylum ~ `Body Site`) : 
Invalid formula of form 'rank ~ factor'

dansmith01 commented 11 months ago

Superseded by taxa_barplot(), which supports pairwise stats.