Если вы хотите загрузить средства по группам, вы можете сделать что-то, следуя приведенным ниже строкам.
Это ваш код, средство по группе Species
, полный набор данных.
res1 <- with(iris, tapply(Sepal.Length, Species, mean))
Теперь функция начальной загрузки и ее применение.
library(boot)
bootMeans <- function(data, indices, val, catg, na.rm = FALSE){
val <- deparse(substitute(val))
catg <- deparse(substitute(catg))
d <- data[indices, ]
tapply(d[[val]], d[[catg]], mean, na.rm = na.rm)
}
set.seed(1234)
b <- boot(iris, bootMeans, R = 1000,
val = Sepal.Length, catg = Species, na.rm = TRUE)
identical(b$t0, res1)
#[1] TRUE
colMeans(b$t)
#[1] 5.006788 5.940504 6.588357