Table du Khi2

[Le chat]

Bonjour,
On a toujours pas 'sous la main' une table du Khi2. Ces quelques lignes de programme vous permettrons de calculer le seuil de confiance pour un degré de liberté et un seuil de tolérance donnés:
Nota: 'SeuilTolerance' doit être saisi sous forme d'un pourcentage (ex: khi2(5,1) ... 1 est sous-entendu 1%)

khi2(DegreLiberte,SeuilTolerance){
local khiX,SeuilConfiance,pasX;
pasX:=1/300;
khiX:=-pasX;
while (SeuilConfiance*100 < 100-SeuilTolerance) {
khiX+=pasX;
SeuilConfiance:=chisquare_cdf(DegreLiberte,khiX);}
return("Le seuil de confiance est :" + round(evalf(khiX),2));
}:;

Merci.

[alb]

bonjour
khi2(5,1) et chisquare_icdf(5,0.99) donne la même réponse c'est bien çà ?

[Le chat]

Bonjour,
Oui, exactement: khi2(5,1)=chisquare_icdf(5,0.99)
Alors pourquoi mon programme ? Uniquement pour inciter les personnes non initiées à la programmation de s'y intéresser.

Nota: dans mon programme, la ligne:
SeuilConfiance:=chisquare_cdf(DegreLiberte,khiX);
pourrait être remplacée par la ligne ci-dessous:

SeuilConfiance:=integration(k*e^(-x/2)*x^(DegreLiberte/2-1),x,0,exact(khiX));
avec k:=evalf(1/(2^(DegreLiberte/2)*Gamma(DegreLiberte/2)));
Toutefois, la vitesse d'éxécution du programme serait ralentie dans ce cas.

[alb]

oui j'ai bien compris que le but était d'écrire un programme
Tu fais ça avec des élèves ou des étudiants ?