重复重采样功能1000次?使用lapply?

时间:2013-12-10 09:17:15

标签: r lm statistics-bootstrap

请我出去!我感谢任何帮助!谢谢!

我无法重复重复采样1000次。 我尝试使用replicate()来做到这一点,但它不起作用。有没有其他方法可以做到这一点?任何人都可以告诉我,如果这可能通过使用lapply完成?  以下是我的代码:

#sampling 1000 betas0 & 1 (coefficients) from the data
get.beta=function(data,indices){ 
  data=data[indices,] #let boot to select sample
  lm.out=lm(y ~ x,data=data)
  return(lm.out$coefficients)
}
n=nrow(data)
get.beta(data,1:n)

bootcoe=boot(data,get.beta,R=1000) #generate 1000 random samples
head(bootcoe$t) #look at the betas

从上面的代码中我可以得到1000个betas0& 1通过随机抽样数据。我想做1000次以获得不同的测试版。除了replicate()之外我该怎么办?

1 个答案:

答案 0 :(得分:22)

这更像是一个扩展评论,我证明replicate应该有效。这是CLT的一个例子。只需更换花括号之间的线条。

x <- replicate(1000, {
  mm <- runif(10)
  mean(mm)
  })
hist(x)

enter image description here

相关问题