geeモデルを交換可能なvar-cov行列に適合させてから、結果のモデルに対して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 Estimator(Empirical)は、EmpiricalOptionを使用したProcMIXEDで簡単に実装できることを知っています。上で定義したすべての参照グループのために、GENMODを使用する必要があります。とにかく、上記のGENMODから取得した残差に基づいてHuberWhiteサンドイッチ推定器を実行するマクロを介して結果を渡すことができますか?
私はあなたの助けに感謝します。-Sepehr