Intitation à RStudio

Introduction

• Nous allons étudier la programmation avec R.
• Nous sommes dans un Notebook.
• Nous pouvons écrire des fomules de mathématiques en LATEX.
• Par exemple, je peux écrire une équation comme :

\[ x^2-12x+14=0. \]

• Des formules un peu plus sophistiquées

\[ \int_{-\infty }^{+\infty }e^{-x^2}dx=\sqrt\pi. \]

Mettre un lien hypertexte Le wikipedia du Cnam.

Les Chunks

• Pour exécuter des lignes de code en R, vous pouvez insérer ces lignes à partir du menu.

1+1

[1] 2

Probability and statistics with R

Vectors

x <- 5 
x

[1] 5

y <- c(7, 3, 5)
y

[1] 7 3 5

z <- c(2, 4, 6, 8)
length(z)

[1] 4

length(x)

[1] 1

length(y)

[1] 3

x + y

[1] 12  8 10

y + z # Opération non souhaitée car pas de même longueur

longer object length is not a multiple of shorter object length

[1]  9  7 11 15

# Supposons que z soient des prix en euros. Pour les convertir en dollars, il suffit de faire
0.87*z

[1] 1.74 3.48 5.22 6.96

LogVec <- (x < z) # logical vector LogVec # 5 < 2, 5 < 4, 5 < 6, 5 < 8 [1] FALSE FALSE TRUE TRUE typeof(LogVec)

LogVec <- (x < z)
LogVec

[1] FALSE FALSE  TRUE  TRUE

typeof(LogVec)

[1] "logical"

typeof(x)

[1] "double"

z

[1] 2 4 6 8

z[2]

[1] 4

typeof(z)

[1] "double"

z<-as.integer(z)
typeof(z)

[1] "integer"

LETTERS ## vecteur contenant les lettres de l'alphabet

 [1] "A" "B" "C" "D" "E" "F" "G" "H" "I" "J" "K" "L" "M" "N"
[15] "O" "P" "Q" "R" "S" "T" "U" "V" "W" "X" "Y" "Z"

LETTERS[10] ## 10ème lettre  

[1] "J"

LETTERS[c(1, 2, 3, 4)] ## Les quatre premières lettres

[1] "A" "B" "C" "D"

LETTERS[1:20] ## Les vingt premières lettres

 [1] "A" "B" "C" "D" "E" "F" "G" "H" "I" "J" "K" "L" "M" "N"
[15] "O" "P" "Q" "R" "S" "T"

1:20 ## crée une liste (ou un vecteur) des nombres de 1 à 20

 [1]  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19
[20] 20

1980:2021     ## Créez un vecteur d'années de 1980 à 2021

 [1] 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990
[12] 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001
[23] 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012
[34] 2013 2014 2015 2016 2017 2018 2019 2020 2021

seq(1980,2021,4)              ##  Créez un vecteur des années bisextiles depuis  1980

 [1] 1980 1984 1988 1992 1996 2000 2004 2008 2012 2016 2020

z

[1] 2 4 6 8

z[1] ## Garder le premier élément

[1] 2

z[-1] ## Enlever le premier élément

[1] 4 6 8

z[c(1,2)]      ## Garder les deux premiers éléments  

[1] 2 4

z[-c(1,2)] ## Enlever les deux premiers éléments

[1] 6 8

z[-(1:2)] ## Idem

[1] 6 8

z[z>5] ## Garder les éléments supérieurs à 5

[1] 6 8

data<-c(rep(c("Nord","Oui"),0.3*30000),
rep(c("Nord","Non"),0.1*30000),
rep(c("Sud","Oui"),0.2*30000),
rep(c("Sud","Non"),0.4*30000))

data<-matrix(data,nrow=30000,ncol=2,byrow=TRUE)
data<-as.data.frame(data)
typeof(data)

[1] "list"

class(data)

[1] "data.frame"

is.data.frame(data)

[1] TRUE

class(x)

[1] "numeric"

typeof(x)

[1] "double"

names(data)<-c("Région","Réponse")
head(data)

te<-table(data) ## Tableau d'effectifs
te          

      Réponse
Région   Non   Oui
  Nord  3000  9000
  Sud  12000  6000

  ## Tableau de proportions
prop.table(te) ## Calcule le tableau en proportions à partir du tableau en effectifs

      Réponse
Région Non Oui
  Nord 0.1 0.3
  Sud  0.4 0.2

addmargins(prop.table(te)) ## Rajoute les marges, sommes en ligne et en colonne

      Réponse
Région Non Oui Sum
  Nord 0.1 0.3 0.4
  Sud  0.4 0.2 0.6
  Sum  0.5 0.5 1.0

addmargins(prop.table(te,1),2) ## Calcule les prop. en ligne (en conditionnant par les lignes) et rajoute les marges qui ont un sens

      Réponse
Région       Non       Oui       Sum
  Nord 0.2500000 0.7500000 1.0000000
  Sud  0.6666667 0.3333333 1.0000000

addmargins(prop.table(te,2),1) ## Calcule les prop. en colonne (conditionnelles) et rajoute les marges qui ont un sens

      Réponse
Région Non Oui
  Nord 0.2 0.6
  Sud  0.8 0.4
  Sum  1.0 1.0

Grades <- c("A", "D", "C", "D", "C", "C", "C", "C", "F", "B") # Crée la variable
Grades # Affiche la variable

 [1] "A" "D" "C" "D" "C" "C" "C" "C" "F" "B"

table(Grades) # Crée une table de fréquence

Grades
A B C D F 
1 1 5 2 1 

xtabs(~Grades) # Idem avec la fonction xtabs

Grades
A B C D F 
1 1 5 2 1 

table(Grades)/length(Grades) # Crée la table en prop.

Grades
  A   B   C   D   F 
0.1 0.1 0.5 0.2 0.1 

prop.table(table(Grades)) # Idem

Grades
  A   B   C   D   F 
0.1 0.1 0.5 0.2 0.1 

prop.table(xtabs(~Grades)) # Idem

Grades
  A   B   C   D   F 
0.1 0.1 0.5 0.2 0.1 

Utiliser des données qui se trouvent dans un Package

addmargins(xtabs(~Age+Sex,data=quine))

     Sex
Age     F   M Sum
  F0   10  17  27
  F1   32  14  46
  F2   19  21  40
  F3   19  14  33
  Sum  80  66 146

Utiliser les pipes

Des graphiques plus sophistiqués

Stem and leaf plot

Statistiques de base

summary(df)

   randomnorm   
 Min.   :17.29  
 1st Qu.:39.07  
 Median :45.74  
 Mean   :46.08  
 3rd Qu.:53.30  
 Max.   :71.25  

Charger des données sur un lien