R devel - Jul 2005 - qgamma with first argument 1

> qgamma(1,1.1)
   [1] Inf

as expected.  But for shape parameter between 0.16 and 1:

   > qgamma(1,.5)
   [1] NaN
   Warning message:
   NaNs produced in: qgamma(p, shape, scale, lower.tail, log.p)

and for shape parameter 0.16 or less:

   > qgamma(1,.1)
   [1] 99.42507

Arguments close to 1 give approximately correct results:

   > qgamma(1-1e-15,.5)
   [1] 32.05055
   > qgamma(1-1e-15,.1)
   [1] 28.8129





R version is the R 2.1.1 binary from CRAN:
   platform i386-pc-mingw32
   arch     i386
   os       mingw32
   system   i386, mingw32
   status
   major    2
   minor    1.1
   year     2005
   month    06
   day      20
   language R

Thank you, Mr Hosking,

yes, this is a buglet.

I had started to fix this (and similar cases) 
--  using a C Macro based approach in src/nmath/dpq.h  --

For this reason, the fix will probably only appear in R-devel
i.e., from R-2.2.0 on, and not yet in [R 2.1.1]-patched.

Regards,
Martin Maechler, ETH Zurich
>>>>> "J" == J Hosking <jh910 at juno.com>
>>>>>     on Sat, 09 Jul 2005 15:27:00 -0400 writes:
    >> qgamma(1,1.1)
    J> [1] Inf

    J> as expected.  But for shape parameter between 0.16 and 1:

    >> qgamma(1,.5)
    J> [1] NaN
    J> Warning message:
    J> NaNs produced in: qgamma(p, shape, scale, lower.tail, log.p)

    J> and for shape parameter 0.16 or less:

    >> qgamma(1,.1)
    J> [1] 99.42507

    J> Arguments close to 1 give approximately correct results:

    >> qgamma(1-1e-15,.5)
    J> [1] 32.05055
    >> qgamma(1-1e-15,.1)
    J> [1] 28.8129