以下代码实现的是
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的值
非常感谢


雷达卡




京公网安备 11010802022788号







