resid_space: resid_space

In ericward-noaa/vista: Visual Inspection of Spatiotemporal Analyses

resid_space

Description

resid_space

Usage

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

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")

ericward-noaa/vista documentation built on June 24, 2022, 12:19 a.m.