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 usesmcycledata. Add partial residuals to the plot. - Make a second plot, making partial residuals more visible by changing the shape using the
pchargument, and size of the residuals using thecexargument. Set both thepchandcexarguments 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
___