Aan de slagGa gratis aan de slag

Onzekerheid in auto-prestaties visualiseren

Betrouwbaarheidsintervallen zijn een heel belangrijke visuele indicator van de modelpassing. Hier oefen je met het aanpassen van de weergave van betrouwbaarheidsintervallen en het transformeren van de schaal van de partial effects-plots.

Deze oefening maakt deel uit van de cursus

Niet-lineaire modellering met Generalized Additive Models (GAM's) in R

Cursus bekijken

Oefeninstructies

  • Plot het model (mod) dat de mpg-gegevens gebruikt en plot alleen het partiële effect van weight. Maak het betrouwbaarheidsinterval gearceerd en geef het de kleur "hotpink".
  • Maak nog een plot van het partiële effect van weight, maar verschuif dit keer de schaal met de waarde van het intercept via het shift-argument, en neem de onzekerheid van het modelintercept mee met het seWithMean-argument.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

library(mgcv)
# Fit the model
mod <- gam(hw.mpg ~ s(weight) + s(rpm) + s(price) + comp.ratio, 
           data = mpg, method = "REML")

# Plot the weight effect with colored shading
plot(mod, select = 1, ___)

# Make another plot adding the intercept value and uncertainty
plot(mod, select = 1, ___)
Code bewerken en uitvoeren