0

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

4

1 に答える 1

0

1つの方法は、procGENMODで使用可能なCOVBオプションを使用して経験的パラメーター共分散行列を使用することです。経験的共分散行列推定量(ロバスト分散推定量、サンドイッチ推定量、Huber-White法とも呼ばれます)を使用するには、procgenmodの繰り返しステートメントにcovbオプションを追加する必要があります。

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