*The program below is for Math 485 HW10.2
* = Multiple Linear and 1D Regression 10.15  
* = Math 583 Applied Robust Statistics 13.21;

options ls = 70;
data cbrain;
infile "a:cbrain.dat";
input obs x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 sex x12;

proc rsquare;
  model sex = x2 x3 x4 x5 x6 x7 x9 x10 / cp;
  
proc reg;
  model sex = x2 x3 x4 x5 x6 x7 x9 x10;
  output out = a p = pred r = resid;
  run;