[R] new dgamma rate argument

Jim Lindsey james.lindsey at luc.ac.be
Fri Jan 11 14:36:25 CET 2002


> 
> >>>>> "Jim" == Jim Lindsey <james.lindsey at luc.ac.be> writes:
> 
>     Jim> Can someone explain to me in what way the new
>     Jim> (dpqr)gamma parameter can be interpreted as a rate
>     Jim> (when shape != 1)? The only gamma rate that I am aware
>     Jim> of is the hazard rate given by dgamma/(1-pgamma), the
>     Jim> log of which is returned by my hgamma function (event
>     Jim> library).  Jim
> 
> NEWS has
> 
>     o	[dpqr]gamma now has third argument `rate' for S-compatibility
> 	(and for compatibility with exponentials).  Calls which use
> 	positional matching may need to be altered.
> 
> i.e. one point of view (close to mine) could be:
> 
>  The authors of R (R&R) called that argument of [dpqr]gamma()
>  `scale' as it should sensibly be called.
>  OTOH, (at least one of) the original S authors used `rate'
>  (for 1/scale) in a loose analogy with the exponential and
>  weibull distribution quite some time before R was born.  Now
>  that there is an increasing drive for S source compatibility
>  between the different S dialects --  whenever it's ``easy'' --
>  the compatible parametrization has been allowed as well.

That does not answer my question. Of course, I read the NEWS but a
misleading error in one system is not a reason for copying it into R.
There is no rate parameter in the Weibull either, only
dweibull/(1-pweibull) as well. Only the exponential has a constant
rate.
  Jim

> 
> Martin
> -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
> r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
> Send "info", "help", or "[un]subscribe"
> (in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
> _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
> 

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list