Data Reg1; Input X Y; 
xdev=X-5; ydev=Y-6; xy=xdev*ydev;xsq=xdev*xdev; ysq=ydev*ydev; 
datalines; 
4 5
7 7
5 6
4 6
6 .
 ; 
proc print; sum xy xsq ysq xdev ydev; 
 proc reg; model Y = X/P clm cli; 
proc gplot; plot Y*X/haxis = 3 to 8 vaxis=4 to 8; 
symbol1 v=dot c=red i=rl;  
run;