我对Stata相当新,所以我的一些问题可能非常基本,但我感谢任何帮助。我的数据示例如下
wage surveyweights Constant Sex
32 14 .56 1
56 45 .96 1
77 88 .25 0
我有调查数据,我正试图引导女性的平均工资差异。我想找出平均结果与[调查权重]和[调查权重*常数]的差异手动计算
mean x [pw=wt] if sex==1
mat x1=e(b)
mean x [pw=wt*constant] if sex==1
mat x2= e(b)
mat dd=x1-x2
**my outcome of interest is the point estimate and bootstrapped SE of dd
我写了以下程序,但我没有得到我想要的结果。特别是,我的SE栏显示为空白,我对平均值和均值差异的估计值与手动计算得出的结果不一致。
program define meandiff, eclass properties (svyb)
args vars
mean `vars'
mat x1=e(b)
mean `vars' [pw=constant]
mat x2=e(b)
mat dd = x1-x2
ereturn scalar dd=e1(dd,1,1)
end
local vars wage
svy bootstrap e(dd), subpop(sex): means `vars'
我已经使用bootstrap权重来设置我的数据。我的问题如下:
当我输入svy:意思是'vars'时,程序似乎开始运行bootstrap复制,当我输入svy bootstrap时:意味着'vars',程序似乎也开始复制。这两个命令有什么区别?
当我mean x
使用常规调查权重时,是否需要执行[pw=wt]
或svy
命令会自动应用调查权重?
如果我必须在第一个意思中写[pw=wt]
,那么当我计算第二个时,我需要创建一个名为gen wtxcons = wt * constant
的变量来做[pw=wtxcons]
意思?
如何计算引导程序SE和我对感兴趣结果的点估计值,即均值的差异。为什么我的点估计与我的手动计算不匹配?