nestedchecker(comm)
nestedn0(comm)
nesteddisc(comm)
nestedtemp(comm, ...)
nestednodf(comm, order = TRUE, weighted = FALSE)
## S3 method for class 'nestedtemp':
plot(x, kind = c("temperature", "incidence"),
col=rev(heat.colors(100)), names = FALSE, ...)
plot
.comm
.
If it is a logical vector of length 2, row and column labels are
returned accordingly.statistic
, but the other components differ among functions. The
functions are constructed so that they can be handled by
oecosimu
.oecosimu
to analyse
the non-randomness of results.
Function nestedchecker
gives the number of checkerboard units,
or 2x2 submatrices where both species occur once but on different
sites (Stone & Roberts 1990). Function nestedn0
implements
nestedness measure N0 which is the number of absences from the sites
which are richer than the most pauperate site species occurs
(Patterson & Atmar 1986).
Function nesteddisc
implements
discrepancy index which is the number of ones that should be shifted
to fill a row with ones in a table arranged by species frequencies
(Brualdi & Sanderson 1999). The original definition arranges species
(columns) by their frequencies, but did not have any method of
handling tied frequencies.
The nesteddisc
function tries to
order tied columns to minimize the discrepancy statistic but this is
rather slow, and with a large number of tied columns there is no
guarantee that the best ordering was found. In that case a warning
of tied columns will be issued.
Function nestedtemp
finds the matrix temperature which is
defined as the sum of nestedtemp
also
has a plot
method which can display either incidences or
temperatures of the surprises. Matrix temperature was rather vaguely
described (Atmar & Patterson 1993), but
vignette
Design decisions and
implementation that you can read using functions
vignette
or vegandocs
. Function
nestedness
in the
Function nestednodf
implements a nestedness metric based on
overlap and decreasing fill (Almeida-Neto et al., 2008). Two basic
properties are required for a matrix to have the maximum degree of
nestedness according to this metric: (1) complete overlap of 1's
from right to left columns and from down to up rows, and (2)
decreasing marginal totals between all pairs of columns and all
pairs of rows. The nestedness statistic is evaluated separately for
columns (N columns
) for rows (N rows
) and combined for
the whole matrix (NODF
). If you set order = FALSE
,
the statistic is evaluated with the current matrix ordering allowing
tests of other meaningful hypothesis of matrix structure than
default ordering by row and column totals (breaking ties by total
abundances when weighted = TRUE
) (see Almeida-Neto et
al. 2008). With weighted = TRUE
, the function finds the
weighted version of the index (Almeida-Neto & Ulrich,
2011). However, this requires quantitative null models for adequate
testing.
Almeida-Neto, M. & Ulrich, W. (2011). A straightforward computational approach for measuring nestedness using quantitative matrices. Env. Mod. Software 26, 173--178. Atmar, W. & Patterson, B.D. (1993). The measurement of order and disorder in the distribution of species in fragmented habitat. Oecologia 96, 373--382.
Brualdi, R.A. & Sanderson, J.G. (1999). Nested species subsets, gaps, and discrepancy. Oecologia 119, 256--264.
Patterson, B.D. & Atmar, W. (1986). Nested subsets and the structure of insular mammalian faunas and archipelagos. Biol. J. Linnean Soc. 28, 65--82.
Stone, L. & Roberts, A. (1990). The checkerboard score and species distributions. Oecologia 85, 74--79.
Wright, D.H., Patterson, B.D., Mikkelson, G.M., Cutler, A. & Atmar, W. (1998). A comparative analysis of nested subset patterns of species composition. Oecologia 113, 1--20.
oecosimu
which generates Null model communities to assess the non-randomness of
nestedness patterns.data(sipoo)
## Matrix temperature
out <- nestedtemp(sipoo)
out
plot(out)
plot(out, kind="incid")
## Use oecosimu to assess the non-randomness of checker board units
nestedchecker(sipoo)
oecosimu(sipoo, nestedchecker, "quasiswap")
## Another Null model and standardized checkerboard score
oecosimu(sipoo, nestedchecker, "r00", statistic = "C.score")
Run the code above in your browser using DataLab