Rstudio TD 1: Clément Rau Laboratoire de Mathématiques de Toulouse Université Paul Sabatier-IUT GEA Ponsan

Prise en main
Quelques fonctions liées aux statistiques descriptives

Importer des données depuis un tableur
Exercices
Pour aller plus (ou moins) loin...
Rstudio TD 1
Clément Rau
Laboratoire de Mathématiques de Toulouse
Université Paul Sabatier-IUT GEA Ponsan
Clément Rau Rstudio TD 1

Prise en main
Exercices
1 Prise en main
Introduction
Opérations élémentaires
Affectation
Structures de contrôle
Autour des séries
2 Quelques fonctions liées aux statistiques descriptives
Quantités cumulées
Calcul des paramètres standards
Représentations graphique
3 Importer des données depuis un tableur
4 Exercices
Exo 1
Exo 2
5 Pour aller plus (ou moins) loin...
Exo 3
Exo 4 Clément Rau Rstudio TD 1
Prise en main Introduction
Quelques fonctions liées aux statistiques descriptives Opérations élémentaires
Importer des données depuis un tableur Affectation
Exercices Structures de contrôle
Pour aller plus (ou moins) loin... Autour des séries
Remarques sur R
R est un logiciel permettant de faire des analyses statistiques et de
produire des graphiques.
R est aussi un langage de programmation complet.
I R est un logiciel gratuit qui se télécharge sur
www.r-project.org (Vidéo tuto ici.)
I Il existe également une version en ligne ici , qui ne
nécessite donc pas d’installation, mais juste la création
d’un compte.
I Enfin, pour les "fans d’ios", sachez qu’il existe un appli
gratuite sur iphone (nommée R compiler) qui fonctionne
très bien. (Probablement qu’une appli similaire sur Android
existe aussi.)
De l’aide sur la fonction bloumbloum s’obtient par
>help(bloumbloum)

Remarques sur R
d’un compte.
existe aussi.)
>help(bloumbloum)

Remarques sur R
d’un compte.
existe aussi.)
>help(bloumbloum)

Remarques sur R
d’un compte.
existe aussi.)
>help(bloumbloum)

Remarques sur R
d’un compte.
existe aussi.)
>help(bloumbloum)

Remarques sur R
d’un compte.
existe aussi.)
>help(bloumbloum)

Interface, Fenêtres dans R

Répertoire de travail
R a par défaut un répertoire de travail. Cela signifie que si
aucun chemin n’est spécifié, R va lire ou écrire dans ce
répertoire.
Pour connaître ce répertoire de travail, saisir dans la console
de commande :
getwd()
Pour changer de répertoire de travail, saisir dans la console de

commande :
setwd("chemin du répertoire souhaité")
Attention, dans les chemins Windows, les \ doivent être

remplacé par des / ou des \\
répertoire.
de commande :
getwd()

commande :

répertoire.
de commande :
getwd()

commande :

répertoire.
de commande :
getwd()

commande :

1 Prise en main
Introduction
Affectation
Autour des séries
4 Exercices
Exo 1
Exo 2
Exo 3
Essayer qques opérations élémentaires : addition,

multiplication etc...
>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))


>2*5
>4+2
>2,5*2
>2.5*2
>2ˆ3
>sqrt(16)
>ln(2)
>log(2)
>log(100,10)
>exp(1)
>log(exp(3))

Nombre de décimales souhaitées
On utilise la fonction :
round(x,n)
où x désigne un réel et n un entier.

La fonction renvoie une valeur approchée du réel x avec
uniquement les n premières décimales.

Nombre de décimales souhaitées
On utilise la fonction :
round(x,n)
où x désigne un réel et n un entier.

La fonction renvoie une valeur approchée du réel x avec
uniquement les n premières décimales.

1 Prise en main
Introduction
Affectation
Autour des séries
4 Exercices
Exo 1
Exo 2
Exo 3
Affectation
Si l’on veut affecter à x la valeur 3, on peut taper :

>x=3 ou bien >x<-3
On peut alors demander que vaut x :
>x
Donner à x et y deux valeurs et faire x+y, 3x-y. Controler.
! ! ! Attention R distingue majuscule et minuscule ! ! !
Pour affecter plusieurs variables, on les sépare par un ;

