有没有办法通过排列测试比较线性回归斜率?

机器算法验证 回归 多重比较 置换检验
2022-03-31 14:39:34

我有一些物种×年矩阵并通过回归年份之间的距离来分析时间动态time lag. 由于数据点不是独立的,因为我得到(n2n)/2时间序列长度的距离n,使用标准方法来确定回归斜率的显着性是有问题的。因此,我应用了蒙特卡洛置换程序,并将随机斜率的分布与观察到的斜率进行了比较。到目前为止没有问题;我的大部分回归斜率在p<0.0001. 我的问题:如何比较不同物种的坡度×年矩阵?由于上述问题,参数化方法仍然不是解决方案,但我不知道基于置换过程的解决方案。

1个回答

原则上,您可以对两组数据的任何函数使用置换检验。要比较回归斜率,您只需选择物种,随机打乱组之间的数据点,但不是比较每个排列组的平均值,而是比较回归斜率。您可能必须在计算回归量之前将数据居中(可以是线性的或自定义构建的,只要它可以在数据居中后由单个参数描述)。我不知道如何处理置换测试中的多重比较——对于一些类别,Bonferroni 校正可以解决问题,但我想有几十个物种会杀死任何影响......关于这个的想法?