Vervuilingsmodellen met multischalige interacties
De meuse-gegevensset bevat enkele voorspellers die op dezelfde schaal liggen (x, y), en enkele die op verschillende schalen liggen (elev, dist, om). In een eerdere oefening paste je een model waarin je cadmiumvervuiling voorspelde als functie van locatie en hoogte:
mod <- gam(cadmium ~ s(x, y) + s(elev),
data = meuse, method = "REML")
In deze oefening bouw je een model waarmee meerdere variabelen kunnen interacteren ondanks deze verschillende schalen, met een tensor smooth, te().
Deze oefening maakt deel uit van de cursus
Niet-lineaire modellering met Generalized Additive Models (GAM's) in R
Oefeninstructies
- Zet dit om naar een model waarin
x,yenelevallemaal met elkaar interacteren in éénte()-term, waarbij ze op hun eigen schaal variëren. - Vat daarna het model samen en visualiseer het met
plot().
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Fit the model
tensor_mod <- ___
# Summarize and plot
___
___