Simulation selon la fonction de répartition

Bonsoir à tous,


Imaginons une grandeur quantitative aléatoire "X" (à valeurs dans R+). Soit l'intervalle [0, 30 000 000] découpé en 30 tranches égales d'1 million, je noterai "x" ce vecteur de taille 30 ( càd : x = (1M , 2M , .....30M)). J'ai aussi un autre vecteur de même taille 30, noté F_X(x) = P(X<x), ce vecteur contient les probabilités que la grandeur "X" soit inférieure à chacune des valeurs de "x".

Je voudrais être capable de simuler des valeurs dans la loi de X, et je n'ai à disposition que les 30 valeurs de "x" et "F_X(x)". Comment faire ?

Je suis à priori dans un contexte où les techniques d'inversion de la fonction de répartition sont inutilisables car je n'ai pas l'expression littérale de F_X(x), il en va de même pour la méthode de rejet-acceptation. Je voyais comme troisième piste d'ajuster des loi usuelles, sauf qu'à ma connaissance les ajustements du type "Maximum de vraisemblance" s'effectuent lorsqu'on dispose d'échantillons de x, et non pas lorsqu'on ne dispose que de F_X(x). Ainsi, que reste-il comme piste ?

Merci d'avance.