Розуміння виводу завантажувальної програми, виконаної в R (tsboot, MannKendall)

У мене виникає запитання щодо інтерпретації виклику tsboot в Р. Я перевірив документацію як Kendall, так і завантажувального пакету, але я не розумніший, ніж раніше.

Коли я запускаю завантажувальний інструмент, використовуючи, наприклад, приклад в пакеті Kendall, де тестовою статистикою є тау Kendall:

library(Kendall)
# Annual precipitation entire Great Lakes
# The Mann-Kendall trend test confirms the upward trend.
data(PrecipGL)
MannKendall(PrecipGL)

що підтверджує тенденцію до зростання:

tau = 0.265, 2-sided pvalue =0.00029206

Потім приклад продовжує використовувати блок завантаження:

#
#Use block bootstrap 
library(boot)
data(PrecipGL)
MKtau<-function(z) MannKendall(z)$tau
tsboot(PrecipGL, MKtau, R=500, l=5, sim="fixed")

Я отримую такий результат:

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

Якщо я правильно розумію, "t1 * original" є оригінальним MKtau, "зміщення" - це середнє значення MKtau від завантаженого часового ряду R = 500, а "std. Помилка" - це стандартне відхилення MKtaus від 500 зразків.

У мене виникають проблеми з розумінням того, що це означає - це, в основному, говорить мені, що всі 500 MKTaus нижче, ніж оригінал, і що вихідний t1 * знаходиться в діапазоні 3 sd завантаженого MKtaus. Ерго це суттєво відрізняється?

Або я б сказав, що MKtau для набору даних становить 0,26 плюс / мінус стандартна помилка?

Вибачте за тривале запитання, але я є початківцем статистикою і навчаюсь через самостійне вивчення, мені не вистачає, щоб хтось підстрибував цю, ймовірно, дуже просту проблему.

r bootstrap

— Марія
джерело

У виході biasє просто різниця між середнім значенням 500 збережених зразків завантажувальної програми та початковою оцінкою. Це std. error- стандартне відхилення 500 зразків завантажувальної камери і є оцінкою стандартної помилки. Результат повідомляє вам, що ваша початкова оцінка вище, ніж середня оцінка 500 завантажених завантажень (тому не всі завантажені MKtaus нижче). Завантажувальний інструмент часто використовується для обчислення стандартних інтервалів помилок / довіри, не роблячи припущень щодо розподілу. Використовуйте boot.ciфункцію для обчислення довірчих інтервалів.

— COOLSerdash

@COOLSerdash, дякую за це! Отже, якщо моя первісна статистика на 3 сд вище середнього показника завантаженої статистики, чи можу я потім зробити що-небудь безпосередньо (скажімо: статистика є значущою 0,99)? Я також використав boot.ciдля обчислення довірчих інтервалів, і знову ж таки, початково обчислена статистика лежить поза цими інтервалами.

— Марія

Ні, ви не порівнюєте завантажену статистику з початковою статистикою з тестом гіпотези. Я просто використовую / повідомляю про завантажений стандартний інтервал помилок та довіри у вашому випадку.

— COOLSerdash

Зіткнувшись із тим самим питанням і дослідивши його за допомогою керованого набору даних - модель y = ax + b з N (0, sig) помилками, я вважаю, що пакет Kendall може не працювати як рекламований. У моєму випадку x1:100 , а y = x , з sig = 100 (дисперсія терміна помилки).

Регресія виглядає добре, як і тау Кендалла. Тут немає ніякої автокореляції, окрім індукованої лінійною моделлю. Запуск тестування Кендалл за рекламою з довжиною блоків 1, 3, 5 та 10 дає дуже великі значення зміщення, і boot.ciне повідомляє про тенденцію.

Згодом я зашифрував завантажувальний файл даних із цією довжиною блоку, і з моїм контрольним рядом я отримую розумні результати щодо середнього значення зразків завантажувальної програми та їх розповсюдження. Отже, можливо, що щось зіпсувалося з пакетом Kendall щодо блочного завантаження.

— user46157
джерело