以下代码实现的是
Draw a 30 observations sample from a uniform distribution, compute the mean.

program onesample, rclass
drop _all
quietly set obs 30
gen x = runiform()
summarize x, meanonly
return scalar meanofonesample = r(mean)
end

现在为了用数据模拟central limit theroem，需要运行以上程序1000次 得出1000个mean然后画出histogram
我想到的是以下修改 蓝色部分为修改部分 红色为不知道如何改的部分

program onesample, rclassdrop _allquietly set obs 30在此处设定一个循环  如此一来就可以得出1000组 30observations的数据forv i = 1(1)1000{gen x`i' = runiform()}summarize, meanonlyreturn scalar meanofonesample = r(mean)end但是如何取出这1000个mean呢？如何将第Xi组数据的第j个数取出或更改？（因为每组数据有30个observations）也就是如何更改一组数据任意一observation的值

非常感谢