Discussion :

[miopie]

bonjour,

voila pour les statisticiens un petit probleme.
j'ai un fichier de mesure qui contient sur 20 cartes environ 200 mesures chacune. mais certaine peut atteindre 4000 mesures
je voudrai recuperer la valeur de P pour chaque ligne afin de savoir si le données par reference suivent une loi normal.
j'ai recuperer une fonction anderson darling mais je ne sais pas comment l'adapter et l'utiliser
voici la fonction

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
RSiteSearch("Anderson-Darling") using R-2.1.0 beta turns up "ad.test" in the
"nortest" package for R on CRAN.  The function looks like:
 
ad.test <- function (x) {
    DNAME <- deparse(substitute(x))
    x <- sort(x[complete.cases(x)])
    n <- length(x)
    if (n < 8) 
        stop("sample size must be greater than 7")
    p <- pnorm((x - mean(x))/sd(x))
    h <- (2 * seq(1:n) - 1) * (log(p) + log(1 - rev(p)))
    A <- -n - mean(h)
    AA <- (1 + 0.75/n + 2.25/n^2) * A
    if (AA < 0.2) {
        pval <- 1 - exp(-13.436 + 101.14 * AA - 223.73 * AA^2)
    }
    else if (AA < 0.34) {
        pval <- 1 - exp(-8.318 + 42.796 * AA - 59.938 * AA^2)
    }
    else if (AA < 0.6) {
        pval <- exp(0.9177 - 4.279 * AA - 1.38 * AA^2)
    }
    else {
        pval <- exp(1.2937 - 5.709 * AA + 0.0186 * AA^2)
    }
    RVAL <- list(statistic = c(A = A), p.value = pval, method =
"Anderson-Darling normality test", 
        data.name = DNAME)
    class(RVAL) <- "htest"
    return(RVAL)
}

je join le fichier ou j'ai recuperer cette fonction avec la facon dont je voudrai arriver a afficher le resultat de p ainsi que le fichier contenant mes mesure a traiter

merci d'avance de votre aide