0

我希望能够用可交换的 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 中使用经验选项实现。由于我在上面定义的所有参考组，我必须使用 GENMOD。无论如何，我可以通过宏传递结果吗？ HuberWhite 三明治估计器基于从上面的 GENMOD 获得的残差？

我感谢您的帮助。-塞佩尔

4

1 回答 1

0

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

repeated subject={subject id} / covb;
于 2012-07-18T21:54:11.877 回答