I am using the phylANOVA()
function in the R package phytools, as follows:
<code>phylANOVA(HoloceneTree, HolDat$Category,
HolDat$Niche_Breadth, nsim = 1000, posthoc = T)
<code>phylANOVA(HoloceneTree, HolDat$Category,
HolDat$Niche_Breadth, nsim = 1000, posthoc = T)
</code>
phylANOVA(HoloceneTree, HolDat$Category,
HolDat$Niche_Breadth, nsim = 1000, posthoc = T)
I keep getting the error message:
<code>Warning: no labels for x. Assuming order of tree$tip.label.
Warning: no labels for y. Assuming order of tree$tip.label.
<code>Warning: no labels for x. Assuming order of tree$tip.label.
Warning: no labels for y. Assuming order of tree$tip.label.
</code>
Warning: no labels for x. Assuming order of tree$tip.label.
Warning: no labels for y. Assuming order of tree$tip.label.
x and y are the categories and the continuous data, respectively. HolDat looks as follows:
<code>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")
<code>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")
</code>
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")
As you can see, I have the species names as row names. I have also tried to rename the first column, with species, as “Species”, “Species_Name”, and “Label”, and yet every time I get the same error that says that there are no labels for x and y. How do I make it so that the function can read the species names as the species labels?