Ну, я чувствую себя глупо.Ответ на мою проблему прост:
Я просто присоединяю дополнительную информацию к возвращаемому вектору в функции загрузки, а позже я просто подставляю матрицу bmat$t
.
Так что ответ можетвыглядеть так:
library(boot)
r=50;c=50
m1 <- (sample(1000,r*c,T))
nboot = 100
boot_fun <- function(m,b){
m <- m[b]
mn <- mean(m)
vr <- var(m)
return(c(mn,vr))
}
bmat <- boot(data=m1,statistic=boot_fun,R=nboot)
mn <- bmat$t[,1]
vr <- bmat$t[,2]
Надеюсь, это может кому-то помочь.