Resource Selection Functions (RSF) with `amt`

knitr::opts_chunk$set(warning = FALSE, message = FALSE)

About

This vignette briefly introduces how one can fit a Resource-Selection Function (RSF) with the amt package. We will be using the example data of one red deer from northern Germany and one covariate: a forest cover map.

Getting the data ready

First we load the required libraries and the relocation data (called deer)

library(amt)
data("deer")
deer

Next, we have to get the environmental covariates. A forest layer is included in the package. Note, that this a regular SpatRast.

sh_forest <- get_sh_forest()

Prepare Data for RSF

Random Points

Before fitting a RSF we have to do some data preparation. We have to generate random points, points that we think the animal could have used. The random points define the availability domain. In amt the function random_points is designed to do just that. The function can be used in 3 different ways, depending to the type of object that is passed to the function call.

  1. A track_* (such as the deer object) can be passed to the function random_points. The function then calculates a home range (the home-range estimator can be controlled with argument hr). Within this home range n random points are generated. The default value of n is ten times the number of present points.
  2. If a hr-object (i.e., the result of a home-range estimation in amt) is passed to random_points, points are generated within the home range. This allows to generate random points within any home range that was previously estimated in amt. Note, that this could be a home range of multiple animals. In this case, the function random_points has one additional argument called presence. This argument takes a trk_* with the presence points and adds these points for convenience to the random points.
  3. A SpatialPolygons*-object or sf-object. The latter must contain POLYGONs or MULTIPOLYGONs as features. This can be useful in situation where a home range needs to be buffered, or when other geographical features are considered as the availability domain. As before, this method for random_points also takes the argument presence to optionally add the observed points to the output.

Lets now illustrate the three different situations. First we take random points from a track_xy

r1 <- random_points(deer)
plot(r1)

With the argument n we can control the number of random points (remember that the default is ten times as many points as we observed points).

r1 <- random_points(deer, n = 100)
plot(r1)

Next, we can create random point within a home range, that we estimated before.

hr <- hr_mcp(deer)
r1 <- random_points(hr, n = 500)
plot(r1)

Here, we can also add the observed points:

hr <- hr_mcp(deer)
r1 <- random_points(hr, n = 500, presence = deer)
plot(r1)

Finally, we can work with the home range and for example a buffer and then generate random points within the this new polygon.

hr <- hr_mcp(deer) |> hr_isopleths() |> 
  sf::st_buffer(dist =3e4) # add a 30km buffer
r1 <- random_points(hr, n = 500)
plot(r1)

And we can also add the observed points.

hr <- hr_mcp(deer) |> hr_isopleths() |> 
  sf::st_buffer(dist =3e4) # add a 30km buffer
r1 <- random_points(hr, n = 500, presence = deer)
plot(r1)

Of course we are not restricted to the sf::st_buffer function. All geometric operations from the sf package can be used to generate arbitrarily complex availability domains.

Extract covariates

As the next step we have to extract the covariates at point. We can do this with extract_covariates.

rsf1 <- deer |> random_points() |> 
  extract_covariates(sh_forest) 

Fitting RSF

Now all pieces are there to fit a RSF. We will use fit_rsf, which is just a wrapper around stats::glm with family = binomial(link = "logit").

rsf1 |> fit_rsf(case_ ~ forest) |> 
  summary()

Session

sessioninfo::session_info()


Try the amt package in your browser

Any scripts or data that you put into this service are public.

amt documentation built on June 25, 2024, 1:14 a.m.