搜索
人大经济论坛 附件下载

附件下载

所在主题:
文件名:  1_Missing20.xlsx
资料下载链接地址: https://bbs.pinggu.org/a-2994928.html
附件大小:
68.58 KB   举报本内容
各位大神:
请教一个问题,我在用stata进行模拟数据的分析。在使用twofold 程序完成多重插补及分析后,输出了多重插补方差信息表和回归分析结果表。我想要在这两个表格里分别提取相应的7个数值,如图红框内。
因为是模拟研究,循环内500次分析,就有500次的结果。我需要将这7个数字的500次结果提取出来求平均值和标准差。
求大神指点,如何将这部分数据提取出来。
附数据:,
附代码:
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)

}





    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载,请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品,拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知,将积极的采取必要措施;同时,本站也将在技术手段和能力范围内,履行版权保护的注意义务。
(如有侵权,欢迎举报)
二维码

扫码加我 拉你入群

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

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

GMT+8, 2025-12-29 19:48