I have a data set where community composition data was collected from 2 sites in a paired design: under the canopy and nearby canopy. We used mvabund package in R to check whether the communities differed. They do differ. The main problem is to find a way to plot the data and show that communities are different, accounting for the paired design. We plotted the data through nMDS, but the communities did not seem to be different at all. I'm trying to use boral package but I haven't found in any example how to control for the pairing (blocks).

  • $\begingroup$ How did you measure "community composition"? $\endgroup$ – AdamO Dec 7 '17 at 13:04
  • $\begingroup$ I have a matrix of species (columns) per samples (rows). Cells are filled with species abundances $\endgroup$ – Raf1987 Dec 7 '17 at 16:36
  • $\begingroup$ Have you considered a simple side-by-side barchart or tornado diagram? $\endgroup$ – AdamO Dec 7 '17 at 16:43
  • $\begingroup$ how would you account for the pairings in such situations? $\endgroup$ – Raf1987 Dec 7 '17 at 19:04

The updated version of BORAL allows for inclusion of site effects, which you can group as needed (replicates, nested, blocked, etc). See examples from https://www.rdocumentation.org/packages/boral/versions/1.7/topics/boral.

You can use results from BORAL in ggplot to create a figure similar to NMDS (using latent variables as axes), to show community differences.


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.