Exemple : a=2 ; b=3; c=0

Affectation

>x=3 ou bien >x<-3
>x


Affectation

>x=3 ou bien >x<-3
>x


Affectation

>x=3 ou bien >x<-3
>x


Affectation

>x=3 ou bien >x<-3
>x


Affectation

>x=3 ou bien >x<-3
>x


Si l’on veut affecter à S1 la série de nombres 3,5,1,12,20.

On tape :
>S1=c(3,5,1,12,20)
Vérifier en tapant >S1
Si l’on veut affecter à S2 la série de caractéres : rouge,
bleu, vert. On tape :
>S2=c("rouge", "bleu","vert")
Vérifier en tapant > S2
En tapant >mode(S2), R nous indique la nature de la série S2.
En tapant >length(S1), R nous indique la taille de la série
S1.


On tape :
>S1=c(3,5,1,12,20)
S1.


On tape :
>S1=c(3,5,1,12,20)
S1.


On tape :
>S1=c(3,5,1,12,20)
S1.


On tape :
>S1=c(3,5,1,12,20)
S1.

Rentrer une série rapidement
En utilisant SCAN la saisie d’une série est moins fastidieuse.

Taper >S=scan()
On rentre successivement les valeurs de la série. Essayer

Pense bête sur les affectations
la fonction >ls.str() permet d’afficher tous les détails

sur les affectations en cours.
la fonction >rm(list=ls()) permet d’effacer toutes les
affectations.

Pense bête sur les affectations
la fonction >ls.str() permet d’afficher tous les détails

sur les affectations en cours.
la fonction >rm(list=ls()) permet d’effacer toutes les
affectations.

1 Prise en main
Introduction
Affectation
Autour des séries
4 Exercices
Exo 1
Exo 2
Exo 3
Maintenant que nous avons fait un tour d’ensemble des

fonctionnalités de R, revenons à quelques rudiments de
programmation.

Structure conditionnelle, Boucle if
Syntaxe :
if (condition) A else B
calcule les instructions A si la condition est vraie, et les

instructions B sinon.

Syntaxe :
if (condition) A else B
calcule les instructions A si la condition est vraie, et les

instructions B sinon.

Exemple 1 :
>a=3
>if (a>0) a=a+1 else a=a-1
>a
[1] 4
Essayer avec a =-5

Exemple 1 :
>a=3
>a
[1] 4
Essayer avec a =-5

Exemple 1 :
>a=3
>a
[1] 4
Essayer avec a =-5

Exemple 2 :
>a=2
>if (a>0) b=log(a) else b=0
>b
[1] 0.6931472
Essayer avec a =-3

Exemple 2 :
>a=2
>b
[1] 0.6931472
Essayer avec a =-3

Exemple 2 :
>a=2
>b
[1] 0.6931472
Essayer avec a =-3

Structure itérative, Boucle for
Syntaxe :
for (var in seq) commandes

Syntaxe :
for (var in seq) commandes

Exemple 1 :
> a=4
> for (i in 1:2) a=a+i
> a
[1] 7

Exemple 1 :
> a=4
> for (i in 1:2) a=a+i
> a
[1] 7

Exemple 2 :
> for (i in 1:5) print(i)
[1] 1
[1] 2
[1] 3
[1] 4
[1] 5

Exemple 2 :
> for (i in 1:5) print(i)
[1] 1
[1] 2
[1] 3
[1] 4
[1] 5


Exemple 3 :
> a=100
> for (i in 1:5) (a=a+1) & print(a) & print(i)
[1] 101
[1] 1
[1] 102
[1] 2
[1] 103
[1] 3
[1] 104
[1] 4
[1] 105
[1] 5

Exemple 3 :
> a=100
> for (i in 1:5) (a=a+1) & print(a) & print(i)
[1] 101
[1] 1
[1] 102
[1] 2
[1] 103
[1] 3
[1] 104
[1] 4
[1] 105
[1] 5
Exemple 4 :
> b=numeric(length(10))
> for (i in 1:10) b[i]=iˆ2
> b
[1] 1 4 9 16 25 36 49 64 81 100
>

Exemple 4 :
> b=numeric(length(10))
> for (i in 1:10) b[i]=iˆ2
> b
[1] 1 4 9 16 25 36 49 64 81 100
>

