Sojabohnenwachstum mit GAM modellieren
In dieser Übung modellierst du das durchschnittliche Blattgewicht einer Sojabohnenpflanze als Funktion der Zeit (nach dem Pflanzen). Wie du sehen wirst, wächst die Sojabohne nicht mit gleichmäßigem Tempo, sondern hat einen „Wachstumsschub“, der später abflacht. Daher lässt sich das Blattgewicht nicht gut mit einem linearen Modell beschreiben.
Denk daran, dass du in einer Formel mit der Funktion s() (Docs) angeben kannst, welche Variable du nichtlinear modellieren möchtest:
y ~ s(x)
Ebenfalls wichtig: gam() (Docs) aus dem Paket mgcv hat die folgende Aufrufsignatur
gam(formula, family, data)
Für Standardregression verwende family = gaussian (Voreinstellung).
Die Trainingsdaten soybean_train sind vorab geladen. Sie haben zwei Spalten: das Ziel weight und die Variable Time. Zum Vergleich wurde auch das lineare Modell model.lin, das mit der Formel weight ~ Time gefittet wurde, bereits geladen.
Diese Übung ist Teil des Kurses
<Kurs>Überwachtes Lernen in R: Regression</Kurs>Interaktive praktische Übung
Versuche dich an dieser Übung, indem du diesen Beispielcode vervollständigst.
# soybean_train is available
summary(soybean_train)
# Plot weight vs Time (Time on x axis)
ggplot(soybean_train, aes(x = ___, y = ___)) +
geom_point()