This function computes a correlogram from taxonomic levels.
correlogram.formula(formula, data = NULL, use = "all.obs")a formula of the type y1+..+yn ~ g1/../gn, where
the y's are the data to analyse and the g's are the
taxonomic levels.
a data frame containing the variables specified in the
formula. If NULL, the variables are sought in the user's
workspace.
a character string specifying how to handle missing
values (i.e., NA). This must be one of "all.obs",
"complete.obs", or "pairwise.complete.obs", or any unambiguous
abbrevation of these. In the first case, the presence of missing
values produces an error. In the second case, all rows with missing
values will be removed before computation. In the last case, missing
values are removed on a case-by-case basis.
An object of class correlogram which is a data frame with three
columns:
the computed Moran's I
the corresponding P-values
the names of each level
or an object of class correlogramList containing a list of objects of class correlogram if several variables are given as response in formula.
See the vignette in R: vignette("MoranI").
# NOT RUN {
data(carnivora)
### Using the formula interface:
co <- correlogram.formula(SW ~ Order/SuperFamily/Family/Genus,
data=carnivora)
co
plot(co)
### Several correlograms on the same plot:
cos <- correlogram.formula(SW + FW ~ Order/SuperFamily/Family/Genus,
data=carnivora)
cos
plot(cos)
# }