我使用下一个程序代码来估计使用Libras 数据移动集的引导程序的标准误差:
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'
您能帮我找出问题所在并给出如何解决的建议吗?提前致谢!