domain              package:adehabitat              R Documentation

_E_s_t_i_m_a_t_i_o_n _o_f _t_h_e _P_o_t_e_n_t_i_a_l _D_i_s_t_r_i_b_u_t_i_o_n _o_f _a _S_p_e_c_i_e_s

_D_e_s_c_r_i_p_t_i_o_n:

     '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.

_U_s_a_g_e:

     domain(kasc, pts, type = c("value", "potential"), thresh = 0.95)

_A_r_g_u_m_e_n_t_s:

    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)

_D_e_t_a_i_l_s:

     This function implements the DOMAIN algorithm described in
     Carpenter et al. (1993).

_V_a_l_u_e:

     Returns a matrix of class 'asc'.

_W_a_r_n_i_n_g:

     'domain' is restricted to maps containing only numerical variables
     (i.e. no factors).

_A_u_t_h_o_r(_s):

     Clment Calenge calenge@biomserv.univ-lyon1.fr

_R_e_f_e_r_e_n_c_e_s:

     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.

_S_e_e _A_l_s_o:

     'as.kasc' for additionnal information on objects of class 'kasc',
     'import.asc' for information on matrices of class 'asc'.

_E_x_a_m_p_l_e_s:

     ## 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)

