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