构造方差分量的精确置信区间
构造设计
若您想要跳过本节中的步骤,请选择帮助 > 样本数据库,然后打开 Design Experiment/
Catalyst Design.jmp
。在 Catalyst Design.jmp 数据表中,点击实验设计模拟脚本旁边的绿色小三角。然后转到拟合模型
1.
选择实验设计 > 定制设计
2.
在“因子”分级显示项中的添加因子数旁边键入 3。
3.
点击添加因子 > 连续
4.
通过双击将这些因子重命名为温度、时间催化剂
5.
对于温度,点击容易并选择极难
这样即可将温度定义为一个整区因子。
6.
对于时间,点击容易并为时间选择困难
这样即可将时间定义为一个子区因子。
7.
点击继续
8.
注意:设置步骤 10中的随机种子和步骤 11中的开始数可重现本例中显示的确切结果。自行构造设计时,这些步骤不是必需的。
10.
(可选)点击“定制设计”红色小三角并选择设置随机种子。键入 12345 并点击确定
11.
(可选)点击“定制设计”红色小三角并选择开始数。键入 1000 并点击确定
12.
点击制作设计
13.
点击制表
注意:YY 模拟列将不同于设计表中所示的列。
设计表
“模拟响应”窗口
设计表和“模拟响应”窗口随即显示。请注意,该设计表包含一个实验设计模拟脚本。您可以随时运行该脚本以指定不同的参数值。
拟合模型
1.
在“分布”面板中(“模拟响应”窗口)的整区 σ 旁边键入 2。
2.
子区 σ 旁边键入 1.5。
3.
点击应用
在数据表中,Y 模拟的公式随之更新以反映您指定的内容。要查看该公式,点击“列”面板中列名右侧的加号。
4.
在数据表中,点击模型脚本旁边的绿色小三角。
5.
点击 Y 按钮旁边的 Y 变量,然后点击删除
6.
点击 Y 模拟,然后点击 Y 按钮。
7.
点击运行
注意:由于 Y 模拟中的值是随机生成的,您报表中的条目将不同于显示 Wald 置信区间的 REML 报表中显示的那些条目。
显示 Wald 置信区间的 REML 报表
探索功效
“模拟”窗口
在您的模拟中,用新的 Y 模拟列(用于为每次模拟生成新的模拟值列)替换用于运行模型的列 Y 模拟。将为“参数估计值”表中列出的每个效应模拟您右击并显示为已选定的方差分量列。
2.
样本数旁边输入 200
3.
(可选)在随机种子旁边,输入 456。
这将重现方差分量模拟结果表(部分视图)中显示的值,只有第 1 行中的值除外。
4.
点击确定
方差分量模拟结果表(部分视图)
“拟合最小二乘法模拟结果(方差分量)”数据表的第一行包含方差分量的初始值,该行已被排除。其余行包含模拟值。
5.
运行分布脚本。
方差分量的分布图(部分视图)