Hi,
how do i get the actual lambda value using fitDiscrete? it's for a binary
trait over a phylogeney. And would i need a p value too, for loglikelihood
ratio test between this lambda model and another model? thanks!
>lambda<-fitDiscrete(tree,trait,treeTransform="lambda")
>lambda
GEIGER-fitted comparative model of discrete data
fitted Q matrix:
1 2
1 -42.0543 42.0543
2 42.0543 -42.0543
model summary:
log-likelihood = -53.756148
AIC = 109.512296
AICc = 109.551900
free parameters = 1
Convergence diagnostics:
optimization iterations = 100
failed iterations = 0
frequency of best fit = 1.00
object summary:
'lik' -- likelihood function
'bnd' -- bounds for likelihood search
'res' -- optimization iteration summary
'opt' -- maximum likelihood parameter estimates
[[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/
