R vegan包中的部分Mantel测试给出了错误

Question

我尝试使用R vegan包进行部分Mantel测试。即使是非常简单的伪数据，我在输出中也会出现以下错误。

我的代码：

  m1 = matrix(c(0, 1, 2,
                1, 0, 3,
                2, 3, 0),nrow = 3)
  m2 = m1 * 2
  m3 = m1 * 3
  mtest = mantel.partial(xdis = m1, ydis = m2,zdis = m3)
  print(mtest)

输出：

Partial Mantel statistic based on Pearson product-moment correlation 

Call:
mantel.partial(xdis = m1, ydis = m2, zdis = m3) 

Mantel statistic r:   NaN 
      Significance: NA 

Error in quantile.default(x$perm, c(0.9, 0.95, 0.975, 0.99)) : 
  missing values and NaN's not allowed if 'na.rm' is FALSE

Answer 1

示例矩阵与线性相关; m2是m1的简单线性变换，因此它们完全相关。与m3一样。一旦考虑了一对矩阵之间的相关性，就没有任何东西可以与之相关联，因此奇怪的NaN和NA以及最终的错误。