Yesterday a user reported that when he tried to reproduce my example analysis for the function evol.rate.mcmc() here, he encountered the following error when pre-processing the posterior sample:
> mcmc<-posterior.evolrate(tree,minsplit,result$mcmc[201:1001,], result$tips[201:1001])
Error in if (where == 0 || where == "root") where <- ROOTx :
missing value where TRUE/FALSE needed
Today I suggested he try updating {ape} - he did, and the problem went away.
I'm not sure what the compatibility issue is here, but I suspect it has something to do with one of drop.tip(), extract.clade(), or bind.tree(), as these are all called extensively in this function. Take home message - update {ape}!
To find out what version of {ape} (or any other package) you are working with, simply type:
> packageDescription("ape")["Version"]
$Version
[1] "2.7-1"
(FYI 2.7-1 is in fact the version I have installed on this machine.)
No comments:
Post a Comment
Note: due to the very large amount of spam, all comments are now automatically submitted for moderation.