stata 循环语句请教

各位大神，现有500的Y，分别命名为Y1、Y2、……Y500。 6个自变量，分别为A、B、C、D、E、F，现在需要分别对500个Y做回归分析，并显示出各个解释变量的回归系数，如何写程序啊？我用的forvalues 语句，不知哪里出错了，老是跑不出来结果，请教大神帮忙，非常感谢。

gen betaa=.       //生成一个变量存储回归系数，这个存储自变量A的回归系数，下面类推
gen betab=.
gen betac=.
gen betad=.
gen betae=.
local i=1
foreach y of varlist Y1-Y500 {
reg `y' A B C D E
replace betaa=_b[A] in `i'
replace betab=_b[B] in `i'
replace betac=_b[C] in `i'
replace betad=_b[D] in `i'
replace betae=_b[E] in `i'
local i=`i'+1
}

要用foreach

1. set matsize 700                                // 设置矩阵大小
   
2. foreach y of varlist Y1-Y500 {
   
3. qui: reg `y' A B C D E
   
4. matrix beta=nullmat(beta) \ e(b)       //将回归结果保存到矩阵beta中
   
5. }
   
6. 
   
7. clear                                               //删去现有的数据，如有需要先保存原数据。
   
8. svmat beta, clear                             //将矩阵beta转化为内存中的数据

复制代码

1. gen betaa=.       //生成一个变量存储回归系数，这个存储自变量A的回归系数，下面类推
   
2. gen betab=.
   
3. gen betac=.
   
4. gen betad=.
   
5. gen betae=.
   
6. local i=1
   
7. foreach y of varlist Y1-Y500 {
   
8. reg `y' A B C D E
   
9. replace betaa=_b[A] in `i'
   
10. replace betab=_b[B] in `i'
    
11. replace betac=_b[C] in `i'
    
12. replace betad=_b[D] in `i'
    
13. replace betae=_b[E] in `i'
    
14. local i=`i'+1
    
15. }

复制代码

506232839 发表于 2017-8-2 00:01

您好。很感谢您的精彩回复，但是我还是遇到几个问题，我采用你的程序进行回归时，总是出现几个错误。我的代码是：
gen betamkt=.   
gen betasmb=.  
Gen betahml=.
local i=1
foreach y of varlist y1-y500 {
reg `y' mkt smb hml
replace betamkt=_b[mkt] in `i'
replace betasmb=_b[smb] in `i'
replace betahml=_b[hml] in `i'
local i=`i'+1
}
这种情况下，就会出现
beta ambiguous abbreviation

如果采用以下程序
gen betamkt=.   
local i=1
foreach y of varlist y1-y500 {
reg `y' mkt
replace betamkt=_b[mkt] in `i'
local i=`i'+1
}
就是出现　　
Obs. nos. out of range
r(198);

能否帮忙指导下呢

应该可以啊，你qq多少？

506232839 发表于 2017-8-2 15:28
应该可以啊，你qq多少？

565986231

506232839 发表于 2017-8-1 22:36
gen betaa=.       //生成一个变量存储回归系数，这个存储自变量A的回归系数，下面类推
gen betab=.
g ...

我根据您的程序改版的
foreach r of varlist r1-r145{
arch `r' L(1).`r', arch(1) garch(1) dist(t) nolog
predict e_`r',res
predict h_`r',variance
}
总是提示我
flat log likelihood encountered, cannot find uphill direction
这是什么情况

506232839 发表于 2017-8-1 22:36
gen betaa=.       //生成一个变量存储回归系数，这个存储自变量A的回归系数，下面类推
gen betab=.
g ...

这样是能够生成的，只是不知道有没有办法讲生成的系数和R平方全部导出成一个序列呢

请据以更改 (请 ssc install asreg)：

1. clear
   
2. set obs 50
   
3. set seed 98034
   
4. 
   
5. gen y1 = runiform()
   
6. gen y2 = runiform()
   
7. gen y3 = runiform()
   
8. gen A = runiform()
   
9. gen B = runiform()
   
10. gen C = runiform()
    
11. 
    
12. gen t = _n
    
13. reshape long y, i(t) j(id)
    
14. sort id t
    
15. bys id: asreg y A B C, fit

复制代码