对于每个对比,t 比是用对比除以 PSE 计算得到的。原假设下这些 t 比的参考分布在计算上很难处理,因此是通过模拟得到的。使用的方法(如下所述)基于 Ye and Hamada (2000) 中的论述。
用 ti 表示第 i 个对比的 t 比:
筛选实验中最重要的是个体误差率,即声明给定效应为活跃的但是其实它不活跃的概率。对于第 i 个效应,当 |ti| 很大、落入其参考分布的上尾时会出现这种情况。
因为平台构造相对较大数量的效应,实验级误差率也很重要。实验级误差率是当没有效应是活跃的时候声明任何效应是活跃的概率。当没有效应是活跃的且绝对 t 比的最大值 max|ti| 很大并落入其参考分布的上尾时,出现实验级误差。
“拟合两水平筛选”平台使用 Monte Carlo 模拟获取两种误差率类型的参考分布。考虑包含 n - 1 个值的一个集合,它们采用均值为 0 且标准差等于 PSE 的正态分布来模拟。这些值表示在无活跃效应的原假设下实验的可能对比值。一共生成了 10,000 组 n - 1 随机对比值。
考虑第 i 个对比。使用这 10,000*(n - 1) 个模拟值中的每一个模拟值构造 Lenth t 比。用这些 t 比的绝对值近似表示个体误差率的参考分布。个体 p 值是在以降序排列的 10,000*(n - 1) 个模拟的绝对 t 比中观测到的绝对 Lenth t 比的小数位置内插值。对于参考分布,这近似于观测的绝对 Lenth t 比的绝对值右侧区域。
若当所有效应都不活跃时任意 t 比导致拒绝原假设,会出现实验级误差。即若绝对t 比的最大值 max|ti| 导致拒绝原假设,会出现实验级误差。
要获取这种情况下的参考分布,在 10,000 个模拟中考虑每个模拟的绝对 t 比的最大值。这 10,000 个最大值构成了参考分布。联合 p 值是在以降序排列的 10,000 个模拟的最大绝对 t 比中观测到的绝对 Lenth t 比的小数位置内插值。对于基于模拟的最大绝对 t 比的参考分布,这近似于绝对 Lenth t 比的绝对值右侧区域。
要更改集合中包含的默认模拟个数 (10,000),您必须为名为 LenthSimN 的全局 JSL 变量赋值。用于生成 p 值的模拟数在报表窗口中报告。
下例将模拟数设置为 50,000。
1. 选择帮助 > 样本数据库,然后打开 Reactor Half Fraction.jmp。
2. 选择实验设计 > 经典 > 两水平筛选 > 拟合两水平筛选。
3. 选择反应百分比并点击 Y。
4. 从进料速度一直选择到浓度,点击 X,然后点击确定。
5. 点击“反应百分比的筛选”旁边的红色小三角菜单,然后选择保存脚本 > 至脚本窗口。
6. 插入 LenthSimN=50000; 作为脚本的第一行。
LenthSimN=50000;
Fit Two Level Screening(
Y( :Percent Reacted),
X( :Feed Rate, :Catalyst, :Stir Rate, :Temperature, :Concentration )
);
7. 在脚本窗口中右击并选择运行脚本。
图 10.10 使用用于生成 p 值的 N = 50,000 次模拟来拟合两水平
模拟数在报表窗口中说明。
注意:若 LenthSimN=0,则使用标准 t 分布且不提供联合 p 值。不推荐使用该方法。