请教一下大侠们怎么做响应分析

shiyiming

研究者为了了解肉雏鸡日粮中不同钙、磷、VD3、锰、锌水平的组合效应，应用5因子部分实施的2次回归正交旋转组合设计，采用玉米-豆粕型基础日粮和432只肉用公雏，研究日粮中锰(Mn)、锌(Zn)、钙(Ca)、有效磷(av.P)、维生素D3(VD3)5因素不同添加水平的组合效应，探讨5因素与指标间的相互关系。试验采用了10个中心点，各因素的最高水平及最低水平如表（1）。

表（1）　因素水平值

Ca(%)(X1) av.P(%)(X2) VD3(X3) Mn (mg/kg) (X4) Zn (mg/kg)
(X5)
最高水平 1.2 0.55 3250 180 180
最低水平 0.8 0.35 1250 80 80

试验结果如表（2）,其中run为该对应的正交组合设计的处理号，y为雏鸡的体重(g)：

表（2） 试验结果

run y run y run y run Y
379.3557 2 359.8637 3 361.3073 4 384.4367
5 357.2201 6 378.6897 7 388.6827 8 366.2155
9 342.1952 10 358.6246 11 363.0592 12 334.5492
13 379.6527 14 369.1787 15 366.3804 16 369.1017
17 380.6571 18 350.2671 19 344.5074 20 361.0806
21 377.8076 22 371.4853 23 372.5273 24 368.0442
25 350.4669 26 391.4660 27 373.3460 28 361.9551
29 362.5661 30 370.6146 31 366.7305 32 371.0900
33 371.5925 34 371.1776 35 373.4910 36 371.0044

（1）做响应面分析，在试验区域内，找出最佳组合及最佳组合时的响应值。
（2）对以上的二次回归模型进行诊断，作出最佳的预测模型。

shiyiming

仿照下面的例子，改一下就可以了
参考书：机械工业出版社 作者：黄燕 书名：SAS统计分析及其应用  Page：240
data source;
        input nbr x1 x2 x3 x4 y ;
        cards;
1 4.5 0.1 40 3 15.9
2 4.5 0.1 50 3 21
3 5.5 0.9 50 3 16.1
4 5.5 0.9 40 3 23.5
5 5.5 0.5 45 4 12.9
6 5.5 0.1 50 5 9.6
7 5 0.5 45 4 21
8 5 0.9 45 4 17.6
9 5 0.5 45 4 21
10 5 0.5 45 3 17.9
11 5 0.1 45 4 23.1
12 4.5 0.9 40 5 10.5
13 4.5 0.5 45 4 26.6
14 5 0.5 50 4 27.4
15 5 0.5 45 4 21
16 5 0.5 40 4 28.9
17 5.5 0.1  40 5 12.2
18 4.5 0.9 50 5 22.4
19 5 0.5 45 5 20.9
;
run;
data source;
        set source;
        x2=x2/100;
        y=y/100;

run;
proc sort data=source;by x1 x2 x3 x4;run;
proc rsreg data=source ;
        model y=x1 x2 x3 x4/lackfit ;
        ridge max;
run;

shiyiming

万分感谢！！！