全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
7786 8
2011-12-08
二维码

扫码加我 拉你入群

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

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

全部回复
2011-12-8 23:35:15
用STATA就可以做Vuong检验.其步骤好像是:

1. 从网上下载一个vuong.ado文件http://personal.anderson.ucla.edu/judson.caskey/data.html1. 将网页保存为vuong.ado [按纯文本格式保存,但是扩展名添加为ado,这一步很重要;或者到上一级目录(页面)直接点击相应连接后“另存为”]
将"vuong.ado" 保存到 STATA目录下

2. 将两个模型各自回归,保存其回归结果

3. 再将结果带入vuong检验

楼主可以参看一些STATA的教程.进行学习.
二维码

扫码加我 拉你入群

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

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

2011-12-8 23:42:37
******************************
* vuong.ado
* Judson A. Caskey
* UCLA
* 26-Apr-2007
*
* Computes Vuong (1989 Econometrica) test
* of two non-nested regressions as implemented
* and described in Dechow (1994 Journal of
* Accounting and Economics)
*
******************************

cap prog drop vuong
prog define vuong, rclass

syntax [anything]

tempname mod1 mod2 rss1 rss2 rsq1 rsq2 n zstat pval nbad
tempvar s1 s2 e1 e2 m

if `"`: word count `anything''"' ~= "2" {
        di as err "You must specify two distinct models"
        exit
        }

est_expand `"`anything'"', min(1) max(2)
local `mod1' : word 1 of `r(names)'
local `mod2' : word 2 of `r(names)'
if "``mod1''"=="``mod2''" {
        di as err "You must specify two distinct models"
        exit
        }

quietly estimates restore ``mod1''
local `rss1'=e(rss)
local `rsq1'=e(r2)
quietly _predict `e1' if e(sample), resid
quietly gen `s1'=e(sample)

quietly estimates restore ``mod2''
local `rss2'=e(rss)
local `rsq2'=e(r2)
quietly _predict `e2' if e(sample), resid
quietly gen `s2'=e(sample)

quietly count if `s1'==1 & `s2'==1 & ~missing(`e1',`e2')
local `n'=r(N)

quietly gen `m'=log(``rss1''/``rss2'')/2 + ``n''*((`e1'^2)/``rss1'' - (`e2'^2)/``rss2'')/2 if `s1'==1 & `s2'==1 & ~missing(`e1',`e2')

quietly reg `m'

local `zstat'=-sqrt((e(N)-1)/e(N))*_b[_cons]/_se[_cons]
local `pval'=(1-normal(abs(``zstat'')))*2

di as result _column(20) "   Model 1"  _column(35) "   Model 2"
di as result "R-Squared" _column(20) %10.4f ``rsq1'' _column(35) %10.4f ``rsq2''
di as result ""
di as result "Vuong Z-Statistic" _column(20) %10.4f ``zstat''
di as result "  p-value" _column(20) %10.4f ``pval''

return scalar r2_1=``rsq1''
return scalar r2_2=``rsq2''
return scalar Z=``zstat''
return scalar p=``pval''

quietly count if (`s1'==1 & `s2'==0) | (`s1'==0 & `s2'==1)
local `nbad'=r(N)
di as result ""
if ``nbad''>0 di as err "``nbad'' observations in one model but not the other"

end

二维码

扫码加我 拉你入群

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

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

2011-12-8 23:42:53
没有单独的软件!
。。
二维码

扫码加我 拉你入群

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

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

2011-12-9 00:02:20
lwfaaa3 发表于 2011-12-8 23:42
******************************
* vuong.ado
* Judson A. Caskey
哇哦
二维码

扫码加我 拉你入群

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

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

2011-12-9 01:23:54
问题已解决。版主真是好人。谢谢!
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群