[R-sig-phylo] threshDIC value error: Inf for DIC Values

[Aja Carter]

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

BMThreshDIC <-threshDIC(scaled.tree, AT.data, ancestralBM, burnin=2000000)

**** 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-***@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-***@r-project.org/

[David Bapst]

Do you have any zero-length branches in your tree, Aja?

Cheers,
-Dave Bapst

Post by Aja Carter
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

BMThreshDIC <-threshDIC(scaled.tree, AT.data, ancestralBM, burnin=2000000)

**** 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]]
_______________________________________________
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo

--
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-***@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/r-sig-***@r-project.org/