响应曲面失拟问题，附上数据和DOE仿真模拟软件，佳能炮

在做试验设计的时候遇到以下问题：
全因子试验：
标准序 运行序 中心点 区组 角度 quantity ball type powder type Fuse type 距离
32 1 0 1 35 7.5 2 2 2 256.2
40 2 0 1 35 7.5 2 2 2 262.5
38 3 0 1 35 7.5 2 1 2 230.8
14 4 1 1 50 5.0 2 2 1 130.7
15 5 1 1 20 10.0 2 2 1 255.4
30 6 0 1 35 7.5 2 1 2 237.3
19 7 0 1 35 7.5 1 2 1 267.3
31 8 0 1 35 7.5 1 2 2 271.9
39 9 0 1 35 7.5 1 2 2 271.0
29 10 0 1 35 7.5 1 1 2 238.8
33 11 0 1 35 7.5 1 1 1 238.3
28 12 0 1 35 7.5 2 2 1 253.9
7 13 1 1 20 10.0 2 1 2 240.5
9 14 1 1 20 5.0 1 2 1 198.0
24 15 0 1 35 7.5 2 2 2 256.7
18 16 0 1 35 7.5 2 1 1 231.0
36 17 0 1 35 7.5 2 2 1 259.5
5 18 1 1 20 5.0 2 1 1 190.0
8 19 1 1 50 10.0 2 1 1 201.2
20 20 0 1 35 7.5 2 2 1 269.3
13 21 1 1 20 5.0 2 2 2 200.3
4 22 1 1 50 10.0 1 1 2 202.7
16 23 1 1 50 10.0 2 2 2 221.3
26 24 0 1 35 7.5 2 1 1 230.8
23 25 0 1 35 7.5 1 2 2 269.6
3 26 1 1 20 10.0 1 1 1 247.1
27 27 0 1 35 7.5 1 2 1 272.0
17 28 0 1 35 7.5 1 1 1 238.5
35 29 0 1 35 7.5 1 2 1 268.5
1 30 1 1 20 5.0 1 1 2 194.1
22 31 0 1 35 7.5 2 1 2 223.8
12 32 1 1 50 10.0 1 2 1 227.1
11 33 1 1 20 10.0 1 2 2 257.1
37 34 0 1 35 7.5 1 1 2 239.0
34 35 0 1 35 7.5 2 1 1 224.3
2 36 1 1 50 5.0 1 1 1 119.8
10 37 1 1 50 5.0 1 2 2 134.9
6 38 1 1 50 5.0 2 1 2 115.0
25 39 0 1 35 7.5 1 1 1 238.5
21 40 0 1 35 7.5 1 1 2 238.9

首先全因子试验删选出主效应 系数标
项 效应 系数 准误 T P
常量 195.95 1.692 115.78 0.000
角度 -53.73 -26.86 1.692 -15.87 0.000
quantity 71.20 35.60 1.692 21.04 0.000
ball type -7.13 -3.56 1.070 -3.33 0.003
powder type 24.14 12.07 1.070 11.28 0.000
Fuse type 0.06 0.03 1.070 0.03 0.978
角度*quantity 16.77 8.39 1.692 4.96 0.000
角度*ball type -0.77 -0.39 1.692 -0.23 0.821
角度*powder type 4.52 2.26 1.692 1.34 0.194
角度*Fuse type -0.80 -0.40 1.692 -0.24 0.815
quantity*ball type -0.60 -0.30 1.692 -0.18 0.861
quantity*powder type 3.05 1.53 1.692 0.90 0.377
quantity*Fuse type -1.88 -0.94 1.692 -0.55 0.585
ball type*powder type -0.03 -0.01 1.070 -0.01 0.989
ball type*Fuse type -0.23 -0.12 1.070 -0.11 0.915
powder type*Fuse type -0.08 -0.04 1.070 -0.04 0.971
Ct Pt 53.57 2.185 24.52 0.000


S = 6.76946 PRESS = 9942.80
R-Sq = 98.45% R-Sq（预测） = 85.38% R-Sq（调整） = 97.37%


距离 的方差分析（已编码单位）

