考虑 Uniformity Trial.jmp 样本数据表。对 8x8 的田块网格执行农事均匀性试验。在均匀性试验中,在没有进行实验处理的一块田地上种植检验作物。测量响应变量(通常为产量)。目标是找出田地变异性特征,以便为在该田地上进行设计实验提供参考。(请参见 Littell 等,2006,pp. 447。)
1.
选择帮助 > 样本数据库,然后打开 Uniformity Trial.jmp
2.
选择分析 > 拟合模型
3.
选择保持对话框打开,以便您可以在下一示例中返回到启动窗口。
4.
选择产量并点击 Y
5.
6.
选择重复结构选项卡。
9.
选择,然后点击重复
显示“重复结构”选项卡的完成的“拟合模型”启动窗口
10.
点击运行
“拟合混合模型”报表 - 空间球体
“拟合混合模型”报表显示在“拟合混合模型”报表 - 空间球体中。 “预测值-实际值”图显示预测产量为单个值。这是因为只拟合了空间协方差。“拟合统计量”报表显示“-2 对数似然”为 227.68,“AICc”为 234.08。
因为拟合了各向同性空间结构,将显示“变差图”。因为是在 8x8 网格设计这些试验的,距离近的地方将比距离很远的地方有更多对点。请参见建议的完整和不完整区组设计的图形生成器图 以了解布局。“变差图”显示球形空间结构对于最大约为 8.4 的距离是极好的拟合。最终距离的距离分类只包括两对对角线点。
拟合独立误差模型
2.
选择重复结构选项卡。
5.
点击运行
您可以在本例中执行似然比检验,因为独立误差模型嵌套在空间模型系列中。独立误差模型是空间相关性 ρ 等于 0 的空间模型。这意味着您可以对这两个模型执行正式的似然比检验。
在本例中,似然比检验统计量为 254.22-227.68=26.54。将它与一个自由度的卡方分布进行比较,当 p 值 < 0.0001 时拒绝无空间相关性这个原假设。您可以得出这些数据包含显著空间变异性的结论。
2.
选择重复结构选项卡。
3.
选择,然后点击重复
4.
从“结构”列表中选择带 Nugget 的空间
6.
点击运行
“拟合混合模型”报表显示在“拟合混合模型”报表 - 带 Nugget 的空间球体中。请注意,该对数似然基本上等于不带 nugget 的球形模型的对数似然值,AICc 略高一点(236.36 对比 234.08)。“重复效应协方差参数估计值”报表显示 Nugget 协方差参数的估计值为零。这说明不存在任何 nugget 效应。
“拟合混合模型”报表 - 带 Nugget 的空间球体
“拟合混合模型”报表 - 变差图
9.
选择重复结构选项卡。
12.
点击运行
“拟合混合模型”报表显示在“拟合混合模型”报表 - 空间各向异性球形中。 拟合统计量指示该模型的拟合效果没有各向同性(空间结构)球形模型的拟合效果好(AICc 240.54 对比 234.08)。“重复效应协方差参数估计值”报表显示“行”(空间球体行)和“列”(空间球体列)协方差的估计值很接近。没有证据表明网格的行内和列内的空间相关性存在差异。
“拟合混合模型”报表 - 空间各向异性球形
确定空间结构的类型
2.
选择重复结构选项卡。
3.
选择,然后点击重复
6.
点击运行
9.
点击运行
12.
点击运行
空间模型拟合的拟合统计量汇总了这些类型以及您执行的其他拟合的观测 AICc 值。
将模型与区组设计比较
建议的完整和不完整区组设计的图形生成器图
2.
选择重复结构选项卡。
4.
从效应中删除
5.
选择随机效应选项卡。
6.
选择季度,然后点击添加
7.
点击运行
2.
选择随机效应选项卡。
3.
选择季度并点击删除
4.
选择子季度并点击添加
5.
点击运行