niche.test: Monte-Carlo Test on Parameters of the Ecological Niche
In adehabitat: Analysis of Habitat Selection by Animals

Description Usage Arguments Details Value Warning Author(s) See Also Examples

niche.test tests for the significance of two parameters of the ecological niche of a species (marginality and tolerance), using Monte-Carlo methods. This is a bivariate test.

1	niche.test(kasc, points, nrep = 999, o.include = TRUE, ...)

`kasc`	a raster map of class `kasc`
`points`	a data frame with two columns, giving the coordinates of the species locations
`nrep`	the number of permutations
`o.include`	logical, passed to `biv.test`. If `TRUE`, the origin is included in the plot
`...`	further arguments passed to `biv.test`

niche.test tests the significance of two parameters describing the ecological niche: the marginality (squared length of the vector linking the average available habitat conditions to the average used habitat conditions in the ecological space defined by the habitat variables), and the tolerance (inertia of the niche in the ecological space, i.e. the sum over all variables of the variance of used pixels).

At each step of the randomisation procedure, the test randomly allocates the n points in the pixels of the map. The marginality and the tolerance are then recomputed on this randomised data set.

Actual values are compared to random values with the help of the function biv.test.

Returns a list containing the following components:

`dfxy`	a data frame with the randomized values of marginality (first column) and tolerance (second column).
`obs`	the actual value of marginality and tolerance.

biv.test uses the function kde2d of the package MASS.

Mathieu Basille basille@ase-research.org
Clement Calenge clement.calenge@oncfs.gouv.fr

biv.test for more details on bivariate tests. histniche for the histograms of the variables of the niche.

## Not run: 
data(lynxjura)

## We keep only "wild" indices.
tmp=lynxjura$loc[,4]!="D"
niche=niche.test(lynxjura$map,
                 lynxjura$locs[tmp, c("X", "Y")],
                 side = "bottom")
names(niche)

## End(Not run)