在数据子集上使用boot

时间:2018-06-13 16:24:44

标签: r statistics-bootstrap

boot文档中,看起来应该可以直接通过命令传递有关用于引导的观察的信息,但我无法弄清楚如何访问索引。

举一个简单的例子,假设我只想使用带自动变速器的汽车作为我的自举。一种方法是直接对数据进行子集化。

# Bootstrap on a subset of mtcars
library(boot)
my_func <- function(data, i) {
  d <- data[i,]
  m <- lm(mpg ~ disp + log(hp), data = d)
  beta <- m$coefficients['log(hp)']
  y_mean <- mean(m$model[, 1])
  # Return estimated elasticity
  return(beta / y_mean)
}
my_boot <- boot(mtcars[mtcars$am == 1, ], my_func, R = 100)

我的问题是是否可以使用boot()中的选项实现相同的效果。在这个简单的例子中,它可能没有多大意义,但想象一下从列表中传递的一组更复杂的条件。

0 个答案:

没有答案