警告:对于包含名义型因子的设计,使用“保存 X 矩阵”选项保存的“模型矩阵”脚本不是拟合线性模型时使用的编码矩阵。通过在运行“模型”脚本时获取的“拟合模型”报表中选择选项“保存列”>“保存编码表”,您可以获取用于拟合模型的编码矩阵。
模型矩阵描述实验设计。在模型矩阵中,每次试验都对应一行,“模型”分级显示项中指定的每一个模型项都对应一列。对于每次试验,模型矩阵的相应行中都会包含模型项的编码值:
由于名义型项经过编码用于构造最优设计的方式不同,所以若设计包含名义型因子,模型矩阵编码不同于拟合线性模型时使用的编码。有关用于拟合线性模型中的名义型项的编码的信息,请参见《拟合线性模型》手册中的名义型效应的编码。
矩量矩阵依赖模型效应但不依赖设计。其定义如下:
其中, 表示与因子向量 的因子组合对应的模型效应,R 表示设计空间。有关矩量和设计矩阵的更多详细信息,请参见 Goos and Jones (2011, pp 88-90) 和 Myers et al. (2009)。请注意,矩量矩阵在 Myers et al. (2009, p. 376) 中称为区域矩量矩阵。
从“定制设计”红色小三角菜单中,选择保存 X 矩阵。创建设计和表之后,在定制设计表中,“矩量矩阵”和“模型矩阵”脚本以及“V 求逆”脚本(若设计为裂区)都将另存为表属性。
•
|
从“矩量矩阵”、“模型矩阵”或“V 求逆”脚本旁边的红色小三角菜单中选择编辑。脚本将显示对应的矩阵。您可以将该矩阵复制到您编写的脚本中。
|
1.
|
2.
|
3.
|
点击交互作用 > 二阶。
|
4.
|
从“定制设计”红色小三角菜单中,选择保存 X 矩阵。
|
5.
|
6.
|
脚本显示在脚本窗口中。该脚本显示称为矩的矩量矩阵。
7.
|
日志中显示的行数为 N Row(Moments)=7。
9.
|
脚本显示在脚本窗口中。该脚本显示称为 X 的模型矩阵。
10.
|
点击运行。
|
日志中显示的行数为 N Row(X)=12。
dt = New Table( "Model Matrix" );
dt << Set Matrix( X );
12.
|
点击运行。
|