부트 스트랩 통계 값

Question

부트 스트랩 통계 값 (원래, 바이어스 및 오류)을 별도의 목록으로 가져오고 싶지만 그 방법을 알 수는 없습니다. 내가 실제 값을 원하는 텍스트 대신, 이제

> library(boot) 
> set.seed(123) 
> mean.fun <- function(data, idx) { mean(data[idx]) } 
> data <- boot(data=rnorm(100), statistic=mean.fun, R=999) 
> names(data) 
[1] "t0"  "t"   "R"   "data"  
[5] "seed"  "statistic" "sim"  "call"  
[9] "stype"  "strata" "weights" 
> data 
ORDINARY NONPARAMETRIC BOOTSTRAP 


Call: 
boot(data = rnorm(100), statistic = mean.fun, R = 999) 


Bootstrap Statistics : 
    original bias std. error 
t1* 0.09040591 0.004751773 0.08823615 

:

다음은 예입니다. 분명히 data$t0은 "원본"이지만 바이어스 및 오류 값을 가져 오는 방법을 알지 못합니다. 함수 이름을 입력하는 것은 당신의 코드를 제공하기 때문에 또한

---

는, 내가 R에 boost를 입력하고, 소스 코드에서 조각을 복사하고, 내 로컬 R 설치에 검색을 시도했다. 그러나 아무것도 찾을 수 없습니다. 왜, R이 로컬 스토리지에서 해당 소스 코드를 가져 오지 않아야합니까?

Answer 1

std.error 및 바이어스는 부트 객체의 일부로 저장되지 않습니다. 그것은 즉시 계산된다 (참조 : https://stat.ethz.ch/pipermail/r-help/2011-July/284660.html를)

를 귀하의 경우, 시도 :

의미가

mean(data$t) - data$t0 
sd(data$t)