Як обчислити стандартну помилку початкового завантаження в r

Бутстрапінг – це метод, який можна використовувати для оцінки стандартної помилки середнього.

Основний процес обчислення початкової стандартної помилки такий:

• Візьміть k повторних зразків із заміною з заданого набору даних.
• Для кожного зразка обчисліть стандартну похибку: s/√ n
• Це призводить до k різних оцінок стандартної помилки. Щоб знайти початкову стандартну помилку, візьміть середнє k стандартних помилок.

У наведених нижче прикладах пояснюється два різні методи, які можна використовувати для обчислення початкової стандартної помилки в R.

Спосіб 1. Використовуйте стартовий пакет

Одним із способів обчислення стандартної помилки завантаження в R є використання функції boot() із бібліотеки завантаження .

Наступний код показує, як обчислити стандартну помилку початкового завантаження для заданого набору даних у R:

 #make this example reproducible
set. seeds (10)

#load boot library
library (boot)

#define dataset
x <- c(12, 14, 14, 15, 18, 21, 25, 29, 32, 35)

#define function to calculate mean
meanFunc <- function (x,i){mean(x[i])}

#calculate standard error using 100 bootstrapped samples
boot(x, meanFunc, 100)

Bootstrap Statistics:
    original bias std. error
t1* 21.5 0.254 2.379263

«Початкове» значення 21,5 показує середнє значення вихідного набору даних. “std. Значення 2,379263 вказує на початкову стандартну помилку середнього значення.

Зауважте, що ми використали 100 початкових зразків для оцінки стандартної помилки середнього в цьому прикладі, але ми могли використати 1000 або 10 000 або будь-яку кількість початкових зразків, яку нам потрібно.

Спосіб 2: Напишіть власну формулу

Інший спосіб обчислити початкову стандартну помилку — це написати нашу власну функцію.

Наступний код показує, як це зробити:

 #make this example reproducible
set. seeds (10)

#load boot library
library (boot)

#define dataset
x <- c(12, 14, 14, 15, 18, 21, 25, 29, 32, 35)

mean(replicate(100, sd( sample (x, replace= T ))/sqrt( length (x))))

[1] 2.497414

Стандартна помилка початкового завантаження виявляється 2,497414 .

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