Hi! Trees are not binary. Do you think this is the problem? Should I convert it?
Thanks! ________________________________ De: Klaus Schliep <[email protected]> Cc: "[email protected]" <[email protected]> Enviadas: Terça-feira, 8 de Abril de 2014 13:20 Assunto: Re: [R-sig-phylo] chronos and gammaStat I am guessing. Are the trees all binary, too? is.binary.tree(chr) On Tue, Apr 8, 2014 at 2:16 PM, Fabricia Nascimento <[email protected]> wrote: Dear Klaus, > > >all trees are rooted. >So I'm not sure why I am getting this warning. > > >Thanks! >Fabricia. > > > >________________________________ > De: Klaus Schliep <[email protected]> >Cc: "[email protected]" <[email protected]> >Enviadas: Terça-feira, 8 de Abril de 2014 13:13 >Assunto: Re: [R-sig-phylo] chronos and gammaStat > > >Hello Frabricia, > >I assume your RAxML tree is not rooted. >is.rooted(chr) >and use the function root if chr is not rooted. > >Regards, >Klaus > > >On Tue, Apr 8, 2014 at 2:03 PM, Fabricia Nascimento < > >> 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]/ >> > > > >-- >Klaus Schliep > > > [[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]/ > > > -- Klaus Schliep [[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]/
