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]/
        [[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]/

Reply via email to