全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
14473 8
2009-03-02
<p>下面是个program,只需注意加粗部分</p><p>program define skewt, rclass<br/>version 10.1<br/>syntax varlist (numeric)[if] [in], [bs saving(string) size(integer 4) reps(integer 1000) replace *]<br/>marksample touse <br/>foreach var of varlist `varlist' { <br/>capture confirm numeric variable `var'<br/> if _rc==0 {<br/>  qui sum `var' if `touse', detail<br/>  di ""<br/>  di in gr _col(5)  "`var'- stats from the sample"<br/>  di ""<br/>  di in gr _col(5) <strong>"N coefficient  = `=sqrt(r(N))'"</strong><br/>  di in gr _col(5) "S-coefficient  = `=r(mean)/r(sd)'"<br/>  di in gr _col(5) "G-coefficient  = `r(skewness)'"<br/>  di in gr _col(5) "Sample mean    = `r(mean)'"<br/>  di ""<br/>  local ratio = (`=sqrt(r(N))')*((`=r(mean)/r(sd)') + ((1/3) * (`r(skewness)') * ((`=r(mean)/r(sd)')^2)) + ((1/(6*((`=sqrt(r(N))')^2)))* (`r(skewness)')))<br/>  return scalar ratio_`var'=`ratio'<br/> }<br/> else {<br/>  di as input "`var'" as text " is not a numeric variable skewness adjusted t-statistic cannot be calculated."<br/> }</p><p>if "`bs'" == "bs" {<br/> if "`saving'" != "" {<br/>  local saving `saving'_`var'<br/> }<br/> bootstrap r(ratio_`var'), saving("`saving'", `replace') reps(`reps') size(`=int(_N/`size')') `options': skewt `var'<br/> }<br/>}<br/>end<br/></p><p>在上述命令中,加粗的部分不理想,因为只能返回no of observation,不管missing还是nonmissing,我想返回no of nonmissing observation,怎么达到,谢谢</p>
二维码

扫码加我 拉你入群

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

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

全部回复
2009-3-2 20:58:00

sum x if x<.

二维码

扫码加我 拉你入群

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

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

2009-3-2 21:20:00

你的意思是改成 qui sum `var' if `touse'&nonmissing, detail?

二维码

扫码加我 拉你入群

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

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

2009-3-2 21:31:00

结果是出现错误r(111)nonmissing not found

我又试了qui sum `var' if `touse'&`var'!=., detail

结果是和没加那个条件一样,返回的还是包括了missing值的数量


二维码

扫码加我 拉你入群

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

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

2009-3-2 21:40:00

sum x等价于sum x if x<.

sum x本来就不会考虑missing吧?

二维码

扫码加我 拉你入群

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

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

2009-3-2 22:02:00

对,sum不考虑missing,谢谢,

我做出来的结果,对每个variable,N coefficient都一样,我怀疑一直用的第一个variable的r(N),可就是查不出程序哪里有问题

对了,这个是用来求Bootstrapped Skeweness adjusted t数值的

二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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