Randomizes a binary matrix speciesData by reshuffling elements
within each row. Sampling weights for each column are proportional to
column sums. Makes a call to the vector_sample function.
Returns a binary presence-absence matrix with the same
dimensions and rowsums as the input matrix.
Details
This algorithm preserves differences among species in occurrence
frequencies, but assumes differences among sites in suitability are
proportional to observed species richness (= column sums).
References
Gotelli, N.J. 2000. Null model analysis of species co-occurrence
patterns. Ecology 81: 2606-2621.