1. proc reg data=one;

model: y=x1 x2 x3/vif;run; quit;

2. proc autoreg

3. 异方差检验

proc model data=hetero1;
      parms a1 b1 b2;
      exp = a1 + b1 * inc + b2 * inc2;
      fit exp / white pagan=(1 inc inc2)
      out=resid1 outresid;
   run;
   quit;