Discussion :

[Djamely]

Assalam alaikom,

Bonjour,

Cher collèges,

Je souhaite avoir la courbe de densité de loi Exponentielle sur la plateforme de ggplot2, j'ai utilisé le suivant:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
library("ggplot2")
p <- ggplot(mpg, aes(displ, hwy))
p+ labs(x = "x", y = "",title = "densité de loi Exponentiell")+
     xlim(0, 5)+
     ylim(0,1)+

ici je connais pas comment ajouter cela:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

curve(dexp(x), xlim=c(0,5),ylim=c(0,1), ylab="", col="red", lwd=2, lty=1)

Merci d'avance pour votre repose.

[VonVelten]

Bonjour,

Vu que sa fonction de densité est connue, il est aussi simple de calculer les points qui la caractérisent selon son paramètre lambda

Code R :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
library("ggplot2")
Lambda<-1
x<-seq(from =0, to=5,length.out = 1000)
y<-1/Lambda*exp(-x/Lambda)
data<-data.frame(Lambda,x,y)
 
p <- ggplot(data, aes(x=x,y=y))+
  labs(x = "x", y = "",title = "densité de loi Exponentiell")+geom_line()
print(p)

VV

[Djamely]

Merci bien,
Monsieur, Est-ce que on peut d'ajouter dans le même graphe la courbe de la densité de la loi paréto la loi normale,....
Si oui comment faire?
encore une fois merci.

VV

[VonVelten]

C'est assez simple, voilà un petit exemple

Code R :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
 
library("ggplot2")
x<-seq(from =0, to=5,length.out = 1000)
 
#Loi Exponentielle
Lambda<-1
y.expo<-1/Lambda*exp(-x/Lambda)
 
#Loi Normale
Moyenne<-0
Sd<-1
y.normale<-1/(pi*sqrt(2*Sd))*exp(-((x-Moyenne)/Sd)^2/2)
 
 
 
 
#Loi de Pareto
library(rmutil)
Pareto<-function(x,xm=2,k=1){
  if(x>=xm){
    return(k*xm^k/x^(k+1))
  }else{
    return(0)
  }
}
y.pareto<-sapply(x,Pareto,k=1,xm=2)
 
 
Loi<-rep(c("Exponentielle","Normale","Pareto"),each=1000)
y<-c(y.expo,y.normale,y.pareto)
data<-data.frame(Lambda,x,y,Loi)
 
p <- ggplot(data, aes(x=x,y=y,color=Loi))+
  labs(x = "x", y = "",title = "densité")+geom_line()
print(p)

[Djamely]

C'est magnifique, Merci beaucoup

Le dernière question : est-ce qu'on peut centrer le titre "densité" ?