pred_resid
pred_resid(df, time = "time", pred = "pred", resid = "resid", by_time = TRUE)
Dataframe, containing locations, time, and predictions
Character string containing the name of time variable
Character string containing the name of predicted variable
Character string containing the name of residual variable
Whether to facet by time (default = TRUE) or not
A ggplot object that can be manipulated further
# \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
pred_resid(d, time = "year")
#> `geom_smooth()` using method = 'loess' and formula 'y ~ x'
# }