Convert the newick format tree to an ancestor matrix
nwk2anc(nwk)the string of newick format tree.
p = c(11, 2, 3, 14, 14, 16, 8, 6, 9, 1, 15, 8, 10, 14, 5, 13, 17)
nwk <- parentVector2nwk(parents=p,n = 16)
nwk2anc(nwk)
#> [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
#> [1,] 1 0 0 0 0 0 0 0 0 0 0 0 0
#> [2,] 0 1 0 0 0 0 0 0 0 1 0 0 0
#> [3,] 0 1 1 0 0 0 0 0 0 0 0 0 0
#> [4,] 0 0 1 1 0 0 0 0 0 0 0 0 0
#> [5,] 0 0 0 0 1 0 0 0 0 0 0 0 0
#> [6,] 0 0 0 0 0 1 0 0 0 0 0 0 0
#> [7,] 0 0 0 0 0 0 1 1 0 0 0 0 0
#> [8,] 0 0 0 0 0 0 0 1 0 0 0 0 0
#> [9,] 0 0 0 0 0 0 1 0 1 0 0 1 0
#> [10,] 0 0 0 0 0 0 0 0 1 1 0 0 0
#> [11,] 0 0 0 0 0 0 0 0 0 0 1 0 1
#> [12,] 1 0 0 0 0 0 0 0 0 0 0 1 0
#> [13,] 0 0 0 0 0 0 0 0 0 0 0 0 1
#> [14,] 0 0 0 0 0 0 0 0 0 0 0 0 0
#> [15,] 0 0 0 1 1 0 0 0 0 0 0 0 0
#> [16,] 0 0 0 0 0 0 0 0 0 0 1 0 0
#> [,14] [,15] [,16]
#> [1,] 0 0 0
#> [2,] 0 0 0
#> [3,] 0 0 0
#> [4,] 0 0 0
#> [5,] 0 0 0
#> [6,] 0 1 0
#> [7,] 0 0 0
#> [8,] 0 0 0
#> [9,] 0 0 0
#> [10,] 0 0 0
#> [11,] 0 0 0
#> [12,] 0 0 0
#> [13,] 0 0 0
#> [14,] 1 0 1
#> [15,] 1 1 0
#> [16,] 0 0 1