Wilcoxon 符号秩检验当我在 JMP 中运行 Wilcoxon Signed-Rank 选项时,我得到的答案与 SAS 中得到的答案截然不同。 有没有办法获取每个选项(SAS 和 JMP)中使用的代码或用于找出结果不同的公式?
我在这里先向您的帮助表示感谢。回复:Wilcoxon符号秩检验如上所述,
大家好 - 我已在下面设置了一些 SAS 代码,按照前面参考文献中的描述进行分析,然后通过 JMP 运行相同的分析...
要在 JMP 中执行相同操作,请下载 James Abbey 在www.personal.psu.edu/jda188/JMP/Signedrank.doc上的以下说明,然后按照那里的说明对数据集的 DIFF 列进行操作。
也就是说,将 FERTILIZER1.SAS7BDAT 或 FERTILIZER3.SAS7BDAT 从您的 WILC 库导入 JMP,然后使用“分析 | 分布”分析 DIFF 列。将 DIFF 列复制到“Y,列”占位符中,然后点击“确定”。接下来,点击第二个红色三角形,选择“检验均值”。在“指定假设均值”中输入 0,并勾选“Wilcoxon 符号秩检验”复选框,然后点击“确定”。在输出结果的“检验均值=值”部分,您应该会发现检验统计量 = 74,P 值为 0.0017——这与 SAS 生成的答案相同。
可以吗?
回复:Wilcoxon符号秩检验很有可能其中一个只是设置不正确,所以最好的办法是先找一个教科书上的例子,然后用 SAS 和 JMP 分析这些数据,找出哪个方法没有产生预期的结果。我个人最喜欢的、通俗易懂的例子是 Sidney Siegel 的《行为科学的非参数统计》,但还有很多其他的例子。
我的第二个想法是,你执行的 SAS 或 JMP 检验可能实际上不是配对样本的符号秩检验,而是独立样本的秩和检验(有时也称为 Mann-Whitney U 检验)。这两个检验分别是配对 t 检验和非配对 t 检验的非参数类似物,并且都与 Wilcoxon 的名字有关:两者的详细信息以及示例都可以在维基百科中找到,链接如下:
http://en.wikipedia.org/wiki/Wilcoxon_signed-rank_test
http://en.wikipedia.org/wiki/Mann-Whitney-Wilcoxon_test回复:威尔科克森符号秩检验以下是来自加州大学洛杉矶分校统计学术技术服务中心的示例数据,该数据与上述观点相符:
回复:威尔科克森符号秩检验感谢二位的反馈。
David Q/RI 用的是书中的例子,所以我才对得到不同的答案感到惊讶。
有没有办法手动运行源代码来解决这个问题,并找出代码在哪一步可能出现不同的行为?
再次感谢。回复:威尔科克森符号秩检验如果有什么方法可以获取实际的源代码,恐怕我也不知道是什么——但我刚才快速浏览了一下MattF在上面帖子中提到的Wilcoxon检验,分别在SAS和JMP中查看了一下,它们看起来都像是对数据对之间的差异进行Mann-Whitney-Wilcoxon检验,这意味着你实际上需要先自己计算出数据对之间的差异,然后再对这些差异进行单因素比较。你实际是这样做的吗?
下面提供另一个参考资料,展示了如何在 SAS 中执行 Wilcoxon 配对符号秩检验,这表明基本检验实际上并不是以简单程序形式提供的(否则为什么要费这么大劲,用这种稍微绕弯的方式呢?):
http://udel.edu/~mcdonald/statsignedrank.html
如果你愿意贴出你正在运行的代码副本,我很乐意检查一下我的判断是否正确。