View source: R/f_phyclust_Pt.r
phyclust.Pt | R Documentation |
This computes transition probabilities of phyclust
given time.
phyclust.Pt(Q, Tt, substitution.model = .substitution.model$model[1])
Q |
a list according to the substitution model. |
Tt |
total evolution time, |
substitution.model |
substitution model. |
The major models for Q
are:
Model | Author and Publication | Parameter |
JC69 | Jukes and Cantor 1969. | t |
K80 | Kimura 1980. | \kappa, t |
F81 | Felsenstein 1981. | \pi, t |
HKY85 | Hasegawa, Kishino, and Yano 1985. | \pi, \kappa, t |
A list of Q
should contains pi
, kappa
based on
substitution models and code types. Tt
may be separately stored.
Depending on identifiers, Q
s can be composite to a QA
,
Q matrix array.
A list with class Pt
will be returned containing several
elements as the following:
'Pt' |
a transition probability matrix. |
'log.Pt' |
a log transition probability matrix. |
'H' |
a negative entropy, |
vectorize Tt
for repeated computation in C.
Wei-Chen Chen wccsnow@gmail.com
Phylogenetic Clustering Website: https://snoweye.github.io/phyclust/
.substitution.model
,
phyclust
,
phyclust.em.step
.
## Not run:
library(phyclust, quiet = TRUE)
Tt <- 0.5
Q <- list(pi = c(0.25, 0.25, 0.25, 0.25), kappa = 0.5)
phyclust.Pt(Q, Tt, "HKY85")
Q <- list(pi = c(0.5, 0.5), kappa = 0.5)
phyclust.Pt(Q, Tt, "SNP_JC69")
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.