options linesize=70 nodate nonumber;
data; input b s;w=1/s**2;s2=s**2;
cards;
4   4
2   5
4   5
8   5
8   5
7   7
7   7
8   8
9   8
11  8
13  8
5   9
5   9
13  9
8   10
14  10
17  10
11  12
19  12
21  12
15  13
18  13
27  13
14  14
16  14
16  15
14  16
19  16
34  16
22  17
29  17
29  18
34  18
47  18
30  19
48  20
39  21
42  21
55  21
56  24
33  25
48  25
56  25
59  25
39  26
41  26
57  27
78  27
64  28
84  28
54  29
68  29
60  30
67  30
101 30
77  31
85  35
107 35
79  36
138 39
110 40
134 40
;
proc reg;model b=s s2/noint;weight w; output out=a p=p r=r;
data c; set a; r=r/s; proc plot; plot r*p/vpos=17;
proc reg; model b=s s2/noint;output out=b p=p1 r=r1;proc plot;
plot r1*p1/vpos=17;
proc print;var p r p1 r1;
run;