web-dev-qa-db-ja.com

リサンプリング機能を1000回繰り返しますか？ lapplyを使用していますか？

私を出してください！助けてくれてありがとう！ありがとう！

1000回リサンプリングを繰り返すのが苦手です。それを行うためにreplicate（）を使用してみましたが、機能しません。それを行う他の方法はありますか？これがラッププライを使用して行われた可能性がある場合、誰かに見せてもらえますか？以下は私のコードです：

#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ベータ0および1を取得できます。そして、私はそれを1000回行って異なるベータ版を入手したいと思っています。 replicate（）以外にどうすればよいですか？

rlmstatistics-bootstrap
10
user2978129

これは、replicateが機能することを示す拡張コメントの詳細です。これはCLTの例です。中かっこの間にある行を置き換えるだけです。

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

enter image description here

24
Roman Luštrik