topic Re: Output dataset with parameter estimates from proc genmod in SAS Procedures

Output dataset with parameter estimates from proc genmod

Dinurik — Fri, 07 Feb 2020 22:56:35 GMT

Hi,

I have been trying to output a dataset from proc genmod that would contain parameter estimates and their associated standard errors. I need to further use them to get pooled estimates from 5 imputed datasets. Here is my code:

proc genmod data=mydata;

weight stweight;

class id var1 var2 var3 var4;

model y = var1 var2 var3 var4 / dist=N link=id type3;

repeated subject=id / type=ind;

ods output parameterestimates=parms parminfo=info;

run;

When I run it, SAS almost crashes and then tells me that the output dataset has not been created. Could anyone please tell me what I am doing wrong?

Thanks!

Re: Output dataset with parameter estimates from proc genmod

PGStats — Fri, 07 Feb 2020 23:20:27 GMT

Please show the messages from the SAS log.

Re: Output dataset with parameter estimates from proc genmod

Dinurik — Mon, 10 Feb 2020 16:21:46 GMT

WARNING: Output 'ParameterEstimates' was not created. Make sure that the output object name,
label, or path is spelled correctly. Also, verify that the appropriate procedure
options are used to produce the requested output object. For example, verify that
the NOPRINT option is not used.

Re: Output dataset with parameter estimates from proc genmod

PaigeMiller — Mon, 10 Feb 2020 16:28:09 GMT

@Dinurik wrote:

WARNING: Output 'ParameterEstimates' was not created.

The PROC GENMOD documentation discusses this explicitly. You need the PRINTMLE option

Re: Output dataset with parameter estimates from proc genmod

Dinurik — Mon, 10 Feb 2020 16:41:20 GMT

Thanks. The SAS explanation lacks a bit of detail. Wasn't clear whether PRINMLE belongs to the model statement or repeated statement. It should be in the repeated statement in my case. This code worked.

proc genmod data=mydata;
class id var1 var2 var2 var4;
model y = var1 var2 var3 var4 ;
repeated subject=id / type=ind printmle ;
ods output ParameterEstimates=gmparms;
run;