Re: [Rd] An extreme quantile of the geometric distribution

2012-06-30 Thread peter dalgaard
On Jun 30, 2012, at 12:13 , wrote: > Hi, > > I'm sorry, I do not clearly understand. > > I'm aware that the source is available at : > > http://svn.r-project.org/R/trunk/src/nmath/qgeom.c Yes, but you were suggesting that non-use of log1p caused the problem, and the source clearly uses it

Re: [Rd] An extreme quantile of the geometric distribution

2012-06-30 Thread michael.baudin
Hi, I'm sorry, I do not clearly understand. I'm aware that the source is available at : http://svn.r-project.org/R/trunk/src/nmath/qgeom.c But a good source does not mean a correct result, because of compilation issues. Moreover, I do not fully understand why the 1e-7 coefficient in the form

Re: [Rd] An extreme quantile of the geometric distribution

2012-06-29 Thread peter dalgaard
On Jun 28, 2012, at 22:49 , wrote: > Hi, > > With R version 2.10.0 (2009-10-26) on Windows, I computed the p=1.e-20 > quantile of the geometric distribution with parameter prob=0.1. > >> qgeom(1.e-20,0.1) > [1] -1 > > But this is not possible, since X=0,1,2,... > > I guess that this might

[Rd] An extreme quantile of the geometric distribution

2012-06-29 Thread michael.baudin
Hi, With R version 2.10.0 (2009-10-26) on Windows, I computed the p=1.e-20 quantile of the geometric distribution with parameter prob=0.1. qgeom(1.e-20,0.1) [1] -1 But this is not possible, since X=0,1,2,... I guess that this might be a bug in the quantile function, which should use the l