phy_rich_mess <- as.data.frame(do.call(cbind, phy_richness)) phy_richness <- data.frame(SampleID = phy_rich_mess$ID_1, phy_rich_mess[,seq(0, length(phy_rich_mess), 2)]) ...