全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
4131 6
2015-11-17
小弟没学过stata,但是需要用stata跑一个程序,数据、代码都有,可否帮看下在怎么操作。
具体就是比较两个回归,第一个回归rate(因变量) credit 与D-T-I Ratio两个自变量,第二个回归是rate(因变量)fico与D-T-I Ratio两个自变量。
两个回归R平方显著性和两者(当自变量为credit时和当自变量为fico时)回归系数差异。
比较两个回归R平方显著性程序已经上传,数据一并上传,希望高人不吝赐教。
程序是ado,只需要把txt后缀改为ado即可,关键求教在command窗口运行什么命令。。。
附件列表

data.xlsx

大小:6.1 KB

 马上下载

数据

data97-03.xls

大小:2.51 MB

 马上下载

数据

vuong test.txt

大小:1.93 KB

 马上下载

程序,改为ado

二维码

扫码加我 拉你入群

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

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

全部回复
2015-11-18 16:09:27
顶一下
二维码

扫码加我 拉你入群

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

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

2015-11-19 10:23:24
太不厚道了,下载还要流量费
二维码

扫码加我 拉你入群

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

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

2015-11-19 18:58:46
陈信研究员 发表于 2015-11-19 10:23
太不厚道了,下载还要流量费
啊?我没设置这个东西。不过不知高人对stata这个有所了解吗?
我把程序贴下来吧
希望能帮我一下,小小论坛币不成敬意
二维码

扫码加我 拉你入群

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

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

2015-11-19 19:00:06
******************************
* 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

二维码

扫码加我 拉你入群

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

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

2015-11-19 20:15:35
我已经查阅了文献然后给你做了具体示范和最后结果,你自己查阅吧!
附件列表

详细步骤及结果.rar

大小:162.42 KB

只需: 200 个论坛币  马上下载

操作步骤与结果

二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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