Converts a tree given as parent vector format to an ancestor matrix

parentVector2anc(parents, n)

Arguments

parents

A integer vector. Parent vector format tree.

n

A single integer. The number of mutations.

Value

a logical matrix;0 and 1

Examples

p = c(11, 2, 3, 14, 14, 16, 8, 6, 9, 1, 15, 8, 10, 14, 5, 13, 17)
parentVector2anc(parents=p,n = 16)
#>       [,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,]    1    1    0    0    0    0    1    1    1     1     0     1     0
#>  [3,]    1    1    1    0    0    0    1    1    1     1     0     1     0
#>  [4,]    1    1    1    1    0    0    1    1    1     1     0     1     0
#>  [5,]    0    0    0    0    1    0    0    0    0     0     0     0     0
#>  [6,]    1    1    1    1    1    1    1    1    1     1     1     1     1
#>  [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,]    1    0    0    0    0    0    1    1    1     0     0     1     0
#> [10,]    1    0    0    0    0    0    1    1    1     1     0     1     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     1     0     1
#> [15,]    1    1    1    1    1    0    1    1    1     1     1     1     1
#> [16,]    0    0    0    0    0    0    0    0    0     0     1     0     1
#>       [,14] [,15] [,16]
#>  [1,]     0     0     0
#>  [2,]     0     0     0
#>  [3,]     0     0     0
#>  [4,]     0     0     0
#>  [5,]     0     0     0
#>  [6,]     1     1     1
#>  [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     1
#> [16,]     0     0     1