使用前缀在Stata中存储二项式置信区间结果
我试图计算每年完成治疗的患者比例的95%二项Wilson置信区间（数据集为每个人列出的行） 我想将结果存储到矩阵中，这样我就可以使用使用前缀在Stata中存储二项式置信区间结果,stata,confidence-interval,binomial-coefficients,Stata,Confidence Interval,Binomial Coefficients,我试图计算每年完成治疗的患者比例的95%二项Wilson置信区间（数据集为每个人列出的行） 我想将结果存储到矩阵中，这样我就可以使用putexcel命令将结果导出到现有的Excel电子表格中，而无需更改工作表的格式。我创建了一个二进制变量dscomplete\u binary，如果治疗未完成，则该变量为0，如果治疗已完成，则为1 我尝试了以下方法： bysort year: ci dscomplete_binary, binomial wilson level(95) 这使得每年的产量具有95
putexceldscomplete\u binarybysort year: ci dscomplete_binary, binomial wilson level(95)
statsby year mean LowerCI UpperCI
r1 2005 .7031588 .69229454 .71379805
r2 2006 .75532377 .74504232 .7653212
r3 2007 .78125924 .77125096 .79094833
r4 2008 .80014324 .79059798 .80935836
r5 2009 .81860977 .80955398 .82732689
r6 2010 .82641232 .81723672 .83522016
r7 2011 .81854123 .80955547 .82719356
r8 2012 .83497983 .82621944 .8433823
r9 2013 .85411799 .84527379 .86253893
r10 2014 .84461939 .83499599 .85377985
svyset id2
bysort year: eststo: ci dscomplete_binary, binomial wilson level(95)
我认为
postfilecipostfilepostfilebybysort/*make sure no postfile is open, in case a previous run did not close the file*/
cap postclose ci_results
/*create the postfile that will store results*/
postfile ci_results year mean lowerCI upperCI using ci_results.dta, replace
/*loop through years*/
forval y = 2004/2014 {
ci dscomplete_binary if year==`y', binomial wilson level(95)
/*store saved results from ci to postfile. Make sure the post statement contains results in the same order stated in postfile command.*/
post (`y') (r(mean)) (r(lb)) (r(ub))
}
/*close the postfile once you've looped through all the cases of interest*/
postclose ci_results
use ci_results.dta, clear
一旦您将ci_results.dta数据加载到内存中，您就可以应用任何您喜欢的Excel导出命令。我认为
postfilecipostfilepostfilebybysort/*make sure no postfile is open, in case a previous run did not close the file*/
cap postclose ci_results
/*create the postfile that will store results*/
postfile ci_results year mean lowerCI upperCI using ci_results.dta, replace
/*loop through years*/
forval y = 2004/2014 {
ci dscomplete_binary if year==`y', binomial wilson level(95)
/*store saved results from ci to postfile. Make sure the post statement contains results in the same order stated in postfile command.*/
post (`y') (r(mean)) (r(lb)) (r(ub))
}
/*close the postfile once you've looped through all the cases of interest*/
postclose ci_results
use ci_results.dta, clear
一旦将ci_results.dta数据加载到内存中，您就可以应用任何您喜欢的Excel导出命令。这是对使用
statsby. sysuse auto, clear
(1978 Automobile Data)
. preserve
. statsby mean=r(mean) ub=r(ub) lb=r(lb), by(rep78) : ci foreign, binomial wilson level(95)
(running ci on estimation sample)
command: ci foreign, binomial wilson
mean: r(mean)
ub: r(ub)
lb: r(lb)
by: rep78
Statsby groups
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.....
. list
+----------------------------------------+
| rep78 mean ub lb |
|----------------------------------------|
1. | 1 0 .6576198 0 |
2. | 2 0 .3244076 0 |
3. | 3 .1 .2562108 .0345999 |
4. | 4 .5 .7096898 .2903102 |
5. | 5 .8181818 .9486323 .5230194 |
+----------------------------------------+
. restore
. describe
descripestatsby. sysuse auto, clear
(1978 Automobile Data)
. preserve
. statsby mean=r(mean) ub=r(ub) lb=r(lb), by(rep78) : ci foreign, binomial wilson level(95)
(running ci on estimation sample)
command: ci foreign, binomial wilson
mean: r(mean)
ub: r(ub)
lb: r(lb)
by: rep78
Statsby groups
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
.....
. list
+----------------------------------------+
| rep78 mean ub lb |
|----------------------------------------|
1. | 1 0 .6576198 0 |
2. | 2 0 .3244076 0 |
3. | 3 .1 .2562108 .0345999 |
4. | 4 .5 .7096898 .2903102 |
5. | 5 .8181818 .9486323 .5230194 |
+----------------------------------------+
. restore
. describe
descripestatsby:help cibysort year:preserverestorepreserverestorepreservestatsbyrestorestatsby:help cibysort year:preserverestorepreserverestorepreservestatsbyrestore