Modellera sojabönsväxt med GAM
I den här övningen modellerar du den genomsnittliga bladmassan på en sojabönsplanta som en funktion av tid (efter plantering). Som du kommer att se växer sojabönsplantan inte i jämn takt, utan har en "tillväxtspurt" som gradvis avtar. Bladmassan beskrivs därför inte väl av en linjär modell.
Kom ihåg att du anger vilken variabel som ska modelleras icke-linjärt i en formel med funktionen s() (docs):
y ~ s(x)
Kom också ihåg att gam() (docs) från paketet mgcv anropas med
gam(formula, family, data)
Vid vanlig regression använder du family = gaussian (standardvärdet).
Träningsdatan för sojaböna, soybean_train, är förinstallerad. Den har två kolumner: utfallet weight och variabeln Time. För jämförelse har även den linjära modellen model.lin, som anpassades med formeln weight ~ Time, redan laddats in.
Den här övningen är en del av kursen
Övervakad inlärning i R: Regression
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# soybean_train is available
summary(soybean_train)
# Plot weight vs Time (Time on x axis)
ggplot(soybean_train, aes(x = ___, y = ___)) +
geom_point()