ozone <- read.table("ozone.txt", header=T, row.names=1, sep=";");

# relation entre température à 12h et température à 15h ? régression linéaire simple
plot(ozone$T12, ozone$T15) # ça a l'air très corrélé
mod.temp <- lm(T15 ~ T12, data=ozone).html
summary(mod.temp) # l'ordonnée à l'origine n'est pas significative, on peut tester que le coeff devant T12 est bien 1


# la vitesse du vent dépend-elle de la direction ? ANOVA

t1 <-ozone[ozone$VENT12=="N", "Vx"];  plot(array(1, length(t1)), t1 , col="red", xlim = c(0,5))
t2 <-ozone[ozone$VENT12=="S", "Vx"];  points(array(2, length(t2)), t2 , col="blue")
t3 <-ozone[ozone$VENT12=="E", "Vx"];  points(array(3, length(t3)), t3 , col="black")
t4 <-ozone[ozone$VENT12=="W", "Vx"];  points(array(4, length(t4)), t4 , col="green")
# on voit bien qu'il y a des différences, on confirme avec le test:
mod.dirVent <- lm(Vx ~ VENT12, data=ozone).html
summary(mod.dirVent)
#le test qu'on veut faire est fait automatiquement à la dernière ligne : l'impression est confirmée, les différences visibles sont significatives

# pour expliquer les maxO3 en gardant les bonnes variables, voir ozone_selection.R
# un modèle raisonnable, obenu par recherche forward, semble garder une ordonnée à l'origine, NE12, maxO3v, T15 et Vx