R素食包中的部分Mantel測試給出了錯誤

Question

我正在嘗試使用R素食主義者包運行部分壁爐架測試。 即使在非常簡單的偽數據上，我在輸出中也會出現以下錯誤。

我的代碼：

  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以及最終誤差。