p3_rsf.R
In amt: Animal Movement Tools

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

## -----------------------------------------------------------------------------
library(amt)
data("deer")
deer

## -----------------------------------------------------------------------------
sh_forest <- get_sh_forest()

## ---- fig.width=4, fig.height=4-----------------------------------------------
r1 <- random_points(deer)
plot(r1)

## ---- fig.width=4, fig.height=4-----------------------------------------------
r1 <- random_points(deer, n = 100)
plot(r1)

## ---- fig.width=4, fig.height=4-----------------------------------------------
hr <- hr_mcp(deer)
r1 <- random_points(hr, n = 500)
plot(r1)

## ---- fig.width=4, fig.height=4-----------------------------------------------
hr <- hr_mcp(deer)
r1 <- random_points(hr, n = 500, presence = deer)
plot(r1)

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

## ---- fig.width=4, fig.height=4-----------------------------------------------
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)

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

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


## -----------------------------------------------------------------------------
sessioninfo::session_info()