全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
1747 1
2019-12-02
各位大神:
请教一个问题,我在用stata进行模拟数据的分析。在使用twofold 程序完成多重插补及分析后,输出了多重插补方差信息表和回归分析结果表。我想要在这两个表格里分别提取相应的7个数值,如图红框内。
因为是模拟研究,循环内500次分析,就有500次的结果。我需要将这7个数字的500次结果提取出来求平均值和标准差。
求大神指点,如何将这部分数据提取出来。
附数据:
1_Missing20.xlsx
大小:(68.58 KB)

 马上下载

2_Missing20.xlsx
大小:(68.79 KB)

 马上下载


附代码:
forvalues i=1/2 {
import excel D:\桌面\`i'_20.xlsx, sheet("Sheet 1") firstrow

******产生五个新变量
generate time1 =1
generate time2 =2
generate time3 =3
generate time4 =4
generate time5 =5

**多重插补-5次,包含插补质量分析
twofold,timein(time1) timeout(time3)  table clear saving(D:\办公\课题\数据\用于方法学习\stata-学习用数据\twofoldtrace`i'.dta,replace) base(time1) depmis(responsevar_) depobs( txvar_ age_) indmis( education_1 income_care_n_1 duration_of_illness_n_1 ) indobs(sex_n_1 ) cat(sex_n_1) m(5) ba(3) bw(5) width(1) im keepoutside savetrace(D:\办公\课题\数据\用于方法学习\stata-学习用数据\responsevartrace`i'.dta,replace)

**重置工作路径
cd "D:\办公\课题\数据\用于方法学习\stata-学习用数据\stata-循环练习用数据"

**wild to long format

mi reshape long NO_ ID_LEAN_ cluster_ sex_n_ age_n_ age_ education_ income_care_n_ duration_of_illness_n_ responsevar_ txvar_ time,i(ID) j(timeperiods)

***GEE estimation

mi xtset ID timeperiods
mi estimate,saving(GEE_`i') vartable: xtgee responsevar_ txvar_ i.timeperiods i.sex_n_ age_ income_care_n_ duration_of_illness_n_, family(normal) link(identity) corr(exchangeable)

}


1ecacd0ad68780aab142d683e491d57.png 5808496c1bd551356688d03c6549a7a.png
二维码

扫码加我 拉你入群

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

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

全部回复
-->

二维码

扫码加我 拉你入群

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

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

栏目导航
热门文章
推荐文章

分享

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