R Documentation |
qq(df, time = "time", resid = "resid", by_time = TRUE)
df |
Dataframe, containing locations, time, and predictions |
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 |
A ggplot object that can be manipulated further
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 qq(d, time = "year")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.