Bonjour à tous, voici mon problème.
Je dispose d'une densité de probabilité (x+1/2 sur [0,1] par exemple). Je dois générer un échantillon de taille n de réalisations indépendantes de cette loi à l'aide de mon brave PC. Comment dois-je m'y prendre ?
Salut,
ça dépend des fonctions dont tu disposes dans ton langage de programmation favori... Il y a quasiment toujours des bibliothèques de génération de nombres aléatoires.
Cordialement.
En général on dispose au moins de la distribution uniforme sur [0,1]. A partir de là il suffit de calculer la densité cumulative,( x(x+1)/2 dans ton cas), de l'inverser (1/2(sqrt(y²+1) - 1) dans ton cas) et d'appliquer à un tirage uniforme sur [0,1] cet inverse.Envoyé par GuYem
Cordialement,
Merci beaucoup mmy et martini, c'est exactement ce que je cherchais
Salut Guyem
J ai passe mon vendredi a me depatouiller dans la librairie C GNU (GSL).
- il y a un generateur de nombre aleatoire de loi uniforme
- puis tu peux utilier ce generateur pour simuler ta loi.
Le probleme est qu il n existe pas de vrai generateur de nombre aleatoire, ce qui est genere est une suite qui, si examiner par quelqu un, semble tres aleatoire. Or, c est une serie recursive, il faut donc l initialiser. Et la, faut faire attention. Utilise par exemple l horloge du PC
Je t ai parle de GSL car ca prends 30 secondes a faire un prgm, en C, qui marche avec.
++
Bonjour,
Cela n'a pas que des inconvénients. Pour le debug ou du test de non-régression, il est plus intéressant de l'initialiser par une valeur fixe, le programme est alors parfaitement répétitif!
Cordialement,
Oups, je suis perdu là (C GNU GSL, horloge, récursif ?)
Mon problème est de tester différentes procédures de choix d'un histogramme approchant une densité inconnue sur [0,1].
Est-ce-que, pour ce problème précis, la manière de générer l'échantillon peut avoir un impact important ?
Pour le moment j'ai simplement utilisé la procédure fournie (dans scilab, le logiciel que j'utilise) pour générer un échantillon uniforme et j'ai transformé cet échantillon en celui de ma densité (x+1/2) avec la technique décrite par mmy.
je ne voudrais pas jouer les rabat-joies, mais c'est un problème bien connu (voir le bouquin de Bosc et Lecoutre chez Economica, ou les articles d'Abou-Jaoude, Deheuvels, etc.)
A priori non; aucun générateur pseudo-aléatoire que l'on trouve usuellement n'est loin de la perfection, même si aucun ne l'atteint!
Cordialement,
Merci ambrosio de ces références, je regarde ce problème depuis un moment maintenant et je l'aborde par les critères d'informations, je ne sais pas si tu connais ces trucs.
Je ne connaissais cependant pas les références que tu cites et j'y jetterai surement un oeil.
mmy, ta réponse me rassure, je n'ai pas envie de chercher trop loin dans le problème de la génération de nombres aléatoires, la technique que tu m'as donnée donne des résultats satisfaisants pour l'instant donc je la garde.
je ne sais pas si cela répond à ton problème, mais il existe aussi des sites Web qui fournissent ce genre de service, par exemple:
générateur quantique:
http://www.randomnumbers.info/
à partir de bruit atmosphérique:
http://www.random.org/
Salut
Scilab fait effectivement les choses bien. Mais je t ai donne ce que j utilise en C (pouraller plus vite). Et hela, il faut aller dans les details pour cette librairie. Dsl si c etait confus
++
beaucoup des travaux des statisticiens français concernent le choix de la taille des intervalles (ou du nombre d'intervalles) de l'histogramme qui assure la convergence vers la densité vraie. Comme il y a pas mal de sortes de convergences, le sujet est assez touffu. Je pense que le Bosc et Lecoutre est une bonne introduction. Autrement, la référence absolue, du moins il y a 20 ans quand je m'intéressais à ça, était le Prakasa Rao: Non Parametric Functional Estimation, Academic Press, Orlando, Florida (1983). xiv + 522 pp.
voir aussi:
http://www.isid.ac.in/~statmath/homepageblsp/index.html
