data life; input age line @@;
label age='Age at Death';
label line = 'Length of Lifeline (cm) ';
  x=line-9.198; y=age-66.66; xsq=x*x; ysq=y*y; xy=x*y;
  Title 'Wilson & Mather JAMA 229 (1974)';
  cards;
19  9.75   40  9.00   42  9.6    42  9.75   47 11.25   49  9.45
50 11.25   54  9.00   56  7.95   56 12.00   57  8.10   57 10.20
58  8.55   61  7.20   62  7.95   62  8.85   65  8.25   65  8.85
65  9.75   66  8.85   66  9.15   66 10.20   67  9.15   68  7.95
68  8.85   68  9.00   69  7.80   69 10.05   70 10.50   71  9.15
71  9.45   71  9.45   72  9.45   73  8.10   74  8.85   74  9.60
75  6.45   75  9.75   75 10.20   76  6.00   77  8.85   80  9.00
82  9.75   82 10.65   82 13.20   83  7.95   86  7.95   88  9.15
88  9.75   94  9.00
;
proc print; sum x y xsq xy ysq; 
proc reg; model Age=line; 
run; quit; 
ods html file='lifeline.html'; 
goptions device=activex; 
proc gplot; plot Age*line/haxis=5 to 15; 
symbol1 v=dot i=rl c=green; 
run; quit; 
ods html close; 
run;