testlist <- list(ances = c(322601176L, -132831231L, 2107072297L, 1069463867L, NA, 1885310409L, 1476627388L, 1692541047L, 497045525L, -1233643725L, 259415119L, 873827053L, 1558959309L, 1247909539L, 1679908808L, 75448353L, 529544917L, NA, -1016829961L, 2045896150L, -591813849L, 439995006L, 1529431647L, -1467580217L, 381066527L, 77936583L, 178175875L, -1885054613L, 1576833655L, 912061841L, -2059669459L, -1032784424L, 647165738L, 925823553L, -860764837L, 1428080562L )) result <- do.call(phylobase:::hasPolytomy,testlist) str(result)