data movafegh;
input r a b y1 y2 y3;
datalines;
1	1	1	3.5 1.21	4.71
1	1	2	2.33	0.92	3.25
1	1	3	1.71	1.45	3.16
1	2	1	3.4	1.63	5.03
1	2	2	0.78	1.86	2.64
1	2	3	0.8	0	0.8
1	3	1	7.45	3.62	11.07
1	3	2	8.93	4.34	13.27
1	3	3	2.5	0.72	3.22
1	4	1	3.27	3.64	6.91
1	4	2	10.7	4.28	14.98
1	4	3	9.95	3.71	13.66
1	5	1	11.24	4	15.24
1	5	2	11.13	4.56	15.69
1	5	3	9.33	2.6	11.93
2	1	1	3.71	1.51	5.22
2	1	2	2.25	0.82	3.07
2	1	3	1.72	1.83	3.55
2	2	1	4.28	1.66	5.94
2	2	2	1.05	1.51	2.56
2	2	3	0.69	0	0.69
2	3	1	8.67	3.46	12.13
2	3	2	8.61	4.6	13.21
2	3	3	3.02	0.76	3.78
2	4	1	2.87	3.66	6.53
2	4	2	9.98	4.45	14.43
2	4	3	10.61	3.43	14.04
2	5	1	12.02	3.87	15.89
2	5	2	12.33	4.5	16.83
2	5	3	8.52	2.47	10.99
3	1	1	3.29	1.63	4.92
3	1	2	2.41	0.84	3.25
3	1	3	2.41	1.96	4.37
3	2	1	4.24	1.59	5.83
3	2	2	0.69	1.68	2.37
3	2	3	0.7	0	0.7
3	3	1	8.56	3.33	11.89
3	3	2	8.77	4.54	13.31
3	3	3	2.66	0.77	3.43
3	4	1	3.2	3.68	6.88
3	4	2	10.69	4.62	15.31
3	4	3	10.7	3.45	14.15
3	5	1	12.3	3.75	16.05
3	5	2	12.5	4.63	17.13
3	5	3	8.78	2.65	11.43
4	1	1	3.5	1.75	5.25
4	1	2	2.33	1.11	3.44
4	1	3	2.78	1.71	4.49
4	2	1	3.54	1.63	5.17
4	2	2	0.87	2.03	2.9
4	2	3	0.58	0.08	0.66
4	3	1	7.6	3.21	10.81
4	3	2	8.69	4.4	13.09
4	3	3	2.84	0.76	3.6
4	4	1	3.53	3.33	6.86
4	4	2	10.33	4.35	14.68
4	4	3	10.79	3.47	14.26
4	5	1	11.7	3.9	15.6
4	5	2	11.54	4.76	16.3
4	5	3	9.11	2.56	11.67
;
proc reg;
model y1=a b/p;
run;