This function bootstraps a star tree from a model fitted by phyclust
.
Each cluster corresponds to a star tree and a center sequence where
sequences will evolve from. This function is called by
bootstrap.star.trees.seq
to generate sequences.
The fitted model can be varied in .identifier
.
bootstrap.star.trees(ret.phyclust, min.n.class = 1)
a phyclust object in Class phyclust
.
minimum number of sequences for a cluster.
Return a list containing \(K\) star trees according to ret.phyclust
.
ret.phyclust
is a phyclust object in Class phyclust
which is usually
fitted by phyclust
, or returned by phyclust.m.step
.
min.n.class
is a boundary condition to avoid degenerate clusters
when some population proportions, \(\eta_k\), are small in the
fitted model.
Phylogenetic Clustering Website: https://snoweye.github.io/phyclust/
# NOT RUN { library(phyclust, quiet = TRUE) set.seed(1234) EMC.1 <- .EMC EMC.1$EM.iter <- 1 # the same as EMC.1 <- .EMControl(EM.iter = 1) ret.1 <- phyclust(seq.data.toy$org, 2, EMC = EMC.1) ret.trees <- bootstrap.star.trees(ret.1) # }