Time until the most recent common ancestor

The function creates a continuous phase-type representation of the time until the most recent common ancestor for a specified number of sequences and following Kingsman's coalescent model.

t_mrca(n)

Arguments

n	positive integer larger than 1.

Examples

ph_mrca <- t_mrca(4)
summary(ph_mrca)
#> 
#> Subintensity matrix:
#>      [,1] [,2] [,3]
#> [1,]   -6    6    0
#> [2,]    0   -3    3
#> [3,]    0    0   -1
#> 
#> Initial probabilities:
#>      [,1] [,2] [,3]
#> [1,]    1    0    0
#> 
#> Defect:
#> [1] 0
#> 
#> Mean: 1.5
#> 
#> Variance: 1.138889
#>