Do you have any zero-length branches in your tree, Aja?
Cheers,
-Dave Bapst
On Mon, Dec 18, 2017 at 3:22 PM, Aja Carter wrote:
> I am using the ancThresh model in the phytools package for 172 taxa. For
> model we have the liabilities, parameters, MCMC and alpha values output.
> When I tried getting DIC values from the threshDIC function I got Infs for
> both the OU model and the Lambda model.The error output I have been
> receiving is posted below:
>
>> BMThreshDIC <-threshDIC(scaled.tree, AT.data, ancestralBM, burnin=200)
> NOTE: no sequence provided, column order in x
>> OUThreshDIC
> Dbar Dhat pD DIC
> 2109.678 Inf -Inf -Inf
>
> Has anyone experienced this error before? Any advice on how to fix this
> error would be greatly appreciated.
>
>
> --
> Aja Carter
> PhD Candidate
>
> Department of Earth and Environmental Sciences
>
> University of Pennsylvania
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-phylo mailing list - R-sig-phylo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
> Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
--
David W. Bapst, PhD
Postdoc, Ecology & Evolutionary Biology, Univ of Tenn Knoxville
Lecturer, Geology & Geophysics, Texas A & M University
https://github.com/dwbapst/paleotree
Google Calender: https://goo.gl/EpiM4J
___
R-sig-phylo mailing list - R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-phylo@r-project.org/
