使用phytools包phylANOVA()时物种标签匹配警告的解决问询

问题描述 #

使用R语言phytools包的phylANOVA()函数执行以下代码时：

phylANOVA(HoloceneTree, HolDat$Category, 
          HolDat$Niche_Breadth, nsim = 1000, posthoc = T)

反复收到如下警告：

Warning: no labels for x. Assuming order of tree$tip.label.

Warning: no labels for y. Assuming order of tree$tip.label.

数据框HolDat已将物种名设为行名，也曾尝试将物种列重命名为"Species"、"Species_Name"、"Label"，但警告仍未消除，其结构如下：

structure(list(label = c("Aegotheles_novazelandiae", "Alauda_arvensis", 
"Anas_chlorotis", "Anas_platyrhynchos", "Anomalopteryx_didiformis", 
"Chen_caerulescens"), Category = c("Flightless", "True Fliers", 
"True Fliers", "True Fliers", "Flightless", "True Fliers"), biovar = c(1L, 
1L, 1L, 1L, 1L, 1L), Niche_Breadth = c(5.259427983, 3.207528743, 4.208010193, 
29.52295142, 6.985928971, 32.54111922)), row.names = c("Aegotheles_novazelandiae", 
"Alauda_arvensis", "Anas_chlorotis", "Anas_platyrhynchos", "Anomalopteryx_didiformis", 
"Chen_caerulescens"), class = "data.frame")

解决方案 #

核心问题是：直接传入的HolDat$Category和HolDat$Niche_Breadth是无名称的向量，phylANOVA()无法将它们与进化树的末端标签对应，只能默认按顺序匹配。以下两种方法可解决：

方法1：为向量添加物种名称标签 #

将分类变量和连续变量转为带物种名（数据框行名）的命名向量：

# 提取带名称的分类变量
x <- HolDat$Category
names(x) <- rownames(HolDat)

# 提取带名称的连续变量
y <- HolDat$Niche_Breadth
names(y) <- rownames(HolDat)

# 运行phylANOVA
phylANOVA(HoloceneTree, x, y, nsim = 1000, posthoc = TRUE)

方法2：直接传入数据框并指定列名 #

phylANOVA()支持直接传入数据框，通过x和y参数指定对应列名，函数会自动读取数据框的行名作为物种标签：

phylANOVA(HoloceneTree, x = "Category", y = "Niche_Breadth", data = HolDat, 
          nsim = 1000, posthoc = TRUE)

关键验证步骤 #

• 确保数据框的行名（或向量名称）与HoloceneTree$tip.label的拼写、大小写、符号完全一致
• 确认进化树的所有末端物种都存在于数据框中，无缺失或多余物种

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