我如何在R中使用示例函数进行引导？替换行问题

Posted 2021-05-03

我正在尝试为数据框'ev_all_clean'创建引导程序

set.seed(1315)
boot_s <- as.data.frame(matrix(NA, ncol = 19, nrow = 1000))
for(k in 1:19){
  for(l in 1:1000){
  boot_s[l,k]<- sample(ev_all_clean[,k], size=1, replace=T)
  }
}

以上代码返回

replacement element 1 has 768 rows to replace 1 rows

非常感谢您的帮助！预先谢谢你。

答案

我刚刚根据您提供的规格创建了一个虚拟ev_all_clean，看看是否可以帮助您入门。

boot_s <- as.data.frame(matrix(NA, ncol = 19, nrow = 1000))
for(k in 1:19){
  for(l in 1:1000){
    boot_s[[l,k]]<- sample(ev_all_clean[,k], size=1, replace=TRUE)
  }
}

仅对脚本的功能进行一些解释：

ev_all_clean[,k]-遍历单个列和sample 1元素进行替换，并且将使用单个元素替换boot_s[[l,k]]处的值。您可以了解[和[[ here之间的区别。

由于要对单个元素进行采样，所以您可能希望替换单个元素，这就是错误消息试图告诉您的内容-我认为。

另一答案

一旦有了原始数据帧boot_s，就可以做第一个引导程序

# first prepare a dataframe of the same shape as your original
samplemat=boot_s 
# then perform bootstrap for each column separately
for(k in 1:19){
  samplemat[,k] = boot_s[sample(1:nrow(boot_s),nrow(boot_s),T),k]
}

样本准备行索引的向量，以保留在引导样本中。您将需要为每个所需的引导操作执行此操作（您将需要很多）。