R语言中boot函数的小小小疑问！！！

求教 关于R语言中的boot函数：
boot(data,statistic,R,...)

statistic= function(dat,ind)
这里的dat是每次从data中抽取一个bootstrap样本吗？
如果R设为100，那么，这100次的dat和ind是都不同的吗？

下面代码是想从bootstrap样本中进行Monte Carlo模拟，请问这样子对吗？？？
set.seed(1234)
theat.acf_95 <- function(dat,ind){
  acfood1 <- dat[ind,1]*0.3
  acfood4 <- dat[ind,4]*0.5
  su <- (acfood1+acfood4)/21.32

  acf1 <- sample(dat[,1],100000,replace = T)*0.3
  acf4 <- sample(dat[,4],100000,replace = T)*0.5
  s <- (acf1+acf4)/21.32

  S <- c(su,s)
  a <- quantile(S,0.95)
  a
}
boot.obj1_95 <- boot(data=ACfood,statistic=theat.acf_95,R=2000)

lion果 发表于 2020-3-21 18:27
请问 “d

这样不对，重抽样是通过d <- dat[ind,]实现的，重抽样样本你该用d，不是用dat

各位大佬看看这里！！！

每次从data中抽取一个bootstrap样本，每次的dat和ind都不一样。

Whig 发表于 2020-3-21 17:27
每次从data中抽取一个bootstrap样本，每次的dat和ind都不一样。

请问 “d <- dat[ind,]”放着但没用可以吗？比如下面这样?
theat.acf_95 <- function(dat,ind){
  d <- dat[ind,]
  acf1 <- sample(dat[,1],1,replace = T)*0.3
  acf4 <- sample(dat[,4],1,replace = T)*0.5
  acf5 <- sample(dat[,5],1,replace = T)*0.3
  acf6 <- sample(dat[,6],1,replace = T)*2.2
  acf8 <- sample(dat[,8],1,replace = T)*21.0
  acf9 <- sample(dat[,9],1,replace = T)*7.9
  acf11 <- sample(dat[,11],1,replace = T)*5.4
  s <- (acf1+acf4+acf5+acf6+acf8+acf9+acf11)/21.32
  a <- quantile(s,0.95)
  a
}

Whig 发表于 2020-3-21 17:27
每次从data中抽取一个bootstrap样本，每次的dat和ind都不一样。

请问每次新抽的dat的样本容量，和data的一样吗？

Whig 发表于 2020-3-22 11:05
这样不对，重抽样是通过d

意思是 d 就是重抽样出来的一个bootstrap样本吗？
请问下面这样对吗？
theat.acf_95 <- function(dat,ind){
  d <- dat[ind,]
  acf1 <- sample(d[,1],1,replace = T)*0.3
  acf4 <- sample(d[,4],1,replace = T)*0.5
  acf5 <- sample(d[,5],1,replace = T)*0.3
  acf6 <- sample(d[,6],1,replace = T)*2.2
  acf8 <- sample(d[,8],1,replace = T)*21.0
  acf9 <- sample(d[,9],1,replace = T)*7.9
  acf11 <- sample(d[,11],1,replace = T)*5.4
  s <- (acf1+acf4+acf5+acf6+acf8+acf9+acf11)/21.32
  a <- quantile(s,0.95)
  a
}

lion果 发表于 2020-3-23 13:22
意思是 d 就是重抽样出来的一个bootstrap样本吗？
请问下面这样对吗？
theat.acf_95

是的，重抽样体现在ind上了

已解决！

Whig 发表于 2020-3-23 18:14
是的，重抽样体现在ind上了

万分感谢！！！