Design Experiment 文件夹下的 Donev Mixture Data.jmp 样本数据表中的数据基于 Atkinson 和 Donev (1992) 的示例。该设计包括三个混料因子和一个非混料因子。
• 响应是丙烯腈粉的电磁阻尼。
• 三个混料成分为硫酸铜 (CuSO4)、亚硫酸钠 (Na2S2O3) 和乙二醛。
• 关注的非混料环境因子为光的波长。
尽管波长在理论上是连续的,研究人员只关注三个离散波长处的预测值。因此,波长被视为具有三个水平的分类因子。
有关使用“定制设计”针对该情形构造设计的详细信息,请参见混料实验。
1. 选择帮助 > 样本数据库,然后打开 Design Experiment/Donev Mixture Data.jmp。
2. 点击列面板中硫酸铜 (CuSO4) 旁边的星号,然后选择混料。
图 A.20 硫酸铜 (CuSO4) 的“混料”列属性面板
请注意以下事项:
‒ “下限”为 0.2,即硫酸铜 (CuSO4) 的最小设计设置。
‒ “上限”为 0.8,即硫酸铜 (CuSO4) 的最大设计设置。
‒ “项总和”设置为 1。这是三个混料因子之和。
‒ 已选择 L 伪成分编码选项。请参见伪成分编码。
3. 点击取消。
4. 点击列面板中乙二醛旁边的星号,然后选择混料。
对于该因子,请注意以下事项:
‒ “下限”为 0,即乙二醛的最小设计设置。
‒ “上限”为 0.6,即乙二醛的最大设计设置。
5. 点击取消。
6. 在 Donev Mixture Data.jmp 数据表中,点击模型脚本旁边的绿色小三角。
模型包含混料因子的主效应和所有四个因子的双因子交互作用。
7. 点击运行。
在“参数估计值”报表中,混料因子以伪成分编码形式出现。混料因子出现在交互作用中时,它们不用编码形式表示。但是,基于伪成分拟合模型。“参数估计值”报表(图 A.21)中的前三项显示混料因子的编码形式。
图 A.21 “参数估计值”报表
8. 点击“响应‘阻尼’”红色小三角并选择估计值 > 显示预测表达式。
“预测表达式”报表显示拟合的模型。请注意混料因子使用 L 伪成分编码进行了变换。
图 A.22 “阻尼”模型的预测表达式
假定您关注波长 L2 处的预测值。还假定亚硫酸钠 (Na2S2O3) 和乙二醛分别设置为其低值 0.2 和 0,硫酸铜 (CuSO4) 设置为其高值 0.8。在这种情况下,预测的阻尼等于硫酸铜 (CuSO4) 的参数估计值 (6.191) 加上硫酸铜 (CuSO4)*波长[L2] 的参数估计值 (1.878)。您可以在“预测刻画器”中对此进行验证。
9. 点击“响应‘阻尼’”红色小三角并选择保存列 > 保存编码表。
图 A.23 显示出编码混料因子的编码表的前三列
对于这个特定设计,“L 伪成分编码”将混料因子变换到 0-1 的范围。请注意这通常不会发生。