Discussion :

[lilly74]

Bonjour,
j'essaie de simuler une distribution asymétrique.
J'ai tout d'abord pensé à la loi log-normale:

Code :

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

X1 <- exp(rnorm(1000, 0,3))

mais il y a beaucoup trop de valeurs extrêmes comme le montre:

Code :

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

 boxplot(X1)

J'ai trouvé un article qui parle de "distribution asymétrique générée par une base 1.4 mise à des exposants tirés aléatoirement d'une distribution normale".
Quelqu'un aurait-il la gentillesse de m'expliquer comment faire cela svp?
Merci d'avance pour votre aide

[lilly74]

Re-bonjour,
en faisant quelque chose du style:

Code :

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

X1 <- (rnorm(1000, 0,30))^(ifelse(runif(1)>0.1,1.4,1)

cela donne:

alors que j'aimerais quelque chose de plus concentré sur un côté comme:

Une idée svp?
Merci d'avance

[ol9245]

et comme ça ?

Code :

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

1
2
x = 1.4^rnorm(1000, 0,1)
boxplot(x)

[lilly74]

et comme ça ?

Code :

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

1
2
x = 1.4^rnorm(1000, 0,1)
boxplot(x)

Merci pour cette réponse
Il me semble que cela correspond à l'idée mentionnée dans l'article, mais je trouve que ce n'est pas si "concentré" que cela
J'aurais vraiment obtenir un boxplot très asymétrique
L'idée me paraît bien mais quand je compare avec la boxplot que j'ai mise en exemple, la vôtre est beaucoup plus centrée...
Peut-être ai-je indiquée une mauvaise piste

[ol9245]

regarde du côté de la loi gamma. C'est une loi qui, selon les paramètres, peut être très asymétrique. De plus c'est une des grandes lois de statistiques donc ton choix sera plus facile à documenter lorsque tu auras à expliquer ton travail.

[lilly74]

regarde du côté de la loi gamma. C'est une loi qui, selon les paramètres, peut être très asymétrique. De plus c'est une des grandes lois de statistiques donc ton choix sera plus facile à documenter lorsque tu auras à expliquer ton travail.

Ok, merci
Je regarde cela et je vous tiens au courant