Modéliser la croissance du soja avec un GAM
Dans cet exercice, vous allez modéliser le poids moyen des feuilles d’un plant de soja en fonction du temps (après la plantation). Comme vous le verrez, le soja ne pousse pas à un rythme constant : il connaît un « pic de croissance » qui finit par ralentir. Le poids des feuilles n’est donc pas bien décrit par un modèle linéaire.
Rappelez-vous que vous pouvez indiquer quelle variable modéliser de façon non linéaire dans une formule avec la fonction s() (docs) :
y ~ s(x)
Souvenez-vous aussi que gam() (docs) du package mgcv s’utilise avec l’interface suivante :
gam(formula, family, data)
Pour une régression standard, utilisez family = gaussian (valeur par défaut).
Les données d’entraînement sur le soja, soybean_train, ont été préchargées. Elles comportent deux colonnes : la cible weight et la variable Time. Pour comparaison, le modèle linéaire model.lin, ajusté avec la formule weight ~ Time, a également été chargé.
Cet exercice fait partie du cours
<cours>Apprentissage supervisé en R : Régression</cours>Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# soybean_train is available
summary(soybean_train)
# Plot weight vs Time (Time on x axis)
ggplot(soybean_train, aes(x = ___, y = ___)) +
geom_point()