[已解决] 求助:可以自动选择确切概率法的函数

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

看到一个程序 写的是 tabauto `var1' `var2'  `in' `if', missing  /*automatically choose exact test*/

我查了help tabauto 可是提示没有这个函数,上谷歌搜索了一下,也没有相关函数,请问各位老师有这个函数吗?  

如果没有,那么我想做四格表频数检验,如何能让stata 根据最小格子数 自动选择pearson卡方或exact test呢?

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：已解决 Automatic missing pearson choose 谷歌搜索 pearson 程序 格子 如何

没有见过你说的这个
你把更多的贴出来
是不是你抄错命令了

我猜测这里执行的是Fisher's exact Probability test。在stata中可以使用tabi命令检验之。详见帮助文件。
如下文献使用了这一方法：
连玉君, 李丹. 财务比率的行业 特征: 差异性及收敛性——基于中国上市公司的非参数检验[J]. 南方经济, 2006, (3): 81-90.

问了一下作者,说是自已写的ado, 难怪找不到

keyzhi 发表于 2009-6-21 16:29
问了一下作者,说是自已写的ado, 难怪找不到

提问要说清楚，不要藏着掖着

eblog 发表于 2009-6-21 18:38

keyzhi 发表于 2009-6-21 16:29
问了一下作者,说是自已写的ado, 难怪找不到

提问要说清楚，不要藏着掖着

不好意思,因为这是别人写的程序,未经同意,不敢擅自发布网上,起初以为只是一个函数而已,想通过自己努力搞懂它,
后来实在忍不住就跟作者交流了一下,他说是他附加的一个ado,希望eblog版主能理解,同时说声抱歉.

征求作者南京医科大学张汝阳同意后,将他写的tabauto 发布,在此感谢计量经济版各位朋友老师,同时感谢张汝阳同志!
program drop _all
program define tabauto,rclass/*this is the sub_program for frept,only for rawdata*/
        syntax varlist(min=2) [in] [if],[missing]       
        token "`varlist'"
        local r_var="`1'"
        local c_var="`2'"
        qui tab `r_var' `c_var' `if' `in' `byoption',chi2

        if "`r(p)'"!=""{
                preserve
                        if "`missing'"!=""{
                                contract `r_var' `c_var' `if' `in' `byoption',freq(fre_q)
                        }
                        else{
                                contract `r_var' `c_var' `if' `in' `byoption',freq(fre_q) nomiss
                        }
       
                        qui sum fre_q
                        local sample=`r(sum)'
                        tempvar rowsum colsum
                        egen `rowsum'=total(fre_q),by(`r_var')
                        egen `colsum'=total(fre_q),by(`c_var')
                        qui sum `rowsum'
                        local rowmin=`r(min)'
                        qui sum `colsum'
                        qui local colmin=`r(min)'
                        local min_t=(`rowmin'*`colmin')/`sample'/*theoretical frequence*/
                restore,preserve

                if (`min_t'<1 | `sample'<40){
                        return local exact_auto="e"
                }
                else{
                        return local exact_auto=""
                }
        }
end

不要在意啊，我只是建议大家提问要说清楚，最好以下面的格式提问：
1.清楚明确的表明想问什么？
2.提问之前你自己做了什么努力？
3.有数据的至少要提供部分数据（Stata中可以很简单的对数据进行抽样）
4.涉及版权而与问题高度相关的东东，不一定要贴到论坛上来，但最好写清来源。

这样子提问可以让别人很快了解你的问题，而不至于在猜测你到底在问什么上花费心思。