Getting started

library("devtools")install_github("cboettig/pmc")

library("pmc")library("geiger")#> Loading required package: ape#> Loading required package: phytools#> Loading required package: mapslibrary("phytools")phy<- sim.bdtree(n=10)dat<- sim.char(rescale(phy,"lambda",.5),1)[,1,]out<- pmc(phy,dat,"BM","lambda",nboot=50)#> Warning in geiger::fitContinuous(phy = tree, dat = data, model = model, :#> Parameter estimates appear at bounds:#>  lambda

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

Citation