I have just submitted a new version of phytools (phytools 0.5-00) to CRAN. Hopefully, it survives the CRAN maintainer scrutiny & posts within the next 24 hours; however it can already be installed from GitHub using devtools:
##  '0.5.0'
This version of phytools has many changes and new functionality when compared to the previous CRAN version from early July, so it will be hard to be fully comprehensive about the updates; however, among the changes in this CRAN version, users will find:
A bug fix in
plot.lttwith a plotted tree overlain that caused the tip labels to be spaced exponentially.
An important bug fix in the phytools function to test for differences in the evolutionary correlation between traits in different parts of the tree, based on Revell & Collar (2009; Evolution).
A very nice update/improvement to the phytols function (
cladelabels) for adding clade labels to a plotted tree (here).
A change in the way object class is checked in phytools that affected the source code of many different phytools function (here).
Some changes to
cladelabelsto permit them to be used more effectively with the function
An update to the function
treeSliceto slice the tree at a particular point & return all subtrees (that, among other things, allows this to be done interactively).
A fix to the default color scheme used by
phenogramwith a mapped discrete character.
An update to the phytools function
phylomorphospacethat allows it to play with ape functions such as
nodelabels(and thus permits us to plot labels at internal or external nodes of a phylogeny projected into 2D morphospace.
A new S3 plotting method for objects of class
"rateshift"created by the function to estimate the temporal position of a shift in the evolutionary rate for a continuous character through time.
A small bug fix in the phytools function
A new S3
An important update to the internal functions used by
make.simmapto compute the likelihood, as well as the conditional likelihoods at internal nodes. This primarily affects users who are using a non-default prior at the root node of the tree.
A new function,
fitMk, to fit the Mk model. (Functionally similar to
A more or less imperceptible update to the function
rerootingMethod. This basically allows the function to handle polytomies “natively” (i.e., without first having to resolve them, perform calculations using the resolved tree, and then match nodes back to the original tree).
An additional estimation method (and a change in the default estimation method) for the phytools function
fitPagel, which fits the Pagel (1994) model for testing whether the evolution of two binary characters is correlated.
An update to the function
rateshiftfor testing for temporal shifts in the rate of evolution for a continuous trait that permits much more effective (if not faster) optimization. (Prior versions really could not maximize the likelihood effectively for more than - at best - two shifts on reasonably sized trees.)
Finally, some modifications to the phytools function
brownieREML, and a REML version of the function
rateshift, which runs much faster than the ML version (here).
Feedback is always welcome on any issues created by these changes. (Of course, I'm always happy to take positive feedback too!)