Hi!
I am using the gammaStat in ape to analyze some trees I generated using
simulations and RaxML.
After using the R code below to convert it to ultrametric tree:
cal <- makeChronosCalib(Tree, age.min=1)
chr <- chronos(Tree, lambda = 0.1, calibration = cal)
gammaStat(chr)
I get the following result:
> gammaStat(chr)
[1] 2.639323
Warning message:
In (2:N) * g :
longer object length is not a multiple of shorter object length
Can I ignore this warning?
Thanks very much!
Fabricia.
[[alternative HTML version deleted]]
_______________________________________________
R-sig-phylo mailing list - [email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/[email protected]/
