在SAS中用quadratic response surface analysis（过程步中用rsreg）方法进行低温、高温与种子萌发率之间的回归分析。具体如下：
title "response surface with guoduxing";
data guoduxing;
input x1 x2 y@;
m=25;
y1=y*4;
output;
label
x1="low temperature"
x2="high temperature";
datalines;
5 5 0
5 5 0
5 5 0
5 5 0
5 10 24
5 10 24
5 10 25
5 10 24
5 15 21
5 15 25
5 15 25
5 15 25
5 20 25
5 20 25
5 20 25
5 20 25
5 25 15
5 25 18
5 25 24
5 25 10
5 30 25
5 30 24
5 30 25
5 30 25
5 35 23
5 35 23
5 35 24
5 35 22
5 40 0
5 40 0
5 40 0
5 40 0
10 10 24
10 10 25
10 10 25
10 10 24
10 15 23
10 15 25
10 15 23
10 15 24
10 20 24
10 20 25
10 20 25
10 20 24
10 25 25
10 25 25
10 25 25
10 25 25
10 30 25
10 30 25
10 30 25
10 30 25
10 35 25
10 35 25
10 35 22
10 35 24
10 40 1
10 40 1
10 40 1
10 40 1
15 15 24
15 15 25
15 15 25
15 15 24
15 20 25
15 20 25
15 20 25
15 20 25
15 25 25
15 25 25
15 25 24
15 25 25
15 30 24
15 30 25
15 30 23
15 30 25
15 35 25
15 35 25
15 35 25
15 35 25
15 40 24
15 40 24
15 40 18
15 40 13
20 20 24
20 20 24
20 20 24
20 20 24
20 25 25
20 25 24
20 25 24
20 25 25
20 30 24
20 30 25
20 30 25
20 30 25
20 35 25
20 35 25
20 35 25
20 35 24
20 40 25
20 40 25
20 40 25
20 40 24
25 25 25
25 25 24
25 25 24
25 25 23
25 30 25
25 30 25
25 30 25
25 30 25
25 35 25
25 35 25
25 35 25
25 35 25
25 40 23
25 40 21
25 40 19
25 40 18
30 30 24
30 30 24
30 30 24
30 30 24
30 35 24
30 35 25
30 35 25
30 35 24
30 40 25
30 40 24
30 40 25
30 40 24
35 35 22
35 35 20
35 35 18
35 35 21
35 40 21
35 40 14
35 40 7
35 40 20
40 40 0
40 40 0
40 40 0
40 40 0
;
proc rsreg;
model y1=x1 x2 /lackfit predict press actual residual u95 l95 u95m l95m;
run;
现在，用以上方法已经得到拟和方程，但是在0.01水平的置信区间得不到（y1值）。
特此请教！谢谢！
会的请恢复联系我sjb263@sohu.com
| 
发表于 2006-2-18 19:15
|