Hi Dave,
The seed is created every time you generate random data (see details in
?.Random.seed):
R> exists(".Random.seed")
[1] FALSE
R> sample(1)
[1] 1
R> exists(".Random.seed")
[1] TRUE
So in your code below the seed is created by calling rtree(). Indeed,
with ape 4.1:
R>
Hi all,
I've been investigating the use of Pagel's lambda and run into a problem
comparing trees that have a low speciation rate and no extinction with
trees that have a high speciation rate and lots of extinction. I'm
comparing phylogenetic signal in two characters. In the first, I simulate
Hi William.
Why would you suspect phytools? ;)
I didn't figure this out for you, but I did eliminate some possibilities.
It does not seem to be due to phytools::fastBM. You can replace fastBM
with phytools:sim.corrs or geiger:sim.char as follows:
Thanks, Emmanuel. I was not familiar with .Random.seed objects until R
CHECK started complaining about such elements in saved workspaces
within .rdata files within packages (requiring rm to remove them), and
then was mainly familiar with them regarding set.seed() (and so I
thought ace must be