# 2005 年 08 月 13 日作成
lambda<-3
no.sample<-10
sample<-rexp(no.sample,lambda)
print(sample)
print(1/mean(sample))
no.rep<-1000
boot<-rep(no.rep,0)
bsample<-rep(no.sample,0)
for (i in 1:no.rep){
for (j in 1:no.sample){
bsample[j]<-sample[trunc(runif(1,0,no.sample)+1)]
}
boot[i]<-1/mean(bsample)
}
plot(density(boot))