一：问题描述
如何在stata程序中进行简单的数值模拟？
1）产生一个观测值个数为 500，服从[0,10]均匀分布的解释变量X;
2）产生一个观测值个数为 500，服从 Normal (0, 36) 的随机扰动项u;
3）根据数据产生机制y𝑖 = 1 + 2𝑥𝑖 + 𝑢𝑖 = 𝛽0 + 𝛽1𝑥𝑖 + 𝑢𝑖获得y𝑖；
4）做回归得到结果；
5）用同样的机制得到一个观测值个数为 500 的样本，重新估计方程，并与步骤
4）中的结果进行对比；
6）验证 OLS 的代数性质：残差的和为 0；残差与自变量乘积的和为 0。
二：解决方案
- 第一步：用OLS估计已知参数来看OLS是否足够接近真实的参数
- 在命令窗口输入命令语代码,会得到如图1所示的结果
图1
- 第二步：比较不同的样本容量（500，5000， standard error是否真的缩小了）
- 在命令窗口输入命令语代码，会得到如图2所示的结果
图2
- 第三步：验证代数性质，残差是否真的为0，以及残差乘以x是否为0
图3
[从结果中可以看出残差无限地接近于0]
三：do-file 命令语
图4
四：存在的问题
- 这只是简单的数值进行模拟，还有待考究......
网友评论