Estpost for summary statistics with complex design survey data

[Abeya Mokhtar]

Hi ,

I am working on quarterly Labor Force Survey from 2008 to 2010 .
I would like to create a table for summary statistics for 4 subpopulation of my data ( Male before, Male After , Female before and Female After ) where After is the variable Post.
My dataset is already svyset , but the estpost tabstat doesn't seem to support svy , only the tabulate but not tabstat.

The code I used is below ;

** generate the column identifiers by flag
gen flag=1 if sex==1 & post==0
replace flag=2 if sex==1 & post==1
replace flag=3 if sex==0 & post==0
replace flag=4 if sex==0 & post==1
estpost tabstat age mar_st reg* edu2 edu3 edu4 edu6 edu8 lfs unemployed [aw=int_wt], by(flag) statistics(mean sd N) columns(statistics) nototal
esttab using table1.rtf, replace main(mean) aux(sd) unstack nomtitle nonumber label nogaps title("Discriptive statistics for sample ")

and the results were as below;

Discriptive statistics for sample

	1	2	3	4
age	27.10	26.42	27.53	26.72
	(19.26)	(19.21)	(19.02)	(19.02)
Marital status	2.049	2.048	2.423	2.381
	(0.986)	(0.984)	(1.403)	(1.409)
region	4.142	4.155	4.127	4.137
	(1.813)	(1.809)	(1.812)	(1.803)
region==Greater Cairo	0.157	0.155	0.159	0.156
	(0.364)	(0.362)	(0.365)	(0.363)
region==Alex & Suez Canal	0.0769	0.0770	0.0774	0.0773
	(0.266)	(0.267)	(0.267)	(0.267)
region==Urban Upper Egypt	0.0996	0.0964	0.102	0.0998
	(0.299)	(0.295)	(0.302)	(0.300)
region==Urban Lower Egypt	0.0873	0.0910	0.0873	0.0910
	(0.282)	(0.288)	(0.282)	(0.288)
region==Rural Upper Egypt	0.309	0.307	0.307	0.308
	(0.462)	(0.461)	(0.461)	(0.462)
region==Rural Lower Egypt	0.250	0.255	0.252	0.252
	(0.433)	(0.436)	(0.434)	(0.434)
region==Frontier	0.0194	0.0186	0.0162	0.0156
	(0.138)	(0.135)	(0.126)	(0.124)
educ_st==Illitrate	0.161	0.149	0.292	0.267
	(0.367)	(0.356)	(0.455)	(0.442)
educ_st==Read and Write	0.139	0.134	0.105	0.0982
	(0.346)	(0.341)	(0.307)	(0.298)
educ_st==Primary or Prep	0.172	0.172	0.144	0.148
	(0.377)	(0.378)	(0.351)	(0.355)
educ_st==Technical/Vocational Diploma	0.172	0.174	0.139	0.148
	(0.377)	(0.379)	(0.346)	(0.355)
educ_st==College Degree	0.0846	0.0880	0.0638	0.0718
	(0.278)	(0.283)	(0.244)	(0.258)
Labor Force	0.729	0.752	0.227	0.228
	(0.445)	(0.432)	(0.419)	(0.419)
All Unemployed	0.0421	0.0383	0.0521	0.0559
	(0.201)	(0.192)	(0.222)	(0.230)
Observations	1008611

mean coefficients; sd in parentheses
^* p < 0.05, ^** p < 0.01, ^*** p < 0.001



My problem is that when I run the following command ;


svy, subpop (if sex==1 & post==0): mean age lfs unemployed reg* ( for any of the 4 subpopulations above ) , I get different results .

The following
Survey: Mean estimation
Number of strata = 154 Number of obs = 852,739
Number of PSUs = 7,427 Population size = 194,685,642
Subpop. no. obs = 273,166
Subpop. size = 62,671,956.6
Design df = 7,273
--------------------------------------------------------------
| Linearized
| Mean Std. Err. [95% Conf. Interval]
-------------+------------------------------------------------
age | 33.70677 .0251904 33.65739 33.75615
lfs | .7654193 .0011021 .7632589 .7675798
unemployed | .0420523 .000682 .0407153 .0433892
region | 4.06441 .0108864 4.04307 4.085751
reg1 | .1657978 .0026677 .1605684 .1710272
reg2 | .0814624 .0011528 .0792026 .0837221
reg3 | .1037393 .0018627 .1000878 .1073908
reg4 | .087979 .0017812 .0844873 .0914708
reg5 | .3117424 .0020266 .3077696 .3157151
reg6 | .2311122 .0018523 .2274812 .2347432
reg7 | .018167 .0015668 .0150955 .0212385
--------------------------------------------------------------
Note: Strata with single sampling unit centered at overall
mean.
.


Any idea why are the results different and how can I get svy like results using estpost ?

Thanks,