对于每个对比,t 比是用对比除以 PSE 计算得到的。原假设下这些 t 比的参考分布在计算上很难处理,因此是通过模拟得到的。下面所述的方法基于 Ye and Hamada (2000) 中的论述。
ti 表示第 i 个对比的 t 比:
筛选实验中最重要的是个体误差率,即声明给定效应为活跃的但是其实它不活跃的概率。对于第 i 个效应,当 |ti| 很大、落入参考分布的上尾时会出现这种情况。
因为平台构造相对较大数量的效应,实验级误差率也很重要。实验级误差率是当没有效应是活跃的时候声明任何效应是活跃的概率。当没有效应是活跃的且绝对 t 比的最大值 max|ti| 很大并落入参考分布的上尾时,出现实验级误差。
“拟合两水平筛选”平台使用 Monte Carlo 模拟获取两种误差率类型的参考分布。考虑包含 n - 1 个值的一个集合,它们采用均值为 0 且标准差等于 PSE 的正态分布来模拟。这些值表示在无活跃效应的原假设下实验的可能对比值。一共生成了 10,000 组对比值。
考虑第 i 个对比。使用每个模拟值构造 Lenth t 比。用这些 t 比的绝对值近似表示个体误差率的参考分布。在报表的个体 p 值列中给定的 p 值是在以降序排列的 10,000*(n - 1) 个模拟的绝对 t 比中观测到的绝对 Lenth t 比的小数位置内插值。对于参考分布,这近似于观测的绝对 Lenth t 比的绝对值右侧区域。
若当所有效应都不活跃时任意 t 比导致拒绝原假设,会出现实验级误差。即若绝对t 比的最大值 max|ti| 导致拒绝原假设,会出现实验级误差。
要获取这种情况下的参考分布,在 10,000 个模拟中计算每个模拟的绝对 t 比的最大值。这 10,000 个最大值构成了参考分布。在报表的联合 p 值列中给定的 p 值是在以降序排列的 10,000 个模拟的最大绝对 t 比中观测到的绝对 Lenth t 比的小数位置内插值。对于基于模拟的最大绝对 t 比的参考分布,这近似于绝对 Lenth t 比的绝对值右侧区域。
要更改集合中包含的默认模拟个数 (10,000),您必须为名为 LenthSimN 的全局 JSL 变量赋值。例如,执行以下操作:
1.
选择帮助 > 样本数据库,然后打开 Reactor Half Fraction.jmp
2.
选择分析 > 专业建模 > 专业实验设计模型 > 拟合两水平筛选
3.
反应百分比选作响应变量 Y
5.
点击确定
8.
突出显示 LenthSimN=50000; 和其余代码。
注意:LenthSimN=0,则使用标准 t 分布且不提供联合 p 值(不推荐)。