Wrapper function that computes the matrices describing Kingsman coalescent ready for reward transformation.

kingsman(n)

Arguments

n

sample size (a positive integer)

Value

A mult_phase_type object containing the subintensity matrix, reward matrix and vector of initial probabilities

Examples

kingsman(4)
#> $subint_mat #> [,1] [,2] [,3] [,4] #> [1,] -6 6 0 0 #> [2,] 0 -3 2 1 #> [3,] 0 0 -1 0 #> [4,] 0 0 0 -1 #> #> $RewardM #> [,1] [,2] [,3] #> [1,] 4 0 0 #> [2,] 2 1 0 #> [3,] 1 0 1 #> [4,] 0 2 0 #> #> $init_probs #> [1] 1 0 0 0 #> #> attr(,"class") #> [1] "mult_phase_type"