adonis, betadiver, and betadisp

HI @devonorourke,

Thanks for tagging me in this discussion! This is an area I'm getting into and getting excited about; the referenced post was awesome.

I think you're looking for two separate functionalities: statistical test and visualization. I think they're closely linked, but slightly different things.

As far as I understand, permadisp is the statistical test that measures whether the dispersion of a group of samples in ordination space is different between groups (essentially a measure of variance). I usually use the scikit-bio implementation, but the result is essentially the same. The permadisp function is statisitcal test that gives a psuedo-F and permutative p-value for whether or not there's dispersion between the two groups or more groups.

Visualization is a second issue and is maybe more subjective, but that's the plot. QIIME (again AFAIK, Im not a developer and do most of my own analysis in python because I like control) doesn't have a way to visually show the centroids for clusters. Personally, I think I'd rather see the centroid as a cloud with a variation limit rather than the center point (like transparent, or something).

I think the steps to getting there for an initial prototype might be the following:

  1. Calculate a median position along each PC using something like beta rarefaction
  2. Somehow merge the median PCoA with the main PCoA using something like bi-plots
  3. Coerce the shiny new biplots into emperor
  4. ????
  5. Profit

My list of steps is probably a vast oversimplification and I'm missing something obvious, but I think that would be my first pass to try it. Probably in the underlying python architecture. Although @yoshiki can probably give limitations of emperor better than I can.

I don't have the bandwidth right now to play around this. Maybe over the summer when everyone here goes off on vacation :airplane: :palm_tree:, Ive got a confusing number of hours of daylight, and theoretically more free time. (Although that feels like a dangerous promise.)

Best,
Justine

4 Likes