[R] Problem with dgamma ?

ecatchpole
Mon, 23 Apr 2007 04:41:03 -0000

dgamma(x=1, shape=nu, rate=nu, log=TRUE)
[1] -2.045951

This is a good example of why you should call parameters by name.

Ted.

Tong Wang wrote on 04/23/2007 01:59 PM:
> Hi All,
>      Here 's what I got using dgamma function :
>
>
>> nu<-.2
>> nu*log(nu)-log(gamma(nu))+(nu-1)*log(1)-nu*(1)
>>
> [1] -2.045951
>
>
>> dgamma(1,nu,nu,1)
>>
> [1] 0.0801333
>
>
>> dgamma(1,nu,nu,0)
>>
> [1] NaN
> Warning message:
> NaNs produced in: dgamma(x, shape, scale, log)
>
> Could anyone tell me what is wrong here ?
> I am using R-2.4.1 on windows XP.
>
> Thanks a lot.
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>


--
 Dr E.A. Catchpole
 Visiting Fellow
 Univ of New South Wales at ADFA, Canberra, Australia
    _	  and University of Kent, Canterbury, England
   'v'	  - www.pems.adfa.edu.au/~ecatchpole
  /   \	  - fax: +61 2 6268 8786
   m m    - ph:  +61 2 6268 8895