Total tree length

This function creates a continuous phase-type representation of the total tree length for a specified number of sequences and following Kingsman's coalescent model.

t_total(n)

Arguments

n	positive integer larger than 1.

Examples

ph_total <- t_total(4)
summary(ph_total)
#> 
#> Subintensity matrix:
#>      [,1] [,2] [,3]
#> [1,] -1.5  1.5  0.0
#> [2,]  0.0 -1.0  1.0
#> [3,]  0.0  0.0 -0.5
#> 
#> Initial probabilities:
#>      [,1] [,2] [,3]
#> [1,]    1    0    0
#> 
#> Defect:
#> [1] 0
#> 
#> Mean: 3.666667
#> 
#> Variance: 5.444444
#>