发布日期: 09/18/2023

置换检验的示例

在本例中,您使用“模拟”执行随机化或置换试验。您正在研究三种药物治疗疼痛的效果,并关注这些药物在效果上是否不同。由于您只有极小的样本并且有些担心违反通常的方差分析假设,所以您可以使用“模拟”来执行置换检验。

首先,您构造了一个公式,该公式在三种药物之间随机排列疼痛测量值。在“无效”这一原假设下,以上任意分配都与其他任何分配相似。由此判定:以这种方式下获得的 F 比近似原假设下的 F 比的分布。最后,您将 F 比的观测值与模拟获得的零分布进行比较。

定义模拟公式

1. 选择帮助 > 样本数据文件夹,然后打开 Analgesics.jmp

2. 选择列 > 新建列

3. 为列名键入随机排列的疼痛

4. 从“列属性”列表中,选择公式

5. 在函数列表中,选择行 > Col Stored Value

6. 在“列”列表中,双击疼痛

7. 点击编辑器面板上方符号列表中的插入键 (^)。

8. 从函数列表中,选择随机 > Col Shuffle

图 10.9 完成的公式 

Completed Formula

该公式随机排列疼痛列中的条目。

9. 在“公式编辑器”窗口中点击确定

10. 在“列信息”窗口中点击确定

执行置换检验

1. 选择分析 > 以 X 拟合 Y

2. 选择疼痛并点击 Y,响应

3. 选择药物并点击 X,因子

4. 点击确定

5. 点击“单因子分析”红色小三角菜单并选择均值/方差分析

图 10.10 “方差分析”报表 

Analysis of Variance Report

请注意,F 比为 6.2780。

6. 在“方差分析”分级显示项中,右击“F 比”列并选择模拟

7. 在“换出列”列表中,点击疼痛

8. 在“换入列”列表中,点击随机排列的疼痛

9. 样本数旁边,输入 1000。

10. (可选)在随机种子旁边,输入 456。

这将重现本例中的值。

图 10.11 完成的“模拟”窗口 

Completed Simulate Window

11. 点击确定

在模拟结果表中,“校正总和”和“误差”列均为空,因为“方差分析”表中的“F 比”值仅应用于药物

12. 在模拟值表中,运行分布脚本。

图 10.12 零分布下 F 比的模拟分布 

Simulated Distribution of F Ratios under the Null Distribution

观测到的 F 比值 6.2780 在直方图中用一条红线表示。该值落在 F 比的模拟零分布的上 0.5% 中。这表示:有很明显的证据表明这三种药物在缓解疼痛的效果上有差异。

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