Mantel检验报错‘Error in m1 * m2 : non-numeric argument to binary operator’的解决求助

这个报错的核心原因很明确：你的mat1和mat2矩阵是字符型（chr），而mantel.test()需要数值型矩阵来执行计算（函数内部会对两个矩阵做元素级乘法等操作，字符无法参与这类运算）。下面是具体的排查和修复步骤：

1. 正确读取数据，分离个体标识与数值数据 #

你提到矩阵包含个体标识（"Ind1"等），推测你在读取CSV时，把个体标识当成了矩阵的一列/一行数据，而非行名/列名。修正读取代码，将个体标识设为行名：

# 读取CSV时指定第一列为行名（假设个体标识在第一列）
mat1 <- read.csv("GroupSinfantprox1.csv", sep=",", header=T, row.names=1)
mat2 <- read.csv("GroupSjuvenileprox1.csv", sep=",", header=T, row.names=1)

2. 转换为数值矩阵并验证类型 #

转换后务必检查矩阵的数据类型，确保是数值型：

mat1 <- as.matrix(mat1)
mat2 <- as.matrix(mat2)

# 查看矩阵结构，确认是numeric类型
str(mat1)
str(mat2)

如果输出显示chr类型，说明你的CSV文件中存在非数值的单元格（比如把NA写成了"na"、空字符串或其他文本），可以在读取时指定识别缺失值：

mat1 <- read.csv("GroupSinfantprox1.csv", sep=",", header=T, row.names=1, na.strings=c("NA", "na", ""))

3. 验证矩阵对称性（可选但推荐） #

Mantel检验通常要求输入的是对称矩阵（如距离矩阵、关联矩阵），可以快速验证：

# 检查矩阵是否对称
all(mat1 == t(mat1))
all(mat2 == t(mat2))

如果返回FALSE，需要检查你的原始数据是否符合对称矩阵的格式要求。

4. 重新运行Mantel检验 #

完成上述步骤后，再执行检验代码：

mantel.test(mat1, mat2, nperm=999, graph=FALSE, alternative="two.sided")

内容的提问来源于stack exchange，提问作者Sae In Lee