我正在尝试使用调查包对马拉维2010年DHS数据集中的家庭牛拥有量(变量hv246a)进行小面积估算。我遇到了svyby函数的问题。
父DHS数据集是家庭数据集(MWHR61FL.DTA)。聚类变量为hv021,家庭体重(* 1,000,000)为hv005,采样层为hv021,牛变量为hv246a。
我的代码是:
dhs.des<-svydesign(id=data$hv021, weights=data$hv005/1000000, strata=data$v022, data=data)
cattle<-svyby(~hv246a, ~hv021, dhs.des, svymean)
此数据集中有849个群集(hv021)和23,020户(单位)。
输出的标准错误是:
> summary(cattle$se)
Min. 1st Qu. Median Mean 3rd Qu. Max.
0.000e+00 0.000e+00 0.000e+00 4.674e-17 3.556e-17 1.082e-15
加权估计(cattle $ cattle)看起来不错,但是这些标准错误对我而言似乎是无意义的。
谢谢, 朱莉安(Julianne)