Partial dependence (PD) plots depict the relationship between the modeled occurrence probability and each of the predictor variables used in the model. Status and Trends provides the data to generate these plots for every stixel.

load_pds(path, ext, return_sf = FALSE)



character; directory that the Status and Trends data for a given species was downloaded to. This path is returned by ebirdst_download() or get_species_path().


ebirdst_extent object; the spatiotemporal extent to filter the data to. The spatial component of the extent object must be provided in unprojected, latitude-longitude coordinates.


logical; whether to return an sf object of spatial points rather then the default data frame.


Data frame, or sf object if return_sf = TRUE, containing PD estimates for each stixel for either the occurrence and relative model. The data frame will have the following columns:

  • stixel_id: unique stixel identifier

  • lat and lon: stixel centroid

  • date: day of year, expressed as a value from 0-1, of the stixel center

  • predictor: name of the predictor that the PD data correspond to, for a full list of predictors consult the ebirdst_predictors data frame

  • predictor_value: value of the predictor variable at which PD is evaluated

  • response: predicted response, occurrence or relative abundance, at the given value of the predictor averaged across all the values of the other predictors


if (FALSE) {
# download example data
path <- ebirdst_download("example_data", tifs_only = FALSE)
# or get the path if you already have the data downloaded
path <- get_species_path("example_data")

# load partial dependence data
pds <- load_pds(path)

# plot the top 15 predictor importances
# define a spatiotemporal extent to plot data from
bb_vec <- c(xmin = -86.6, xmax = -82.2, ymin = 41.5, ymax = 43.5)
e <- ebirdst_extent(bb_vec, t = c("05-01", "05-31"))
plot_pds(pds, "solar_noon_diff", ext = e, n_bs = 5)