proc sgplot data=faithful;
scatter x=eruptions y=waiting;
run;

proc reg data=faithful;
model waiting=eruptions / dw;
run; quit;

proc reg data=faithful;
model waiting=eruptions / clm cli clb;
run; quit;