如何为Stata回归程序添加参数以支持y1、y2双因变量回归

核心改动思路 #

给reg_and_save程序新增因变量参数，让它能接收y1或y2作为输入，同时调整结果文件名避免不同因变量的结果互相覆盖。

修改后的完整代码 #

global xlist1 "x1"
global xlist2 "x2"
global xlist3 "x1 x2"
global xlist4 "x1 x3"
global xlist5 "x1 x4"
global xlist6 "x1 x5"

capture program drop reg_and_save
program define reg_and_save
    // 新增第一个参数depvar，接收因变量名称
    args depvar j
    // 使用传入的因变量执行回归
    reg `depvar' ${xlist`j'}
    // 文件名加入因变量标识，防止结果覆盖
    est save results_`depvar'_`j'.ster, replace
end

// 循环两个因变量，再分别执行6组自变量组合的回归
foreach depvar in y1 y2 {
    forvalues j=1(1)6 {
        reg_and_save `depvar' `j'
    }
}

关键改动说明 #

• 程序参数扩展：把原args j改成args depvar j，让程序能接收因变量和模型编号两个参数，灵活支持不同因变量。
• 回归命令适配：将固定的y1替换为传入的参数depvar，实现因变量动态切换。
• 结果文件区分：保存结果时用results_depvar'_j'.ster命名，比如results_y1_1.ster、results_y2_6.ster，避免不同因变量的同编号模型结果互相覆盖。
• 循环逻辑优化：用foreach遍历y1和y2，再嵌套forvalues循环模型编号，一次性完成所有回归任务。

内容的提问来源于stack exchange，提问作者Z.D