求助求助，关于环境不确定性的测量，每个公司运用过去5年的数据，采用普通最小二乘法(OLS)运行如下模型，分别估计过去5年的非正常销售收入: “Sale= Φ0+Φ1Year + ε“中，Sale为销售收入；Year为年度变量，如果观测值是过去第4年的，则Year=1；如果观测值是过去第3年的，则Year=2；依次类推，如果观测值是当前年度的，则Year= 5。模型的残差即为非正常销售收入；计算公司过去5年非正常销售收入的标准差，再除以过去5年销售收入的平均值，从而得到当年未经行业调整的环境不确定性；同一年度同一行业内所有公司的未经行业调整的环境不确定性的中位数，即为行业环境不确定性，各公司未经行业调整的环境不确定性除以行业环境不确定性，即为公司经行业后的环境不确定性。一个code代表同一家公司，现在就是不知道每年每家公司的残差怎么搞，以下面的数据为例，大神们能不能帮忙写一个命令（stata软件用），不胜感激。

code	sale	year
1	287670026.6	5
1	120454422.5	4
1	80608820	3
1	72784567.16	2
1	97363301.61	1
2	481863433.3	5
2	84125970.69	4
2	52776994.29	3
2	52061702.69	2
2	99421227.75	1