[R-sig-phylo] Pagel's lambda for single trait

2015-10-22 Thread xinxin xue
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

Re: [R-sig-phylo] Pagel's lambda for single trait

2015-10-22 Thread Liam J. Revell
Hi Xinxin. In the current CRAN version of geiger treeTransform is not a valid argument of fitDiscrete, instead you probably want to do: fit<-fitDiscrete(tree,trait,transform="lambda") Then the fitted value of lambda will be: fit$opt$lambda All the best, Liam Liam J. Revell, Associate