📦 Phylogenetic Monte Carlo: Quantifying Uncertainty in Phylogenetic Methods (now on CRAN)
Switch branches/tags
Nothing to show
Clone or download
Fetching latest commit…
Cannot retrieve the latest commit at this time.
Type Name Latest commit message Commit time
Failed to load latest commit information.
R minor edits for compatibility with new tidyr & dplyr Jun 10, 2016
data add the original anoles data from pmc, (bimac data with additional tr… Jun 23, 2015
inst use as.person compatible author string Aug 22, 2017
man re-knit readme May 16, 2018
vignettes vignette title set Jun 11, 2016
.Rbuildignore use new location for README figures Aug 23, 2017
.travis.yml vignettes build during build! May 16, 2018
CONTRIBUTING.md push basic skeleton Jun 11, 2015
DESCRIPTION drop keywords for CRAN May 17, 2018
NAMESPACE update docs Jun 10, 2016
NEWS re-knit readme May 16, 2018
README.md re-knit readme May 16, 2018
codemeta.json add codemeta.json, bump version Aug 22, 2017
cran-comments.md update comments May 17, 2018
pmc.Rproj avoid partial matching Sep 13, 2015


Travis-CI Build Status Coverage Status CRAN status DOI

This is a lightweight implementation of my pmc package focusing on what I think are the more common use cases (e.g. it will no longer support comparisons of a geiger model against an ouch model). Further, it does not cover many of the newer model fitting that have been implemented since pmc was first released.

The goal of this release is mostly to provide compatibility with current versions of geiger.

Getting started

Install the package:


A trivial example with data simulated from the lambda model.

#> Loading required package: ape
phy <- sim.bdtree(n=10)
dat <- sim.char(rescale(phy, "lambda", .5), 1)[,1,]
out <- pmc(phy, dat, "BM", "lambda", nboot = 50)

Plot the results:

dists <- data.frame(null = out$null, test = out$test)
#> Attaching package: 'dplyr'
#> The following object is masked from 'package:ggplot2':
#>     vars
#> The following objects are masked from 'package:stats':
#>     filter, lag
#> The following objects are masked from 'package:base':
#>     intersect, setdiff, setequal, union
dists %>% 
  gather(dist, value) %>%
  ggplot(aes(value, fill = dist)) + 
  geom_density(alpha = 0.5) + 
  geom_vline(xintercept = out$lr)


Carl Boettiger, Graham Coop, Peter Ralph (2012) Is your phylogeny informative? Measuring the power of comparative methods, Evolution 66 (7) 2240-51. http://doi.org/10.1111/j.1558-5646.2011.01574.x