Rでbootstrap

Rでbootstrap推定をする場合の基本的な使い方:
require(boot)
boot(data, statistic=function(x,i)mean(x[i]), R=100, stype="i")

この例は平均値を計算するという単純なもので、bootstrapを使う必要はないが、functionの部分を書き変えれば基本的にはどのような推定量でもよい。なお、R=100がbootstrapの繰り返し回数、iはデフォルトでインデックスとして扱われるので、stype="i"はなくてもよい。