Gamma funtion(s) bug

On 30-Mar-09 18:40:03, Kjetil Halvorsen wrote:

With R 2.8.1 on ubuntu I get:

gamma(-1)

[1] NaN
Warning message:
In gamma(-1) : NaNs produced

lgamma(-1)

[1] Inf
Warning message:
value out of range in 'lgamma'

Is'nt the first one right, and the second one (lgamma)
should also be NaN?
Kjetil

That is surely correct! Since lim[x->(-1)+] gamma(x) = +Inf,
while lim[x->(-1)-] gamma(x) = -Inf, at gamma(-1) one cannot
choose between +Inf and -Inf, so surely is is NaN.