2012-07-12 2 views
0

저는 교환 가능한 var-cov 행렬로 gee 모델을 적합하게하고 편향된 결과를 방지하기 위해 결과 모델에 Huber-White 샌드위치 평가기를 실행하고 싶습니다. 내 GEE 모델 내 코드는 다음과 같습니다 :후버 - 화이트 샌드위치 추정기 SAS?

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; 

나는 후버 - 화이트 샌드위치 추정 (실증)는 쉽게 실증 옵션과 혼합 PROC에서 구현 될 수 있다는 것을 알고있다. 위에서 정의한 모든 참조 그룹 때문에 GENMOD를 사용해야합니다. 어쨌든 위의 GENMOD에서 얻은 잔차를 기반으로 한 HuberWhite 샌드위치 측정기를 수행하는 매크로를 통해 결과를 전달할 수 있습니까?

감사합니다. -Sepehr

답변

0

proc GENMOD에서 사용 가능한 COVB 옵션을 사용하여 경험적 매개 변수 공분산 행렬을 사용하는 것이 한 가지 방법입니다. 경험적 공분산 행렬 추정기 (견고한 분산 추정기 또는 샌드위치 추산 기 또는 Huber-White 방법이라고도 함)를 사용하려면 proc genmod에 반복 명령문에 covb 옵션을 추가해야합니다.

repeated subject={subject id}/covb;