Re: [R-sig-phylo] Unexpected behavior in birth death simulations

2011-05-18 Thread Matthew Vavrek
Carl, was the other machine this was working on a Linux box? I just ran into this exact problem two days ago, and it seems to have something to do with the tolerance level in Linux, as it came up in a course and I was surrounded by non-Linux computers that seemed to have no problem with the

Re: [R-sig-phylo] Unexpected behavior in birth death simulations

2011-05-18 Thread Carl Boettiger
Hi Matt, all, The other machine is also linux, but 64 bit. I've also tested it on a third linux machine, which is running R-2.10 instead of R2.13, (but is 64 bit), and it works fine there too. Only on my 32 bit linux have I observed the error. I don't seem to have another 32 bit machine lying

Re: [R-sig-phylo] Unexpected behavior in birth death simulations

2011-05-18 Thread Emmanuel Paradis
09:08:23 To: Matthew Vavrekmatt...@matthewvavrek.com Cc: r-sig-phylo@r-project.org Subject: Re: [R-sig-phylo] Unexpected behavior in birth death simulations Hi Matt, all, The other machine is also linux, but 64 bit. I've also tested it on a third linux machine, which is running R-2.10 instead

Re: [R-sig-phylo] Unexpected behavior in birth death simulations

2011-05-18 Thread Carl Boettiger
...@matthewvavrek.com Cc: r-sig-phylo@r-project.org Subject: Re: [R-sig-phylo] Unexpected behavior in birth death simulations Hi Matt, all, The other machine is also linux, but 64 bit. I've also tested it on a third linux machine, which is running R-2.10 instead of R2.13, (but is 64 bit

Re: [R-sig-phylo] Unexpected behavior in birth death simulations

2011-05-18 Thread Carl Boettiger
Hi Tanja,list Yes, everything is going extinct, but as you say, that shouldn't happen when mu = 0. sim.bd.age(age=2,numbsim=3,lambda=2,mu=0,frac,mrca=FALSE,complete=FALSE) [[1]] [1] 0 [[2]] [1] 0 [[3]] [1] 0 It's a rather malicious error in that, as you point out, it's a possible outcome

Re: [R-sig-phylo] Unexpected behavior in birth death simulations

2011-05-18 Thread Brian O'Meara
As a general approach, when this sort of thing happens, I'll modify and load a copy of the function that has been changed to produce a lot of debugging info (lots of dput and print statements, for example). Then what happens at each simulation step to identify where the problem occurs. [Better

Re: [R-sig-phylo] Unexpected behavior in birth death simulations

2011-05-18 Thread Carl Boettiger
Hi list, Thanks all for the help and suggestions, think I've fixed the issue! Traced the behavior can to the locale used by R. You can replicate the error by using Sys.setlocale(locale=en_US.UTF-8) and fix the issue using Sys.setlocale(locale=C) It's not clear why having certain locales as