#SÉANCE PRATIQUE DE R***********************************************************

#Si le texte n'est pas bien affiche allez-sur File>Reopen With Encoding et 
#choissisez UTF-8 ou Windows-1252

#Partie 1_______________________________________________________________________
#1.0. Importer le jeu de données "FDAdata.csv" avec:
#File>Import Dataset>From text(base)

#1.1. Affecter les valeurs de la colonne Conc.As à l'objet conc.as
conc.as <- FDAdata$Conc.As
conc.as

#1.2. Sélectionner données des échantillons de l'Arkansas
as.ak <- FDAdata$Conc.As[1:27]

#1.3. Pour créer un nouveau tableau de données avec les données de la Californie
data.ca.df <- FDAdata[28:55,]

#1.4. Pour réinitialiser la numérotation des lignes
row.names(data.ca.df) <- NULL

#1.4b. Façons alternatives de sélectionner les données
data.ca.df2 <- subset(FDAdata,Origin=="California")
data.ca.tx.df <- subset(FDAdata,Origin!="Arkansas")

#1.5. Créer une fonction en R
calcul <-function(y,w) {y^2/w}

#1.6. Pour executer la fonction calcul
calcul(1,2)

#1.7. Créer un tableau de données
tab.df <- data.frame(Origin=c("Arkansas", "California", "Texas"),
                     Conc.As=c(5.4,4.2,7.5))

#1.8. Fonction d'aide. Aide sur la fonction "mean"
help(mean)
mean(as.ak)
sd(as.ak)

#1.9. Sommaire de statistiques descriptives
summary(conc.as)

#EXERCICE 1 ******************************************************************

#1. Créez les vecteurs as.ca et as.tx pour stocker, respectivement, les données de 
#concentration d’arsenique dans les échantillons de riz provenant de la Californie 
#et le Texas. 
as.ca <- conc.as[28:55]
as.ca2 <- FDAdata$Conc.As[28:55] 
as.ca3 <- FDAdata[28:55,4]
as.tx <- conc.as[56:71]


#2. Créez une fonction p.etr pour déterminer le pourcentage écart-type relatif 
#des concentrations d’arsenique des échantillons.


#3. Créez un tableau de données tab.riz.df avec les moyennes, médianes et 
#écart-type des données de l’Arkansas, Texas et Californie.

tab.riz.df <- data.frame(Origine=c("Arkansas", "California", "Texas"),
                         Moyenne=c(mean(as.ak), mean(as.ca), mean(as.tx)),
                         Mediane=c(median(as.ak),median(as.ca),median(as.tx)),
                         EcartType=c(sd(as.ak), sd(as.ca), sd(as.tx)))


#******************************************************************************

par(mfrow=c(1,3))

par(mfrow=c(1,1))