来源 自由度 Seq SS Adj SS Adj MS F P
主效应 5 38159.1 38159.1 7631.8 166.54 0.000
角度 1 11545.5 11545.5 11545.5 251.94 0.000
quantity 1 20277.8 20277.8 20277.8 442.50 0.000
ball type 1 508.4 508.4 508.4 11.09 0.003
powder type 1 5827.4 5827.4 5827.4 127.16 0.000
Fuse type 1 0.0 0.0 0.0 0.00 0.978
2因子交互作用 10 1265.8 1265.8 126.6 2.76 0.021
角度*quantity 1 1125.6 1125.6 1125.6 24.56 0.000
角度*ball type 1 2.4 2.4 2.4 0.05 0.821
角度*powder type 1 81.9 81.9 81.9 1.79 0.194
角度*Fuse type 1 2.6 2.6 2.6 0.06 0.815
quantity*ball type 1 1.4 1.4 1.4 0.03 0.861
quantity*powder type 1 37.2 37.2 37.2 0.81 0.377
quantity*Fuse type 1 14.1 14.1 14.1 0.31 0.585
ball type*powder type 1 0.0 0.0 0.0 0.00 0.989
ball type*Fuse type 1 0.5 0.5 0.5 0.01 0.915
powder type*Fuse type 1 0.1 0.1 0.1 0.00 0.971
弯曲 1 27546.1 27546.1 27546.1 601.11 0.000
残差误差 23 1054.0 1054.0 45.8
失拟 7 773.1 773.1 110.4 6.29 0.001
纯误差 16 280.9 280.9 17.6
合计 39 68025.0

quantity的交互作用响应曲面试验：
标准序 运行序 点类型 区组 角度 quantity 距离
8 1 -1 1 35 10.0 275.0
1 2 1 1 20 5.0 200.8
10 3 0 1 35 7.5 268.8
7 4 -1 1 35 5.0 191.9
12 5 0 1 35 7.5 268.5
3 6 1 1 20 10.0 258.2
4 7 1 1 50 10.0 231.2
6 8 -1 1 56 7.5 186.0
2 9 1 1 50 5.0 133.7
9 10 0 1 35 7.5 270.3
5 11 -1 1 13 7.5 223.3
11 12 0 1 35 7.5 270.7
13 13 0 1 35 7.5 271.9

中心复合设计

因子: 2 仿行: 1
基础次数: 13 总试验数: 13
基础区组: 1 合计区组数: 1

两水平因子：全因子

立方点: 4
立方体的中心点: 5
轴点: 4
轴点的中心点: 0

Alpha: 1.41421


响应曲面回归:距离 与 角度, quantity

分析是使用已编码单位进行的。

距离 的估计回归系数

系数标
项 系数 准误 T P
常量 270.04 4.592 58.801 0.000
角度 -18.36 3.631 -5.056 0.001
quantity 34.05 3.631 9.379 0.000
角度*角度 -35.96 3.893 -9.237 0.000
quantity*quantity -21.56 3.893 -5.539 0.001
角度*quantity 10.03 5.134 1.952 0.092


S = 10.2689 PRESS = 5205.44
R-Sq = 96.94% R-Sq（预测） = 78.42% R-Sq（调整） = 94.75%


距离 的方差分析

来源 自由度 Seq SS Adj SS Adj MS F P
回归 5 23386.5 23386.5 4677.31 44.36 0.000
线性 2 11972.3 11972.3 5986.14 56.77 0.000
角度 1 2695.6 2695.6 2695.62 25.56 0.001
quantity 1 9276.7 9276.7 9276.66 87.97 0.000
平方 2 11012.3 11012.3 5506.13 52.22 0.000
角度*角度 1 7777.5 8997.5 8997.50 85.32 0.000
quantity*quantity 1 3234.8 3234.8 3234.75 30.68 0.001
交互作用 1 402.0 402.0 402.00 3.81 0.092
角度*quantity 1 402.0 402.0 402.00 3.81 0.092
残差误差 7 738.2 738.2 105.45
失拟 3 730.3 730.3 243.43 123.69 0.000
纯误差 4 7.9 7.9 1.97
合计 12 24124.7


R-Sq = 98.45% 和 R-Sq（调整） = 97.37% 值都已经在比较好的水平，显著项也已经识别，但是一直失拟，软件中也没有其他因子作用了，请教大家为什么？如何分析这种情况谢谢！







用于仿真模拟的软件也分享给大家。