SAS中的Huber-White Sandwich Estimator？

Question

我希望能够使用可交换的var-cov矩阵拟合gee模型，然后在结果模型上运行Huber-White三明治估算器以防止偏差结果。我的GEE模型代码如下：

Proc GENMOD data = Cohort1ONLY;
class SSID SCHIID0809 Ethnicity(ref = "500") ELLbaseline GENDER freeLunch failedInd
GRADE0809(ref = "3")/param = ref;
Model SSMATH0809 = TRT0809 SSMATH0708 SSENG0708 GRADE0809 ELLbaseline GENDER freeLunch 
ethnicity failedInd; 
repeated subject = SCHIID0809/ type = exch /*corrw: to print the varcov matrix*/; 
Run;

我知道Huber-White Sandwich估算器（经验）可以在Proc MIXED中使用Empirical Option轻松实现。我必须使用GENMOD，因为我已经在上面定义了所有参考组。无论如何，我是否可以通过一个宏来传递结果，该宏基于从上面的GENMOD获得的残差来执行HuberWhite三明治估算器？

感谢您的帮助。 -Sepehr

Answer 1

一种方法是使用proc GENMOD中可用的COVB选项使用经验参数协方差矩阵。为了使用经验协方差矩阵估计器（也称为鲁棒方差估计器，或夹心估计器或Huber-White方法），我们应该将covb选项添加到proc genmod中的重复语句中：

repeated subject={subject id} / covb;