[Rd] An extreme quantile of the geometric distribution

michael.baudin at contrib.scilab.org michael.baudin at contrib.scilab.org
Thu Jun 28 22:49:34 CEST 2012


 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 log1p function, instead of the naive formula.

 Am I correct ?

 Best regards,


More information about the R-devel mailing list