R TD1

Télécharger au format pptx, pdf ou txt
Télécharger au format pptx, pdf ou txt
Vous êtes sur la page 1sur 19

Chapitre 2

Introduction à la statistique avec R

Représentations Graphiques

1 Pr. Bruno Falissard


Plan
Introduction à la statistique avec R > Représentations Graphiques

• Variables qualitatives
• Variables quantitatives
• Diagramme en bâtons, camemberts,
histogrammes, boîtes à moustaches,
diagrammes cartésiens, diagrammes en fagot

2 Pr. Bruno Falissard


Le fichier smp.c
Introduction à la statistique avec R > Représentations Graphiques

• L’étude santé mentale en prison (smp)


– 799 détenus tirés ausort
– Âge
– Profession
– Dépression, schizophrénie (diagnostic issudu
consensus de deux cliniciens)
– Gravité de la pathologie éventuelle
– Nombre d’enfants

3 Pr. Bruno Falissard


Le fichier smp.c
Introduction à la statistique avec R > Représentations Graphiques

• L’étude santé mentale en prison (smp),variables


évaluant la personnalité des détenus
– Recherche de sensation (rs) : curiosité,attrait pour
le risque et lanouveauté
– Évitement du danger (ed) : timidité,précautionneux
– Dépendance à la récompense (dr) : sensibilitéaux
relations sociales, influençable

4 Pr. Bruno Falissard


Le fichier smp.c
Introduction à la statistique avec R > Représentations Graphiques

> smp.c <- read.csv2("D:/MOOC/Data/smp1.csv")


> str(smp.c)
'data.frame': 799 obs. of 9 variables:
$ age : int 31 49 50 47 23 34 24 52 42 45 ...
$ prof : Factor w/ 8 levels "agriculteur",..: 3 NA 7 6...
$ dep.cons : int 0 0 0 0 1 0 1 0 1 0 ...
$ scz.cons : int 0 0 0 0 0 0 0 0 0 0 ...
$ grav.cons: int 1 2 2 1 2 1 5 1 5 5 ...
$ n.enfant : int 2 7 2 0 1 3 5 2 1 2 ...
$ rs : int 2 2 2 2 2 1 3 2 3 2 ...
$ ed : int 1 2 3 2 2 2 3 2 3 2 ...
$ dr : int 1 1 2 2 2 1 2 2 1 2 ...

5 Pr. Bruno Falissard


Diagramme en bâton
Introduction à la statistique avec R > Représentations Graphiques

> barplot(table(smp.c$prof))
200
150
100
50
0

agriculteur artisan cadre prof.intermédiaireemployé ouvrier autre sans emploi

6 Pr. Bruno Falissard


Diagramme en bâton
Introduction à la statistique avec R > Représentations Graphiques

> barplot(table(smp.c$prof))

> str(smp.c$prof)
Factor w/ 8 levels "agriculteur",..: 3 NA 7 6 8 6 3 2 6 6 ...
> table(smp.c$prof)

agriculteur artisan autre


6 90 31
cadre employé ouvrier
24 135 227
prof.intermédiaire sans emploi
58 222

7 Pr. Bruno Falissard


Diagramme en bâton
Introduction à la statistique avec R > Représentations Graphiques

> barplot(table(smp.c$prof))
200
150
100
50
0

agriculteur artisan cadre prof.intermédiaireemployé ouvrier autre sans emploi

8 Pr. Bruno Falissard


Camembert
Introduction à la statistique avec R > Représentations Graphiques

> pie(table(smp.c$prof))

employé
cadre
autre

artisan

agriculteur
ouvrier

sans emploi
prof.intermédiaire

9 Pr. Bruno Falissard


Histogramme
Introduction à la statistique avec R > Représentations Graphiques

> hist(smp.c$age)

120
100 Histogram of smp.c$age
Frequency

80
60
40
20
0

20 30 40 50 60 70 80

