| domain {adehabitat} | R Documentation |
domain uses the DOMAIN algorithm to estimate the potential
distribution of a species based on a list of species occurrences and on
maps of the area.
domain(kasc, pts, type = c("value", "potential"), thresh = 0.95)
kasc |
an object of class kasc |
pts |
a data frame giving the x and y coordinates of the species occurrences. |
type |
a character string. The "value" of the suitability
may be returned
or the "potential" area of distribution |
thresh |
if value = "potential", a threshold value should be
supplied for the suitability (by default 0.95) |
This function implements the DOMAIN algorithm described in Carpenter et al. (1993).
Returns a matrix of class asc.
domain is restricted to maps
containing only numerical variables (i.e. no factors).
Clement Calenge clement.calenge@oncfs.gouv.fr
Carpenter, G., Gillison, A.N. and Winter, J. (1993) DOMAIN: a flexible modelling procedure for mapping potential distributions of plants and animals. Biodiversity and conservation, 2, 667–680.
kasc for additionnal information on objects
of class kasc, asc for additionnal information on
matrices of class asc.
## Preparation of the data data(puechabon) kasc <- puechabon$kasc kasc$Aspect <- NULL pts <- puechabon$locs[puechabon$locs$Name == "Brock", 4:5] ## View of the data elevation <- getkasc(kasc, "Elevation") image(elevation) points(pts, col = "red", pch = 16) ## Estimation of habitat suitability map hsm <- domain(kasc, pts) image(hsm, col = grey((1:256)/256)) contour(hsm, add = TRUE) ## Lighter areas are the most preferred areas ## Potential distribution hsm <- domain(kasc, pts, type = "potential") image(elevation, main = "Habitat suitability map") image(hsm, add = TRUE, col = "orange") points(pts, col = "red", pch = 16)