Structure itérative, Boucle while
Syntaxe :
while (condition) commandes

Structure itérative, Boucle while
Syntaxe :
while (condition) commandes

Structure itérative, while
Exemple 1 :
> n=1
> i<-1
>
> while (i<6) {
+ print(i)
+ n<-i*n
+ i<-i+1
+ print(n) }
Remarque : shift + entrée, vous permet d’aller à la ligne dans
les instructions sous R.
Vous venez de fabriquer factoriel de 5, noté 5!

Exemple 1 :
> n=1
> i<-1
>
> while (i<6) {
+ print(i)
+ n<-i*n
+ i<-i+1
+ print(n) }

Exemple 1 :
> n=1
> i<-1
>
> while (i<6) {
+ print(i)
+ n<-i*n
+ i<-i+1
+ print(n) }

Exemple 1 :
> n=1
> i<-1
>
> while (i<6) {
+ print(i)
+ n<-i*n
+ i<-i+1
+ print(n) }

1 Prise en main
Introduction
Affectation
Autour des séries
4 Exercices
Exo 1
Exo 2
Exo 3
Accès à une ou des valeurs de la série
Soit S une série, par exemple S = (12, 10, 11, 19, 9, 10). On
souhaite parfois avoir accès à la ième valeur, on utilise alors :
>S[3]
Essayer !
Si l’on souhaite avoir uniquement les 3, 4 et 5ième valeurs de
la série, on tape :
>S[3:5]

>S[3]
Essayer !
>S[3:5]

>S[3]
Essayer !
>S[3:5]

Concaténation de 2 séries
Soit S1, et S2 deux séries. Si l’on veut rassembler les donner

en une série S3, on tape :
>S3=c(S1,S2)
Essayer !

Concaténation de 2 séries
Soit S1, et S2 deux séries. Si l’on veut rassembler les donner

en une série S3, on tape :
>S3=c(S1,S2)
Essayer !

Opération sur les séries
Soit S une série, par exemple S = (12, 10, 11, 19, 9, 10). Si l’on
souhaite multiplier par 2 chaque valeurs puis addtionner 1 à
chaque valeur, on tape alors :
>S= 2*S+1
Ici on a écrasé les anciennes valeurs de S. Essayer

>S= 2*S+1

>S= 2*S+1

Remplacement d’une valeur d’une série
Soit S une série, par exemple S = (12, 10, 11, 19, 9, 10).
Si l’on souhaite remplacer la 2ème valeur par un 5, on
tape : >S[2]=5
Essayer
Si l’on souhaite remplacer tous les 10 par 100, on tape :
>S[S==10]=100
Essayer
Si l’on souhaite remplacer toutes les valeurs supérieures à
10 par 15, on tape : >S[S>10]=15
Essayer

tape : >S[2]=5
Essayer
>S[S==10]=100
Essayer
10 par 15, on tape : >S[S>10]=15
Essayer

tape : >S[2]=5
Essayer
>S[S==10]=100
Essayer
10 par 15, on tape : >S[S>10]=15
Essayer

tape : >S[2]=5
Essayer
>S[S==10]=100
Essayer
10 par 15, on tape : >S[S>10]=15
Essayer

tape : >S[2]=5
Essayer
>S[S==10]=100
Essayer
10 par 15, on tape : >S[S>10]=15
Essayer

tape : >S[2]=5
Essayer
>S[S==10]=100
Essayer
10 par 15, on tape : >S[S>10]=15
Essayer

tape : >S[2]=5
Essayer
>S[S==10]=100
Essayer
10 par 15, on tape : >S[S>10]=15
Essayer

Tri
Si l’on souhaite trier les valeurs de la série par ordre
croissant, on tape :
>sort(S)
Essayer !
Essayer ensuite
>rev(sort(S))

Tri
>sort(S)
Essayer !
Essayer ensuite
>rev(sort(S))

Tri
>sort(S)
Essayer !
Essayer ensuite
>rev(sort(S))

Regroupement par caractère d’une série
La fonction table donne les effectifs de chaque modalité de la

variable. Rentrer une série qualitative ou quantitative avec
quelques redondances, notée S. Puis taper :
>table(S)

