problème de statistique

[thomas]

Bonjour,

Mon problème est le suivant:
Sur Excel existe une fonction de statistique nommée "KHIDEUX.INVERSE()" qui, à partir de la probabilité et du nombre de degrés de liberté n du système, renvoie la valeur d'une variable aléatoire suivant une loi du khi² à n degré de liberté.

Or, je dois programmer sur Visual Basic une fonction analogue. Mais, devant faire face à des degrés de liberté très grand (500), et cette fonction impliquant des exponentielles, mes variables n'ont pas assez de capacité pour contenir de tels nombres.

La fonction d'Excel citée ci-dessus fonctionne, en revanche, parfaitement dans ce cas de figure.
Bref, je cherche désespérément l'algorithme de cette fonction d'Excel ou, ne serait-ce que le calcul sur lequel il est basé.

Tout renseignement à ce sujet est susceptible de m'aider.

Merci.

 

[Zon]

Salut,

En VBA tu peux appeler une formule comme ceci:

function Khi2%(Proba%,Degre%)
khi2=application.chiinv(proba,degre)
end function

A+++

 

[albert]

bonjour thomas, Zon, le forum,

un peu tardivement, mais peut-être à temps, j'apporte ma contribution :
<http://wwwpsy.univ-bpclermont.fr/~meot/FicExcelZippes&pdf/TestAdeqIndEtNP.pdf>

il y a une référence à
http://perso.club-internet.fr/hdelboy/Nonparam.htm

les outils aideront peut-être à la programmation

 

[thomas]

Salut messieurs,

Je vous remercie pour votre participation.

En fait, ce problème (je m'en suis rendu compte au bout de quatre jour de boulot), posé tel quel, n'était pas possible à résoudre par rapport aux équations brutes définissant le test du chi².

Par chance, j'ai trouvé un code source en javascript de ce calcul, mais basé sur des approximations polynomiales. J'avoue que je n'aurais pas pu trouvé ça tout seul, d'ailleurs, même mes profs (je suis en école d'ingénieur) étaient largués à ce sujet.

Si ça vous intéresse, l'adresse du site dont j'ai pompé le code est le suivant :

http://www.fourmilab.ch/rpkp/experiments/analysis/chiCalc.html

Ensuite, je n'ai eu qu'à traduire du javascript au langage VB.

Merci encore et à la prochaine,

Tom.

 

[albert]

merci thomas, pour le code en java ,bonne occasion pour s'y mettre
...pour ma part, je ne saurais traduire en vba

albert

 

[thomas]

Salut, Albert,

Je t'envoie ma traduction du code source en javascript en vb.
Ca vaut ce que ça vaut, je ne suis pas spécialiste de ces deux langages, mais je l'ai testé et ça marche.

De plus, cette page contient mon codage pour la loi normale centrée réduite.

Amuse-toi bien!

Tom.