使用evolvability包计算演化度均值报错：list无法转为double类型

我尝试用R语言的evolvability包计算14个物种协方差矩阵的EvolvabilityMeans。先准备了包含17个颅部距离数据的协方差矩阵列表和均值列表，为测试演化度对矩阵做了均值调整（所有矩阵均为正定矩阵），但调用evolvabilityMeans()时报错：

'list' object cannot be coerced to type 'double'

代码示例 #

cov.matrices <- list(matriz.cov.1, matriz.cov.2, matriz.cov.3,
                     matriz.cov.4, matriz.cov.5, matriz.cov.6, matriz.cov.7,
                     matriz.cov.8, matriz.cov.9, matriz.cov.10, matriz.cov.11, matriz.cov.12,
                     matriz.cov.13, matriz.cov.14)

means.list<- list(Mean_1, Mean_2, Mean_3, Mean_4, Mean_5, Mean_6, Mean_7, Mean_8, Mean_9, Mean_10,  Mean_11,  Mean_12, Mean_13, Mean_14)

mean.std.matrices = function (cov.matrix, means){
  return(cov.matrix/(means%*%t(means)))
}

mean.std.cov.matrices <- list()
for(i in 1:length(cov.matrices)){
  mean.std.cov.matrices[[i]] <- mean.std.matrices(cov.matrices[[i]], means.list[[i]])
}
names(mean.std.cov.matrices) <- names(cov.matrices)

mean.std.cov.matrices

evolvabilityMeans(mean.std.cov.matrices, means = means.list)

问题原因与解决方案 #

报错的核心原因是evolvabilityMeans()的means参数不接受列表格式，该参数要求传入数值向量（单矩阵场景）或矩阵（多矩阵场景，每行对应一个协方差矩阵的均值向量）。

修正步骤：

1. 将均值列表转换为矩阵，每行对应一个物种的均值向量：

# 把列表中的均值向量绑定为矩阵
means_matrix <- do.call(rbind, means.list)

2. 使用转换后的矩阵调用evolvabilityMeans()：

evolvabilityMeans(mean.std.cov.matrices, means = means_matrix)

额外检查

确认均值标准化函数的逻辑是否符合你的分析需求：means%*%t(means)生成的是外积矩阵，协方差矩阵与该矩阵的逐元素除法是否是你需要的标准化方式。若需按变量维度单独标准化，可能需要调整为cov.matrix / (means %o% means)（与外积等价）或其他标准化逻辑，但当前报错与该步骤无关。

内容的提问来源于stack exchange，提问作者Carolline Raidan