MATH_LoiKhiDeux    [Statistiques]

Par Teddy Linet, Philip Burns (Math4D v2)
Nouvelle recherche
Si (Faux)
   ` MATH_LoiKhiDeux [ Teddy LINET 02/03 et Philip BURNS 1985]
   `-------------------------------------------------
   `ATTENTION : nécessite MATH_IntegraleGamma
   `-------------------------------------------------
   ` RAPPELS MATHEMATIQUES :
   ` Evalue la probabilité d'une distribution selon la fonction KhiDeux
   ` à N degrés de liberté
   ` 
   ` La fonction Khi2 à ddl degrés de liberté est
   `            x^((ddl-2)/2)*exp(-x/2)
   ` f(x)=----------------------------------
   `             2^(ddl/2)*gamma(nu/2)
   ` ou gamma est l'intégrale sur x de 0 à + l'infini de la fonction
   ` gamma(x)=x^(z-1)*exp(-x) (cf fonction gamma)
   ` 
   ` La distribution du Khi2 sert pour le test statistique du même nom
   `-------------------------------------------------
   ` MATH_LoiKhiDeux ($vX_F;ddl)-->Réel
   ` $vX_F (Reel) valeur du Khi Carré
   ` Ddl (Reel)  Degrés de liberté
   ` Du fait de l'utilisation de la fonction gamma on peu utiliser des 
   ` degrés de libertés non entier (pas banal non )
   ` ---------------------------------------------
   ` MATHERROR: idem renvoyé par la fonction MATH_IntegraleGamma
   ` +1 -> Ddl non entier (alors qu'il aurait du)
Fin de si 

  ` Déclaration des "constantes"
C_ENTIER($vMaxIter_I;$vDprec_I)
C_ENTIER LONG(MATHERROR)
C_ENTIER($vIter_I)
C_REEL($vCprec_F;$p_F;$vX_F)
C_REEL($vddl_F;$0;$1;$2)
$vMaxIter_I:=200
$vDprec_I:=12
$vX_F:=$1
$vddl_F:=$2


$p_F:=1-MATH_IntegraleGamma ($vX_F/2;$vddl_F/2;$vDprec_I;$vMaxIter_I)

Au cas ou 
 : (MATHERROR<0)
 $p_F:=-1
 : ((Ent($vddl_F)-$vddl_F)#0)
 MATHERROR:=1
Fin de cas 

$0:=$p_F