# ==================================================
# Exemple de régression linéaire avec la commande lm
# ==================================================

# on stocke les données de data.txt dans l’objet « klein »
klein = read.table("data.txt", header=TRUE)
klein                           # pour voir…
fix(klein)                      # pour éditer (refermer la fenêtre d'édition pour poursuivre)

summary(klein)                  # pour des détails

summary(wtot)                   # appel nominatif impossible
attach(klein)
summary(wtot)                   # wtot est à présent connu

# la commande lm

help(lm)                        # si on veut lire la documentation

lm(cons~p+plag+wtot)            # la régression
reg = lm(cons~p+plag+wtot)
reg                             # reproduit l’affichage minimal des résultats
summary(reg)                    # affiche les résultats détaillés

reg$residuals                   # les résidus
reg$fitted                      # la série ajustée
plot(reg$residuals)             # on figure les résidus
lines(reg$residuals)            # idem avec la ligne qui les joint
plot(wtot, reg$residuals)       # un graphique croisé

# la "stepwise regression"
step(reg)