доверительный интервал пропорции начальной загрузки

Я хотел бы создать доверительные интервалы для пропорций, используя загрузочный пакет, если это возможно. У меня есть вектор, и я хотел бы установить порог, а затем рассчитать пропорции ниже указанного уровня. После этого я хотел бы использовать функцию начальной загрузки в загрузочном пакете для расчета доверительных интервалов для пропорций. Простой пример того, что у меня есть до сих пор:

library(boot)

vec <- abs(rnorm(1000)*10) #generate example vector

data_to_tb <- vec

tb <- function(data) {
  sum(data < 10, na.rm = FALSE)/length(data) #function for generating the proportion
}

tb(data_to_tb)

boot(data = data_to_tb, statistic = tb, R = 999)
quantile(boot.out$t, c(.025,.975))

Однако я получаю это сообщение об ошибке:

> boot(data = data_to_tb, statistic = tb, R = 999)
Error in statistic(data, original, ...) : unused argument (original)

Я не могу заставить его работать, хотя помощь приветствуется


person jonas    schedule 10.02.2015    source источник


Ответы (1)


Ваша проблема в вашей функции tb - ей нужны два аргумента. Из файла справки ?boot

статистика Функция, которая при применении к данным возвращает вектор, содержащий интересующую(ие) статистику(и). Когда sim = "параметрический", первым аргументом статистики должны быть данные. Для каждой реплики будет передан смоделированный набор данных, возвращенный ran.gen. Во всех остальных случаях статистика должна принимать как минимум два аргумента.

person csgillespie    schedule 10.02.2015
comment
Посмотрите файл справки. В частности, запись statistic и примеры внизу. - person csgillespie; 10.02.2015