William Zhou 问答专帖

Tiaoti11

周老师，

您好，我来自SAE，曾当面听过您非常清晰精彩的Jmp功能介绍。在此请教一个问题：

Jmp的DOE中，2k因子设计+中心点试验结果的弯曲检验问题。
Jmp的“拟合模型”分析结果中只有“失拟”项的检验结果，没有弯曲项的检验结果，请问我们只需参考“失拟”吗？或者把中心点作为一个因子（两个水平：中心点或非中心点）进行回归分析？

在Minitab中是有单独的“弯曲”项检验结果的。
根据黑带指南，失拟主要是看高次项的误差与试验误差对比，而弯曲是看消除线性后的中心点的误差（2次项）与试验误差对比，失拟是不能完全表现弯曲的。
即使我们没有设计中心点，Jmp的正常分析结果也是有失拟选项的。
我搜索过论坛的帖子，有人建议增加2次项检验，但作为2k因子设计，这样做是没有意义的。

谢谢

zijuanchen

周老师，你好！我最近做DOE实验，总是遇见R2较高，失拟显著的情况？请问原因可能是什么？如何从数据中找到原因？该模型可否用于预测


最小二乘法拟合
响应“acid peak”


“预测值-实际值”图



拟合汇总

         
R 方        0.97691
调整 R 方        0.936503
均方根误差        0.741912
响应均值        22.59
观测数（或权重和）        23




方差分析

源        自由度        平方和        均方        F 比率
模型        14        186.30614        13.3076        24.1766
误差        8        4.40346        0.5504                概率>F
校正合计        22        190.70960                <.0001       




失拟

源        自由度        平方和        均方        F 比率
失拟        6        4.3813955        0.730233        66.1842
纯误差        2        0.0220667        0.011033                概率>F
总误差        8        4.4034622                0.0150       
                                        最大 R 方
                                0.9999       






排序参数估计值

项                 估计值        标准误差        t 比率                t 比率                概率>|t|
capacity(20,34)                3.0796296        0.225757        13.64                  ++++++++++        <.0001
load acid(23,28)                2.0964286        0.198284        10.57                  ++++++++          <.0001
pH(5.7,6.2)                0.5607407        0.225757        2.48                  ++                0.0379
concentration(242,302)                0.527963        0.225757        2.34                  ++                0.0475
capacity*capacity                -1.178333        0.67727        -1.74                 -                  0.1201
concentration*capacity                -2.415        1.679574        -1.44                 -                  0.1884
pH*capacity                -2.3675        1.679574        -1.41                 -                  0.1963
concentration*concentration                -0.693333        0.67727        -1.02                 -                  0.3359
pH*pH                -0.658333        0.67727        -0.97                 -                  0.3595
pH*concentration                -1.48        1.679574        -0.88                 -                  0.4039
load acid*capacity                0.22        0.262305        0.84                  +                 0.4260
pH*load acid                -0.2025        0.262305        -0.77                 -                  0.4623
concentration*load acid                -0.175        0.262305        -0.67                                    0.5234
load acid*load acid                -0.053333        0.67727        -0.08                                    0.9392

admin

Tiaoti11 发表于 2015-12-18 18:15
周老师，

您好，我来自SAE，曾当面听过您非常清晰精彩的Jmp功能介绍。在此请教一个问题：

你好，同意“失拟是不能完全表现弯曲的”，因为弯曲项只代表平方项、立方项等，失拟不仅包含平方项、立方项等，还包含一次项、交互作用项等等，所以看失拟检验是比较齐全的。
如果是2k因子设计，确实不用增加2次项检验，因为只有2个水平，加了也检验不出是否弯曲。

admin

zijuanchen 发表于 2016-2-16 13:37
周老师，你好！我最近做DOE实验，总是遇见R2较高，失拟显著的情况？请问原因可能是什么？如何从数据中找到 ...

你好，看到您另开了一贴，问题相似，还是建议您最好Case by case地来看，一般来讲，说明这个模型可以通过增加高次项来改进，但如果数据中因子水平数量少，想改进也改进不了。
这时最好能增加几个不同水平值的实验，然后再分析。如果其他模型诊断指标（包括残差）也还不错的话，做些简单的预测还是可以的。

zijuanchen

admin 发表于 2016-2-18 10:27
你好，看到您另开了一贴，问题相似，还是建议您最好Case by case地来看，一般来讲，说明这个模型可以通过 ...

数据如下，问题还是失拟显著，R2高
X1        X2        X3        X4        Y1        Y2
6.2        302        28        20        23.42        67.77
5.7        242        23        20               
6.2        302        23        20        19.47        70.98
6.2        242        23        20        17.54        74.11
5.7        302        23        20        17.38        71.6
6.2        242        28        20        20.71        71.23
5.7        302        28        20        20.8        69.89
5.95        272        25.5        20        18.66        71.51
5.7        242        28        20               
5.95        302        25.5        27        23.75        66.23
5.7        272        25.5        27        21.27        69.34
6.2        272        25.5        27        23.56        66.66
5.95        272        25.5        27        23.07        67.5
5.95        272        25.5        27        23.18        67.48
5.95        272        25.5        27        22.97        67.7
5.95        242        25.5        27        21.01        71.15
5.95        272        23        27        20.96        69.41
5.95        272        28        27        25.08        65.86
5.7        242        28        34        29.31        66.64
5.7        242        23        34        22.97        71.99
6.2        302        23        34               
5.95        272        25.5        34        25.13        66
5.7        302        28        34        26.8        63.22
5.7        302        23        34        22.64        66.31
6.2        302        28        34        27.79        62.13
6.2        242        23        34        22.85        67.44

admin

zijuanchen 发表于 2016-2-18 16:18
数据如下，问题还是失拟显著，R2高
X1        X2        X3        X4        Y1        Y2
6.2        302        28        20        23.42        67.77

你好，我仔细研究了一下，如果要彻底解决这个问题，还是要从数据源头看起，至少有以下3个疑点：
1 根据4因子Box-Behnken设计应该生成27次实验，但实际我只看到25次实验记录，最好补做一下
2 在已有的25次实验中，有3次是缺失实验结果的，最好也补充一下；
3 有三次重复实验（5.95,272,25.5,27），有2次误差非常小，建议再做2次，确认重复误差是否真的这么小。

zijuanchen

admin 发表于 2016-2-22 17:06
你好，我仔细研究了一下，如果要彻底解决这个问题，还是要从数据源头看起，至少有以下3个疑点：
1 根据4 ...

周老师，你好！确实是有27组实验，有两组由于收率为0，无纯度数据，所以这两组数据空着。重复的误差小是否会影响失拟值的计算，失拟是否就是残差/纯误差？
老师是否方便留下联系方式，想跟您再请教一下？谢谢！

zijuanchen

zijuanchen 发表于 2016-2-23 12:19
周老师，你好！确实是有27组实验，有两组由于收率为0，无纯度数据，所以这两组数据空着。重复的误差小是 ...

周老师，我将我这组数据中心点减少了一个，发现只有自由度从2降到1，原本失拟P值为0.01变为0.06，失拟显著变为不显著，我这样处理数据可以吗？请问这样的变化原因是什么？

shipl1106

周老师您好，请问哪里能找到关于使用JMP神经网络功能的参考资料

admin

shipl1106 发表于 2016-7-20 17:18
周老师您好，请问哪里能找到关于使用JMP神经网络功能的参考资料

你好，可以参见JMP帮助文档~