resid_space

resid_space

resid_space(
  df,
  X = "X",
  Y = "Y",
  time = "time",
  resid = "resid",
  by_time = TRUE,
  outlier_sd = 3
)

Arguments

df

Dataframe, containing locations, time, and predictions

X

Character string containing the name of X variable

Y

Character string containing the name of Y variable

time

Character string containing the name of time variable

resid

Character string containing the name of residual variable

by_time

Whether to facet by time (default = TRUE) or not

outlier_sd

Threshold for how residuals are classified as outliers; defaults to 3 (e.g. a residual value more than 3 SDs is classified as an outlier)

Value

A ggplot object that can be manipulated further

Examples

# \donttest{
set.seed(2021)
d <- data.frame(
  X = runif(1000), Y = runif(1000),
  year = sample(1:10, size = 1000, replace = TRUE)
)
d$density <- rnorm(0.01 * d$X - 0.001 * d$X * d$X + d$Y * 0.02 - 0.005 * d$Y * d$Y, 0, 0.1)
m <- mgcv::gam(density ~ 0 + as.factor(year) + s(X, Y), data = d)
d$resid <- resid(m)
# the default names match, with the exception of year -- so change it
resid_space(d, time = "year")

# }