Importation des données 

fileName="https://pmbo-sd.fr/STID/Outils_pilotage_2/paca_demog_aire.txt"
X=read.delim(fileName,header=T,sep=";",dec=",")
attach(X)
head(X)
##   code_dep       nom_departement population_totale population_0_19_ans
## 1        4 Alpes-de-Hte-Provence            161664               34245
## 2        5       Alpes (Hautes-)            141576               31059
## 3        6       Alpes-Maritimes           1080614              229665
## 4       13      Bouches-du-Rhone           2043941              489732
## 5       83                   Var           1071832              225121
## 6       84              Vaucluse            566096              135393
##   naissances_2004_2015 superficie
## 1                18577    6925.22
## 2                17107    5548.68
## 3               143486    4298.58
## 4               307980    5087.49
## 5               131517    5972.54
## 6                82920    3567.13

Courbes à 1 variable

bases

 options                               Détail           
 
 pch                                  forme du point
 cex                                  taille du point (*constant expand*), 1 par défaut
 col                                  couleur du point
 type                                 type de courbe
 lwd                                  épaisseur de courbe
par(bg="slategray3",fg="cyan")
plot(population_0_19_ans,pch=15,cex=2,main="Population - de 20 ans",col.main="darkblue")

text(1:6,population_0_19_ans,nom_departement,pos=3,xpd=NA,col="yellow")

pour voir les formes …

plot(1:25,rep(1,25), pch=1:25,cex=2,xaxt="n",yaxt="n",xlab="",ylab="")
text(1:25,rep(1,25), label=1:25,pos="3",cex=0.8)

Avec des formes diverses et des effets …

par(mfrow=c(2,3),bg="slategray3",fg="cyan")
plot(population_0_19_ans,pch=16,cex=2,xpd=NA,main="Population - de 20 ans",col.main="darkblue")
plot(population_0_19_ans,pch=17,cex=2,xpd=NA,main="Population - de 20 ans",col.main="darkblue")
plot(population_0_19_ans,pch=18,cex=2,xpd=NA,main="Population - de 20 ans",col.main="darkblue")
plot(population_0_19_ans,pch=19,cex=2,xpd=NA,main="Population - de 20 ans",col.main="darkblue")
plot(population_0_19_ans,pch="a",cex=2,xpd=NA,main="Population - de 20 ans",col.main="darkblue")
plot(population_0_19_ans,pch="S",cex=2,xpd=NA,main="Population - de 20 ans",col.main="darkblue")

# plot de la population 0-19_ans (type lignes, et couleur)
par(mfrow=c(1,3),bg="darkblue",fg="cyan",col.axis="white")
plot(population_0_19_ans,pch=15,cex=2,type="l",col="yellow",
     main="Population - de 20 ans",col.main="blue")
plot(population_0_19_ans,pch=15,cex=2,type="l",col="yellow",lwd=2,
     main="Population - de 20 ans",col.main="blue")
# plot de la population_0_19_ans (lignes en pointille)
plot(population_0_19_ans,pch=15,cex=2,type="l",col="yellow",lwd=2,lty=2,
     main="Population - de 20 ans",col.main="blue")

les types p , l , h , b , s sur un même ecran (utilisation d’une boucle)

par(bg="slategray1",fg="cyan",mfrow=c(2,3))
for (t in c("p","l","h","b","s"))
{
plot(population_0_19_ans,cex=2,type=t,col="red",lwd=3,xpd=NA,xlab="population_0_19_ans",main=t)
text(1:6,population_0_19_ans,nom_departement,pos=3,xpd=NA,col="darkblue")
}

plot de la population 0-19 ans avec une couleur différente par point

rainbow(6)
## [1] "#FF0000FF" "#FFFF00FF" "#00FF00FF" "#00FFFFFF" "#0000FFFF" "#FF00FFFF"
par(bg="slategray3",fg="cyan",mfrow=c(1,5))
par(mfrow=c(1,1))
plot(population_0_19_ans,pch=15,cex=2,col=rainbow(6),lwd=2,lty=2,
     main="Population - de 20 ans",col.main="blue")
text(1:6,population_0_19_ans,nom_departement,pos=3,xpd=NA,col=rainbow(6))

ajout de points

# un plot des naissances ...
par(bg="slategray3",fg="cyan")
plot(naissances_2004_2015,pch=19,
     main="Naissances 2004 à 2015")

# ... l'ajout de points
points(naissances_2004_2015 + 20000,pch=15,col="yellow",xpd=NA)


