use*_*743 5 for-loop regression stata output
我使用forval循环运行3,500次回归,每组一次.然后我需要总结一下结果.通常,当我使用循环来运行回归时,我使用估计存储函数,然后使用estout.下面是一个示例代码.但我相信此代码可以处理的限制为300.如果有人能让我知道如何自动完成3,500次回归的过程,我将非常感激.
示例代码:
forval j = 1/3500 {
regress y x if group == `j'
estimates store m`j', title(Model `j')
}
estout m* using "Results.csv", cells(b t) ///
legend label varlabels(_cons constant) ///
stats(r2 df_r N, fmt(3 0 1) label(R-sqr dfres N)) replace
Run Code Online (Sandbox Code Playgroud)
下面是一个示例,statsby我对由 rep78 变量定义的 5 个组中的每组运行 mpg 价格回归,并将结果存储在名为 my_regs 的 Stata 数据集中:
sysuse auto, clear
statsby _b _se, by(rep78) saving(my_regs): reg price mpg
use my_regs.dta
Run Code Online (Sandbox Code Playgroud)
如果您愿意,可以省略该saving()选项,然后您的数据集将在内存中被回归结果替换,因此您无需直接使用use.