Köra en modell
Den mjuka trendlinjen du såg i graferna över skörd över tid använder en generaliserad additiv modell (GAM) för att avgöra var linjen ska ligga. Den här typen av modell lämpar sig väl för att anpassa icke-linjära kurvor. För att kunna göra förutsägelser om framtida skördar ska vi nu köra modellen explicit. Syntaxen för att köra en GAM ser ut så här.
gam(response ~ s(explanatory_var1) + explanatory_var2, data = dataset)
Här betyder s() "gör variabeln mjuk", där mjuk ungefär sett innebär icke-linjär.
mgcv och dplyr är inlästa. Datamängderna corn och wheat finns tillgängliga och har utökats med metriska enheter och information om folkräkningsregioner.
Den här övningen är en del av kursen
Introduktion till att skriva funktioner i R
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# Run a generalized additive model of yield vs. smoothed year and census region
___(___ ~ s(___) + ___, data = ___)