R 使用引导程序计算标准误差，is.data.frame(x) 中的错误：

Question

我使用下一个程序代码使用bootstrap for Libras Data movement set 估计标准误差：

mydata<-read.table('C:/Users/Desktop/libra.txt', sep=',', header=TRUE)
head(data)
custom.boot <- function(times, data=mydata) {
  boots <- rep(NA, times)
  for (i in 1:times) {
    boots[i] <- sd(sample(data, length(data), replace=TRUE))/sqrt(length(data))  
  }
  boots
}
# Mean standard error
mean(custom.boot(times=1000))

但我得到了下一个错误：

Error in is.data.frame(x) : 
  (list) object cannot be coerced to type 'double'

您能帮我找出问题并给出解决方法吗？提前致谢！

Answer 1

mydata 对象被读取为 395 x 91 data.frame，custom.boot 正在努力解决这个问题。当我添加更改为data = mydata[, 1] 时，该函数运行时没有出现错误。如果你想保持你的函数不变，我要么循环遍历每一列，要么将所有列堆叠成一个长列。

编辑：

如果您想遍历此 data.frame 的所有列，我会编写一个类似于您所拥有的循环，但按每一列切片：

for(i in 1:ncol(mydata)){
  custom.boot(times=1000, data=mydata[, i])
  print(mean(custom.boot(times=1000)))
}