topic Re: Reproduce Stata GLM is SAS in SAS Procedures

Reproduce Stata GLM is SAS

JonKetchup — Wed, 07 Sep 2022 22:01:16 GMT

I am trying to reproduce a GLM in SAS that I originally created in Stata. Although my estimates are exactly the same, the standard errors are not, and I'm not quite sure why. Note that cost is continuous while treatment, var1, and var2 are all categorical. Here is both codes:

SAS

proc genmod data=data;
 class clinician treatment var1 var2 / param=ref;
 model cost = treatment var1 var2/ link=LOG dist=gamma type3;
repeated SUBJECT=clinician/ PRINTMLE TYPE=ind;
run;

Stata

glm cost i.treatment i.var1 i.var2, f(gam) l(log) cluster(clinician)
testparm i.treatment

Appreciate the help

Re: Reproduce Stata GLM is SAS

PaigeMiller — Wed, 07 Sep 2022 13:50:27 GMT

Are the degrees of freedom the same? Show us.

Re: Reproduce Stata GLM is SAS

JonKetchup — Wed, 07 Sep 2022 13:59:50 GMT

Yes, DFs are all the same along with the estimates. Only the standard errors (and thus the p-values) are different)

Re: Reproduce Stata GLM is SAS

PaigeMiller — Wed, 07 Sep 2022 14:22:23 GMT

How far apart are these estimates of standard errors? Give us a specific case, what are the actual numbers from SAS and what is the actual number in STATA?

What is the estimate of Scale from each program?

Re: Reproduce Stata GLM is SAS

JonKetchup — Wed, 07 Sep 2022 14:44:12 GMT

The estimate for treatment is -0.0356. In SAS, the standard error is 0.1301 while in Stata, the standard error is 0.1469.

I'm not sure how to find the scale in Stata, but the scale is SAS is 0.5432.

Re: Reproduce Stata GLM is SAS

PaigeMiller — Wed, 07 Sep 2022 14:41:43 GMT

I don't understand, it seems as if you have not answered my question, or you did answer my question but worded it confusingly.

The standard error for treatment in SAS is _________________

The standard error for treatment in STATA is _________________

Re: Reproduce Stata GLM is SAS

JonKetchup — Wed, 07 Sep 2022 14:51:13 GMT

My apologies. My last reply had a few errors. It is rewritten, but to again answer your question:

The standard error for treatment in SAS is 0.1301.

The standard error for treatment in STATA is 0.1469.

Re: Reproduce Stata GLM is SAS

PaigeMiller — Wed, 07 Sep 2022 15:00:22 GMT

Ok, so at this point, I have to say I don't really have enough familiarity with STATA to help further, but the problem (in my opinion) is either

the SAS algorithm estimates scale and standard errors differently than STATA algorithm does; or
the models specified are not identical (despite the matching estimates) because of the REPEATED command in PROC GENMOD produces different estimates of scale and different standard errors than STATA model specification does

Re: Reproduce Stata GLM is SAS

JonKetchup — Wed, 07 Sep 2022 21:59:58 GMT

Figured it out. The answer is twofold.

1) Genmod and Stata estimate the scale parameter and correlation structure slightly differently. Details are given on xtgee-related FAQ.

2) Genmod by default produces score type III tests while Stata produces Wald type III tests. Wald tests can be requested for SAS, but score tests cannot be produced by Stata.