如何一次性获得分年度分行业回归的残差？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

这两天处理数据，希望能够得到分年度、分行业的回归残差。也就是与一次性求得分年度、分行业的可操纵性应计利润一样。
不懂编程，结果只能一次一次进行回归，并保留残差，实在太繁琐了！请哪位帮忙一下，谢谢！

有点担心没有说清楚，举例如下：
比如y＝a＋bx＋年度＋行业
假设年度有：2003、2004、2005、2006、2007、2008、2009
假设行业有：A、B、C、D、E、F、G.。。。。。
要求：分别对：2003年的A行业回归，并将残差保留到一个series中
                          2004年的A。。。。。
                          2005年的A。。。。。
                          。。。。。。。。。
                          2003年的B行业。。。。
                          2004年的B行业。。。。
                          。。。。。。。
并要求，将所有回归的残差保留在一个同一个series中。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏14 回帖

关键词：分行业 一次性 Series Serie seri 残差 年度 分行业

*-假设行业变量为 sic        

egen t = group(year)
qui sum t
local Nt = r(max)
egen s = group(sic)
qui sum s
local Ns = r(max)

gen res = .

forvalues t = 1/`Nt'{
   forvalues s = 1/`Ns'{
      cap qui reg tl size if (t==`t' & s==`s')
          cap qui predict e if e(sample)
          cap qui replace res = e if e(sample)
   }
}

非常感谢版主帮忙呀！

arlionn 发表于 2010-7-8 09:22
*-假设行业变量为 sic        

egen t = group(year)

如果删除了一些样本值，导致年份不连续，这些操作也适用吗？为什么我进行了这些操作，得到的残差全是缺漏值呀？

请arlionn老师帮忙，我按照这个循环执行，只有t=1,s=1这组的拟合值，stata提示e已经有定义，就停止循环了，找不出原因？

arlionn老师，问题解决了，老师的循环里predict后面应该加res,replace 后面应该加一条drop e

lsyjy 发表于 2013-5-9 10:26
请arlionn老师帮忙，我按照这个循环执行，只有t=1,s=1这组的拟合值，stata提示e已经有定义，就停止循环了， ...

老师，问题解决了。呵呵

lsyjy 发表于 2013-5-9 11:09
老师，问题解决了。呵呵

你好，能不能把正确的代码贴出来呢？谢谢了

在老师指令的后面加一条drop
set more off
egen g=group(year industry)
sum g
lacal Ng=r(max)
gen aq=.
gen wc_hat=.  
forvalues i=1/'Ng` {
   reg wc_delta cfo_l cfo cfo_f  sale_delta PPE if g==`i'
   predict aq1 if g==`i',res
   predict WC_pre1 if g==`i'
   replace aq=aq1 if g==`i'
   replace wc_hat=WC_pre1  if g==`i'
   drop aq1 WC_pre1
}

想问一下，按年度分行业进行回归并预测因变量的时候为什么会出现invalid syntax是怎么回事