Calcul p-valeur test de Shapiro-Wilk

Bonjour,

J'essaie d'adapter l'algorithme AS 181.2 disponible ici en python: https://sci2s.ugr.es/keel/pdf/algorithm/articulo/royston1982.pdf ou encore là en C: https://fossies.org/dox/qtiplot-0.9.8.9/ShapiroWilkTest_8cpp_source.html et qui permet de déterminer le polynôme algébrique d'ordre d-1 et basé sur les coefficients de Shapiro-Wilk.

Voilà ma fonction:

Code:


poly = function(cc,nord,x) {
 ret_val = cc[1]
 if (nord == 2) {ret_val = ret_val + x * cc[2]}
 if (nord > 2) {
 p = x * cc[nord]
 for (j in (nord-1):2) {
 p = (p + cc[j]) * x
 ret_val = ret_val + p
 }
 }
 return(ret_val)
}


Et ici un exemple pour lequel je devrais trouver la même solution:

Code:


n = 100
c5 = c(-1.5861,-0.31082,-0.083751,0.0038915)
c6 = c(-0.4803,-0.082676,0.0030302)
T = shapiro.test(runif(n,0,1))
# La solution selon R
W = T$statistic
y = log(1 - W)
m = poly(c5,4,log(n))
s = exp(poly(c6,3,log(n)))
# Le calcul manuel
pnorm(y,mean=m,sd=s,lower.tail=FALSE,log.p=FALSE)


Sauf que non, bon on est plus ou moins sur le même ordre mais c'est pas ça.

Si quelqu'un a déjà codé la fonction poly je suis preneur afin de voir si j'arrive enfin aux mêmes résultats que R.