Plotting the motorcycle crash model and data
For our first plotting exercise, we'll add partial residuals to the partial effect plot of a GAM so as to compare the model to the data.
This exercise is part of the course
Nonlinear Modeling with Generalized Additive Models (GAMs) in R
Exercise instructions
- Plot the model (
mod
) that usesmcycle
data. Add partial residuals to the plot. - Make a second plot, making partial residuals more visible by changing the shape using the
pch
argument, and size of the residuals using thecex
argument. Set both thepch
andcex
arguments to 1.
Hands-on interactive exercise
Have a go at this exercise by completing this sample code.
library(mgcv)
# Fit the model
mod <- gam(accel ~ s(times), data = mcycle, method = "REML")
# Make the plot with residuals
plot(mod, ___)
# Change shape of residuals
___