Stata模拟中心极限定理

关于中心极限定理的模拟

do文件如下：

set seed 52
set obs 500000  // 设置总样本数

* 如果变量 group_id 已经存在，删除旧变量
capture drop group_id
generate group_id = ceil(_n / 20)  // 将样本分成组，每组20个样本

* 一次性生成随机数据
generate x = runiformint(1, 6)

* 按组计算均值
collapse (mean) group_means = x, by(group_id)

* 绘制直方图叠加正态分布曲线
histogram group_means, start(1) bin(40) frequency normal ///
    xscale(range(1 6)) xlabel(1(1)6) ///
    title("中心极限定理模拟")

如图所示：