求助一个关于DOE实例的问题

本帖最后由 hell55 于 2011-8-14 17:39 编辑

RT，我求出的结果与书中的不符，例子是部分析因实验的例子，响应值是3个，Mean,BetWfStd,WithUnif,做2的5-1次试验，共16次，其中前两项的显著项分别是A,D,E；A，B，C，D，AC，BD；最后一项我的结果与书中不符，书中结果是B，D，BD；我的结果是A,B,D,BD。
请大家帮忙看一下数据，我分析的哪里有毛病？
数据如下
中心点 区组 BT TF O2 PB SP Mean BetWfStd WithUnif
1 1 30 5.15 3.12 46.6 1 0.84080 0.075141 0.027930
1 1 30 5.15 3.72 40.6 1 0.88710 0.007119 0.069623
1 1 30 6.60 3.12 40.6 1 0.81227 0.075197 0.018262
1 1 30 6.60 3.72 46.6 1 0.85647 0.064485 0.039342
1 1 40 5.15 3.12 40.6 1 0.94543 0.006124 0.084862
1 1 40 5.15 3.72 46.6 1 1.11563 0.095045 0.040396
1 1 40 6.60 3.12 46.6 1 1.20910 0.063213 0.051896
1 1 40 6.60 3.72 40.6 1 1.01517 0.113340 0.031026
1 1 30 5.15 3.12 40.6 2 0.84990 0.029511 0.071349
1 1 30 5.15 3.72 46.6 2 0.92790 0.040018 0.031014
1 1 30 6.60 3.12 46.6 2 0.99281 0.041149 0.053680
1 1 30 6.60 3.72 40.6 2 0.77937 0.001044 0.023782
1 1 40 5.15 3.12 46.6 2 1.25300 0.075592 0.033034
1 1 40 5.15 3.72 40.6 2 1.05070 0.063943 0.112929
1 1 40 6.60 3.12 40.6 2 1.15447 0.083123 0.027263
1 1 40 6.60 3.72 46.6 2 1.25223 0.059693 0.047204

#########################################################################

我的结果是
拟合因子: WithUnif 与 BT, TF, PB

WithUnif 的估计效应和系数（已编码单位）

项 效应 系数 系数标准误 T P
常量 0.04772 0.002324 20.54 0.000
BT 0.01170 0.00585 0.002324 2.52 0.029
TF -0.02234 -0.01117 0.002324 -4.81 0.001
PB -0.01433 -0.00716 0.002324 -3.08 0.010
TF*PB 0.03727 0.01864 0.002324 8.02 0.000


S = 0.00929427 PRESS = 0.00201038
R-Sq = 90.37% R-Sq（预测） = 79.63% R-Sq（调整） = 86.87%


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

来源 自由度 Seq SS Adj SS Adj MS F P
主效应 3 0.0033642 0.0033642 0.00112139 12.98 0.001
BT 1 0.0005479 0.0005479 0.00054789 6.34 0.029
TF 1 0.0019955 0.0019955 0.00199545 23.10 0.001
PB 1 0.0008208 0.0008208 0.00082082 9.50 0.010
2因子交互作用 1 0.0055569 0.0055569 0.00555688 64.33 0.000
TF*PB 1 0.0055569 0.0055569 0.00555688 64.33 0.000
残差误差 11 0.0009502 0.0009502 0.00008638
失拟 3 0.0003869 0.0003869 0.00012896 1.83 0.220
纯误差 8 0.0005633 0.0005633 0.00007042
合计 15 0.0098713


WithUnif 的异常观测值

标准化
观测值 标准序 WithUnif 拟合值 拟合值标准误 残差 残差
14 8 0.112929 0.090542 0.005196 0.022387 2.90R

R 表示此观测值含有大的标准化残差


WithUnif 的估计系数（使用未编码单位的数据）

项 系数
常量 2.39613
BT 0.00117035
TF -0.388983
PB -0.0527265
TF*PB 0.00856833


WithUnif 的效应图


WithUnif 的效应 Pareto


别名结构
I
BT
TF
PB
TF*PB