smp.c$age

10 Pr. Bruno Falissard


Histogramme
Introduction à la statistique avec R > Représentations Graphiques

> hist(smp.c$age,col="grey",main="",
xlab="age")

120
100
Frequency

80
60
40
20
0

20 30 40 50 60 70 80

age

11 Pr. Bruno Falissard


Boîte à moustache
Introduction à la statistique avec R > Représentations Graphiques

> boxplot(smp.c$age,xlab="age")

80
70

25% desdonnées
60
50

25% desdonnées
40

25% desdonnées
30

25% desdonnées
20

age

12 Pr. Bruno Falissard


Boîte à moustache
Introduction à la statistique avec R > Représentations Graphiques

> boxplot(smp.c$age~smp.c$rs,ylab="age",
xlab="Recherche de sensation")

80
70
60
age

50
40
30
20

1 2 3

Recherche de sensation

13 Pr. Bruno Falissard


Diagramme cartésien (« ou diagramme en x,y »)
Introduction à la statistique avec R > Représentations Graphiques

> plot(smp.c$age,smp.c$n.enfant)

14 Pr. Bruno Falissard


Diagramme cartésien (« ou diagramme en x,y »)
Introduction à la statistique avec R > Représentations Graphiques

> plot(jitter(smp.c$age),
jitter(smp.c$n.enfant))

15 Pr. Bruno Falissard


Diagramme temporel
Introduction à la statistique avec R > Représentations Graphiques

> repdat <- read.csv2("C:/Users/Bruno/outils hdrs.csv")


> str(repdat)
'data.frame': 1053 obs. of 3 variables:
$ NUMERO: int 96 96 96 96 96 96 96 96 157 157 ...
$ VISIT : int 0 4 7 14 21 28 42 56 0 4 ...
25
$ HDRS : int 34 26 12 7 5 1 1 1 27 19 ...
repdat$HDRS

20
15
10

n=146 n=139 n=137 n=133 n=129 n=127 n=122 n=120

0 4 7 14 21 28 42 56

repdat$VISIT

16 Pr. Bruno Falissard


Diagramme temporel
Introduction à la statistique avec R > Représentations Graphiques

> library(gplots)
> plotmeans(repdat$HDRS~repdat$VISIT,gap=0,
barcol="black")
25
repdat$HDRS

20
15
10

n=146 n=139 n=137 n=133 n=129 n=127 n=122 n=120

0 4 7 14 21 28 42 56

repdat$VISIT

17 Pr. Bruno Falissard


Diagramme en fagot
Introduction à la statistique avec R > Représentations Graphiques

> interaction.plot(repdat$VISIT,repdat$NUMERO,
repdat$HDRS,lty=1,legend=FALSE)
mean of repdat$HDRS

40
30
20
10
0

0 4 7 14 21 28 42 56

repdat$VISIT

18 Pr. Bruno Falissard


Conclusion
Introduction à la statistique avec R > Représentations Graphiques

smp.c <- read.csv2("D:/MOOC/Data/smp1.csv")


str(smp.c)
barplot(table(smp.c$prof))
pie(table(smp.c$prof))
hist(smp.c$age)
hist(smp.c$age,col="grey", main="",xlab="age")
boxplot(smp.c$age,xlab="age")
boxplot(smp.c$age~smp.c$rs,ylab="age",xlab=
"Recherche de sensation")
plot(smp.c$age,smp.c$n.enfant)
plot(jitter(smp.c$age),jitter(smp.c$n.enfant))
repdat <- read.csv2("D:/MOOC/Data/hdrs.csv")
str(repdat)
library(gplots)
plotmeans(repdat$HDRS~repdat$VISIT,gap=0,
barcol="black")
interaction.plot(repdat$VISIT,repdat$NUMERO,
repdat$HDRS,lty=1,legend=FALSE)

19 Pr. Bruno Falissard

Vous aimerez peut-être aussi