##  automatisation des traitements usuels d'une régression linéaire simple
##  on utilise des fonctions de statgh.r pour les titres

if (!exists("cats")) {  void <- capture.output( source("http://forge.info.univ-angers.fr/~gh/statgh.r",encoding="latin1") ) }

#############################################################################

reglin <- function(x,y,titre) {

#############################################################################

# il serait bon de rajouter les unités,
# de commenter le résultat du test de Fisher et de l'anova

   catss(titre)

   cats("Résumé des données")
   print(summary(data.frame(x,y)))

   # calcul du modèle linéaire
   ml <- lm( y ~ x )

   cats("Coefficients")
   print(ml)

   cats("Modèle complet")
   print( summary(ml) )

   cats("ANOVA du modèle")
   print( anova(ml) )

   idx <- order(x)
   plot(y[idx]~x[idx],type="b",main=titre)
   abline( ml, col="red",lwd=2)

} # fin fonction reglin

## exemples d'utilisation :

# km

km <- lit.dar("http://forge.info.univ-angers.fr/~gh/wstat/Eda/km.dar")
reglin(km$distance,km$essence,"Essence en fonction de distance")

# iris (1)

data(iris)
reglin(iris$Sepal.Width,iris$Sepal.Length ,"Length en fonction Width pour les sépales d'iris")

# iris (2)

reglin(iris$Petal.Width,iris$Sepal.Length ,"Sepal.Length en fonction de Petal.Width pour les données d'iris")