Tengo una pregunta con respecto a la interpretación de la tsboot llamada en R. he comprobado la documentación de Kendall y el inicio de paquete, pero no soy más inteligente que antes.
Cuando ejecuto un bootstrap utilizando, por ejemplo, en el ejemplo en el Kendall paquete, donde el estadístico de prueba es la tau de Kendall:
library(Kendall)
# Annual precipitation entire Great Lakes
# The Mann-Kendall trend test confirms the upward trend.
data(PrecipGL)
MannKendall(PrecipGL)
lo que confirma la tendencia al alza:
tau = 0.265, 2-sided pvalue =0.00029206
El ejemplo a continuación, sigue el uso de un bloque de bootstrap:
#
#Use block bootstrap
library(boot)
data(PrecipGL)
MKtau<-function(z) MannKendall(z)$tau
tsboot(PrecipGL, MKtau, R=500, l=5, sim="fixed")
Me aparece el siguiente resultado:
BLOCK BOOTSTRAP FOR TIME SERIES
Fixed Block Length of 5
Call:
tsboot(tseries = PrecipGL, statistic = MKtau, R = 500, l = 5,
sim = "fixed")
Bootstrap Statistics :
original bias std. error
t1* 0.2645801 -0.2670514 0.09270585
Si entiendo correctamente, el "t1* original" es el original de MKtau, el "sesgo" es la media de las MKtau de R=500 bootstrap de series de tiempo, y la "enfermedad de transmisión sexual. error" es la desviación estándar de la MKtaus de las 500 muestras.
Tengo problemas para entender lo que esto significa - básicamente, esto me dice que todos los 500 MKTaus son inferiores a la original, y que el original t1* está en el rango de 3 sd de el bootstrap MKtaus. Ergo es significativamente diferente?
O le voy a decir la MKtau para el conjunto de datos es de 0,26 más/menos el error estándar?
Lo siento por la larga pregunta, pero soy de las estadísticas de novato y estoy aprendiendo a través de auto-estudio, que carecen de alguien para platicar de esto probablemente es muy simple problema de ida y vuelta.