2022年stata建模中的各种小问题

1、保存估计结果的命令：est store 名称使用保存结果的命令： ，estimates（名称）如果你把那个显示你用过的命令的窗口：窗口操作： windows review 如果你把那个显示变量的窗口：窗口操作： windows variables 时间序列填充和扩展时间区间: 命令： tsappend ,add(n) 增加 n 个观测值窗口操作：在上面找data edit 即像一个表格一样的图标点开即可编辑数据时间序列存在间断点问题，需要补齐处理：命令： tsfill 信息准则赤池信息准则（aic ） 判断判断模型的最大滞后阶数stata 命令：1先回归2estat ic 如何看 aic 统计

2、量：breusch-pagan,cook-weisberg 异方差检验stata 命令：1先回归2estat hettest varlist 或者在statistics postestimation(倒数第二个) reports and statistics(倒数第二个 )在里面选择（hettest）如何看统计量：white 异方差检验 : stata 命令：3先回归4estat imtest,white varlist 或者在statistics postestimation( 倒数第二个 ) reports and statistics(倒数第二个 )在里面选择（imtest）如何看统计量

3、：ramsey 回归设定误差检验：stata 命令：1先回归2estat ovtest 或者在 statistics postestimation(倒数第二个 ) reports and statistics( 倒数第二个 )在里面选择（ovtest）如何看统计量：多重共线性方差膨胀因子检验：精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 1 页，共 4 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 1 页，共 4 页 - - - - - - - - -1先回归2es

4、tat vif,uncentered 或者在statistics postestimation(倒数第二个) reports and statistics(倒数第二个 )在里面选择（vif ）如何看统计量：一般的当最大的方差膨胀因子超过10（相对保守的临界值定位30）后者平均方差膨胀因子超过 1 表示模型存在多重共线性的问题。uncentered 用于当模型没有常数项时的未中心化的方差膨胀因子。多重共线性的其他侦查方法：值高而显著的t 比率小：多重共线性的“ 经典 ” 征兆克里安经验法则：仅当来自一个辅助回归的大于得自y 对全部回归元中的总时，多重共线性才算是一个麻烦的问题。做拟合图 （前提是

5、先回归）stata 命令：1解释变量对成分残差图用于考察模型形式是否设定准确。cprplot 被解释变量acprplot 被解释变量2增加变量图用于考察数据是否存在异常值avplotd 被解释变量3拟合值对残差图的散点图用于考察残差是否满足经典的假设条件rvfplot 4解释变量对残差的散点图rvpplot 被解释变量stata 对于数据的储存与重现est 命令的用法：（1）储存回归结果：reg y x1 x2 x3（不限于 reg，也可储存 ivreg 、mvreg、reg3）est store a （2）重现回归结果：est replay a （3）对回归结果进行进一步分析est for

6、a:sum（对 a回归结果中的各个变量运行sum命令）在非时间序列的数据的情况下，异方差的修正用gls 具体的方法如下：1.quietly regress y x 做回归2.predict u，residual 取残差stdp 表示样本内预测的标准差stdr 表示样本外预测的标准差2r2r2r精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 2 页，共 4 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 2 页，共 4 页 - - - - - - - - -3.predi

7、ct yf ，xb (xb 表示拟合值 ) 将拟合值取出放到yf 里4.gen lnu2=ln(u2) 将残差做平方且取对数的处理5.gen yf2=yf2( 将 yf 这个拟合值同上面的残差做相同的处理) 6.quietly regress lnu2 yf yf2 对处理过的残差对拟合值以及 处理过的拟合值做回归7.predictnl u2f=exp(xb() 再将回归后的拟合值取出并作对数处理放到u2f 里8.gen sd=sqrt(u2f) 将 u2f 做平方处理然后，利用vwls 进行加权估计vwls y x , sd(sd) gls 也可以通过regress命令中的weight 选项

8、来实现。存在自相关的修正用广义差分自相关的修正用广义差分具体的方法如下：1一阶自相关的修正prais y x，rhotype(regress) prais y x，corc rhotype(regress) 2高阶自相关的修正以二阶 自回归为例quietly regress d.y x 对被解释变量取差分并且做回归predict u,resid 取出残差quietly reg u l( 2).u，noconstant 令 u 对其 二阶 滞后期做自回归（无截距）matrix mat=e(b) 生成矩阵mat 将回归的结果放到矩阵里（如果是更高阶可能有多个自相关系数）gen m=d.y-el(m

9、at,1,1)*l( 2)d.y 对 y 的一阶差分与y 的一阶差分的滞后期与权重的乘积做差分，这个权重就是mat矩阵里的第一行第一列的系数，刚好使我们刚刚回归出来的自相关系数（如果是高阶可能不只做一个差分，会更为复杂）gen n=d.x-el(mat,1,1)*l( 2)d.x 对 x 进行和 y 一样的处理方法reg m n 然后令 m 对 n 做回归关于参数约束的模型估计问题(p178 张晓峒 ) stata 命令：cnsreg 被解释变量解释变量条件 if in weight ，constraints(constraints) options 其中 constraints(constr

10、aints)表示线性约束例如：约束规模报酬不变的估计模型命令如下：1.constraint define 1 lnk+lnl=1 将等于 1 的线性约束 (不可以为不等号)条件定义为1 sreg lny lnk lnl ，constraints(1) 约束为 lnk=lnl=0 的估计模型：命令如下 : 1.constraint define 2 lnk lnl 等于 0 的条件可以不写出来sreg lny lnk lnl ， constraints(2) 约束条件为lnk=0.6, lnl=0.4的估计模型：命令如下：1.constraint define 3 lnk=0.6 predict

11、nl 表示模型估计后的非线性预测，比如指数预测xb 表示线性预测exp 表示指数预测pr 表示概率预测se 表示线性预测 （prediction）的标准差stdf 表示线性预测（ forecast）的标准差取对数用ln(var)函数取平方用sqrt(var)函数el(mat,i,j) 矩阵的第i行第 j 列精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 3 页，共 4 页 - - - - - - - - -精品学习资料 可选择p d f - - - - - - - - - - - - - - 第 3 页，共 4 页 - - - - - - - - -2.constraint define 3 lnl=0.4 sreg lny lnk lnl ，constraints