Simulation d'une variable de pareto
bedalfin
-
karatetiger -
karatetiger -
Bonjour,
Est-ce que quelqu'un sait comment simuler une loi de pareto sous le logiciel R
merci.
Est-ce que quelqu'un sait comment simuler une loi de pareto sous le logiciel R
merci.
Configuration: Windows XP Internet Explorer 7.0
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! -
/**
* 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;
} -
##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)