topic Re: proc surveyreg in New SAS User

proc surveyreg

knighsson — Sat, 30 Jan 2021 22:08:23 GMT

Hello,

When we try to use linear regression, we need to do the diagnostics first. So I wonder when I use "proc surveyreg" to investigate the linear regression, should I use "proc reg" statement to do the diagnostics test (residual plot ...) or should I use "proc surveyreg" statement to do the diagnostics? If I should use "proc surveyreg" statement, how can I get the residual plot?

Thank you

Re: proc surveyreg

ballardw — Sun, 31 Jan 2021 01:16:50 GMT

Which specific diagnostics are you wanting?

If the data comes from a complex sample and requires Proc Surveyreg then Reg is very likely not the tool to use for diagnostics. Proc Reg would not properly apply the weights from a complex sample and so the residuals would be extremely likely to be incorrect.

Perhaps Surveymeans or Surveyfreq would come into play.

Re: proc surveyreg

knighsson — Sun, 31 Jan 2021 16:17:40 GMT

Thank you so much!

I want to check the assumption of linear regression including residual normality, outliers, linear relationship between dependent variable and independent variables, homogeneity, and multicollinearity. could you please tell me how to do those test in proc surveyreg statement?

Thank you!

Re: proc surveyreg

PaigeMiller — Sun, 31 Jan 2021 16:56:30 GMT

You can output the residuals from PROC SURVEYREG using the OUTPUT statement, and then you can plot them to take care of "residual normality, outliers, linear relationship between dependent variable and independent variables, homogeneity". I think the COVB option in the MODEL statement would address multicollinearity.

Re: proc surveyreg

knighsson — Sun, 31 Jan 2021 17:49:10 GMT

Thank you so much, that would be very helpful!

Re: proc surveyreg

knighsson — Sun, 31 Jan 2021 18:45:01 GMT

Hello, one more question.

I used following code to try to get the residual:

ods graphics on;

PROC SURVEYREG DATA= nh.outcomes nomcar;

STRATA sdmvstra;

CLUSTER sdmvpsu;

CLASS alpha16 age RIAGENDR PIR SDDSRVYR RIDRETH1;

WEIGHT glucwt4yr;

DOMAIN eligible;

model BMXBMI= age RIAGENDR PIR SDDSRVYR RIDRETH1 EIEER totalcounts alpha16/ stb adjrsq clparm solution vadjust=none COVB ;

lsmeans alpha16/ lines adjust=tukey;

output out=bmi p= predict r = residual ;

run;

quit;

ods graphics off;

But could you please tell me how to plot the residual? I tried proc forecast and proc sgplot statements, but they are not working

proc forecast data=bmi

out=pred outfull outresid;

id seqn;

var age RIAGENDR PIR SDDSRVYR RIDRETH1 EIEER totalcounts alpha16;

run;

proc sgplot data=pred;

where _type_='RESIDUAL';

needle x=seqn y=age RIAGENDR PIR SDDSRVYR RIDRETH1 EIEER totalcounts alpha16 / markers;

run;

Re: proc surveyreg

PaigeMiller — Sun, 31 Jan 2021 22:07:30 GMT

@knighsson wrote:

proc sgplot data=pred;

where _type_='RESIDUAL';

needle x=seqn y=age RIAGENDR PIR SDDSRVYR RIDRETH1 EIEER totalcounts alpha16 / markers;

run;

The NEEDLE statement allows only one variable after Y=. But I doubt you really want a NEEDLE here, looking at residuals is usually done via scatter plots, so you can use the SCATTER statement.

Regarding _type_='RESIDUAL', you need to look (with your own eyes) inside the data set that is created named BMI (it is not named PRED) and see how the data set is structured, that will identify if you need a WHERE statement and what the WHERE statement should say; and it will identify the variable names you can use. Essentially, if you look at BMI with your own eyes, you will see everything you need to code some sort of residual plot.