发布日期: 11/15/2021

保存 X 矩阵

该选项将名为“矩量矩阵”和“模型矩阵”的脚本(它们分别包含矩量矩阵和模型矩阵)保存至设计数据表。矩量矩阵和模型矩阵用于计算“预测平均方差”,该方差显示在“设计评估”分级显示项的“设计诊断”部分中。请参见 Goos and Jones (2011)。若设计属于裂区设计,则还会保存“V 求逆”脚本。“V 求逆”脚本包含响应的协方差矩阵的逆协方。

警告:对于包含名义型因子的设计,使用“保存 X 矩阵”选项保存的“模型矩阵”不是拟合线性模型时使用的编码矩阵。通过在运行“模型”脚本时获取的“拟合模型”报表中选择选项“保存列”>“保存编码表”,您可以获取用于拟合模型的编码矩阵。

注意:您可以通过设置首选项始终保存矩阵脚本。选择文件 > 首选项 > 平台 > 实验设计。选中“保存 X 矩阵”。

模型矩阵

模型矩阵描述实验设计。在模型矩阵中,每次试验都对应一行,“模型”分级显示项中指定的每一个模型项都对应一列。对于每次试验,模型矩阵的相应行中都会包含模型项的编码值:

连续项按照从 -1 到 1 的范围来编码。

对在线性模型拟合中使用到的编码向量,通过应用 Gram-Schmidt 正交过程对名义项编码。

有关 Gram-Schmidt 正交过程的详细信息,请参见 Horn and Johnson (2012)。

注意:功效分析的编码与拟合线性模型的编码相匹配。有关用于拟合线性模型的名义型项的编码的信息,请参见《拟合线性模型》中的名义型效应的编码

矩量矩阵

矩量矩阵依赖模型效应但不依赖设计。其定义如下:

Equation shown here

其中,f(x) 表示与因子向量 x 的因子组合对应的模型效应,R 表示设计空间。有关矩量和设计矩阵的更多详细信息,请参见 Goos and Jones (2011) 和 Myers et al. (2009)。请注意,矩量矩阵在 Myers et al. (2009)。

脚本

从“定制设计”红色小三角菜单中,选择保存 X 矩阵。创建设计和表之后,在定制设计表中,“矩量矩阵”和“模型矩阵”脚本以及“V 求逆”脚本(若设计为裂区)都将另存为表脚本。

右击并从“矩量矩阵”、“模型矩阵”或“V 求逆”脚本选择编辑。脚本将显示对应的矩阵。您可以将该矩阵复制到您编写的脚本中。

运行“矩量矩阵”脚本时,日志中显示矩量矩阵中的行数(称为“矩”)。

运行“模型矩阵”脚本时,日志中显示模型矩阵中的行数(称为 X)。

运行“V 求逆”脚本时,日志中显示逆协方差矩阵中的行数(称为“V 求逆”)。

示例

本例演示模型矩阵脚本的用法:

提示:要查看日志,请选择视图 > 日志(在 macOS 上为窗口 > 日志)。

1. 选择实验设计 > 定制设计

2. 添加 3 个连续因子并点击继续

3. 点击交互作用 > 二次

4. 点击“定制设计”红色小三角并选择保存 X 矩阵

5. 使用默认试验次数 (12),依次点击制作设计制表

6. 在“表”面板中,右击“矩量矩阵”脚本并选择编辑

脚本显示在脚本窗口中。该脚本显示称为的矩量矩阵。

图 4.28 “矩量矩阵”脚本 

Image shown here

7. 若该脚本尚未打开,请选择视图 > 日志(在 macOS 上为窗口 > 日志)。

8. 在“表”面板中,点击“矩量矩阵”脚本旁边的绿色小三角。

日志中显示的行数为 N Row(Moments)=7

9. 在“表”面板中,右击“模型矩阵”脚本并选择编辑

脚本显示在脚本窗口中。该脚本显示称为 X 的模型矩阵。

10. 点击运行

日志中显示的行数为 N Row(X)=12。

11. 要以数据表形式查看模型矩阵,请向脚本中添加以下行:

dt = New Table( "Model Matrix" );
dt << Set Matrix( X );

12. 点击运行

需要更多信息?有问题?从 JMP 用户社区得到解答 (community.jmp.com).