Bonjour,

Un de mes algorithmes essaie différentes valeurs flottantes entre x et y. Il s'arrête dès qu'une valeur tirée au sort est acceptée, sachant que l'on a aucun a priori sur les valeurs en question, on les choisit donc aléatoirement. S'il échoue avec une valeur tirée au sort, il faut essayer une autre valeur dans cet intervalle, en diminuant la probabilité des valeurs environnantes (par exemple selon une gaussienne).

Ainsi j'obtiens une densité à trous, par exemple si j'ai déjà testé les valeurs z et t sans succès :

Il y a-t-il une fonction simple de scipy ou autre pour tirer au sort de cette manière ?

L'idéal serait une fonction du genre 

rand(x, y, [[moyenne1, variance1], [moyenne2, variance2], ...])

qui renvoie un flottant antre x et y en évitant les gaussiennes fournies en paramètre par exemple.

Je précise que je voudrais éviter de discrétiser mon intervalle [x, y] et que je suis pas un spécialiste des probas, je cherche juste comment tirer au sort mes valeurs de manière efficace

Merci

-
Edité par courrier 21 janvier 2015 à 19:16:23

Salut,

Je pense que c'est impossible avec SciPy. Mais tu peux d'inspirer de cet ancien message : le tirage est uniforme sur plusieurs intervalles. A toi de l'adapter avec des gaussiennes.

Salut,

Je précise que je voudrais éviter de discrétiser mon intervalle [x, y]

Dans ce cas là, laisse tomber les ordi. Le calcul continu, c'est pas vraiment leur fort.

Je ne sais pas si ça peut répondre à ton besoin mais tu peux éliminer des valeurs comprises entre des intervalles en utilisant random.choice().

nombres = range(a,b) + range(c,d)
x = random.choice(nombres)

Merci pour vos suggestions, mais j'ai fini par remplacer le tirage aléatoire par un parcours avec un pas qui diminue, du style, pour a=0 et b= 100 : 0, 100, 50, 25, 75, 12.5, 37.5, 62.5, 87.5, ....

J'ai donc contourné le problème mais c'est satisfaisant pour mon utilisation.