The output is plots of cumulative residuals.
Arguments
- model
Fitted model object (glm or gam)
- varlist
Vector of covariate names (continous covariates only)
- label
Label printed at the end of the plot name to identify it if
save=TRUE.- save
(
default=FALSE). Logical stating whether plot should be saved into working directory.
Value
Cumulative residual plots are returned for residuals ordered by each covariate in varlist, predicted value and index of observations (temporally).
The blue dots are the residuals
The black line is the line of cumulative residual.
On the covariate plots (those in varlist) the grey line indicates what we would expect from a well fitted covariate. i.e. one that is fitted with excessive knots.
Note: if the covariate is discrete in nature (like the example below), there will be a lot of overplotting of residuals.
Examples
# load data
data(ns.data.re)
model<-gamMRSea(birds ~ observationhour + as.factor(floodebb) + as.factor(impact),
family='quasipoisson', data=ns.data.re)
plotCumRes(model, varlist=c('observationhour'))
#> [1] "Calculating cumulative residuals"
#> Error in mutate(., Predicted = fitted(model), Index = 1:n()): could not find function "mutate"