用保留残差命令后，全部残差显示缺失？

20324

收藏 2011-05-02

亲爱的版主：请问我需要提取回归的残差项，用以下命令后，结果是所有值全部是缺失值，是何原因，该如何处理？谢谢！
statsby _b _se,by( code year): reg growth population dum1 dum2 dum3 dum4

扫码加我拉你入群

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

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

全部回复

sungmoo

2011-5-2 23:27:55

sharlagao 发表于 2011-5-2 21:33 请问我需要提取回归的残差项，用以下命令后，结果是所有值全部是缺失值，是何原因，该如何处理？谢谢！
statsby _b _se,by( code year): reg growth population dum1 dum2 dum3 dum4

*你首先看一下它运行的结果是否可能得到残差。

egen g=group(code year)
su g
loc n=r(max)
g e=.
forv i=1/`n'{
reg growth population dum1 dum2 dum3 dum4 if g==`i'
predict r if e(sample),r
replace e=r if e(sample)
drop r
}
drop g

扫码加我拉你入群

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

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

sharlagao

2011-5-3 00:03:56

当用到这个：
drop r
}
drop g
系统提示：invalid numlist has too many elements
e 全部是缺失。
我是否可以直接用predict residuals,r ? 但如何在数据中显示出来？谢谢！

扫码加我拉你入群

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

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

sungmoo

2011-5-3 00:11:13

sharlagao 发表于 2011-5-3 00:03 系统提示：invalid numlist has too many elements

修正了一下。code、population对应的unique观测值组太多。

扫码加我拉你入群

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

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

sharlagao

2011-5-3 00:14:11

我直接用predict res,r 出来了。这二个命令难道结果不一样吗，为何用statsby _b _se,by( code year): reg growth population dum1 dum2 dum3 dum4,就不行？

扫码加我拉你入群

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

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

sharlagao

2011-5-3 01:50:34

问题解决了，是我的数据出问题。但我还有一问：用statsby _b _se命令
得到的结果：_b_popultaion _b_cons _se_growth _se_cons, 代表的分别是什么？哪个是模型的残差？其他的变量是什么?谢谢！

扫码加我拉你入群

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

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

点击查看更多内容…

sharlagao

2011-5-3 01:57:10

知道了这样的结果，只能得到每年每个行业的回归系数和方差。谢谢版主！

扫码加我拉你入群

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

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

sharlagao

2011-5-3 02:04:20

sungmoo: 我试用您的命令，有用了版上讨论的命令predict da, res ，结果不一样。请问：da表示什么？Nda表示什么？谢谢！

扫码加我拉你入群

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

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

sharlagao

2011-5-3 22:52:54

连老师能帮我解答上面的问题吗？谢谢！

扫码加我拉你入群

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

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

sungmoo

2011-5-4 00:06:03

sharlagao 发表于 2011-5-3 02:04 sungmoo: 我试用您的命令，有用了版上讨论的命令predict da, res ，结果不一样。请问：da表示什么？Nda表示什么？

predict a ,r

a可以随便命名，生成残差变量a。

扫码加我拉你入群

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

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

sungmoo

2011-5-4 00:15:28

sharlagao 发表于 2011-5-3 00:14 我直接用predict res,r 出来了

首先你要在理论上清楚，“残差”是一个向量，在stata数据库中应对应一个变量——最起码地，你也应该清楚：残差这个变量的观测值要与自变量及因变量的观测值对应。

statsby已经改变了原数据库，不可能再生成残差（变量）。

其次，若你只用一步predict,r，则只会根据最后得到的系数估计计算残差，而你本想分类别得到系数估计，这样每类别的残差本应根据本类别系数估计计算。

简单说，你有的各种困惑，首先不是是否会使用stata的问题，而是理论上还不清楚相关概念的问题（这些问题不解决，使用各种软件都仍有困惑）。

其次，你要自己去尝试了解stata命令的结果，每一命令都有自己的逻辑线条，当然，理解这些逻辑线条，首先还是要把基本的统计概念与原理弄清楚。

很多原理上的东西不清楚，别人再告诉应该使用某某命令，你也仍然不理解该命令得到的各种结果的意义。

扫码加我拉你入群

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

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

sharlagao

2011-5-4 00:39:08

谢谢Sungmoo的耐心解答。翻书后还是不知道da表示什么？Nda表示什么？能否再提醒一下？谢谢！

扫码加我拉你入群

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

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

sungmoo

2011-5-4 01:06:29

sungmoo 发表于 2011-5-4 00:06
sharlagao 发表于 2011-5-3 02:04 sungmoo: 我试用您的命令，有用了版上讨论的命令predict da, res ，结果不一样。请问：da表示什么？Nda表示什么？
predict a ,r
a可以随便命名，生成残差变量a。

