Добро пожаловать на курс!

Статистический вывод для числовых данных в R

Mine Cetinkaya-Rundel

Associate Professor of the Practice, Duke University

Аренда в Манхэттене

В один из дней на Craigslist Manhattan были случайно выбраны двадцать однокомнатных квартир, размещённых в разделе «от владельца» (в отличие от объявлений агентств).

Что лучше отражает типичную арендную плату в Манхэттене — среднее или медиана?

Аренда однокомнатных квартир в Манхэттене

Методы бутстрэпа

Предполагаем, что данные репрезентативны
Метод бутстрэпа

Наблюдаемая выборка

медиана выборки = $2 350

Бутстрэп-выборка и бутстрэп-генеральная совокупность

Бутстрэп-генеральная совокупность

Бутстрэп-выборка и бутстрэп-генеральная совокупность

Схема бутстрэпа

Возьмите бутстрэп-выборку — случайную выборку с возвращением из исходной выборки того же размера.
Вычислите бутстрэп-статистику — например, среднее, медиану, долю и т. д. — по бутстрэп-выборкам.
Повторите шаги (1) и (2) много раз, чтобы построить бутстрэп-распределение статистик.

Схема бутстрэпа в R

library(infer)

___ %>%                               # start with data frame
  specify(response = ___) %>%         # specify the variable of interest

Схема бутстрэпа в R

library(infer)

___ %>%                               # start with data frame
  specify(response = ___) %>%         # specify the variable of interest
  generate(reps = ___, type = "bootstrap") %>%  # generate bootstrap samples

Схема бутстрэпа в R

library(infer)

___ %>%                               # start with data frame
  specify(response = ___) %>%         # specify the variable of interest
  generate(reps = ___, type = "bootstrap") %>%  # generate bootstrap samples
  calculate(stat = "___")             # calculate bootstrap statistic

Построение бутстрэп-интервала

library(infer)

___ %>%                               # start with data frame
  specify(response = ___) %>%         # specify the variable of interest
  generate(reps = ___, type = "bootstrap") %>%  # generate bootstrap samples
  calculate(stat = "___")             # calculate bootstrap statistic

Бутстрэп-распределение без осей

Построение бутстрэп-интервала

library(infer)

___ %>%                               # start with data frame
  specify(response = ___) %>%         # specify the variable of interest
  generate(reps = ___, type = "bootstrap") %>%  # generate bootstrap samples
  calculate(stat = "___")             # calculate bootstrap statistic

Бутстрэп-распределение без осей со скобками

Давайте потренируемся!

Статистический вывод для числовых данных в R