数据如下所示

目前需要按照申慧慧老师的环境不确定性的表征方法，数据处理要求是：

其中，Sales 为企业的销售收入，Year为年份，从计算当年开始倒退四年，分别赋值为 5-1。例如：计算2017 年的环境不确定性，要以 2007-2011 年这五年的销售收入作为原始数据，将 2011 年赋值为 5，2010年赋值为 4，2009 年赋值为3，2008年赋值为 2，2007 年赋值为 1。同理，计算 2012 年的环境不确定性，要以 2008-2012年的销售收入为基础，以此类推。由该模型计算得出的残差项，即为企业的非正常销售收入。将该残差的标准差除以过去五年销售收入的平均值，即为初始的环境不确定性。由于不同行业受外界环境因素的影响程度不同，本文再将初始的环境不确定性除以同一行业的中位数，得到最终求出的经过行业调整的环境不确定性 EU。为了进一步验证本文提出的假设，增添了对环境不确定性新的衡量方式，根据行业中位数进行调整，将环境不确定性 EU 高于行业中位数的赋值为 1，否则为 0，最终得到环境不确定性高低 EUhigh 哑变量。
我参考之前的关于环境不确定性的发帖，某个回帖给的命令如下：

gen msale=.
gen dev=.
forvalues i = 2006(1)2011{
          forvalues j = 1(1)2{
                  gen YEAR= year-`i'
                  qui reg sale YEAR if YEAR <6 &YEAR>0 & stkcd ==`j'
                  qui predict p if YEAR<6 & YEAR>0& stkcd ==`j',r
                  
                  qui sum p if stkcd==`j',de
                  qui dis "i=`i' j=`j'"
                  qui replace dev=r(sd) if stkcd==`j' & YEAR==5
                  qui sum sale if stkcd==`j' & YEAR<6 & YEAR>0 ,de
          qui replace  msale=r(mean) if stkcd==`j' & YEAR==5
                  
                  drop YEAR p  
        
                  }
}
dis "loop is over"
gen eu=dev/msale

但是我在执行倒数第三步的时候，stata显示“no observations r(2000);”
目前，我想请教各位大佬，上述命令是否正确，如果正确，为何没办法完全执行，望不吝赐教（我的这个问题貌似有点复杂，但我有在尽量简单解释，若我表达有问题，望告知，感恩）