nestedness(m, null.models = TRUE, n.nulls = 100, popsize = 30,
n.ind = 7, n.gen = 2000, binmatnestout=FALSE)
m
is the matrix object for which the temperature is calculated. m
will be converted to a binary matrix as temperature is only based on binary datanull.models
=TRUE.n.nulls=100
.popsize
, default is 30n.ind
is 7.n.gen
is 2000.Here we use BINMATNEST by calling a tweaked version of the C++ program binmatnest. In principle nestedness temperature is calculated by using a line of perfect order (using a genetic algorithm) to determine the reordering of rows and columns that leads to minimum matrix temperature of given size and fills. The deviation from this minimun temperature is the matrix temperature. In addition nestedness uses different null models to check for statistical significance of the matrix temperature. For details on what BINMATNEST does different, and better, than the original NTC see reference below.
Notice also that the original software BINMATNEST is available as a stand-alone application, too. Check out Miguel's homepage:
data(vazarr)
nestedness(vazarr) # null models are calculated
# no null models, much faster for bigger matrices:
nestedness(vazarr, null.models=FALSE)
nestedness(vazarr, n.nulls=300, n.gen=300)
Run the code above in your browser using DataLab