The ecostats group are happy to introduce a new R package called boral -- Bayesian Ordination and Regression AnaLysis, for analysis of multivariate data (community composition data especially) in ecology!!!
Boral uses Bayesian MCMC estimation via JAGS (Just Another Gibbs Sampler) to fit three types of models:
1) GLMs fitted independently to each species (like in another R package mvabund, developed by us)
2) Purely latent variable models for model-based unconstrained ordination (see Hui et al., 2014, http://onlinelibrary.wiley.com/doi/10.1111/2041-210X.12236/abstract for details)
3) GLMs fitted to each species while accounting for correlation between species e.g., due to species interaction.
Check it out at: http://cran.r-project.org/web/packages/boral/index.html
And even better, check out the video (which we promise has no dancing and no singing unlike another certain video...): https://www.youtube.com/watch?v=vyMsgyytcUI