Got to the point of testing the model on communities. Got an error in rmultinom, too few positive probabilities. This must be because none of the consumer species can colonize a resource. Will fix now.