Calcul p-valeur test de Shapiro-Wilk

Aller en bas

test - Calcul p-valeur test de Shapiro-Wilk Empty Calcul p-valeur test de Shapiro-Wilk

Message par joyeux_lapin13 le Lun 5 Aoû 2019 - 11:49

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.
joyeux_lapin13
joyeux_lapin13

Nombre de messages : 1925
Age : 36
Localisation : Mayotte
Date d'inscription : 21/04/2010

https://lemakistatheux.wordpress.com/

Revenir en haut Aller en bas

Revenir en haut

- Sujets similaires

 
Permission de ce forum:
Vous ne pouvez pas répondre aux sujets dans ce forum