2012-12-03 17 views
5

Chcę uzyskać wartości statystyki ładowania początkowego (oryginał, odchylenie i błąd) na osobnej liście - ale nie mogę wymyślić, jak to zrobić.Wartości statystyki ładowania początkowego

Oto przykład:

> 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 

Teraz, zamiast tekstu chcę rzeczywiste wartości. Podobno data$t0 jest "oryginalny", ale nie widzę sposobu, aby uzyskać wartości błędu i błędu.


Ponadto, ponieważ wpisując nazwę funkcji daje swój kod, wpisałem boost w badania i skopiowany fragment z kodu źródłowego, i starał się szukać go na mojej instalacji lokalnej R. Ale nie mogłem nic znaleźć. Dlaczego R nie powinien pobierać tego kodu źródłowego z lokalnego magazynu?

Odpowiedz

1

Std.error i stronniczość nie są przechowywane jako część obiektu rozruchowego. Jest ona obliczana na bieżąco (patrz: https://stat.ethz.ch/pipermail/r-help/2011-July/284660.html)

W twoim przypadku, spróbuj:

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

Ma to sens, dziękuję. Czy możesz również powiedzieć, dlaczego nie udało mi się znaleźć kodu R dla 'boot' na moim komputerze? Wpisałem "boot" w R, podniosłem linię z wyjścia (np. 'Lapply (seq_len (RR), fn)') i szukałem go - nic nie znaleziono w mojej lokalnej instalacji R. – andreister

Powiązane problemy