¿Podemos calcular el valor P de la prueba de Anderson Darling con el estadístico AD sin utilizar las tablas dadas?

Question

Estoy ajustando mis datos experimentales con diferentes distribuciones, estoy calculando el estadístico Anderson Darling para mis datos y las distribuciones teóricas. Quiero calcular el valor P a partir del estadístico de Anderson Darling sin utilizar las tablas, ¿cómo puedo calcular el valor P?

Answer 1

Si no quiere interpolar a partir de tablas precalculadas, puede hacer una simulación basada en bootstrap. Pruebe con adSim ( https://cran.r-project.org/package=qualityTools )

Utilicemos la distribución normal en nuestro ejemplo. Tendrás que cambiar la cadena por otra cosa si no quieres probar la normalidad.

Interpolación

adSim(x, "normal", NA)

$distribution
[1] "normal"

$parameter_estimation
 mean        sd 
32.450662  1.717755 

$Anderson_Darling
[1] 0.5028619

$p_value
[1] 0.1865325

Los valores críticos de la tabla son los percentiles 75%, 90%, 95%, 97,5% y 99% (disponibles en el código fuente https://github.com/cran/qualityTools/blob/master/R/adSim.R ).

Simulación de Boostrap

adSim(x, "normal", 10000)

... simulating the Anderson-Darling distribution by 10000 bootstraps for normal distribution... 

$distribution
[1] "normal"

$parameter_estimation
     mean        sd 
32.450662  1.717755 

$Anderson_Darling
[1] 0.5028619

$p_value
[1] 0.1917