topic proc glm in Statistical Procedures

proc glm

Lkahr — Wed, 31 Mar 2021 12:11:02 GMT

Hi,

My sample includes 779 observations and many different variables.

I have fittet a generel linear regression model.

Due to missing data in some variables the model only use 702 observations.

proc glm data=myapaf;
class jobkatreg genderkat ;
model ra=durationew age bmi family_status mg_adl_total isi_total genderkat jobkatreg /solution clparm;
run;

I need to describe these 702 persons (mean age etc.)

Does anyone know how to generate a dataset of these 702 persons? Preferable in the code above.

Thank you !

Re: proc glm

PaigeMiller — Wed, 31 Mar 2021 12:40:51 GMT

proc glm data=myapaf;
    class jobkatreg genderkat ;
    model ra=durationew age bmi family_status mg_adl_total isi_total genderkat jobkatreg  /solution clparm;
    output out=predicteds p=predicted;
run;
data want;
    set predicteds(where=(not missing(predicted)));
run;

The output data set from PROC GLM contains all of your original variables, and a new column named PREDICTED which contains the predicted values. The variable PREDICTED will be missing if any one of the terms in the model are missing, so those are the records you don't want, or to say it the opposite way, you do want the records where PREDICTED is not missing.

Re: proc glm

jiltao — Wed, 31 Mar 2021 13:53:50 GMT

Like @PaigeMiller pointed out, you can use the OUTPUT statement to get what you wanted. But I would use the residuals rather than the predicted values to get the observations that are used in the analysis. This is because if the response variable is missing and none of the independent variables is missing, then this observation will not be used in the analysis but you would still get the predicted value for this observation. If you use the residuals, (output out=out r=residuals;) you would correctly identify those observations included in the analysis (residuals not missing).