बूटस्ट्रैप आंकड़ों के मान

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 "मूल" है लेकिन मुझे नहीं लगता कि पूर्वाग्रह और त्रुटि के मान कैसे प्राप्त करें।

---

इसके अतिरिक्त, क्योंकि एक समारोह का नाम लिखकर आप अपने कोड देता है, मैं आर में boost टाइप किया और स्रोत कोड से एक टुकड़ा की नकल की है, और अपने स्थानीय आर स्थापना पर यह खोज करने के लिए कोशिश की। लेकिन कुछ भी नहीं मिला। क्यों, स्थानीय स्रोत से उस स्रोत कोड को आर नहीं लेना चाहिए?

Answer 1

std.error और पूर्वाग्रह बूट ऑब्जेक्ट के हिस्से के रूप में संग्रहीत नहीं हैं। यह मक्खी पर गणना की जाती है (देखें: https://stat.ethz.ch/pipermail/r-help/2011-July/284660.html)

आपके मामले में, कोशिश:

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