1 Prise en main
Introduction
Affectation
Autour des séries
4 Exercices
Exo 1
Exo 2
Prise en main
Quelques fonctions liées aux statistiques descriptives Quantités cumulées
Importer des données depuis un tableur Calcul des paramètres standards
Exercices Représentations graphique
Exo 3
Exo 4

Prise en main
Effectifs cumulés
Soit S une série. Pour obtenir la série des valeurs cumulées

croissantes, on tape :
>cumsum(S)
Rentrer une série S, puis taper >cumsum(S).
Question : Soit S une série. Comment obtenir les effectifs

cumulés croissants ?

Prise en main
Effectifs cumulés

>cumsum(S)


Prise en main
Effectifs cumulés

>cumsum(S)


Prise en main
Effectifs cumulés
Prenons par exemple S = (12, 10, 10, 12, 11, 19, 9, 10). Taper
>EC=cumsum(table(S)) Demander >EC, vous obtenez :
> EC
9 10 11 12 19
1 4 5 7 8
Ce qui correspond bien aux effectifs cumulés croissants de la
série S regroupée par "valeur" :
> table(S)
9 10 11 12 19
1 3 1 2 1

Prise en main
Effectifs cumulés
> EC
9 10 11 12 19
1 4 5 7 8
> table(S)
9 10 11 12 19
1 3 1 2 1

Prise en main
Effectifs cumulés
> EC
9 10 11 12 19
1 4 5 7 8
> table(S)
9 10 11 12 19
1 3 1 2 1

Prise en main
Effectifs cumulés
> EC
9 10 11 12 19
1 4 5 7 8
> table(S)
9 10 11 12 19
1 3 1 2 1

Prise en main
Effectifs cumulés
> EC
9 10 11 12 19
1 4 5 7 8
> table(S)
9 10 11 12 19
1 3 1 2 1

Prise en main
Fréquences cumulées
Question : Soit S une série. Fabriquez les fréquences

cumulées croissantes ?
Prenons la même série S = (12, 10, 10, 12, 11, 19, 9, 10). Taper
>FC=cumsum(table(S))/8 Demander >FC, vous obtenez :
> FC
9 10 11 12 19
0.125 0.500 0.625 0.875 1.000

Prise en main

> FC
9 10 11 12 19
0.125 0.500 0.625 0.875 1.000

Prise en main

> FC
9 10 11 12 19
0.125 0.500 0.625 0.875 1.000

Prise en main

> FC
9 10 11 12 19
0.125 0.500 0.625 0.875 1.000

Prise en main
Max, Min
Rentrer une série, que l’on nommera "notes". Par exemple vos
notes en maths !
Taper >max(notes) puis >min(notes)
Fabriquer alors l’étendue.

Prise en main
Max, Min
notes en maths !

Prise en main
Max, Min
notes en maths !

Prise en main
Moyenne
Calculer à la main la moyenne de votre série notes

Taper >sum(notes), puis >length(notes) et enfin
>sum(notes)/length(notes)
Taper >mean(notes)
Conclusion ?

Prise en main
Moyenne

Taper >mean(notes)
Conclusion ?

Prise en main
Moyenne

Taper >mean(notes)
Conclusion ?

Prise en main
Moyenne

Taper >mean(notes)
Conclusion ?

Prise en main
Moyenne avec données manquantes
Parfois, certaines données manquantes d’une série sont

rentrées comme "NA". Pour calculer la moyenne en occultant
ces données, on tape :
>mean(notes, na.rm = TRUE)
Essayer en rentrant une série contenant des NA.

Prise en main


Prise en main


Prise en main
Moyenne de zone de série
Si l’on ne veut calculer la moyenne que des 3 premières

notes, on tapera :
>mean(notes[1:3])
Exercices : Calculer la moyenne des 5 notes les plus
faibles.
Indication : utiliser mean(sort(notes)[1 :5])

Prise en main
Moyenne de zone de série
Si l’on ne veut calculer la moyenne que des 3 premières

notes, on tapera :
>mean(notes[1:3])
Exercices : Calculer la moyenne des 5 notes les plus
faibles.
Indication : utiliser mean(sort(notes)[1 :5])

Prise en main
Médiane
Taper
>median(notes)

