Dane o wypadku motocyklowym: podejście nieliniowe

Teraz dopasuj nieliniowy model do tych samych danych mcycle, korzystając z funkcji gam() z pakietu mgcv.

To ćwiczenie jest częścią kursu

Nieliniowe modelowanie z użyciem uogólnionych modeli addytywnych (GAM) w R

Wczytaj pakiet mgcv.
Dopasuj model do danych mcycle, w którym accel ma gładką, nieliniową zależność od times – użyj do tego funkcji gam().
Zwizualizuj dopasowanie modelu, korzystając z dostarczonego wywołania funkcji plot().

Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.

# Load mgcv
___

# Fit the model
gam_mod <- gam(___ ~ s(___), data = mcycle)

# Plot the results
plot(gam_mod, residuals = TRUE, pch = 1)