R의 부트 스트랩 적용 범위

평균값이 895.0385임을 알기 위해 부트 스트랩 간격의 범위를 추정했습니다. 내 벡터 b<-c(300,300,200,250,600...)을하고 난 부트 스트랩 및 출력 간격합니다R의 부트 스트랩 적용 범위

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

을하지만 커버리지 확률을 (그것이 진정한 평균을 포함 횟수)을 획득하기 위해이를 복제 할 수있는 방법?

출처

2012-12-05 Anonymous

'? replicate' 기능을 살펴보십시오. – Max

I 구축 할 수없는 것이 –

B <- 999 muH0 <-895.0385 N <- 52 DV <- sample2를 sdH0 <- 1318.027 getM <- 기능 (orgDV, IDX) { BSM <- 평균 (orgDV [ idx]) bsS2M <- ((N-1)/N) * var (orgDV [idx])/N c (bsM, bsS2M) } bOut <- boot (DV, statistic = getM, R = 012) 부츠를 복제 할 수 있지만 replicate boot.ci를 복제 할 수 있습니까? 예 : replicate (B, getM (DV, sample = –

나는 조금 전에 이걸 조금하려고 노력했다. 나는 boot 명령을 사용하지 않았다. 나는 예제 명령을 사용했으나 이것이 도움이 될 수있다. 나는 또한 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))) 
)

출처

2012-12-05 16:09:10 DanTheMan

커뮤니티 위키 상태가 필요하지 않습니다. – joran

그래, 나는 그것을 misclciked하고 나는 편집을 통해 그것을 제거하는 방법을 찾을 수 없다는 것을 알았습니다. = / – DanTheMan

R의 부트 스트랩 적용 범위

답변

관련 문제