sre: Surface Range Envelope
In BIOMOD: Ensemble platform for species distribution modeling

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

Run a rectilinear surface range envelop (equivalent to BIOCLIM) using the extreme percentiles as recommended by Nix or Busby. The SRE performs a simple analysis of within which range of each variable the data is recorded and renders predictions.

1	sre(Response = NULL, Explanatory = NULL, NewData = NULL, Quant=0.025)

`Response`	a vector giving your species' presences and absences data
`Explanatory`	a matrix containing the environmental variables for the sites given in Response. It must have as many rows as there are elements in Response.
`NewData`	The data for which you want to render predictions with the sre. It must be a matrix of the same type as the one given in Explanatory and with precisely the same variable names.
`Quant`	the value defines the most extreme values for each variable not to be taken into account for determining the tolerance boundaries for the considered species.

The more variables you put in, the more restrictive your model will be (if non-colinear variables).

This method is very much influenced by the data input, and more specifically by the extremes. Where a linear model can discriminate the extreme values from the main tendency, the SRE considers it equal as any other data point which leads to notable changes in predictions. Note that, as a consequence of its functionning, the predictions are directly given in binary, a site being either potentially suitable for all the variables, either out of bounds for at least one variable and therefore considered unsuitable.

The quants argument determines the threshold at which the data will be taken into account for calibration : the default of 0.05 induces that the 5% most extreme values will be avoided for each variable on each side of its distribution along the gradient. So it in fact takes 5% away at each end of the variables distribution, giving a total of 10% of data not considered.

A vector of the same length as there are rows in NewData giving the prediction in binary (1=presence, 0=absence)

Wilfried Thuiller, Bruno Lafourcade

Models

data(Sp.Env)
data(CoorXY)

g <- sre(Sp.Env[,11], Sp.Env[,4:10], NewData=Sp.Env[,4:10])
gg <- sre(Sp.Env[,11], Sp.Env[,4:10], NewData=Sp.Env[,4:10], Quant=0.025)
ggg <- sre(Sp.Env[,11], Sp.Env[,4:10], NewData=Sp.Env[,4:10], Quant=0.05)


par(mfrow=c(1,4))
level.plot(Sp.Env[,11], XY=CoorXY[,1:2], show.scale=FALSE, title="original data")
level.plot(g, XY=CoorXY[,1:2], show.scale=FALSE, title="full data calibration")
level.plot(gg, XY=CoorXY[,1:2], show.scale=FALSE, title="Perc025=TRUE")
level.plot(ggg, XY=CoorXY[,1:2], show.scale=FALSE, title="Perc05=TRUE")