Prise en main
Quantiles
Taper
>quantile(notes)
Pour obtenir uniquement le 1er quantile, on tape :
>quantile(notes, probs = .25)
Et si l’on veut uniquement la valeur, on tape :
> quantile(notes, probs = .25, names =
FALSE)
Exercice : Fabriquer l’écart inter quantile

Prise en main
Quantiles
Taper
>quantile(notes)
FALSE)

Prise en main
Quantiles
Taper
>quantile(notes)
FALSE)

Prise en main
Quantiles
Taper
>quantile(notes)
FALSE)

Prise en main
Déciles
Taper
> quantile(notes, probs = .1)
Pour obtenir tous les déciles, on tape :
>quantile(notes, probs = 0:10/10)
Taper >0:10/10 pour comprendre...

Prise en main
Déciles
Taper

Prise en main
Déciles
Taper

Prise en main
Variance
Taper
>var(notes)
Taper >mean( (notes-mean(notes))ˆ2)
Comparer vos résultats ? Explication ?

Prise en main
Variance
Taper
>var(notes)

Prise en main
Variance
Taper
>var(notes)

Prise en main
Variance, explication
La fonction >var(notes) de R, renvoie l’expression
suivante :
1 X
ni [xi − E(X )]2
n−1
i
alors que la "vraie" variance correspond à :

1X
ni [xi − E(X )]2
n
i
⇒ Ainsi pour obtenir la "vraie" variance, vous avez 2 options :

-Taper >mean( (notes-mean(note))ˆ2)
-Taper
>((length(notes)-1) /length(notes))*var(notes)
Prise en main
suivante :
1 X
ni [xi − E(X )]2
n−1
i

1X
ni [xi − E(X )]2
n
i

-Taper
Prise en main
suivante :
1 X
ni [xi − E(X )]2
n−1
i

1X
ni [xi − E(X )]2
n
i

-Taper
Prise en main
Remarque : la variance "modifiée" de R, est sans biais et sera

utile dans les estimations et tests statistiques.

Prise en main
Ecart type
Taper >sd(notes)
Là encore, pour obtenir le "vrai" écart-type, vous pouvez :
-soit le calculer manuellement en prenant la racine carré de la
"vraie" variance
-soit taper
>sqrt((length(notes)-1)/length(notes))
*sd(notes)

Prise en main
Ecart type
Taper >sd(notes)
Là encore, pour obtenir le "vrai" écart-type, vous pouvez :
-soit le calculer manuellement en prenant la racine carré de la
"vraie" variance
-soit taper
>sqrt((length(notes)-1)/length(notes))
*sd(notes)

Prise en main
Paramètres en un seul clic
Taper
>summary(notes)

Prise en main
1 Prise en main
Introduction
Affectation
Autour des séries
4 Exercices
Exo 1
Exo 2
Prise en main
Exo 3
Exo 4

Prise en main
Graphe basique
Rentrer une série S de quelques valeurs, puis Taper
>plot(S)

