cboettig/pmc

error in mixing ouch + geiger trees with ouch-formatted data

Closed this issue · 1 comments

Error filed originally by Danny Rojas.

The following example reproduces the error with version 0.0-5:

library(pmc)
library(ouch)
data(bimac)
tree <- with(bimac,ouchtree(node,ancestor,time/max(time),species))
eb_v_ou <- pmc(tree, log(bimac['size']), modelA="brown", modelB="lambda", nboot=5)

Problem actually arises because ouch handles data frames nicely for data, while passing to geiger confuses things. Added code to convert data frame to labeled numeric in commit f53237d