Bootstrap-покрытие в R

Я бы оценил покрытие интервала начальной загрузки для среднего значения, зная, что истинное среднее значение равно 895.0385. У меня есть вектор b<-c(300,300,200,250,600...), и я делаю начальную загрузку и интервал вывода:

mean.fun <- function(dat, idx) mean(dat[idx], na.rm = TRUE)
boot.out <- boot(b, mean.fun, R=999)
boot.ci(boot.out)

Но как я могу воспроизвести это, чтобы получить вероятность покрытия (сколько раз оно содержало истинное среднее значение)?


r statistics-bootstrap
person Community    schedule 05.12.2012    source источник
comment
Посмотрите на функцию?replicate.   -  person assumednormal    schedule 05.12.2012
comment
B ‹- 999 мкH0 ‹-895,0385 N ‹- 52 DV‹- образец2 sdH0 ‹- 1318,027 getM ‹- function(orgDV, idx) { bsM ‹- среднее (orgDV[idx]) bsS2M ‹- (((N-1) / N) * var(orgDV[idx])) / N c(bsM, bsS2M) } bOut ‹- boot(DV, statistic=getM, R=B) boots ‹- t(replicate(B, getM(DV, sample (seq(along=DV), replace=TRUE)))) Я могу реплицировать загрузку, но для репликации boot.ci?   -  person    schedule 05.12.2012

Ответы (1)


arrow_upward
1
arrow_downward

Я пытался сделать что-то вроде этого немного назад. Я не использовал команду загрузки, я использовал примерную команду, но это может помочь. Я также могу ошибаться на 100%, я еще не очень хорошо разбираюсь в R.

mean.fun <- function(dat, idx) mean(dat[idx], na.rm = TRUE)
bootoutput <- data.frame(
  bootoutput = replicate(10000, boot.ci(boot(b, mean.fun, R=999)))
)
person Community    schedule 05.12.2012
comment
Да, я ошибся, и к тому времени, когда я увидел, я не мог найти способ удалить это с помощью редактирования. знак равно - person DanTheMan; 05.12.2012