各位大神:
请教一个问题,我在用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)
}