Prise en main
Diagramme en bâtons
Rentrer une série S de quelques valeurs (avec des

redondances), puis Taper
>plot(table(S))
On peut améliorer la présentation en tapant :
>plot(table(S), main="titre", ylab="axe des

y")

Prise en main
Diagramme en bâtons
Rentrer une série S de quelques valeurs (avec des

redondances), puis Taper
>plot(table(S))
On peut améliorer la présentation en tapant :
>plot(table(S), main="titre", ylab="axe des

y")

Prise en main
histogramme
Avec la même série S de quelques valeurs, taper :

>hist(S)
On constate que les classes sont imposées...
Pour choisir les classes, on tape :
>hist(S,breaks=c(0,4,15,20))
si l’on souhaite que les bornes des classes soient 0, 4, 15
et 20.

Prise en main
histogramme

>hist(S)
et 20.

Prise en main
histogramme

>hist(S)
et 20.

Prise en main
Histogramme

Prise en main
Camembert
Soit S une série qualitative, on tape :
>pie(table(S))

Prise en main
Boxplot-Boîte à moustache

>boxplot(S)
Pour obtenir la boxplot horizontalement, taper :
>boxplot(S,horizontal=TRUE)

Prise en main

>boxplot(S)
Pour obtenir la boxplot horizontalement, taper :
>boxplot(S,horizontal=TRUE)

Prise en main

Prise en main
Exercices
1 Prise en main
Introduction
Affectation
Autour des séries
4 Exercices
Exo 1
Exo 2
Prise en main
Exercices
Exo 3
Exo 4

Prise en main
Exercices
Il est assez courant de vouloir importer des données saisies

avec un tableur du type Open Office ou Excel. Rstudio permet
d’importer (entre autres) :
des fichiers .csv
des fichiers .xlsx
Il suffit d’aller dans la fenêtre en bas à droite de Rstudio, de
selectionner l’emplacement du fichier de données. Puis par un
clic gauche sur la souris, selectionner sur "import dataset".
Le tableau de données apparaît alors dans la fenêtre en haut à
gauche de Rstudio.

Prise en main
Exercices

des fichiers .csv
des fichiers .xlsx
gauche de Rstudio.

Prise en main
Exercices

des fichiers .csv
des fichiers .xlsx
gauche de Rstudio.

Prise en main
Exercices

des fichiers .csv
des fichiers .xlsx
gauche de Rstudio.

Prise en main
Exercices
On peut visualiser ce tableau de données en tapant

>View(nom du fichier)
Si l’on souhaite accéder seulement aux effectifs associé à
un caractére X (une colonne), on tapera
>nom du fichier$X

Prise en main
Exercices
On peut visualiser ce tableau de données en tapant

>View(nom du fichier)
Si l’on souhaite accéder seulement aux effectifs associé à
un caractére X (une colonne), on tapera
>nom du fichier$X

Prise en main
Exercices
1 Prise en main
Introduction
Affectation
Autour des séries
4 Exercices
Exo 1
Exo 2
Prise en main
Exo 1
Exo 2
Exercices
Exo 3
Exo 4

Prise en main
Exo 1
Exo 2
Exercices
Exo 1
A l’aide du fichier formulaire.xlxs ou csv (formulaire
sur les moyens de transport) à télécharger sur la page web,
répondre aux questions suivantes.
1 Faire apparaître ce tableau de données dans Rstudio.
2 Combien de personnes utilisent le bus comme moyen de
transport ?
3 Quel est le budget moyen par mois pour le transport ?
Calculer l’écart type. Tracer la boxplot relative au budget.
4 Après avoir déterminer des classes raisonnables
concernant le budget, tracer l’ histogramme des effectifs
associés à vos classes.
5 Tracer le diagramme en bâton relatif aux marques de
véhicules, puis tracer le camembert associé. (commencer
par faire un regroupement par marque.)
Prise en main
Exo 1
Exo 2
Exercices
Exo 1 suite
6 Calculer les pourcentages de chaque marque de

véhicules.
7 Combien de personnes ont une voiture personnelle ?
8 Calculer les quantiles des distances parcourues par les
personnes de cette étude. Tracer la boxplot.
9 Combien de personnes utilisent uniquement le bus ?

Prise en main
Exo 1
Exo 2
Exercices
Exo 1 indication pour question 7
On pourra taper :
>b=0
> for (i in 1:length(Formulaire$VOITUREPERSO))
(if (Formulaire$VOITUREPERSO[i]=="Oui") b=b+1)
> b
[1] 27

Prise en main
Exo 1
Exo 2
Exercices
Exo 1 indication pour question 7
On pourra taper :
>b=0
> for (i in 1:length(Formulaire$VOITUREPERSO))
(if (Formulaire$VOITUREPERSO[i]=="Oui") b=b+1)
> b
[1] 27

Prise en main
Exo 1
Exo 2
Exercices
Exo 2
On a relevé les temps de retard (en min) des employés dans
une entreprise :
5, 7, 12, 1, 5, 9, 4, 0, 0, 0, 14, 4, 22, 2, 6, 10, 5
1 Calculer la moyenne, l’étendue, la variance, l’écart type

2 Calculer les quantiles, ainsi que l’écart inter quantile.
3 Donner la valeur du 3ième décile.
4 Représenter la boxplot
5 Tracer un diagramme en batons.
6 Y a t’il une (ou des) valeur(s) aberrante(s) ?
(Rappel : une valeur est dite aberrante si elle est supérieure à
Q3 + 1, 5(Q3 − Q1 ) ou inférieure à Q1 − 1, 5(Q3 − Q1 ))
Prise en main
Exo 1
Exo 2
Exercices
Exo 2
On a relevé les temps de retard (en min) des employés dans
une entreprise :
5, 7, 12, 1, 5, 9, 4, 0, 0, 0, 14, 4, 22, 2, 6, 10, 5
1 Calculer la moyenne, l’étendue, la variance, l’écart type

2 Calculer les quantiles, ainsi que l’écart inter quantile.
3 Donner la valeur du 3ième décile.
4 Représenter la boxplot
5 Tracer un diagramme en batons.
6 Y a t’il une (ou des) valeur(s) aberrante(s) ?
(Rappel : une valeur est dite aberrante si elle est supérieure à
Q3 + 1, 5(Q3 − Q1 ) ou inférieure à Q1 − 1, 5(Q3 − Q1 ))
Prise en main
Exo 3
Exo 4
Exercices
1 Prise en main
Introduction
Affectation
Autour des séries
4 Exercices
Exo 1
Exo 2
Exo 3
Prise en main
Exo 3
Exo 4
Exercices
Pour définir une fonction en R
Supposons que l’on veuille définit la fonction f définie par
f (x) = x 2 − 5x + 1.
On peut taper en R :
f<-function(x) { xˆ2-5? x +1}
Si l’on veut avoir la valeur de f (1). Il suffit de taper :
f(1)
[1] -3

Prise en main
Exo 3
Exo 4
Exercices
Pour définir une fonction en R
Supposons que l’on veuille définit la fonction f définie par
f (x) = x 2 − 5x + 1.
On peut taper en R :
f<-function(x) { xˆ2-5? x +1}
Si l’on veut avoir la valeur de f (1). Il suffit de taper :
f(1)
[1] -3

Prise en main
Exo 3
Exo 4
Exercices
Exo 3
1 Programmer la fonction factorielle avec une boucle "while".

On la nommera fact
2 Programmer la fonction factorielle de manière récursive.
On l’appelera factrec.
3 A l’aide d’une
des 2 fonctions, en déduire un moyen
n 5
d’afficher . Tester en calculant par ex .
k 2

Prise en main
Exo 3
Exo 4
Exercices
Eléments de réponse (exo 3)
En utilisant "while"
fact<-function(x){ res=1 ;
while(x>1) { res=res ? x; x= x-1 }
res }
Vérifier la valeur de fact(5)

Prise en main
Exo 3
Exo 4
Exercices
En utilisant "while"
fact<-function(x){ res=1 ;
while(x>1) { res=res ? x; x= x-1 }
res }
Vérifier la valeur de fact(5)

Prise en main
Exo 3
Exo 4
Exercices
Forme récursive
factrec<-function(x){tmp<-1
if(x>0){ tmp<-factrec(x-1)*x }
return(tmp)}
Vérifier la valeur de factrec(5)

Prise en main
Exo 3
Exo 4
Exercices
Forme récursive
factrec<-function(x){tmp<-1
if(x>0){ tmp<-factrec(x-1)*x }
return(tmp)}
Vérifier la valeur de factrec(5)

Prise en main
Exo 3
Exo 4
Exercices
Calcul des coefficients binomiaux
n=6; k=3;
print(paste(k, "parmi", n ,"vaut"));
fact(n)/(fact(k)*fact(n-k))

Prise en main
Exo 3
Exo 4
Exercices
Tours de Hanoi

Prise en main
Exo 3
Exo 4
Exercices
Tours de Hanoi
Le jeu consiste à déplacer les n palets de la tour 1 sur la tour 3.

Le jeu n’aurait evidemment guère d’intêret si on ne rajoutait pas
la contrainte suivante : à chaque étape, aucun palet ne doit être
posé sur un palet de diamètre inférieur.
Ecrire une fonction hanoi(n,a,b,c) qui prend comme
arguments : n (nombre de palets sur la tour 1), et a, b, c qui
sont les noms de tours, et qui renvoie la suite des instructions
pour résoudre le problème.
Prise en main
Exo 3
Exo 4
Exercices
Tours de Hanoi
Le jeu consiste à déplacer les n palets de la tour 1 sur la tour 3.

Le jeu n’aurait evidemment guère d’intêret si on ne rajoutait pas
la contrainte suivante : à chaque étape, aucun palet ne doit être
posé sur un palet de diamètre inférieur.
Ecrire une fonction hanoi(n,a,b,c) qui prend comme
arguments : n (nombre de palets sur la tour 1), et a, b, c qui
sont les noms de tours, et qui renvoie la suite des instructions
pour résoudre le problème.
Prise en main
Exo 3
Exo 4
Exercices
Eléments de réponse : Tours de Hanoi
hanoi<-function(n,a,b,c){ if (n==3)
{print(paste("mettre palet de la tour", a, "a la tour",c))
print(paste("mettre palet de la tour", a, "à la tour",b))
print(paste("mettre palet de la tour", c, "à la tour",b))
print(paste("mettre palet de la tour", a, "à la tour",c))
print(paste("mettre palet de la tour", b, "à la tour",a))
print(paste("mettre palet de la tour", b, "à la tour",c))
} else hanoi(n-1,a,c,b)
hanoi(n-1,b,a,c)
}
Tester avec hanoi(3,1,2,3)

Prise en main
Exo 3
Exo 4
Exercices
Eléments de réponse : Tours de Hanoi
hanoi<-function(n,a,b,c){ if (n==3)
{print(paste("mettre palet de la tour", a, "a la tour",c))
print(paste("mettre palet de la tour", a, "à la tour",b))
print(paste("mettre palet de la tour", c, "à la tour",b))
print(paste("mettre palet de la tour", b, "à la tour",a))
print(paste("mettre palet de la tour", b, "à la tour",c))
} else hanoi(n-1,a,c,b)
hanoi(n-1,b,a,c)
}
Tester avec hanoi(3,1,2,3)

Prise en main
Exo 3
Exo 4
Exercices
>hanoi(4,1,2,3)
[1] "mettre palet de la tour 1 à la tour 2"

Rstudio TD 1: Clément Rau Laboratoire de Mathématiques de Toulouse Université Paul Sabatier-IUT GEA Ponsan

Transféré par

Droits d'auteur :

Formats disponibles

Rstudio TD 1: Clément Rau Laboratoire de Mathématiques de Toulouse Université Paul Sabatier-IUT GEA Ponsan

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Rstudio TD 1: Clément Rau Laboratoire de Mathématiques de Toulouse Université Paul Sabatier-IUT GEA Ponsan

Transféré par

Droits d'auteur :

Formats disponibles

Prise en main

Quelques fonctions liées aux statistiques descriptives

Clément Rau Rstudio TD 1

Clément Rau Rstudio TD 1

Clément Rau Rstudio TD 1

Clément Rau Rstudio TD 1

Clément Rau Rstudio TD 1

Clément Rau Rstudio TD 1

Clément Rau Rstudio TD 1

Interface, Fenêtres dans R

Clément Rau Rstudio TD 1

Pour changer de répertoire de travail, saisir dans la console de

setwd("chemin du répertoire souhaité")

Attention, dans les chemins Windows, les \ doivent être

Pour changer de répertoire de travail, saisir dans la console de

setwd("chemin du répertoire souhaité")

Attention, dans les chemins Windows, les \ doivent être

Pour changer de répertoire de travail, saisir dans la console de

setwd("chemin du répertoire souhaité")

Attention, dans les chemins Windows, les \ doivent être

Pour changer de répertoire de travail, saisir dans la console de

setwd("chemin du répertoire souhaité")

Attention, dans les chemins Windows, les \ doivent être

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Essayer qques opérations élémentaires : addition,

Clément Rau Rstudio TD 1

Nombre de décimales souhaitées

où x désigne un réel et n un entier.

Clément Rau Rstudio TD 1

Nombre de décimales souhaitées

où x désigne un réel et n un entier.

Clément Rau Rstudio TD 1

Si l’on veut affecter à x la valeur 3, on peut taper :

! ! ! Attention R distingue majuscule et minuscule ! ! !

Pour affecter plusieurs variables, on les sépare par un ;

Clément Rau Rstudio TD 1

Si l’on veut affecter à x la valeur 3, on peut taper :

! ! ! Attention R distingue majuscule et minuscule ! ! !