Как рассчитать стандартную ошибку начальной загрузки в 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 показывает среднее значение исходного набора данных. «Станд. Значение 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 .

Обратите внимание, что эта стандартная ошибка очень похожа на рассчитанную в предыдущем примере.