sharlagao 发表于 2011-5-4 00:39 翻书后还是不知道da表示什么？Nda表示什么？能否再提醒一下？

前面不是说过了吗？predict后面的变量名你可以随便换成别的（用a、b、c、d、da都无所谓），生成的都是一样内容的东西。

扫码加我拉你入群

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

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

sharlagao

2011-5-4 12:43:50

我不知道的是，除了残差da外，还能预测个什么，那个nda是预测的什么？真不好意思，让您反复回答这个简单的问题。

扫码加我拉你入群

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

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

sharlagao

2011-5-4 12:46:58

我只知道因变量才有残差，难道难道自变量也能有残差？

扫码加我拉你入群

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

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

sungmoo

2011-5-4 13:28:44

sharlagao 发表于 2011-5-4 12:46 我只知道因变量才有残差，难道难道自变量也能有残差？

因变量-自变量*系数估计=残差

你的“因变量才有残差”是什么意义呢？

扫码加我拉你入群

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

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

sungmoo

2011-5-4 13:30:20

sharlagao 发表于 2011-5-4 12:43 我不知道的是，除了残差da外，还能预测个什么，那个nda是预测的什么？真不好意思，让您反复回答这个简单的问题。

da与nda只是你起的变量名。

predict能预测什么，可以用help predict。

比如，",r"、",xb"等等选项的作用。

扫码加我拉你入群

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

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

sharlagao

2011-5-5 22:05:09

请问下面这个predict 的da是残差，但nda,是什么？非常感谢！
drop g nda da
egen g=group(year ind)
gen nda=.
gen da=.
forvalue i=1/220 {
reg a b c d if g==`i'
predict nda1 if g==`i'
predict da1 if g==`i', res
replace nda=nda1 if g==`i'
replace da=da1 if g==`i'
drop nda1 da1
}

扫码加我拉你入群

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

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

sungmoo

2011-5-6 15:18:25

sharlagao 发表于 2011-5-5 22:05 predict nda1 if g==`i'

nda1是因变量的预测值。

扫码加我拉你入群

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

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

llg79

2011-5-8 20:45:46

看看，学学，谢谢

扫码加我拉你入群

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

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

wangying1778

2011-5-8 22:56:41

xiexie,tingxiangjie

扫码加我拉你入群

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

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

sharlagao

2011-5-23 10:20:53

谢谢sungmoo的耐心

扫码加我拉你入群

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

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

309

2012-4-7 22:01:11

那只是个名字
关键是逗号后面不一样
你做一次
browse 下
对比就知道了

扫码加我拉你入群

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

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

quanyi88888

2012-5-22 10:29:45

sungmoo 发表于 2011-5-2 23:27
*你首先看一下它运行的结果是否可能得到残差。

egen g=group(code year)

请问一下，如果年份不连续，这个命令也凑效吗？

扫码加我拉你入群

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

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

水蓝色qing

2012-6-15 17:19:15

弱弱地问一下，loc n=r(max)是什么意思呢？

扫码加我拉你入群

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

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

水蓝色qing

2012-6-16 09:37:53

sungmoo 发表于 2011-5-2 23:27
*你首先看一下它运行的结果是否可能得到残差。

egen g=group(code year)

请教版主，当输入完循环命令之后，提示是invalid syntax，不知道这是哪里出错了呢/

扫码加我拉你入群

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

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

水蓝色qing

2012-6-16 10:30:58

sungmoo 发表于 2011-5-2 23:27
*你首先看一下它运行的结果是否可能得到残差。

egen g=group(code year)

我的程序是

use "G:\15\dta\total.dta", clear
egen g = group(Stkcd)
loc n=r(max)
g e=.
forv i=1/`n'{
reg Dretwd Dretwdtl if g==`i'
predict a if e(sample),r
replace e=a if e(sample)
drop a
}

结果就会提示invalid syntax。
不知道是哪里出错了希望版主帮下忙呀，嘿嘿

扫码加我拉你入群

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

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

sungmoo

2012-6-16 23:05:36

*少了一步
egen g = group(Stkcd)
su g
loc n=r(max)

扫码加我拉你入群

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

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

llouy

2012-11-10 18:03:07

版主您好!看了您的这版帖子相当受用。不过还得请教版主，如果我要获得的不是回归的残值，而是回归的R SQUARE，不知道命令应该怎么改呢，，先谢谢您的热心帮助阿。。。

扫码加我拉你入群

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

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

sungmoo

2012-11-10 19:43:34

如果我要获得的不是回归的残值，而是回归的R SQUARE

https://bbs.pinggu.org/thread-1018129-1-1.html

https://bbs.pinggu.org/thread-605094-1-1.html

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群