Rthreshml: R interface for threshml
Description
This function is an R interface for threshml in the PHYLIP package (Felsenstein 1989, 2013). threshml fits the threshold model of Felsenstein (2005; 2012). Note that threshml is new & not in the currently released version of PHYLIP (as of December 2013). It can be downloaded from its webpage here: http://evolution.gs.washington.edu/phylip/download/threshml/. If not specifying path, the executable file for threshml (e.g., threshml.exe in Windows) should be placed in the folder containing all other executable files for PHYLIP (e.g., C:/Program Files/phylip=3.695/exe in Windows).
Obviously, use of any of the functions of this package requires that PHYLIP (Felsenstein 2013) should first be installed. More information about installing PHYLIP can be found on the PHYLIP webpage: http://evolution.genetics.washington.edu/phylip.html.Usage
Rthreshml(tree, X, types=NULL, path=NULL, ...)
Value
This function returns a list containing the results from threshml.Details
Optional arguments include the following: quiet suppress some output to R console (defaults to quiet = FALSE); burnin burnin generations for the MCMC; nchain number of chains of the MCMC; ngen number of generations in each chain; proposal variance on the proposal distribution for the MCMC; lrtest logical value indicating whether to conduct a likelihood-ratio test of the hypothesis that some correlations are zero (does not appear to work in the present version); and cleanup remove PHYLIP input/output files after the analysis is completed (defaults to cleanup = TRUE).References
Felsenstein, J. (1989) PHYLIP--Phylogeny Inference Package (Version 3.2). Cladistics, 5, 164-166.
Felsenstein, J. (2005) Using the quantitative genetic threshold model for inferences between and within species Philosophical Transactions of the Royal Society London B, 360, 1427-1434.
Felsenstein, J. (2012) A comparative method for both discrete and continuous characters using the threshold model. American Naturalist, 179, 145-156.
Felsenstein, J. (2013) PHYLIP (Phylogeny Inference Package) version 3.695. Distributed by the author. Department of Genome Sciences, University of Washington, Seattle.