#aides
help(q) # Affiche l'aide sur la fonction q
?q # Comme help(q)
help.search("quit") # Cherche le mot quit dans l'aide
??quit # Comme help.search("quit")
q() # Quitte le mode interactif

#charger un fichier
rouge <- read.table("/home/dgm2261a/tp_R/data/WineRed",header=TRUE,sep=";")
blanc <- read.table("/home/dgm2261a/tp_R/data/WineWhite",header=TRUE,sep=";")

#type
typeof(blanc)
typeof(rouge)

#nom des colonnes
names(blanc)
names(rouge)

#nombre de lignes et de colonnes
dim(blanc)
dim(rouge)

#liste colonnes, premieres valeurs
str(blanc)
str(rouge)

#acceder à des valeurs
pH_rouge <- rouge$pH # Avec le nom de la colonne
pH_rouge <- rouge[[9]] # Avec le numéro de la colonne

blanc_density <- blanc$density # Avec le nom de la colonne
blanc_density <- blanc[[8]] # Avec le numéro de la colonne

pH_rouge <- rouge[,9] # Avec le numéro de la colonne
observationRouge1 <- rouge[1,] # Avec le numéro de la ligne
observationRouge1[4]
rouge[1,4]

#moyenne
moy_pH_rouge <- mean(rouge$pH)
moy_pH_blanc <- mean(blanc$pH)

moy_pH_blanc
moy_pH_rouge

pH_blanc_rouge <- c(rouge$pH, blanc$pH) #vecteurs contenant toutes le svaleurs pH
mean(pH_blanc_rouge)

poid <- c(4898, 1599)/(1599+4898)
moy <- c(moy_pH_blanc, moy_pH_rouge)

weighted.mean(moy, poid) #moyenne pondérée

v <- 1:10 #liste de 1 à 10
quantile(v) # Calcul des quantiles de base
quantile(v, probs=seq(0, 1, 0.1)) # Calcul des déciles

quantile(v, probs=seq(0, 1, 0.1), type=7) #je sais pas à quoi type correspond !

quantile(blanc$alcohol)
quantile(rouge$alcohol)

#boxplot
boxplot(rouge$pH) # Diagramme pour le pH des vins rouges
boxplot(blanc$pH) # Diagramme pour le pH des vins blancs
boxplot(rouge$pH, blanc$pH) # Les deux simultanément
boxplot(rouge) # Pour toutes les variables de la liste
boxplot(scale(rouge)) #Pour toutes les variables centrées réduites

boxplot(rouge$residual.sugar)
f <- log(rouge$residual.sugar)/(rouge$residual.sugar**(1/4))
boxplot(f) #plus lisible

#toutes les fonctions possible de boxplot
boxplot(x, range = 1.5, width = NULL, varwidth = FALSE,
        notch = FALSE, outline = TRUE, names, plot = TRUE,
        border = par("fg"), col = NULL, log = "",
        pars = list(boxwex = 0.8, staplewex = 0.5, outwex = 0.5),
        horizontal = FALSE, add = FALSE, at = NULL)

#graphique
qqplot(rouge$alcohol, blanc$alcohol, type="l")
abline(0, 1) #trace une droite sur le graphique, droite de coefficient directeur 1, bx+a
qqplot(rouge$pH, blanc$pH, type="l")

help(abline)

#fonction
variance <- function(v) { #sans poids
  moy <- mean(v)
  var <- 0
  for (i in 1:length(v)) {
    var <- var + (v[i]-moy)^2
  }
  return(var / length(v))
}

variance <- function(v, poids=rep(1, length(v))) { #avec poids qui a une valeur pas defaut
  moy <- mean(v)
  var <- 0
  for (i in 1:length(v)) {
    var <- var + poids[i]*(v[i]-moy)^2
  }
  return(var/sum(poids))
}

variance2 <- function(v, poids=rep(1, length(v))) { #sans boucle
  moy <- mean(v)
  var <- sum(poids*(v-moy)^2)
  return(var/ sum(poids))
}

#calcul le temps d'execution des 2 fonctions variances
source("/home/dgm2261a/tp_R/functions.R") #importe les fonctions tic et toc
tic(); variance(1:100000); toc()
tic(); variance2(1:100000); toc()

variance2(rouge$pH)
variance2(blanc$pH)
variance2(pH_blanc_rouge)