optim.phylo.ls: Phylogeny inference using the least squares method
Description
This function performs phylogeny inference using least-squares.Usage
optim.phylo.ls(D, stree=NULL, set.neg.to.zero=TRUE, fixed=FALSE,
tol=1e-10, collapse=TRUE)Arguments
stree
an optional starting tree for the optimization.
set.neg.to.zero
a logical value indicating whether to set negative branch lengths to zero (default TRUE).
fixed
a logical value indicating whether to estimate the topology - if TRUE only the branch lengths will be computed.
tol
a tolerance value used to assess whether the optimization has converged.
collapse
a logical indicating whether to collapse branches with zero length.
Value
- a
"phylo" that (may be) the least-squares tree with branch lengths; also returns the sum of squares in attr(tree,"Q-score").
Details
Function uses nni from the "phangorn" package (Schliep 2011) to conduct NNIs for topology estimation.References
Cavalli-Sforza, L. L., and A. W. F. Edwards. 1967. Phylogenetic analysis: Modesl and estimation procedures. American Journal of Human Genetics, 19, 233-257.
Felsenstein, J. 2004. Inferring Phylogenies. Sinauer.
Paradis, E., J. Claude, and K. Strimmer. 2004. APE: Analyses of phylogenetics and evolution in R language. Bioinformatics, 20, 289--290.
Schliep, K. P. 2011. phangorn: phylogenetic analysis in R. Bioinformatics, 27, 592-593.