秩变换的 2 路 ANOVA

机器算法验证 方差分析 连续数据 序数数据 排名
2022-04-09 18:02:39

我知道这些天我问了很多问题!对此感到抱歉,但我正在努力研究过去 4 年收集的研究生论文数据,并且反复绊倒我的初学者对统计数据的掌握。

背景:我的基本问题与本网站另一个问题 中描述的相同简而言之,我试图找出一种合理的方法来比较两个地点之间的好年份和差年份,以获得一堆变量(鸟类的繁殖指标)。我只对比较站点之间的好年份和站点之间的差年份感兴趣。我对一个站点的好年份是否比另一个站点的糟糕年份不感兴趣(所以我想这意味着我对交互不感兴趣?虽然我想知道的是一年内的站点是否-type 是 sig. 不同的,所以也许交互重要?)。尽管这些指标是相关的,因为它们都来自相同的 2 个总体,但我正在查看的变量是独立的。

这个问题是关于比较这 2 个地点和 2 年类型之间的筑巢日期(转换为序数日,因此是序数数据)。从直方图和箱线图一目了然,数据是倾斜的,并且方差是异质的。我的样本量范围(对于良好的站点内和较差的站点内)从 n = 69 到 864。我认为我应该使用非参数 2 路 ANOVA 并对其进行一些搜索。我已经看到 Friedmans 的测试作为 2 路 ANOVA 的非参数替代方案提供,但据我所知,它仅适用于重复测量数据。我发现一些人建议您可以对数据进行排名转换并在其上运行 2-way ANOVA。

为此,我只是通过使第一天(例如,154)= 1、第 155 天 = 2 等将每个站点和每个年份类型中的天数转换为排名,但这并没有改变我的任何东西数据,因为无论如何我都有序数日期。

那么,我能做些什么呢?我已经在等级和原始中心数据上运行了 2-way ANOVA。虽然每次测试的数字都不相同,但最终结果是:Holm-Sidak 多重比较表明,在每一年的好坏类型中,地点存在显着差异 (p<0.001)。我可以接受这些结果吗?或者有没有更好的方法来分析这些数据?

我目前正在使用 Sigmaplot 11(具有内置 Sigmastat 功能)来处理基本内容,并在它无法处理我需要做的事情时转向 R。

感谢您的阅读,以及您可以提供的任何提示!

编辑:我也对如何处理非正态(偏斜)且方差不等的连续数据感兴趣。等级转换也不适合吗?

1个回答

比例优势 (PO) 序数逻辑模型是 Wilcoxon 和 Kruskal-Wallis 检验的推广,允许协变量、交互作用以及您可以在回归模型中对单变量响应执行的任何其他操作。排名上的双向 ANOVA 不是基于强大的统计原理。

PO 模型的许多计算工具之一是lrmR rms包中的函数。rms 包的contrastanovasummaryplotnomogram函数可以提供帮助。