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, model = c("occurrence", "count"), 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.


character; whether to make estimates for the occurrence or count model.


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 or count model. The data frame will have the following columns:

  • stixel_id: unique stixel identifier

  • latitude and longitude: stixel centroid

  • day_of_year: center day of year for stixel

  • 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 count, 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 for occurrence model
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_mid", ext = e, n_bs = 5)