View source: R/ecospat.boyce.R

ecospat.boyce | R Documentation |

Calculate the Boyce index as in Hirzel et al. (2006). The Boyce index is used to assess model performance.

```
ecospat.boyce (fit, obs, nclass=0, window.w="default", res=100, PEplot = TRUE,
rm.duplicate = TRUE, method = 'spearman')
```

`fit` |
A vector or a SpatRaster containing the predicted suitability values |

`obs` |
A vector containing the predicted suitability values or xy-coordinates (if "fit" is a SpatRaster) of the validation points (presence records) |

`nclass` |
The number of classes or vector with class thresholds. If |

`window.w` |
The width of the moving window (by default 1/10 of the suitability range) |

`res` |
The resolution of the moving window (by default 100 focals) |

`PEplot` |
If |

`rm.duplicate` |
If |

`method` |
Method used to compute the correlation. The default is |

rm.duplicate = TRUE, method = 'spearman'

The Boyce index only requires presences and measures how much model predictions differ from random distribution of the observed presences across the prediction gradients (Boyce et al. 2002). It is thus an appropriate metric in the case of presence-only models. It is continuous and varies between -1 and +1. Positive values indicate a model which present predictions are consistent with the distribution of presences in the evaluation dataset, values close to zero mean that the model is not different from a random model, negative values indicate counter predictions, i.e., predicting poor quality areas where presences are more frequent (Hirzel et al. 2006).

Removing the successive duplicated P/E values (`rm.duplicate = TRUE`

) focuses more on the discriminative aspect of the predictions, lowering the assessment of the evaluation of the model resolution (sensu Hirzel et al. 2006 p. 150). However, it seems that in the initial version, dupplicated values were not removed.

In the initial publication on the continuous Boyce index, the correlation was set to `method = 'spearman'`

. However, using `method = 'kendall'`

or `method = 'pearson'`

might be more informative about the accuracy of the predictions.

The function returns a list that contains a vector F.ratio (the predicted-to-expected ratio for each class-interval) and a numeric Spearman.cor (the Boyce index value)

Blaise Petitpierre bpetitpierre@gmail.com and Frank Breiner frank.breiner@unil.ch with the updates of Flavien Collart

Boyce, M.S., P.R. Vernier, S.E. Nielsen and F.K.A. Schmiegelow. 2002. Evaluating resource selection functions. *Ecol. Model.*, **157**, 281-300.

Hirzel, A.H., G. Le Lay, V. Helfer, C. Randin and A. Guisan. 2006. Evaluating the ability of habitat suitability models to predict species presences. *Ecol. Model.*, **199**, 142-152.

```
obs <- (ecospat.testData$glm_Saxifraga_oppositifolia
[which(ecospat.testData$Saxifraga_oppositifolia==1)])
ecospat.boyce (fit = ecospat.testData$glm_Saxifraga_oppositifolia , obs, nclass=0,
window.w="default", res=100, PEplot = TRUE)
```

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.