大家好：

新年快樂！

目前遇到的問題，請問該用什麼方法解決呢？謝謝！
原資料：應變數(y)；自變數有4個(x1、x2、x3、x4)
使用stata的向前逐步迴歸的指令，跑出來的模式可能如下：
y=x2+x4
接下來，把逐步回歸的結果所選出的兩個變數x2與x4，再去估計一條迴歸模式如：z=x2+x4，
第二家公司估計程序也是一樣，
使用stata的向前逐步迴歸的指令，跑出來模式如下：
y=x1+x3
再把x1及x3變數選出，再去估計z=x1+x3

請教大家：
上述步驟，共有100間公司，估計程序都相同，但是選擇變數會不同，
如何使用stata來執行這個步驟呢？謝謝大家！