Я внедряю mantel.test() из пакета ape несколько раз, когда я перебираю различные матрицы, и диапазон значений p для всех проведенных мною тестов составляет от 0,001 до 1,999.Почему этот тест дает мне р-значения выше 1?Ниже приведен пример двух матриц, которые приводят к значению p выше 1:
library(ape)
m1 <- matrix(c(0, 12656.18,54236.03, 55983.36, 55983.36, 12656.18, 0,63570.61,
54487.93, 54487.93,54236.03, 63570.61, 0, 106487.00, 106487.00,
55983.36, 54487.93, 106487.00, 0, 0, 55983.36, 54487.93, 106487.00,
0, 0), nrow = 5, ncol = 5, byrow = TRUE)
m2 <- matrix(c(0, 0.05252174, 0.05252174, 1, 0, 0.05252174, 0, 0, 1, 0.052521740,
.05252174, 0, 0, 1, 0.05252174, 1, 1, 1, 0, 1, 0, 0.05252174,
0.05252174, 1, 0), nrow = 5, ncol = 5, byrow = TRUE)
set.seed(123)
mantel.test(m1, m2)
Вывод на моем конце выглядит следующим образом:
$z.stat
[1] 228926.3
$p
[1] 1.071
$alternative
[1] "two.sided"