我有一个指数分布的样本,比方说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))
但似乎它不起作用。