Simulation d'une variable de pareto

Fermé
bedalfin - 9 avril 2009 à 01:59
 karatetiger - 19 janv. 2011 à 11:37
Bonjour,
Est-ce que quelqu'un sait comment simuler une loi de pareto sous le logiciel R
merci.
A voir également:

3 réponses

Salut

Je me pose exactement la même question mais en Matlab. Enfin le principe est le même. Tu peux me le dire si tu as la réponse. Merci!
0
/**
* Return a real number with a Pareto distribution with parameter alpha.
*/
public static double pareto(double alpha) {
return Math.pow(1 - uniform(), -1.0/alpha) - 1.0;
}
0
##inverse de la fonction de répartition d'une pareto:

invfdrpareto=function(y,l,a)
{x=l*(1/(1-y)^(1/a)-1)
x}


##génération de nombres suivant la loi de Pareto :

rpareto=function(l,a)
{f=function(n=1)
{U=runif(n)
Y=sapply(U,invfdrpareto,l,a)
Y}
f}

##Exemple d'exécution:

V=rpareto(1.4,1.5e6)(1000)
plot(V)
hist(V,prob=TRUE)
0
Autre méthode qui va plus vite, en remarquant que Y=T*exp(X)-T suit une pareto (a,T) si X suit une exponentielle(a):

rparexp=function(a,T)
{f=function(N=1)
{V=T*exp(rexp(N,a))-T
V}
f}
0
Bonjour je voudrais savoir si je dois simuler des nombre suivant une loi de pareto comment choisir mes paramètres? Par exmple toi tu prend 1.4 et 1.5e6 au hasard ou cela vient de quelques part?

Merci
0