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

From: peter dalgaard <pdalgd_at_gmail.com>
Date: Fri, 29 Jun 2012 14:21:50 +0200

On Jun 28, 2012, at 22:49 , <michael.baudin_at_contrib.scilab.org> <michael.baudin_at_contrib.scilab.org> 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 be a bug in the quantile function, which should use the log1p function, instead of the naive formula.
>
> Am I correct ?

Nope. (The source is availably, you know....).

The problem is that a slight fuzz is subtracted inside ceil(....), but there's no check that the result is positive.

qnbinom(...., size=1) is equivalent and does get right, by the way.

-pd

>
> Best regards,
>
> MichaŽl
>
> ______________________________________________
> R-devel_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel

-- 
Peter Dalgaard, Professor,
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Email: pd.mes_at_cbs.dk  Priv: PDalgd_at_gmail.com

______________________________________________
R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
Received on Fri 29 Jun 2012 - 12:25:25 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 30 Jun 2012 - 14:20:30 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive