[R-sig-phylo] Estimating discrete ancestral states with ARD model in ace

Hi, I am interested in estimating discrete ancestral states in ace using the calls: fitER<-ace(tipRegion,tree,model="ER",type="discrete") fitARD<-ace(tipRegion,tree,model="ARD",type="discrete", use.expm=TRUE) When I fit the ARD model (with either use.expm=TRUE or use.eigen=TRUE), I get the

Re: [R-sig-phylo] Estimating discrete ancestral states with ARD model in ace

Hi Katharine. You can also do ancestral state estimation for the model "ARD" using the packages diversitree. There is a worked example on my blog here: http://blog.phytools.org/2016/09/possible-bug-in-acemodelard-and-other.html (also highlighting a possible bug in ace for model="ARD" & other