Bootstrap和样本均值

时间:2014-10-14 00:16:43

标签: r statistics statistics-bootstrap

我有一个指数分布的样本,比方说x<-rexp(30,0.2)

我用替换重新取样5000次并计算样本均值:

resample<-replicate(5000,mean(sample(x,30,replace=TRUE)))

我执行以下直方图来查看T(X *) - T(X)的分布: hist(resample-mean(x),freq=FALSE)

我知道,因为我有一系列的iid指数,所以这个序列的总和有一个Gamma分布(用我考虑的指数rv的数量来缩放(即30))。

如何将此Gamma分布叠加到之前的直方图? 我正在尝试使用以下res.range<-seq(min(resample),max(resample),.001) lines(res.range, dgamma(res.range,shape=1,rate=0.2/30)) 但似乎它不起作用。

0 个答案:

没有答案