From bolker at ufl.edu  Tue Mar 31 17:08:45 2009
From: bolker at ufl.edu (Ben Bolker)
Date: Tue, 31 Mar 2009 15:08:45 +0000 (UTC)
Subject: [Rd] Gamma funtion(s) bug
References: <49D12D9F.1000203@stats.uwo.ca>
	<XFMail.090330222854.Ted.Harding@manchester.ac.uk>
	<18897.49007.275534.467989@lynne.math.ethz.ch>
Message-ID: <loom.20090331T150244-162@post.gmane.org>

Martin Maechler <maechler <at> stat.math.ethz.ch> writes:

>     >> But lgamma(x) is log(abs(gamma(x))), so it looks okay to me.
>     >> 
>     >> Duncan Murdoch
> 
>     TH> Oops, yes! That's what comes of talking off the top of my head
>     TH> (I don't think I've ever had occasion to evaluate lgamma(x)
>     TH> for negative x, so never consciously checked in ?lgamma).
> 
>     TH> Thanks, Duncan!
> 
> Indeed.... as we all know, a picture can be worth a thousand words,
> and a simple R call such as
>        plot(lgamma, -7, 0, n=1000)
> would have saved many words, and notably spared us from
> yet-another erroneous non-bug report.
> 
> Martin

 In Kjetil's defense, he didn't submit an actual bug report --
and although his subject line does contain the word "bug",
I read his "bug report" as asking a question.  People are
allowed to make mistakes ...

  While I was reading ?lgamma I noticed that the "See Also"
section refers to gammaCody(), which is now defunct.  Perhaps
remove the sentence?

  Ben Bolker



[Rd] Gamma funtion(s) bug
Ben Bolker