| phyclust.se {phyclust} | R Documentation |
The phyclust.se is an advanced function of phyclust.
The phyclust.se implements finite mixture models for sequence data
with sequencing errors. The same as phyclust,
the mutation processes are modeled by evolution processes based on
Continuous Time Markov Chain theory.
phyclust.se(X, K, EMC = .EMC, manual.id = NULL, byrow = TRUE)
X |
nid/sid matrix with N rows/sequences and L columns/sites. |
K |
number of clusters. |
EMC |
EM control. |
manual.id |
manually input class ids. |
byrow |
advanced option for |
phyclust.se considers mutations with sequencing error, but only
for NUCLEOTIDE data and only the EM algorithm is
implemented. Currently, phyclust.se implements the following steps:
| 1 | assume non-sequencing errors as the phyclust does. |
| 2 | use the initialization as the phyclust does. |
| 3 | run the phyclust to find the non-sequencing error MLE's. |
| 4 | initial by the results of phyclust. |
| 5 | update all parameters including probabilities of sequencing errors. |
See the help page of phyclust for the explanations of
input arguments since both functions share the same arguments. Note that
the underling model assumptions are very different of both functions.
A list with class phyclust will be returned containing
several elements as described in phyclust. But, the
phyclust.se returns extra parameters for sequencing errors, and
they are shown in the following:
'SE' |
a list returning parameters of sequencing error models, including:
|
Wei-Chen Chen wccsnow@gmail.com
Phylogenetic Clustering Website: http://snoweye.github.io/phyclust/
.EMC,
.EMControl,
phyclust.se,
phyclust.se.update.
## Not run: library(phyclust, quiet = TRUE) X <- seq.data.toy$org set.seed(1234) (ret.1 <- phyclust.se(X, 3)) ## End(Not run)