nombre degres libertes

[alb]

salut,
on ne peut pas aller au dela de 99 ddl avec student et les autres fonctions stat ?
student_icdf(100,0.95) renvoie undef
les tableurs n'ont pas ce genre de limite.

[parisse]

Peut-etre parce que student_icdf resoud l'equation student_cdf(...)=... alors que dans ce cas il doit etre plus judicieux d'approcher par une loi normale ?

[alb]

c'est une question de precision
sur LibreOfficeCalc la commande TINV(0,1;200) renvoie 1,6525081009
sur Xcas, normal_icdf(0.95) renvoie 1.64485362695

[parisse]

je peux renvoyer l'approximation initiale utilisee dans la methode de Newton avec student_cdf, elle n'est sans doute pas tres bonne, ca donnerait:
student_icdf(100,0.95)
Low accuracy
1.66023333873

student_icdf(200,0.95)
Low accuracy
1.65250797776

[alb]

c'est bizarre qu'on ne puisse pas atteindre la precision des tableurs ?
student_cdf ne marche pas non plus avec plus de 99 ddl

int(student(100,x),x,-inf,1.66023333873) renvoie 0.949999868764 // 1.66023333873 avec la nouvelle approx de Xcas
int(student(100,x),x,-inf,1.6602343261) renvoie 0.949999968296 // 1.6602343261 avec un tableur

int(student(200,x),x,-inf,1.65250797776) renvoie 0.949999987443 // 1.65250797776 avec la nouvelle approx de Xcas
int(student(200,x),x,-inf,1.6525081009) renvoie 0.949999999999 // 1.6525081009 avec un tableur

[parisse]

C'est student_cdf le coupable. Il faut dire que j'utilise des recettes de cuisine que je n'ai jamais regardees en details pour le calculer (vu que ca ne m'interesse guere...), une d'elle faisait un cas particulier pour dof>=100, et cette regle semble foireuse. En la commentant ca a l'air de marcher (au sens ou student_icdf(200,0.95) renvoie un resultat correct).

[alb]

je pense que les limites 30 ou 100 etaient valables quand on utilisait des tables.
Les fonctions stat se trouvent dans quel fichier du source ?

[parisse]

dans moyal.cc
j'ai copie la methode de calcul sur les calculatrices HP, donc c'est surement assez obscur... Il faudrait peut-etre rajeunir tout ca a coup de Abramowitz-Stegun...