ici, on a ajouté des points d’ordonnée naissances + 20 000 (d’une autre couleur)

Nuages de points (courbes à 2 variables)

superficie en abscisses, population en ordonnées

par(bg="slategray3",fg="yellow")
plot(superficie,population_totale)
text(superficie,population_totale,nom_departement,pos=3,xpd=NA)

Courbes mathématiques


Tracés successifs de :

  • la courbe de la fonction x -> x^2-3x-4
  • la meme sur l’intervalle [-2,5]
  • avec les axes des abscisses et des ordonnées
  • avec la droite x=3 en couleur et y=4 en couleur et pointillés
par(bg="royalblue3",fg="yellow",mfrow=c(2,2),col.axis="white")
curve(x^2-3*x-4)

curve(x^2-3*x-4 , -2 , 5 )

curve(x^2-3*x-4 , -2 , 5 )
abline(h=0);abline(v=0)

curve(x^2-3*x-4 , -2 , 5 )
abline(h=3,col="orange")
abline(v=4,col="orange",lty=2)


Successivement :

  • courbe de l’exponentielle entre -3 et 6 (avec axes abscisses / ordonnées)
  • logarithme de 0 a 10
  • courbe de y=x, de -2 a 2 sur repere normé (mêmes échelles abscisses / ordonnées)
par(bg="royalblue3",fg="yellow",mfrow=c(1,3),col.axis="white")
curve(exp(x),-3,6,main="exp(x)",col.main="white")
abline(h=0);abline(v=0)

curve(log(x),0,10,main="log(x)",col.main="white")
abline(h=0);abline(v=0)

curve(x^1,-2,2,asp=1,main="f(x)=x",col.main="white")
abline(h=0);abline(v=0)


Les courbes de y = ax, pour a dans {-2 , -1.5 , -1 , -0.5 , 0 , 0.5 , 1 , 1.5 , 2} et, sur l’axe d’abscisse 1, les coefficients (pentes) des droites :

par(bg="royalblue3",fg="yellow",col.axis="white")
curve(x^1,-2,2,asp=1)
abline(h=0);abline(v=0)
for (a in seq(-2,2,0.5)) curve(a*x,add=T)
text(1,seq(-2,2,0.5),seq(-2,2,0.5))


les courbes des puissances 2 a 5 de x ainsi que des racines carrée a cinquième

par(bg="royalblue3",fg="yellow",col.axis="white")
curve(x^1,0,5);abline(h=0);abline(v=0)
for (i in 2:5) curve(x^i,add=T)
for (i in 2:5) curve(x^(1/i),add=T)


Réafficher les courbes des puissances positives et negatives précédentes en différentes couleurs

par(bg="slategray4",fg="black",col.axis="white")
curve(x^1,0,5,col="white",lwd=2)
abline(h=0);abline(v=0)

for (i in 2:5) curve(x^i,add=T,col=rainbow(10)[i],lwd=2)
for (i in 2:5) curve(x^(1/i),add=T,col=rainbow(10)[i+4],lwd=2)


Courbes paramétrées : en fonction de t (comme “temps”), un point se déplace dans le plan. On obtient une trajectoire …

par(bg="slategray3",fg="white",col.axis="white")
par(mfrow=c(1,2))
t=seq(0,20,1)
plot(t*cos(t),t*sin(t),type="b",asp=1)

t=seq(0,20,0.1)
plot(t*cos(t),t*sin(t),type="b",asp=1)


afficher la courbe de la loi normale de parametres moyenne m = 0 et écart-type s = 1, entre -5 et 5

curve(dnorm(x),-5,5)
abline(v=0,col="red",lwd=2)


idem pour m=7 et s= 2, entre m-4s et m+4s

m=7;s= 2
curve(dnorm(x,m,s),m-4*s,m+4*s,lwd=2)
abline(v=m,col="red",lwd=2)


ajouter les courbes des normales de moyenne 5 a 9 (même écart-type s) en couleurs distinctes

m=7;s= 2
curve(dnorm(x,m,s),m-4*s,m+4*s,lwd=2)
abline(v=m,col="red",lwd=2)
for (m in 5:9) curve(dnorm(x,m,s),add=T,col=rainbow(5)[m-4])


afficher les courbes des normales de moyenne 0 et ecart-type 1 a 2 par pas de 0.1 entre -4 et 4

par(bg="royalblue3",fg="yellow",col.axis="white")
curve(dnorm(x),-4,4)
for (i in 1:10) curve(dnorm(x,0,1+i/10),add=T,col=rainbow(10)[i])