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/