I am investigating patterns of phenotypic overdispersion/clustering in a
reginal scale. I am using a set of traits to calculate phenotypic
diversity. I need to simulate this set of correlated traits to generate a
null model.
I would like to know whether the following approach is appropriate:
- I used fitContinuous::Geiger to verify which model of evolution best fits
to my traits.The models used were White Noise, Ornstein Uhlenbeck, Brownian
Motion, Early Bust and lambda.
- I compared the models according to its AICc values. Most of my traits
fitted better to lambda model with values next to 1 (0.9).
Can I transform the original tree based on lambda values for each trait
separately (using transform.phylo) and then to use a function to simulate
the evolution of each song trait?
Which is the appropriate function to simulate ( "RtraitCont", "sim.corrs"
or fastBM)?
Thank you,
Viviane
[[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]/
