Generator of subintensity matrix for the special case of segregating sites, ie. summary of all frequency counts
segsites(n, theta)
| n | sample size (positive integer) |
|---|---|
| theta | mutation parameter (positive) |
A `disc_phase_type` object containing subintensity matrix (P), vector of initial probabilities (alpha) and defect (probability of not entering any transient state prior to absorption)
segsites(n = 4, theta = 2)#> $subint_mat #> [,1] [,2] [,3] #> [1,] 0.4 0.3 0.2000000 #> [2,] 0.0 0.5 0.3333333 #> [3,] 0.0 0.0 0.6666667 #> #> $init_probs #> [,1] [,2] [,3] #> [1,] 1 0 0 #> #> $defect #> [1] 0 #> #> attr(,"class") #> [1] "disc_phase_type"