stata + 循环语句 死结

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

          
         keep if quantity!=. & price !=.
         qui sum hscode_i_link
          global N = r(max)
          cap dropvars q_hat e
          gen q_hat = .
          forvalues i = 1/$N{
              qui reg quantity price  if (hscode_i_link==`i')
              qui predict e if e(sample), res
              qui replace q_hat = e if e(sample)
              drop e
            }
这样总是出现“insufficient observations  r(2001);”
r(2001):
[P]     error . . . . . . . . . . . . . . . . . . . . . . . . Return code 2001
        insufficient observations;
        You have requested some statistical calculation, and while there
        are some observations, the number is not sufficient to carry out
        your request.
即使看懂这个字面意思，还是完全不知道哪里出错了，请各位指点迷津啊！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Stata tata 循环语句 observations Insufficient

写得不错                                                                        
                                       
                                                     
                                                     
                                             
                                                                 
                                                                                 

小迷彩魂 发表于 2014-10-2 21:00
写得不错                                                                        
                  ...

呀，真是个坏小孩！！

把样本数据上传来看看

ywh19860616 发表于 2014-10-3 08:17
把样本数据上传来看看

循环——样本数据.xls (11 KB)

2014-10-3 08:43:19 上传

循环——样本数据

回归的时候是按hscode_i_link分类回归的。因为加了  cap noisily
   forvalues i = 1/$N{
              qui capture noisily  reg quantity price if (hscode_i_link==`i')
              qui predict e if e(sample)
              qui replace q_hat = e if e(sample)
              drop e
            }

循环能进行下去了，但是结果让我非常迷惑：一、明明都是hscode_i_link==1,却会有断的呢？
                                                                      二、既然有中断，为什么q_hat又有取值呢？
                                                                      三、q_hat都有了，为什么系数又没有呢？
                                                                     

yuquanguihua 发表于 2014-10-3 08:49
回归的时候是按hscode_i_link分类回归的。因为加了  cap noisily
   forvalues i = 1/$N{
              qui ...

抱歉，刚看到
你为什么说hscode_i_link是都为1，明显hscode_i_link有等于2的。

ywh19860616 发表于 2014-10-6 09:16
抱歉，刚看到
你为什么说hscode_i_link是都为1，明显hscode_i_link有等于2的。

嗯，我说的是在hscode_i_link==1的情况下，有样本参与回归，但是有的却没有

虽然有关这个问题的细节还是不清楚，但是最大的问题应该是出在分类中有些类别的样本量太少了，把样本量少于5的类别剔除后，重新回归，现在已经得到想要的系数、拟合值和预测值了。谢谢各位的热心回答！