Abstract
A neutral model of community dynamics has been built using the hierarchical Bayesian framework and fitted with Markov Chain Monte Carlo methods to three community datasets. To fit the data well, the model would need parameter values that are impossible. This suggests that variation in species abundances cannot be explained solely by random drift between species as suggested by the neutral model.