小弟没学过stata，但是需要用stata跑一个程序，数据、代码都有，可否帮看下在怎么操作

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

小弟没学过stata，但是需要用stata跑一个程序，数据、代码都有，可否帮看下在怎么操作。
具体就是比较两个回归，第一个回归rate（因变量） credit 与D-T-I Ratio两个自变量，第二个回归是rate（因变量）fico与D-T-I Ratio两个自变量。
两个回归R平方显著性和两者（当自变量为credit时和当自变量为fico时）回归系数差异。
比较两个回归R平方显著性程序已经上传，数据一并上传，希望高人不吝赐教。
程序是ado，只需要把txt后缀改为ado即可，关键求教在command窗口运行什么命令。。。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Stata tata 怎么操作 Command Credit command credit 因变量 自变量 程序

顶一下

太不厚道了，下载还要流量费

陈信研究员 发表于 2015-11-19 10:23
太不厚道了，下载还要流量费

啊？我没设置这个东西。不过不知高人对stata这个有所了解吗？
我把程序贴下来吧
希望能帮我一下，小小论坛币不成敬意

******************************
* 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
我已经查阅了文献然后给你做了具体示范和最后结果，你自己查阅吧！

多谢啦！谢谢